Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*11,451,1240CA38.9% 44.6 / 25.6 36intergenic (+118602/+150676)unknown/unknowntRNA‑Met‑CAT/tRNA‑Val‑GAC
Reads supporting (aligned to +/- strand):  ref base C (21/1);  new base A (0/14);  total (21/15)
Fisher's exact test for biased strand distribution p-value = 3.95e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TATTAATTTGGGTGGTGGAATCAATCTCGGCAATTAAATCGCCTTTTTTCACCTGTTGTCCCAGTTTGACATAAAGTTTGGTGACTTT‑CCCTGACACCTGCGCGCCCACATCAACAGTATTGATACTTTCGAT‑GGAACCGGTAGCAACCACGTTCTTTTCAATATTGCCGCGTTTGACCTCTTCC‑GTTAAGTAAGTGGTTT‑GTGCGCTGGAATCGCGTACAAAAAAATAACCGCCGGTAATCACAGCAGCTGCAATAGCAGTAATAATAAGAGGTTTTTTCTTCATCTTTTTCTCTCATAAGTCCACAATGAATCGCCACACTATACAACAAAACTTAAGAAATTATTAAAATTTTTAGCGTGTTGCCATTATTCTCCCAATGATACTCAATTAAAAATAAGGGTAGTTTGCACATTTATTTCAACAATAAGCCGCATTCCTCGCTTTTTATGCTGAAG  >  1/1450905‑1451363
                                                                                                                                                                                                                               |                                                                                                                                                                                                                                               
tATTAATTTGGGTGGTGGAATCAATCTCGGCAATTAAATCGCCTTTTTTCACCTGTTGTCCCAGTTTGACATAAAGTTTGGTGACTTT‑CCCTGACACCTGCGCGCCCACATCAACAGTATTGATACTTTCGAT‑GGAACCGGTAGCAACCACGTTCTTTTCAATATTGCCGCGTTTGACCTCTTCC‑GTTAAGTAAGTGGTTT‑GTGCGCTGGAATCGCGTACAAAAAAATAACCGCCGGTAATCAcagcag                                                                                                                                                                                                                    >  1:735312/1‑249 (MQ=255)
    aaTTTGGGTGGTGGAATCAATCTCGGCAATTAAATCGACTTTTTTCACCTGTTGTCCCAGTTTGACATAAAGTTTGGTGACTTT‑CCCTGACACCTGCGCGCCCACATCAACAGTATTGATACTTTCGAT‑GGAACCGGTAGCAACCACGTTCTTTTCTATATTGCCGCGTTTGACCTCTTCC‑GTTAAGTAAGTGGTTT‑GTGCGCTGGAATCGCGTACAAAAAAATAACCGCCGGTAATCACAGCAGCTGCaa                                                                                                                                                                                                              >  1:887223/1‑251 (MQ=255)
           gtggtgGAATCATCCTCGGCAATTAATTCCCTTTTTTTCCCTTGTTGCCCCATTTTAACATAATGTTTGGTG‑CTTTCCCCTGCCCCCGGCGCCCCCACATCAACATTTTTATTTTTTTGGAG‑GGACCCGGTAGCACCCCCTTTTTTTTAAATTTGGCCGCTTTTGCCCTCTCCC‑TTTAAGTAGGTGGTTT‑GGGCGCGGAAACCGCGTAAAAAAAAATAACCGCCGGTAATCACAGCAGCTGCAATAGCAg                                                                                                                                                                                                        <  1:175128/250‑1 (MQ=255)
                       ctctcGGAAATTAAATCGCTTTTTTTCCCCTTTTGCCCCATTTTGACAAAAATTTGGGTTACTTT‑CCCTGCCCCCGGCGCCCCCACATCAACAGTATTGATCCTTTCGAT‑GGAACCGGTAGCACCCACTTTTTTTTAAATATTGCCGCGTTTGCCCTCTTCC‑GTTAAGTAGGTGGTTT‑GTGCGCGGGAATCGCGTAAAAAAAAATAACCGCCGGTAATCACAGCAGCTGCAATAGCAGTAATAATAAGAgg                                                                                                                                                                                           <  1:535260/250‑1 (MQ=255)
                           cGGCAATTAAATCGCCTTTTTTCACCTGTTGTCCCAGTTTGACATAAAGTTTGGTGACTTT‑CCCTGACACCTGCGCGCCCACATCAACAGTATTGATACTTTCGAT‑GGAACCGGTAGCAACCACGTTCTTTTCAATATTGCCGCGTTTGACCTCTTCC‑GTTAAGTAAGTGGTTT‑GTGCGCTGGAATCGCGTACAAAAAAATAACCGCCGGTAATCACAGCAGCTGCAATAGCAGtaataata                                                                                                                                                                                                >  1:256640/1‑242 (MQ=255)
                                aTTAAATCGCCTTTTTTCACCTGTTGTCCCAGTTTGACATAAAGTTTGGTGACTTT‑CCCTGACACCTGCGCGCCCACATCAACAGTATTGATACTTTCGAT‑GGAACCGGTAGCAACCACGTTCTTTTCAATATTGCCGCGTTTGACCTCTTCC‑GTTAAGTAAGTGGTTT‑GTGCGCTGGAATCGCGTACAAAAAAATAACCGCCGGTAATCACAGCAGCTGCAATAGCAGTAATAATAAGAGGTTTTttctt                                                                                                                                                                                  >  1:702675/1‑251 (MQ=255)
                                         ccTTTTTTCACCTGTTGTCCCAGTTTGACATAAAGTTTGGTGACTTT‑CCCTGACACCTGCGCGCCCACATCAACAGTATTGATACTTTCGAT‑GGAACCGGTAGCAACCACGTTCTTTTCAATATTGCCGCGTTTGACCTCTTCC‑GTTAAGTAAGTGGTTT‑GTGCGCTGGAATCGCGTACAAAAAAATAACCGCCGGTAATCACAGCAGCTGCAATAGCAGTAATAATAAGAGGTTTTTTCTTCATCttttt                                                                                                                                                                         >  1:82538/1‑251 (MQ=255)
                                              tttCACCTGTTGTCCCAGTTTGACATAAAGTTTGGTGACTTT‑CCCTGACACCTGCGCGCCCACATCAACAGTATTGATACTTTCGAT‑GGAACCGGTAGCGACCACGTTCTTTTCAATATTGCCGCGTTTGACCTCTTCC‑GTTAAGTAAGTGGTTT‑GTGCGCTGGAATCGCGTACAAAAAAATAACCGCCGGTAATCACAGCAGCTGCAATAGCAGTAATAATAAGAGGTTTTTTCTTCATCTTTTtctct                                                                                                                                                                     >  1:785732/1‑250 (MQ=255)
                                                   cttgttgcCCCATTTTAACATAATGTTTGGTTTCTTT‑CCCTGCCCCCTGCGCCCCCAAATAAAAATTTTTGATTTTTTGGAT‑GGACCGGGTAGCAACCTCTTTTTTTTAAATTTTGCCGTTTTTGCCCT‑TTCCTGTTAAGTAGGGGGTTT‑GTGGGGTGAAAGCGCGTAAAAAAAAATAACCGCCGGTAATCACAGCAGCTGCAATAGCAGTAATAATAAGAGGTTTTTTCTTCATCTTTTTCTCTCATAAg                                                                                                                                                               <  1:15885/249‑1 (MQ=255)
                                                               gTTTGACATAAAGTTTGGTGACTTT‑CCCTGACACCTGCGCGCCCACATCAACAGTATTGATACTTTCGAT‑GGAACCGGTAGCAACCACGTTCTTTTCAATATTGCCGCGTTTGACCTCTTCC‑GTTAAGTAAGTGGTTT‑GTGCGCTGGAATCGCGTACAAAAAAATAACCGCCGGTAATCACAGCAGCTGCAATAGCAGTAATAATAAGAGGTTTTTTCTTCATCTTTTTCTCTCATAAGTCCACAATGAAt                                                                                                                                                   >  1:528935/1‑251 (MQ=255)
                                                                                 ttcccttccccTGACCCCTGCCCCCCCACATCAACAGTATTGATACTTTGGAT‑GGACCGGGTAGAACCCTCTTTCTTTTAAATATTCCCGCGTTTGCCCTCTTCC‑TTTAAGTAGGGGGTTG‑GTGCGCTGAAAGCGCGAAAAAAAAAAAATCCGCCGGTAATCACAGCAGCTGCAATAGCAGTAATAATAAGAGGTTTTTTCTTCATCTTTTTCTCTCATAAGTCCACAATGAATCGCCACACTATacaac                                                                                                                                   <  1:403679/242‑1 (MQ=255)
                                                                                 ttaCTTT‑CCCTGACCCCGGCGCCCCAACAACAACAGTATTGATATTTTCGAT‑GGAACGGGTAGCAACCGCTTTTTTTTCAATATTGCCGCGTTTGCCCTCTTCC‑GTTAAGTAAGTGTTTG‑GTGCGGTGGAATTGTGAACAAAAAAAAAGCCGCCGGTAATCACAGCAGCTGCAATAGCAGTAATAATAAGAGGTTTTTTCTTCATCTTTTTCTCTCATAAGTCCACAATGAATCGCCACACTATAcaacaa                                                                                                                                 <  1:406803/249‑1 (MQ=255)
                                                                                 gtgccttccccTGACCCCTGCGCCCCCACAAAAACAGTATTGATTTTTTAGAT‑GAACCCGGAAGCACCCACTTTTTTTAAAATATGGCCGTGTTTGCCCTCTCCC‑GTAAAGGAAG‑GGTTTGGGGCGCGGGAATCGCGAAAAAAAAAAAAACCGCCGGTAATCACAGCAGCTGCAATAGCAGTAATAATAAGAGGTTTTTTCTTCATCTTTTTCTCTCATAAGTCCACAATGAATCGCCACACTATacaac                                                                                                                                   <  1:437108/242‑1 (MQ=255)
                                                                                         cccTGACCCCTGCGCGCCCACATCAACATTTTTGATATTTTCGAT‑GGACCCGGTAGCACCCACTTTTTTTTAAATATTGCCGCTTTTGCCCTCTTCC‑GTTAAGTAGGTGGTTT‑GTGCGCTGGAATCGCGTAAAAAAAAATAACCGCCGGTAATCACAGCAGCTGCAATAGCAGTAATAATAAGAGGTTTTTTCTTCATCTTTTTCTCTCATAAGTCCACAATGAATCGCCacac                                                                                                                                           <  1:210317/234‑1 (MQ=255)
                                                                                             cccccctgcccccccaaataaaCAGTATTGATTTTTTCGAT‑GGACCCGGTAGCACCCTTTTTTTTTTAAATTTTGCCGCTTTTGCCCTCTTCC‑TTAAAGTAGGTGGTTT‑GTGCGCGGGAATCGCGTAAAAAAAAATAACCGCCGGTAATCACAGCAGCTGCAATAGCAGTAATAATAAGAGGTTTTTTCTTCATCTTTTTCTCTCATAAGTCCACAATGAATCGCCACACTATACAACAAAACTTAAGaaa                                                                                                                      <  1:727897/231‑1 (MQ=255)
                                                                                                         cccACATCAACAGTATTGATACTTTCGAT‑GGAACCGGTAGCAACCACGTTCTTTTCAATATTGCCGCGTTTGACCTCTTCC‑GTTAAGTAAGTGGTTT‑GTGCGCTGGAATCGCGTACAAAAAAATAACCGCCGGTAATCACAGCAGCTGCAATAGCAGTAATAATAAGAGGTTTTTTCTTCATCTTTTTCTCTCATAAGTCCACAATGAATCGCCACACTATacaac                                                                                                                                   >  1:754424/1‑226 (MQ=255)
                                                                                                           cacaTCAACAGTATTGATACTTTCGAT‑GGAACCGGTAGCAACCACGTTCTTTTCAATATTGCCGCGTTTGACCTCTTCC‑GTTAAGTAAGTGGTTT‑GTGCGCTGGAATCGCGTACAAAAAAATAACCGCCGGTAATCACAGCAGCTGCAATAGCAGTAATAATAagag                                                                                                                                                                                            >  1:18176/1‑167 (MQ=255)
                                                                                                             caTCAACAGTATTGATACTTTCGAT‑GGAACCGGTAGCAACCACGTTCTTTTCAATATTGCCGCGTTTGACCTCTTCC‑GTTAAGTAAGTGGTTT‑GTGCGCTGGAATCGCGTACAAAAAAATAACCGCCGGTAATc                                                                                                                                                                                                                           >  1:832924/1‑134 (MQ=255)
                                                                                                                     gTATTGATACTTTCGAT‑GGAACCGGTAGCAACCACGTTCTTTTCAATATTGCCGCGTTTGACCTCTTCC‑GTTAAGTAAGTGGTTT‑GTGCGCTGGAATCGCGTACAAAAAAATAACCGCCGGTAATCACAGCAGCTGCAATAGCAGTAATAATAAGAGGTTTTTTCTTCATCTTTTTCTCTCATAAGTCCACAATGAATCGCCACACTATACAACAAAACTTAAGAAATTATTAAAATTTTTAGCGTGTTg                                                                                               >  1:594922/1‑250 (MQ=255)
                                                                                                                              tcttCGATGGGACCCGGAAGCAACCTCTTTCTTTTAAATATTGCCGGTTTTCCCCTCTCCG‑GTTAAGTAGGTGTTTT‑GTGCGCTGGAAGCGCTTAAAAAAAAAAAGCCGCCGGTAATCACAGCAGCTGCAATAGCAGTAATAATAAGAGGTTTTTTCTTCATCTTTTTCTCTCATAAGTCCACAATGAATCGCCACACTATACAACAAAACTTAAGAAATTATTAAAATTTTTAGCGTGTTGCCATTATtc                                                                                      <  1:224867/249‑1 (MQ=255)
                                                                                                                              cTTTCGAT‑GGAACCGGTAGCAACCACGTTCTTTTCAATATTGCCGCGTTTGACCTCTTCC‑GTTAAGTAAGTGGTTT‑GTGCGCTGGAATCGCGTACAAAAAAATAACCGCCGGTAATCACAGCAGCTGCAATAGCAGtaat                                                                                                                                                                                                    >  1:272191/1‑140 (MQ=255)
                                                                                                                                         accgggTAGCAACCACTTTTTTTTAAATATGGCCGTGTTTGCCCTCTTCC‑GTTAAGTAGGTGGTTT‑GTGCGCTGAAATCGCGAAAAAAAAAATAACCGCCGGTAATCACAGCAGCTGCAATAGCAGTAATAATAAGAGGTTTTTTCTTCATCTTTTTCTCTCATAAGTCCACAATGAATCGCCACACTATacaac                                                                                                                                   <  1:50827/191‑1 (MQ=255)
                                                                                                                                         aaCCGGTAGCAACCACGTTCTTTTCAATATTGCCGCGTTTGACCTCTTCC‑GTTAAGTAAGTGGTTT‑GTGCGCTGGAATCGCGTACAAAAAAATAACCGCCGGTAATCACAGCAGCTGCAATAGCAGTAATAATAAGAGGTTTTTTCTTCATCTTTTTCTCTCATAAGTCCACAATGAATCGCCacac                                                                                                                                           >  1:878531/1‑187 (MQ=255)
                                                                                                                                                              tttCAATATTGCCGCGTTTGACCTCTTCC‑GTTAAGTAAGTGGTTT‑GTGCGCTGGAATCGCGTACAAAAAAATAACCGCCGGTAATCACAGCAGCTGCAATAGc                                                                                                                                                                                                          >  1:332401/1‑103 (MQ=255)
                                                                                                                                                                 cAATATTGCCGCGTTTGACCTCTTCC‑GTTAAGTAAGTGGTTT‑GTGCGCTGGAATCGCGTACAAAAAAATAACCGCCGGTAATCACAGCAGCTGCAATAGCAGTAATAATAAGAGGTTTTTTCTTCATCTTTTTCTCTCATAAGTCCACAATGAATCGCCACACTATAcaacaa                                                                                                                                 >  1:40478/1‑173 (MQ=255)
                                                                                                                                                                                    ctctTCC‑GTAAAGTAAGTGTTTT‑GTGCGCGGGAACCGCGTAAAAAAAAATAACCGCCGGTAATCACAGCAGCTGCAATAGCAGTAATAATAAGAGGTTTTTTCTTCATCTTTTTCTCTCATAAGTCCACAATGAATCGCCACACTATACAACAAAACTTAAGAAATTATTAAAATTTTTAGCGTGTTGCCATTATTCTCCCAATGa                                                                             <  1:56299/206‑1 (MQ=255)
                                                                                                                                                                                           ctttAAGTAGGTGTTTT‑TTGCGCTGAAATCGCGTAAAAAAAAATAACCGCCGGTAATCACAGCAGCTGCAATAGCAGTAATAATAAGAGGTTTTTTCTTCATCTTTTTCTCTCATAAGTCCACAATGAATCGCCACACTATACAACAAAACTTAAGAAATTATTAAAATTTTTAGCGTGTTGCCATTATTCTCCCAATGAt                                                                            <  1:685730/199‑1 (MQ=255)
                                                                                                                                                                                            gTTAAGTAAGTGGTTT‑GTGCGCTGAAATCGCGAAAAAAAAAATAACCGCCGGTAATCACAGCAGCTGCAATAGCAGTAATAATAAGAGGTTTTTTCTTCATCTTTTTCTCTCATAAGTCCACAATGAATCGCCACACTATacaac                                                                                                                                   <  1:855885/145‑1 (MQ=255)
                                                                                                                                                                                             ttAAGTAAGTGGTTT‑GTGCGCTGGAATCGCGTACAAAAAAATAACCGCCGGTAATCACAGCAGCTGCAATAGCAGTAATAATAAGAGGTTTTTTCTTCATCTTTTTCTCTCATAAGTCCACAATGAATCGCCACACTATACAACAAAACTTAAGAAATTATTAAAATTTTTAGCGTGTTGCCATTATTCTCCCAATGATACTCAATTAAAAATAAGGGTAGTTTGCACATTTATTTcaacaa                                 >  1:710502/1‑242 (MQ=255)
                                                                                                                                                                                                agtaagtGGTTT‑GTGCGCTGGAATCGCGTACAAAAAAATAACCGCCGGTAATCACAGCAGCTGCAATAGCAGTAATAATAAGAGGTTTTTTCTTCATCTTTTTCTCTCATAAGTCCACAATGAATCGCCACACTATACAACAAAACTTAAGAAATTATTAAAATTTTTAGCGTGTTGCCATTATTCTCCCAATGATACTCAATTAAAAATAAAGGTAGTTTGCACATTTATTTCAACAATAAGCCGCAtt                      >  1:515189/1‑250 (MQ=255)
                                                                                                                                                                                                 gtaagtGGTTT‑GTGCGCTGGAATCGCGTACAAAAAAATAACCGCCGGTAATCACAGCAGCTGCAATAGCAGTAATAATAAGAGGTTTTTTCTTCATCTTTTTCTCTCATAAGTCCACAATGAATCGCCACACTATACAACAAAACTTAAGaa                                                                                                                       >  1:703714/1‑152 (MQ=255)
                                                                                                                                                                                                       ggTTT‑GTGCGCTGGAATCGCGTACAAAAAAATAACCGCCGGTAATCACAGCAGCTGCAATAGCAGTAATAATAAGAGGTTTTTTCTTCATCTTTTTCTCTCATAAGTCCACAATGAATCGCCACACTATACAACAAAACTTAAGAAATTATTAAAATTTTTAGCGTGTTGCCATTATTCTCCCAATGATACTCAATTAAAAATAAGGGTAGTTTGCACATTTATTTCAACAATAAGCCGCATTCCTCGCtt              >  1:326357/1‑251 (MQ=255)
                                                                                                                                                                                                        tttt‑GTGTGCTGGAATCGCAAAAAAAAAAAAAACCGCCGGTAATCACAGCAGCTGCAATAGCAGTAATAATAAGAGGTTTTTTCTTCATCTTTTTCTCTCATAAGTCCACAATGAATCGCCACACTATACAACAAAACTTAAGAAATTATTAAAATTTTTAGCGTGTTGCCATTATTCTCCCAATGATACTCAATTAAAAATAAGGGTAGTTTGCACATTTATTTCAACAATAAGCCGCATTCCTCGCttt             <  1:218815/250‑1 (MQ=255)
                                                                                                                                                                                                             gtgCGCTGGAATCGCGTACAAAAAAATAACCGCCGGTAATCACAGCAGCTGCAATAGCAGTAATAATAAGAGGTTTTTTCTTCATCTTTTTCTCTCATAAGTCCACAATGAATCGCCACACTATACAACAAAACTTAAGAAATTATTAAAATTTTTAGCGTGTTGCCATTATTCTCCCAATGAt                                                                            >  1:872817/1‑184 (MQ=255)
                                                                                                                                                                                                               gcgcTGGAACCGCGAAAAAAAAAATAACCGCCGGTAATCACAGCAGCTGCAATAGCAGTAATAATAAGAGGTTTTTTCTTCATCTTTTTCTCTCATAAGTCCACAATGAATCGCCACACTATACAACAAAACTTAAGAAATTATTAAAATTTTTAGCGTGTTGCCATTATTCTCCCAATGATACTCAATTAAAAATAAGGGTAGTTTGCACATTTATTTCAACAATAAGCCGCATTCCTCGCTTTTTATg        <  1:782645/250‑1 (MQ=255)
                                                                                                                                                                                                                    ggAATCGCGTACAAAAAAATAACCGCCGGTAATCACAGCAGCTGCAATAGCAGTAATAATAAGAGGTTTTTTCTTCATCTTTTTCTCTCATAAGTCCACAATGAATCGCCACACTATACAACAAAACTTAAGAAATTATTAAAATTTTTAGCGTGTTGCCATTATTCTCCCAATGATACTCAATTAAAAATAAGGGTAGTTTGCACATTTATTTCAACAATAAGCCGCATTCCTCGCTTTTTATGCTGAAg  >  1:511983/1‑251 (MQ=255)
                                                                                                                                                                                                                               |                                                                                                                                                                                                                                               
TATTAATTTGGGTGGTGGAATCAATCTCGGCAATTAAATCGCCTTTTTTCACCTGTTGTCCCAGTTTGACATAAAGTTTGGTGACTTT‑CCCTGACACCTGCGCGCCCACATCAACAGTATTGATACTTTCGAT‑GGAACCGGTAGCAACCACGTTCTTTTCAATATTGCCGCGTTTGACCTCTTCC‑GTTAAGTAAGTGGTTT‑GTGCGCTGGAATCGCGTACAAAAAAATAACCGCCGGTAATCACAGCAGCTGCAATAGCAGTAATAATAAGAGGTTTTTTCTTCATCTTTTTCTCTCATAAGTCCACAATGAATCGCCACACTATACAACAAAACTTAAGAAATTATTAAAATTTTTAGCGTGTTGCCATTATTCTCCCAATGATACTCAATTAAAAATAAGGGTAGTTTGCACATTTATTTCAACAATAAGCCGCATTCCTCGCTTTTTATGCTGAAG  >  1/1450905‑1451363

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑