Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,395,4440AC39.7% 102.1 / 18.2 69G151G (GGT→GGG) nuoJNADH:quinone oxidoreductase subunit J
Reads supporting (aligned to +/- strand):  ref base A (31/10);  new base C (0/27);  total (31/38)
Fisher's exact test for biased strand distribution p-value = 5.21e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CCGGCGACCACGAAGGCCAGCGCGGAGGCGTTAATCATGATTTCCAGACCAATCAACATAAACAGCAGATTGCGACGGATAACCA‑GACCGGTTAAGCCAA‑GAACGAATAAGATTGCCGCGAGGATCAGTCCATGTT‑GTAAGGGG‑ATC‑ATGCGTGCTCCTCCGTTTTTCTTT‑T‑CGCGCTGT‑CGTCTTTACGATTGCTCAGCACTTCACCCGCACGCTCTTCACGACCGACGTGGAAGGCCACAACCAGACCTG‑CGAGCAGCAGCATAGAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACTGCTTTAGCACTGATTGGCGTACCGTCGATACCCTGATCGTTAACACCGAGGATGGCGTAAACAATCACCACCAGCATGATGGCCGACAAAATTGCCGGACCAATCCACACCTGCGGTTTCAGCCACTGGCGTTCCTGTTCGATTTCTGAACCGCCCAGGTTGAGCATCATC  >  NC_000913/2395197‑2395684
                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                 
ccGGCGACCACGAAGGCCAGCGCGGAGGCGTTAATCATGATTTCCAGACCAATCAACATAAACAGCAGATTGCGACGGATAACCA‑GACCGGTTAAGCCAA‑GAACGAATAAGATTGCCGCGAGGATCAGTCCATGTT‑GTAAGGGG‑ATC‑ATGCGTGCTCCTCCGTTTTTCTTT‑T‑CGCGCTGT‑CGTCTTTACGATTGCTCAGCACTTCACCCGCACGCTCTTCACGACCGACGTGGAAGGCCACAACCAGACCt                                                                                                                                                                                                                                                <  1:91935/251‑1 (MQ=255)
             aGGCCAGCGCGGAGGCGTTAATCATGATTTCCAGACCAATCAACATAAACAGCAGATTGCGACGGATAACCA‑GACCGGTTAAGCCAA‑GAACGAATAAGATTGCCGCGAGGATCAGTCCATGTT‑GTAAGGGG‑ATC‑ATGCGTGCTCCTCCGTTTTTCTTT‑T‑CGCGCTGT‑CGTCTTTACGATTGCTCAGCACTTCACCCGCACGCTCTTCACGACCGACGTGGACGGCCACAACCAGACCTG‑CGagcagcagc                                                                                                                                                                                                                                   >  2:584199/1‑250 (MQ=255)
              ggCCAGCGCGGAGGCGTTAATCATGATTTCCAGACCAATCAACATAAACAGCAGATTGCGACGGATAACCA‑GACCGGTTAAGCCAA‑GAACGAATAAGATTGCCGCGAGGATCAGTCCATGTT‑GTAAGGGG‑ATC‑ATGCGTGCTCCTCCGTTTTTCTTT‑T‑CGCGCTGT‑CGTCTTTACGATTGCTCAGCACTTCACCCGCACGCTCTTCACGACCGACGTGGAAGGCCACAACCAGACCTG‑CGAGCAGCAGCAt                                                                                                                                                                                                                                 <  1:228279/251‑1 (MQ=255)
               gCCAGCGCGGAGGCGTTAATCATGATTTCCAGACCAATCAACATAAACAGCAGATTGCGACGGATAACCA‑GACCGGTTAAGCCAA‑GAACGAATAAGATTGCCGCGAGGATCAGTCCATGTT‑GTAAGGGG‑ATC‑ATGCGTGCTCCTCCGTTTTTCTTT‑T‑CGCGCTGT‑CGTCTTTACGATTGCTCAGCACTTCACCCGCACGCTCTTCACGACCGACGTGGAAGGCCACAACCAGACCTG‑CGAGCAGCAGCAt                                                                                                                                                                                                                                 >  1:477482/1‑250 (MQ=255)
                 cAGCGCGGAGGCGTTAATCATGATTTCCAGACCAATCAACATAAACAGCAGATTGCGACGGATAACCA‑GACCGGTTAAGCCAA‑GAACGAATAAGATTGCCGCGAGGATCAGTCCATGTT‑GTAAGGGG‑ATC‑ATGCGTGCTCCTCCGTTTTTCTTT‑T‑CGCGCTGT‑CGTCTTTACGATTGCTCAGCACTTCACCCGCACGCTCTTCACGACCGACGTGGAAGGCCACAACCAGACCTG‑CGAGCAGCAGCAt                                                                                                                                                                                                                                 >  2:568323/1‑248 (MQ=255)
                 cAGCGCGGAGGCGTTAATCATGATTTCCAGACCAATCAACATAAACAGCAGATTGCGACGGATAACCA‑GACCGGTTAAGCCAA‑GAACGAATAAGATTGCCGCGAGGATCAGTCCATGTT‑GTAAGGGG‑ATC‑ATGCGTGCTCCTCCGTTTTTCTTT‑T‑CGCGCTGT‑CGTCTTTACGATTGCTCAGCACTTCACCCGCACGCTCTTCACGACCGACGTGGAAGGCCACAACCAGACCTG‑CGAGCAGCAGCAt                                                                                                                                                                                                                                 >  1:174003/1‑248 (MQ=255)
                    cgcgGAGGCGTTAATCATGATTTCCAGACCAATCAACATAAACAGCAGATTGCGACGGATAACCA‑GACCGGTTAAGCCAA‑GAACGAATAAGATTGCCGCGAGGATCAGTCCATGTT‑GTAAGGGG‑ATC‑ATGCGTGCTCCTCCGTTTTTCTTT‑T‑CGCGCTGT‑CGTCTTTACGATTGCTCAGCACTTCACCCGCACGCTCTTCACGACCGACGTGGAAGGCCACAACCAGACCTG‑CGAGCAGCAGCATAGAAGc                                                                                                                                                                                                                           <  2:307204/251‑1 (MQ=255)
                                   cATGATTTCCAGACCAATCAACATAAACAGCAGATTGCGACGGATAACCA‑GACCGGTTAAGCCAA‑GAACGAATAAGATTGCCGCGAGGATCAGTCCATGTT‑GTAAGGGG‑ATC‑ATGCGTGCTCCTCCGTTTTTCTTT‑T‑CGCGCTGT‑CGTCTTTACGATTGCTCAGCACTTCACCCGCACGCTCTTCACGACCGACGTGGAAGGCCACAACCAGACCTG‑CGAGCAGCAGCATAGAAGCCAGTTCCACCGcc                                                                                                                                                                                                              >  1:249113/1‑249 (MQ=255)
                                     tGATTTCCAGACCAATCAACATAAACAGCAGATTGCGACGGATAACCA‑GACCGGTTAAGCCAA‑GAACGAATAAGATTGCCGCGAGGATCAGTCCATGTT‑GTAAGGGG‑ATC‑ATGCGTGCTCCTCCGTTTTTCTTT‑T‑CGCGCTGT‑CGTCTTTACGATTGCTCAGCACTTCACCCGCACGCTCTTCACGACCGACGTGGAAGGCCACAACCAGACCTG‑CGAGCAGCAGCATAGAAGCCAGTTCCACCGCCAGt                                                                                                                                                                                                           >  2:363018/1‑250 (MQ=255)
                                      gATTTCCAGACTAATCAACATAAACAGCAGATTGCGCCGGATAACCA‑GACCGGTTAAGCCAA‑GAACGAATAAGCTTGCAGCGAGGATCAGTCCATGTT‑GTAAGGGG‑ATC‑ATGCGTGCTCCTCCGTTTTTATTT‑T‑CGCGCTGG‑CGTCTATACGATTGCTCAGCACTTCACCCGCACGCTCTTCACGACCGACGTGGAAGGCAACAACCAGATCTG‑CGAGCAGAAGCATAGAAGCCAGTTCAACCGCCAGTAc                                                                                                                                                                                                         >  1:179175/1‑251 (MQ=255)
                                        tttCAAGACCAAACAACATAAACATCAGATTGCGAAGAATCACCA‑GACCGGTAACGCCAA‑GAACGAAGAAGATTGCCGCGGGGATCATGCCATGTT‑GAAAGGTG‑ATC‑ATGCGTGCCCCTCCGTTTTTCTTT‑T‑CGCTCTGT‑CTTCTTTACGATCGCTCAGCCCTTCACCCGCACCCTCTTCGCGACCGAGGTGGAAGCCCACACCCAG‑CCTGGCGAGCAGCAGCATAGAAGCCAGTTCCACCGCCAGTACg                                                                                                                                                                                                        <  2:477482/250‑1 (MQ=255)
                                         ttCCAGACCAATCAACATAAACAGCAGATTGCGACGGATAACCA‑GACCGGTTAAGCCAA‑GAACGAATAAGATTGCCGCGAGGATCAGTCCATGTT‑GTAAGGGG‑ATC‑ATGCGTGCTCCTCCGTTTTTCTTT‑T‑CGCGCTGT‑CGTCTTTACGATTGCTCAGCACTTCACCCGCACGCTCTTCACGACCGACGTGGAAGGCCACAACCAGACCTG‑CGAGCAGCAGCATAGAAGCCAGTTCCACCGCCAGTACGt                                                                                                                                                                                                       >  1:488964/1‑250 (MQ=255)
                                           ccaacccaacaaaaaaaaaaagcagTTTGCGACGGAAACCCG‑GCCCGGTTAAGCAAA‑GAAAAAATAAGATTGCCGCGAGGATCAGTCC‑TTTT‑GGAAGGGGGATC‑TTGCGTGCCCCCCCTTTTTTCTTT‑T‑CGGGCTGT‑CGTCTTTACGATTGCTCACCACTCCCCCCGCACGCTCTCCACGCCGGAGGTGGAGGGCCACAACCGGCCCTG‑GGGGCAGCAGAATAGAAGCCAGTTCCACCGCCAGTACGTAAg                                                                                                                                                                                                    <  2:473772/231‑1 (MQ=255)
                                                ccaaaaaaaaTAAACAACAGATTGCGACGGATAACAA‑GCCCGTTAAAGCCAA‑GAACGAATAATTTTGCCGCGAGGATCAGTCCATGTT‑GAAGGGGG‑ATC‑ATGCTTGCCCCTCCTTTTTTCTTT‑T‑CGCGCTGT‑CGTCTTTACGATTCCTCACCCCTTCCCCCGCACGCTCTCCCCGACCGACGTGGAAGGCCACAACCAGCCCGG‑CGAGCAGCAGCATAGAAGCCATTTCCACCGCCAGTACGTAATGcccn                                                                                                                                                                                               <  2:285815/242‑2 (MQ=255)
                                                     cAACATAAACAGCAGATTGCGACGGATAACCA‑GACCGGTTAAGCAAA‑GAACGAATAAGATTGCCGCGAGGATCAGTCCATGTT‑GTAAGGGG‑ATC‑ATGCGTGCTCCTCCGTTTTTCTTT‑T‑CGCGCTGT‑CGTCTTTACGATTGCTCAGCACTTCCCCCGCACGCTCTTCACGACCGACGTGGAAGGCCACAACCAGCCCTG‑CGAGCAGCAGCATAGAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAg                                                                                                                                                                                          <  1:461830/251‑1 (MQ=255)
                                                       aCATAAACAGAAGATTGCGACGGATAACCA‑GACCGGTTAAGCCAA‑GAACGAATAAGATTGCCGCGAGGATCAGTCCATGTT‑GTAAGGGG‑ATC‑ATGCGTGCTCCTCCGTTTTTCTTT‑T‑CGCGCTGT‑CGTCTTTACGATTGCTCAGCACTTCACCCGCACGCTCTTCACGACCGACGCGGAGGGCAACCACCAGACCTG‑CGAGCAGCAGCATAGAAGCCAGTTCCACCGCCACTACGTCAGGCCCGAACAGCg                                                                                                                                                                                        >  2:183813/1‑251 (MQ=255)
                                                          aaaaaagcagATGGCGACGGA‑AACCAAGCCCGGTTAAGCAAA‑GAACGAATAAGATTGCCGGGAGGATCAGTCCATGTT‑GTAAGGGA‑ATC‑GTGCGTCCCCCTCCTTTTTCCTTT‑C‑CGCTCTGC‑CGTCTTTACGATTGCCCAGCACTTCCCCCGCCCCCTCTACCCGCCCGAGGTGGAAGGCCACAACCGGCCCTG‑CGAGCAGCAGCATAGAAGCCAGTTCCACCGCCAGTACGGAAGGCCCGAACAGCGTaa                                                                                                                                                                                     <  2:488964/246‑1 (MQ=255)
                                                          aaaaCAGCAGATTGCGACGGATAACCA‑GACCGGTTAAGCCAA‑GAACGAATAAGATTGCCGCGAGGATCAGTCCATGTT‑GTAAGGGG‑ATC‑ATGCGTGCTCCTCCTTTTTTCTTT‑C‑CGCGCTGT‑CGTCTTTACGATTGCTCAGCCCTTCCCCCGCACGCTCTTCACGACCGACGTGGAAGGCCACAACCAGCCCTG‑CGAGCAGCAGCATAGAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTaa                                                                                                                                                                                     <  1:29477/250‑1 (MQ=255)
                                                                 aagATTGCGACGAATAACCA‑GACCGGTTAAGCCAA‑GAACGAATAAGATTGCCGCGAGGATCAGTCCATGTT‑GTAAGGGG‑ATC‑ATGCGTGCTCCTCCTTTTTTCTTT‑T‑CGCGCTGT‑CGTCTTTACGATTCCTCAGCCCTTCCCCCGCACGCTCTTCACGACCGACGGGGAAGGCCACAACCAGCCCTG‑CGAGCAGCAGCATAGAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCAc                                                                                                                                                                              <  1:175423/250‑1 (MQ=255)
                                                                  agATTGCGAGGGATAACCA‑G‑CCCGTTAAACCAAAGAACGAAAAAGATGGCCGGGAGGATCA‑TCCCTTTT‑TGAAGGGG‑GTCCATGCTTGCCCCCCCTTTTTTTTTT‑C‑CGCGCGGC‑CGTTTTTACGATTGCCCAGCACTCCCCCCGCCCGCTCTTCCCGCCCGACGGGGAGGGCCACACCCAGCCCGG‑CGAGCAGCAGCATAAAAGCCAGTCCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACt                                                                                                                                                                             <  2:238578/251‑1 (MQ=255)
                                                                  agATTGCGACGGATACCCA‑GACCGGTTAAGCCAA‑GAACGAATAAGATTGCCGCGAGGATCAGTCCATGTT‑GTAAGGGG‑ATC‑ATGCGTGCTCCTCCGTTTTTCTTT‑T‑CGCGCTGT‑CGTCTTTACGATTGCTCAGCCCTCCCCCCGCACGCTCTTCACGACCGACGTGGAAGGCCACACCCAGCCCTG‑CGAGCAGCAGCATAGAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACt                                                                                                                                                                             <  1:568323/251‑1 (MQ=255)
                                                                  agATTGCGACGGATAACCA‑GACCGGTTAAGCCAA‑GAACGAATAAGATTGCCGCGAGGATCAGTCCATGTT‑GTAAGGGG‑ATC‑ATGCGTGCTCCTCCGTTTTTCTTT‑T‑CGCGCTGT‑CGTCTTTACGATTGCTCAGCACTTCACCCGCACGCTCTTCACGACCGACGTGGAAGGCCACAACCAGACCTG‑CGAGCAGCAGCATAGAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGCAATACCCACt                                                                                                                                                                             >  2:357384/1‑251 (MQ=255)
                                                                          ccgaaaaaCCA‑GACCGGTTAATCAAA‑GAACAAATAAGATTGCCGCGAG‑ATCAATTCCTGTTTGAATGGGC‑ATC‑TTGCGTCCCCCTC‑GTTTTTTTTTCT‑CGCGCTGT‑CGTCTTTACGATTGCTCAGCACTCCCCCCGCACGCTCTTCACGCCCGACGTGGAAGGCCACACCCAGCCCTG‑CGACCAGCAGAATAAAAGCCATTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACTGCTTTAGc                                                                                                                                                                     <  1:377835/245‑1 (MQ=255)
                                                                            gaaTAACCA‑ACCCGGTAAAGCAAA‑GAACGAATAAGATGGCCGCGAGGACCAGTCCATGTT‑GTAAGGGG‑ATC‑ATGCGTCCCCCTCCGTTTTTCTTT‑C‑CGCGCTGT‑CGTCTTTCCGATTGCTCACCCCTCCCCCCGCACGCTCTTCCCGACCGCCGTGGAAGGCCACAACCAGACCCG‑CGGGCAGCAGCAAAAAACCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACTGCTTTAGCAc                                                                                                                                                                   <  1:542505/249‑1 (MQ=255)
                                                                                aaaccagcccttttAAGCCAA‑GAACAAATAAAGTTGCCGAGGGAATCAGTCCATGTT‑GTGGGGGG‑ATC‑ATGCGTGCTCCTCCGTTTTTCTTT‑T‑CGCG‑TGTTCGTCTTTCCGATGGCTCACCCCTCCCCCCGCACGCTCTCCACGACCGGGGTGGAAGCCCACAACCGGCCCGG‑CGACCAGCAGCAGAGAAGCCTGTCCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACTGCTTTAGCACTGa                                                                                                                                                                <  2:183243/239‑1 (MQ=255)
                                                                                  ccA‑GACCGGTTAAGCCAA‑GAACGAATAAGATTGCCGCGAGGATCAGTCCATGTT‑GTAAGGGG‑ATC‑ATGCGTGCTCCTCCGTTTTTCTTT‑T‑CGCGCTGT‑CGTCTTTACGATTGCTCAGCACTTCACCCGCACGCTCTTCACGACCGACGTGGAAGGCCACAACCAGACCTG‑CGAGCAGCAGCATAGAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACTGCTTTAGCACTGATTg                                                                                                                                                             >  1:253217/1‑251 (MQ=255)
                                                                                           ggtAAGCCAAAGAACAAAAAAGATGCCCGGGAGGATCAGTCCGTGTT‑GTAAGGGG‑ATC‑GTGCGTGCTCCCCCTTTTTTCTTT‑T‑CGCGCTGT‑CGTCTTTACGATTGCTCAGCCCTCCACCCGCACGCTCTTCACGCCCGACGTGGAAGGCCACAACCAGCCCTG‑CGAGCAGCACAATAGAACCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACTGCTTTAGCACTGATTGGCGTAcc                                                                                                                                                      <  1:140595/249‑1 (MQ=255)
                                                                                             aaaGCCAA‑GAACGAAAAAGATTGCCGCGAGGATCAGTCCATGTT‑GTAGGGGG‑ATC‑ATGCGTGCCCCCCCGTTTTTCTTT‑T‑CGCGCTGT‑CGTTTTTACGATTGCTAAGCACTCCCCCCGCCCGCTCTTCACGCCCGACGGGGAAGGCCACACCCAGCCCTG‑CGAGCAGCAGAATAGAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACTGCTTTAGCACTGATTGGCGTACCGTc                                                                                                                                                   <  1:239428/250‑1 (MQ=255)
                                                                                                 gcaaagaaaaaaTAAGATGGCCGCGAGGATCAGTCCATGTG‑GTAAGGGG‑ATC‑ATGCGTGCTCCTCCTTTTTTCTTT‑T‑CGCGCTGT‑CGTTTTTACGATTGCTCAGCACTTCCCCCGCACGCTCTTCCCGACCGACGTGGAAGGCCACACCCAGACCTG‑CGAGCAGCAGCATAGAAGCCAGTCCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACTGCTTTAGCACTGATTGGCGTACCGTCGAt                                                                                                                                                <  1:189441/241‑1 (MQ=255)
                                                                                                          gAATAAGATTGCCGCGAGGATCAGTCCATGTT‑GTAAGGGG‑ATC‑ATGCGTGCTCCTCCGTTTTTCTTT‑T‑CGCGCTGT‑CGTCTTTACGATTGCTCAGCACTTCACCCGCACGCTCTTCACGACCGACGTGGAAGGCCACAACCAGACCTG‑CGAGCAGCAGCATAGAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACTGCTTTAGCACTGATTGGCGTACCGTCGATACCCTGATc                                                                                                                                       >  2:539320/1‑251 (MQ=255)
                                                                                                          aaaaaaGATCGCCGCGGGTATCAGTCCATGTT‑GTAGGGGG‑TTC‑AGGCGGTCTCCTCCGTTTT‑CTTT‑TGCCCGCTGT‑CGTCTTTACTTTGGCTCACCCCTCCACCCGCACGCCCTTCACGACCGAGGTGAAAGGCCCAAACCAGCCCTG‑CGAGCAGCAGCATAGAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACTGCTTTAGCACTGATTGGCGTACCGTCGATACCCTGATc                                                                                                                                       <  2:338768/247‑1 (MQ=255)
                                                                                                                gATTGCCGCGAGGATCAGTCCATGTT‑GTAAGGGG‑ATC‑ATGCGTGCTCCTCCGTTTTTCTTT‑T‑CGCGCTGT‑CGTCTTTACGATTGCTCAGCACTTCACCCGCACGCTCTTCACGACCGACGTGGAAGGCCACAACCAGACCTG‑CGAGCAGCAGCATAGAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACTGCTTTAGCACTGATTGGCGTACCGTCGATACCCTGATCGTTAac                                                                                                                                 >  2:16550/1‑251 (MQ=255)
                                                                                                                   tGCCGCGAGGATCAGTCCATGTT‑GTAAGGGG‑ATC‑ATGCGTGCTCCTCCTTTTTTCTTT‑T‑CGCGCTGT‑CGTCTTTACGATTGCTCAGCACTTCCCCCGCCCGCTCTTCACGACCGACGTGGAAGGCCACAACCAGCCCTG‑CGAGCAGCAGCATAGAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACTGCTTTAGCACTGATTGGCGTACCGTCGATACCCTGATCGTTAacac                                                                                                                               <  1:127271/250‑1 (MQ=255)
                                                                                                                    gCCGCGAGGATCAGTCCATGTT‑GTAAGGGG‑ATC‑ATGCGTGCTCCTCCGTTTTTCTTT‑T‑CGCGCTGT‑CGTCTTTACGATTGCTCAGCACTTCACCCGCACGCTCTTCACGACCGACGTGGAAGGCCACAACCAGACCTG‑CGAGCAGCAGCATAGAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACTGCTTTAGCACTGATTGGCGTACCGTCGATACCCTGATCGTTAACACCg                                                                                                                             >  1:425699/1‑251 (MQ=255)
                                                                                                                              tCAGTCCATGTT‑GTAAGGGG‑ATC‑ATGCGTGCTCCTCCGTTTTTCTTT‑T‑CGCGCTGT‑CGTCTTTACGATTGCTCAGCACTTCACCCGCACGCTCTTCACGACCGACGTGGAAGGCCACAACCAGACCTG‑CGAGCAGCAGCATAGAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACTGCTTTAGCACTGATTGGCGTACCGTCGATACCCTGATCGTTAACACCGAGGATGGCGt                                                                                                                   >  2:168860/1‑251 (MQ=255)
                                                                                                                               cAGTCAATTTT‑GTAGGGGG‑ATC‑AGGCGTCCTCCTCCTTTTTTCTTT‑C‑CGCGCTGT‑CGTCTTTACATTTGCTCAGCCCTCCCCCCGCCCGCTCTCCCCGACCGACGGGGAGGGCCACACCCAGCCCGG‑CGAGCAGCAGCATAGAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACTGCTTTAGCACTGATTGGCGTACCGTCGATACCCTGATCGTTAACACCGAGGATGGCGTa                                                                                                                  <  1:621020/251‑1 (MQ=255)
                                                                                                                                      tgtt‑gtAAGGGG‑ATC‑ATGCGGGCTCCTCCGTTTTTCTTT‑C‑CGCGCTGT‑CGTCTTTACGATTGCTCAGCACTCCCCCCGCACGCTCTCCCCGACCGACGTGGAAGCCCAAAACCAGACCGG‑CGAGCAGCAGCAAAGAAGCCATTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACTGCTTTAGCACTGATTGGCGTACCGTCGATACCCTGATCGTTAACACCGAGGATGGCGTAAACAAt                                                                                                            <  1:153230/250‑1 (MQ=255)
                                                                                                                                        gttgaaggggg‑ATC‑ATGCGTGCTCCTCCGTTTTTCTTT‑T‑CGCGCTGT‑CGTTTTTACGTTTGCTCAGCACTTCCCCCGCCCGCTCTTCACGCCCGACGTGGAAGGCCACAACCAGCCCTG‑CGAGCAGCAGCATAGAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACTGCTTTAGCACTGATTGGCGTACCGTCGATACCCTGATCGTTAACACCGAGGATGGCGTAAACAATc                                                                                                           <  1:619125/243‑1 (MQ=255)
                                                                                                                                           gtAAGGGG‑ATC‑ATGCGTGCTCCTCCGTTTTTCTTT‑C‑CGCGCTGT‑CGTCTTTACGATTGCTCAGCCCTCCACCCGCACGCTCTTCACGCCCGACGTGGAAGGCCACAACCAGCCCGG‑CGAGCAGCAGCATAGAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACTGCTTTAGCACTGATTGGCGTACCGTCGATACCCTGATCGTTAACACCGAGGATGGCGTAAACAATcacc                                                                                                        <  1:173542/250‑1 (MQ=255)
                                                                                                                                              aGGGG‑ATC‑ATGCTTGCTCCTCCGTTTTTCTTT‑T‑CGCGCTGT‑CGTCTTTACGATTGCTCAGCACTCCCCCCGCACGCTCTTCCCGACCGACGTGGAAGGCCACAACCAGCCCTG‑CGAGCAGCAGCATAGAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACTGCTTTAGCACTGATTGGCGTACCGTCGATACCCTGATCGTTAACACCGAGGATGGCGTAAACAATCAccacca                                                                                                    <  1:581218/251‑1 (MQ=255)
                                                                                                                                              aGGGG‑ATC‑ATGCGTGCTCCTCCGTTTTTCTTT‑T‑CGCGCTGT‑CGTCTTTACGATTGCTCAGCACTTCACCCGCACGCTCTTCACGACCGACGTGGAAGGCCACAACCAGACCTG‑CGAGCAGCAGCATAGAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACTGCTTTAGCACTGATTGGCGTACCGTCGATACCCTGATCGTTAACACCGAGGATGGCGTAAACAATCaccacc                                                                                                     >  1:41916/1‑250 (MQ=255)
                                                                                                                                               gggg‑ATC‑ATGCGTGCTCCTCCGTTTTTCTTT‑T‑CGCGCTGT‑CGTCTTTACGATTGCTCAGCACTTCACCCGCACGCTCTTCACGACCGACGTGGAAGGCCACAACCAGACCTG‑CGAGCAGCAGCATAGAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACTGCTTTAGCACTGATTGGCGTACCGTCGATACCCTGATCGTTAACACCGAGGATGGCGTAAACAATCaccacc                                                                                                     >  2:311789/1‑249 (MQ=255)
                                                                                                                                                ggg‑ATC‑ATGCGTGCTCCTCCGTTTTTCTTT‑T‑CGCGCTGT‑CGTCTTTACGATTGCTCAGCACTTCACCCGCACGCTCTTCACGACCGACGTGGAAGGCCACAACCAGACCTG‑CGAGCAGCAGCATAGAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACTGCTTTAGCACTGATTGGCGTACCGTCGATACCCTGATCGTTAACACCGAGGATGGCGTAAACAATCACCACCAGc                                                                                                  >  1:70241/1‑251 (MQ=255)
                                                                                                                                                    aTC‑ATGCGTGCTCCTCCGTTTTTCTTT‑T‑CGCGCTGT‑CGTCTTTACGATTGCTCAGCACTTCACCCGCACGCTCTTCACGACCGACGTGGAAGGCCACAACCAGACCTG‑CGAGCAGCAGCATAGAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACTGCTTTAGCACTGATTGGCGTACCGTCGATACCCTGATCGTTAACACCGAGGATGGCGTAAACAATCACCACCAGCatg                                                                                               >  2:596926/1‑251 (MQ=255)
                                                                                                                                                          gCGTGCTCCTCATTTTTTCTTT‑C‑CGCGCTGT‑CGTCTTTCCGATGGCTCAGCCCTTCCCCCGCCCGCTCTCCACGCCCAAGGGGGAAGGCCACAACCAGCCCTG‑CGAGCAGCAGCAAAGAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACTGCTTTAGCACTGATTGGCGTACCGTCGATACCCTGATCGTTAACACCGAGGATGGCGTAAACAATCACCACCAGCatgatg                                                                                            <  1:85404/249‑1 (MQ=255)
                                                                                                                                                            gTGCTCCTCCGTTTTTCTTT‑T‑CGCGCTGT‑CGTCTTTACGATTGCTCAGCACTTCACCCGCACGCTCTTCACGACCGACGTGGAAGGCCACAACCAGACCTG‑CGAGCAGCAGCATAGAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACTGCTTTAGCACTGATTGGCGTACCGTCGATACCCTGATCGTTAACACCGAGGATGGCGTAAACAATCACCACCAGCATGATGGCCg                                                                                        >  2:616931/1‑251 (MQ=255)
                                                                                                                                                                      gTTTTTCTTT‑T‑CGCGCTGT‑CGTCTTTACGATTGCTCAGCACTTCACCCGCACGCTCTTCACGACCGACGTGGAAGGCCACAACCAGACCTG‑CGAGCAGCAGCATAGAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACTGCTTTAGCACTGATTGGCGTACCGTCGATACCCTGATCGTTAACACCGAGGATGGCGTAAACAATCACCACCAGCATGATGGCCGACAAAATTGc                                                                              >  2:507900/1‑251 (MQ=255)
                                                                                                                                                                            cttt‑t‑cGCGCTGT‑CGTCTTTACGATTGCTCAGCACTTCACCCGCACGCTCTTCACGACCGACGTGGAAGGCCACAACCAGACCTG‑CGAGCAGCAGCATAGACGCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACTGCTTTAGCACTGATTGGCGTACCGTCGATACCCTGATCGTTAACACCGAGGATGGCGTAAACAATCACCACCAGCATGATGGCCGACAAAATTGCCGGAcc                                                                        >  2:561275/1‑251 (MQ=255)
                                                                                                                                                                              tt‑t‑cGCGCTGT‑CTTCTTTGCGATGGCCCAGCACTTCCCCCGCACCCTCTTCCCGCCCGACGGGGAAGCCCCCAACCAGCCCTG‑CGAGCAGCAGCATAGAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACTGCTTTAGCACTGATTGGCGTACCGTCGATACCCTGATCGTTAACACCGAGGATGGCGTAAACAATCACCACCAGCATGATGGCCGACAAAATTGCCGGAcc                                                                        <  2:174003/249‑1 (MQ=255)
                                                                                                                                                                                ttccgcgCTGT‑CGTCTTTACGATTGCTCAGCACTCCCCCCGCCCGCTCTTCCCGCCCGACGTGAAGGGCCACACCCAGCCCTG‑CGAGCAGCAGCATAGAAGCCATTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACTGCTTTAGCACTGATTGGCGTACCGTCGATACCCTGATCGTTAACACCGAGGATGGCGTAAACAATCACCACCAGCATGATGGCCGACAAAATTGCCGGAcc                                                                        <  1:507900/246‑1 (MQ=255)
                                                                                                                                                                                               cTTTACGATTGCTCAGCACTTCACCCGCACGCTCTTCACGACCGACGTGGAAGGCCACAACCAGACCTG‑CGAGCAGCAGCATAGAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACTGCTTTAGCACTGATTGGCGTACCGTCGATACCCTGATCGTTAACACCGAGGATGGCGTAAACAATCACCACCAGCATGATGGCCGACAAAATTGCCGGACCAATCCACACCTGCGGt                                                        >  2:6346/1‑251 (MQ=255)
                                                                                                                                                                                                 ttACGATTGCTCAGCCCTTCCCCCGCACGCTCTTCACGACCGACGGGGAAGGCCACACCCAGCCCGG‑CGAGCAGCAGCATAGAAGCAGGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACTGCTTTAGCACTGATTGGCGTACCGTCGATACCCTGATCGTTAACACCGAGGATGGCGTAAACAATCACCACCAGCATGATGGCCGACAAAATTGCCGGACCAATCCACACCTGCGGtt                                                       <  1:592865/250‑1 (MQ=255)
                                                                                                                                                                                                 ttACGATTGCCCAGCACTTCCCCCGCCCGCTCTTCACGACCGACGTGGAAGGCCACACCCAGACCGG‑CGAGCAGCAGAATAGAAGCCAGTCCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACTGCTTTAGCACTGATTGGCGTACCGTCGATACCCTGATCGTTAACACCGAGGATGGCGTAAACAATCACCACCAGCATGATGGCCGACAAAATTGCCGGACCAATCCACACCTGCGGtt                                                       <  1:168860/250‑1 (MQ=255)
                                                                                                                                                                                                  tccgtttGCTCACCACTCCACCCGCACGCTCTCCACGACCGACGGGGAAGGCCACACCCAGCCCTG‑CGAGCAGCAGCATAGAAGCCATTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACTGCTTTGGCACTGATTGGCGTACCGTCGATACCCTGATCGTTAACACCGAGGATGGCGTAAACAATCACCACCAGCATGATGGCCGACAAAATTGCCGGACCAATCCACACCTGCGGttt                                                      <  1:363018/245‑1 (MQ=255)
                                                                                                                                                                                                  gaCGATTGCTCACCACTCCCCCCGCCCGCCCTTCCCGACCGACGGGGAAGGCCACACCAAGCCCAG‑GGGGCAGCAGCATAGAAGCCAGTTCCACCGCCAGTACGTAAGGCCCAAACAGCGTAATACCCACTGCTTTAGCACTGATTGGCGTACCGTCGATACCCTGATCGTTAACACCGAGGATGGCGTAAACAATCACCACCAGCATGATGGCCGACAAAATTGCCGGACCAATCCACACCTGCGGttt                                                      <  2:249113/249‑1 (MQ=255)
                                                                                                                                                                                                    cGATTGCTCAGCACTTCACCCGCACGCTCTTCACGACCGACGTGGAAGGCCACAACCAGACCTG‑CGAGCAGCAGCATAGAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACTGCTTTAGCACTGATTGGCGTACCGTCGATACCCTGATCGTTAACACCGAGGATGGCGTAAACAATCACCACCAGCATGATGGCCGACAAAATTGCCGGACCAATCCACACCTGCGGTTTc                                                     >  1:83037/1‑249 (MQ=255)
                                                                                                                                                                                                         gCTCAGCCCTCCCCCCGCACGCTCTCCACGACCGACGTGGAAGGCCACAACCAGACCCG‑CGAGCAGCAGCATAGAAGCCATTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACTGCTTTAGCACTGATTGGCGTACCGTCGATACCCTGATCGTTAACACCGAGGATGGCGTAAACAATcaccac                                                                                                      <  1:341790/195‑1 (MQ=255)
                                                                                                                                                                                                         gCTCAGCACTTCACCCGCACGCTCTTCACGACCGACGTGGAAGGCCACAACCAGACCTG‑CGAGCAGCAGCATAGAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACTGCTTTAGCACTGATTGGCGTACCGTCGATACCCTGATCGTTAACACCGAGGATGGCGTAAACAATcaccac                                                                                                      >  2:341790/1‑195 (MQ=255)
                                                                                                                                                                                                                   ccccccgcccgCTCTCCCCGACCGACGTGGAAGGCCACACCCAGCCCGG‑CGAGCAGCAGCATAGAGCCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACTGCTTTAGCACTGATTGGCGTACCGTCGATACCCTGATCGTTAACACCGAGGATGGCGTAAACAATCACCACCAGCATGATGGCCGACAAAATTGCCGGACCAATCCACACCTGCGGTTTCAGCCACTGGCGTTCCTg                                    <  1:357384/248‑1 (MQ=255)
                                                                                                                                                                                                                           ccgcccctCACGACCCGCGTGGAAGCCCACAACCAGACCTG‑CGAGCAGCAGCATAGAAGCCAATTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACTGCTTTAGCACTGATTGGCGTACCGTCGATACCCTGATCGTTAACACCGAGGATGGCGTAAACAATCACCACCAGCATGATGGCCGACAAAATTGCCGGACCAATCCACACCTGCGGTTTCAGCCACTGGCGTTCCTGTTCGAttt                            <  1:152323/244‑1 (MQ=255)
                                                                                                                                                                                                                           aCGCTCTTCACGACCGACGTGGAAGGCCACAACCAGACCTG‑CGAGCAGCAGCATAGAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACTGCTTTAGCACTGATTGGCGTACCGTCGATACCCTGATCGTTAACACCGAGGATGGCGTAAACAATCACCACCAGCATGATGGCCGACAAAATTGCCGGACCAATCCACACCTGCGGTTTCAGCCACTGGCGTTCCTGTTCGAttt                            >  1:90884/1‑251 (MQ=255)
                                                                                                                                                                                                                             gCTCTTCACGACCGACGTGGAAGGCCACAACCAGACCTG‑CGAGCAGCAGCATAGAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACTGCTTTAGCACTGATTGGCGTACCGTCGATACCCTGATCGTTAACACCGAGGATGGCGTAAACAATCACCACCAGCATGATGGCCGACAAAATTGCCGGACCAATCCACACCTGCGGTTTCAGCCACTGGCGTTCCTGTTCGATTTCt                          >  2:498795/1‑251 (MQ=255)
                                                                                                                                                                                                                              ctctTCACGACCGACGTGGAAGGCCACAACCAGACCTG‑CGAGCAGCAGCATAGAAGCCAGATCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACTGCTTTAGCACTGATTGGCGTACCGTCGATACCCTGATCGTTAACACCGAGGATGGCGTAAACAATCACCACCAGCATGATGGCCGACAAAATTGCCGGACCAATCCACACCTGCGGTTTCAGCCACTGGCGTTCCTGTTCGATTTCTg                         >  1:373232/1‑251 (MQ=255)
                                                                                                                                                                                                                                 ttCACGACCGACGTGGAAGGCCACAACCAGACCTG‑CGAGCAGCAGCATAGAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACTGCTTTAGCACTGATTGGCGTACCGTCGATACCCTGATCGTTAACACCGAGGATGGCGTAAACAATCACCACCAGCATGATGGCCGACAAAATTGCCGGACCAATCCACACCTGCGGTTTCAGCCACTGGCGTTCCTGTTCGATTTCTGAAc                      >  2:518301/1‑251 (MQ=255)
                                                                                                                                                                                                                                      gcccgacGTGGAAGCCCACACCCAGCCCTG‑CGAGCAGCACCATAGAACCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACTGCTTTAGCACTGATTGGCGTACCGTCGATACCCTGATCGTTAACACCGAGGATGGCGTAAACAATCACCACCAGCATGATGGCCGACAAAATTGCCGGACCAATCCACACCTGCGGTTTCAGCCACTGGCGTTCCTGTTCGATTTCTGAACCGccc                 <  1:611258/249‑1 (MQ=255)
                                                                                                                                                                                                                                      gaccgacGTGGCAGGCCACAACCAGACCTG‑CGAGCAGCAGCATAGAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACTGCTTTAGCACTGATTGGCGTACCGTCGATACCCTGATCGTTAACACCGAGGATGGCGTAAACAATCACCACCAGCATGATGGCCGACAAAATTGCCGGACCAATCCACACCTGCGGTTTCAGCCACTGGCGTTCCTGTTCGATTTCTGCACCGccc                 >  2:161513/1‑251 (MQ=255)
                                                                                                                                                                                                                                           ccgtgaacggccccacccAGCCCTG‑CGAGCAGCGGCATAGAAGCTATTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACTGCTTTAGCACTGATTGGCGTACCGTCGATACCCTGATCGTTAACACCGAGGATGGCGTAAACAATCACCACCAGCATGATGGCCGACAAAATTGCCGGACCAATCCACACCTGCGGTTTCAGCCACTGGCGTTCCTGTTCGATTTCTGAACCGCCCAGGtt            <  1:236157/235‑1 (MQ=255)
                                                                                                                                                                                                                                           aggtgaaaggcacaaaCCAGACCTG‑CGAGAAGAAACAAAGAGGCCAGTCCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACTGCTTTAGCACTGATTGGCGTACCGTCGATACCCTGATCGTTAACACCGAGGATGGCGTAAACAATCACCACCAGCATGATGGCCGACAAAATTGCCGGACCAATCCACACCTGCGGTTTCAGCCACTGGCGTTCCTGTTCGATTTCTGAACCGCCCAGGtt            <  1:39210/237‑1 (MQ=255)
                                                                                                                                                                                                                                                      cacaACCAGACCTG‑CGAGCAGCAGCATAGAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACTGCTTTAGCACTGATTGGCGTACCGTCGATACCCTGATCGTTAACACCGAGGATGGCGTAAACAATCACCACCAGCATGATGGCCGACAAAATTGCCGGACCAATCCACACCTGCGGTTTCAGCCACTGGCGTTCCTGTTCGATTTCTGAACCGCCCAGGTTGAGCatcatc  >  1:191459/1‑250 (MQ=255)
                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                 
CCGGCGACCACGAAGGCCAGCGCGGAGGCGTTAATCATGATTTCCAGACCAATCAACATAAACAGCAGATTGCGACGGATAACCA‑GACCGGTTAAGCCAA‑GAACGAATAAGATTGCCGCGAGGATCAGTCCATGTT‑GTAAGGGG‑ATC‑ATGCGTGCTCCTCCGTTTTTCTTT‑T‑CGCGCTGT‑CGTCTTTACGATTGCTCAGCACTTCACCCGCACGCTCTTCACGACCGACGTGGAAGGCCACAACCAGACCTG‑CGAGCAGCAGCATAGAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACTGCTTTAGCACTGATTGGCGTACCGTCGATACCCTGATCGTTAACACCGAGGATGGCGTAAACAATCACCACCAGCATGATGGCCGACAAAATTGCCGGACCAATCCACACCTGCGGTTTCAGCCACTGGCGTTCCTGTTCGATTTCTGAACCGCCCAGGTTGAGCATCATC  >  NC_000913/2395197‑2395684

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