Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,228,2750TC44.8% 54.1 / 13.3 58pseudogene (335/1287 nt)mdtQpseudogene; putative channel/filament proteins
Reads supporting (aligned to +/- strand):  ref base T (25/7);  new base C (12/14);  total (37/21)
Fisher's exact test for biased strand distribution p-value = 1.51e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCG‑CG‑TTCCGCCGCCCGTGC‑TTTAACCGTACCCAGGCGGGCAG‑T‑AACCTCCG‑CCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCT  >  NC_000913/2228128‑2228422
                                                                                                                                                        |                                                                                                                                                    
tctTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTAC‑CGTTTCCCCCCCCCGTGT‑TTTAACCTTACCC‑CGCGGGGAG‑T‑AAACTc                                                                                                                                                     >  2:597123/1‑149 (MQ=255)
    tttCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCG‑CG‑TTCCGCCCCCCCCGC‑TTTTACCCCACCCACGCGGGCAG‑G‑AAAccccc‑ccc                                                                                                                                               >  2:685731/1‑150 (MQ=255)
    tttCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCG‑CG‑TTCCCCCCCCCCTGC‑TTTAAACGTACCCCAGCGGGCAG‑T‑AACCCCCG‑ccc                                                                                                                                               >  1:111230/1‑150 (MQ=255)
     ttCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCG‑CG‑TTCCGCCGCCCGTGC‑TTTAACCGTACCCAGGCGGGCAG‑T‑AACCTCCg                                                                                                                                                   <  2:334183/146‑1 (MQ=255)
       cTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTTGCGGGTTTTCTCC‑CG‑GTCCCCCGCCCGTGC‑TTTTACCCGACCCAGGCGGGCGG‑T‑AACCTccc‑ccc                                                                                                                                               >  2:299607/1‑147 (MQ=255)
        tATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCC‑CG‑GTTCCCCGCCCGTGC‑TTTTACCCGAACCCGGCGGGGAG‑G‑AAACCCCG‑ccccctt                                                                                                                                           >  1:330131/1‑146 (MQ=255)
           ttGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCG‑CG‑TTCCCCCGCCCCCGC‑TTTTACCGCACCCCGGCGGGGCG‑T‑AAACTCCC‑CCCCATTCttt                                                                                                                                       >  2:567680/1‑151 (MQ=255)
              ctgcAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTTCTCG‑CC‑TTTCCCCCCCCCCGC‑TTTTAACGCACCCCGGCGGGGAG‑G‑AAACCCCC‑CCCCATTCTTTtcc                                                                                                                                    >  2:122102/1‑151 (MQ=255)
              ctgcAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCACCCAGCAACTGGCGGGCTTGCCCGCCG‑TTCC‑CCCCCCCTTC‑TTTAACCCTACCCACACGGGGAGTT‑ATCCTCCC‑CCGGGTTCTTTcc                                                                                                                                     >  2:59692/1‑151 (MQ=255)
                    gACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCG‑CG‑GTCCGCCCCCCCGGC‑TTTAAACGTACCCCGGCGGGGAG‑T‑AAACCCCG‑CCCCGTTTTTTTCCCCGat                                                                                                                               >  1:620972/1‑150 (MQ=255)
                     aCCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCG‑CG‑TTCCCCCGCCCGTGC‑TTTAAACGTAACCACGCGGGCAG‑G‑TAACTCCG‑CCCGATTCTTTtccccaat                                                                                                                               >  1:43837/1‑145 (MQ=255)
                     aCCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCG‑CG‑TTTCGCCGCCCGGGC‑TTTAACCGTACCCAGGCGGGGAG‑G‑AAACCCCG‑CCCCCTTCTTTTCCCCGata                                                                                                                              >  1:121768/1‑150 (MQ=255)
                      ccGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCG‑CG‑TTCCGCCGCCCGTGC‑TTTAACCGTACCCAGGCGGGCAG‑T‑AACCTCCG‑CCCGATTCTTTCCCCAGatat                                                                                                                             <  2:32857/150‑1 (MQ=255)
                      ccGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCG‑CG‑GTCCCCCCCCCCTGC‑TTTTACCCTACCCCCGCGGGCCG‑T‑AAACCCCC‑CCCCCTTCTTTtccc                                                                                                                                   >  2:955780/1‑144 (MQ=255)
                       cGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCG‑CG‑TTTCCCCGCCCCCGC‑CTTTACCCGACCCCAGCGGGCCG‑TAAACCCCCC‑CCCCATCTTTTcccccgt                                                                                                                                >  2:668717/1‑144 (MQ=255)
                              aaCGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCG‑CG‑TTCCGCCGCCCGGGC‑TTTAAACGTACCCAGGCGGGCAG‑T‑AACCTCCC‑CCCGGTTTTTTCCCCCGATATCGAGATGcc                                                                                                                    >  1:1001181/1‑151 (MQ=255)
                               aCGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTTCTCC‑CG‑GTCCCCCCCCCCTGC‑TTTTACCCTACCCCAGCGGGGAG‑G‑AAACCCCGCCCCCATTTTTTCCCCA‑ATATCCAGATGcc                                                                                                                    >  2:962062/1‑150 (MQ=255)
                                cgccgcCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCG‑CG‑TTCCGCCGCCCGTGC‑TTTAACCGTACCCAGGCGGGCAG‑T‑AACCTCCG‑CCCGATTCTTTCCCCAGATATCGAGATGCCAg                                                                                                                  <  2:73533/151‑1 (MQ=255)
                                 gccgccTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCG‑CG‑TTCCGCCGCCCGTGC‑TTTAACCGTACCCAGGCGGGCAG‑T‑AACCTCCG‑CCCGATTCTTTCCCCAGATATCGAGATGCCAGc                                                                                                                 <  2:182746/151‑1 (MQ=255)
                                    gccTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCCCC‑CG‑TTTCGCCCCCCCTGC‑TTTTACCGTACCCAAGCGGGCAG‑T‑AAACTCCC‑CCCCCTTTTTTTCCCCGATATCGAGATGCCAGccc                                                                                                               >  1:268025/1‑150 (MQ=255)
                                            ttGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCG‑CG‑TTCCGCCCCCCGCGC‑TTTTAACGTACCCAGGCGGGCAG‑G‑AAACTCCC‑CCCCATTCTTTCCCCCGATATCGAGAAGCCAACccgcccgtaa                                                                                                       >  1:815432/1‑145 (MQ=255)
                                            ttGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCC‑CG‑TTCCCCCCCCCCCGCTTTTAACCG‑ACCCCCGCGGGCAG‑G‑AAACCCCC‑CCCCCTTCTTTTCCCCGATATCCAGAAGCCCGcccccccgt                                                                                                         >  2:192409/1‑145 (MQ=255)
                                             tGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTTCTCC‑CG‑TTTCGCCGCCCCCGC‑TTTTACCCTACCCCGGGGGGGCG‑G‑AACCCCCC‑CCCGGTTTTTTTCCCCCATATCGAGATGCCCGCCCCCCCTTaa                                                                                                       >  2:240944/1‑149 (MQ=255)
                                                  cTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTTCTCC‑CG‑GTCCGCCGCCCCTGC‑TTTAAACCGACCCCGGCGGGCCG‑T‑AAACCCCG‑CCCCGATCTTTTCCCCGATATCGAGATGCCAACCCCCCCTTTAGCCaa                                                                                                  >  1:823239/1‑149 (MQ=255)
                                                         tACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTGTCC‑CG‑GTCCGCCGGCCGTGG‑TTTTACCGGAACCCGGCGGGGAG‑G‑AAACCCCC‑CCCCAGTCTTTTCCCCGATAATCAGATACCCGCCCCCCGATTAACCaaaaagaac                                                                                           >  2:118309/1‑144 (MQ=255)
                                                            aGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCG‑CG‑TTCCGCCGCCCGTGC‑TTTAACCGTACCCCGGCGGGCAG‑T‑AACCTCCG‑CCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCCAAAGTACCgct                                                                                        >  2:379400/1‑147 (MQ=255)
                                                            aGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCG‑CG‑TTCCGCCGCCCGTGC‑TTTAACCGTACCCAGGCGGGCAG‑T‑AACCTCCG‑CCCGATTCTTTCCCCCGATATCGAGATGGCAGCCCGCCGGTTAGCCCAAAGTACCGTTgg                                                                                      >  2:513663/1‑151 (MQ=255)
                                                              gCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCG‑CG‑TTCCGCCGCCCGTGC‑TTTAACCGTACCCAGGCGGGCAG‑T‑AACCTCCG‑CCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGgtg                                                                                    >  1:719256/1‑151 (MQ=255)
                                                                 ggCTACGCTGCCAGCCAGCAATGGGGGGGTTTGCTCG‑CG‑TCCCCCCGCCCGTGC‑TTTACCCGTCCCCGGGGGGGAAG‑T‑ACCCCCCG‑CCCGTTTCTTCCCCAAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTAc                                                                                 <  2:849786/151‑1 (MQ=255)
                                                                       gCTGCCAGCCAGCAATTGGCGGGTTTGCTCG‑CG‑TTCCGCCGCCCGTGC‑TTTAACCGTCCCCGGGCGGGCAG‑T‑ACCCTCCG‑CCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACgg                                                                            <  1:561624/150‑1 (MQ=255)
                                                                            cagccagcAATTGGCGGGTTTGCTCG‑CG‑TTCCGCCGCCCGTGC‑TTTAACCGTACCCAGGCGGGCAG‑T‑AACCTCCG‑CCCGATTCTTTTCCCAGATATCGAGATGCCAGCCCGCCGTTAAACCAAGAGTACCGTTGGGGTACCACGGGCCgg                                                                       >  2:674957/1‑150 (MQ=255)
                                                                                 agcAATTGGCGGGGTTGCTCG‑CG‑TTCCGCCGCCCGTGC‑TTTAACCGTACCCAGGCGGGCAG‑T‑AACCTCCG‑CCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCg                                                                        >  1:546549/1‑144 (MQ=255)
                                                                                         gCGGGTTTGCTCG‑CG‑TTCCGCCGCCCGTGC‑TTTAACCGTACCCAGGCGGGCAG‑T‑AACCTCCG‑CCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGcc                                                         >  1:712419/1‑151 (MQ=255)
                                                                                           ggttttGCCGG‑CG‑TCCCGCCCCCGGTGT‑TTAACCCGTCCCCCGGCGGACAG‑A‑ACCCCCCC‑CCCAATTTTTTCCCCAGTTATCGAGACGCCGCCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCgg                                                           <  2:18026/144‑1 (MQ=255)
                                                                                                       cg‑TTCCGCCGCCCGTGC‑TTTAACCGTACCCAGGCGGGCAG‑T‑AACCTCCG‑CCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCag                                              >  2:775745/1‑149 (MQ=255)
                                                                                                        g‑TTCCGCCGCCCGTGC‑TTTAACCGTACCCAGGCGGGCAG‑T‑AACCTCCG‑CCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGc                                           >  1:729540/1‑151 (MQ=255)
                                                                                                          ttCCGCCGCCCGTGC‑TTAACCCGTCCCCAGGCGGGAAG‑T‑ACCCCCCC‑CCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCag                                              <  2:79412/147‑1 (MQ=255)
                                                                                                            ccgccgCCCGTGC‑TTTACCCGTACCCAGGCGGGCAG‑T‑ACCCCCCG‑CCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCa                                          <  1:514250/149‑1 (MQ=255)
                                                                                                               ccgccCGTGC‑TTTAACCGTCCCCAGGCGGGCAG‑T‑ACCCCCCG‑CCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACgg                                     <  1:269819/151‑1 (MQ=255)
                                                                                                               ccgccCATGC‑TTTAACCGTACCCAGGCGGGCAG‑T‑AACCTCCG‑CCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACgg                                     >  1:403388/1‑151 (MQ=255)
                                                                                                                cgccCGTGC‑TTTAACCGTACCCAGGCGGGCAG‑T‑AACCTCCG‑CCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGc                                    >  1:89913/1‑151 (MQ=255)
                                                                                                                 gccCGTGC‑TTTAACCGTACCCAGGCGGGCAG‑T‑AACCTCCG‑CCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGcc                                   >  1:333951/1‑151 (MQ=255)
                                                                                                                 gccCGTGC‑TTTAACCGTACCCAGGCGGGCAG‑T‑AACCTCCG‑CCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGcc                                   >  1:435750/1‑151 (MQ=255)
                                                                                                                      tGC‑TTTACCCGTCCCCAGGGGGGCAG‑T‑ACCCCCCC‑CCCGATTCTTCCCCCAGATAGCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCa                                <  1:675589/149‑1 (MQ=255)
                                                                                                                          tttACCCGTACCCTGCCGGGATT‑T‑ACCCCCCG‑CCCATTTTTTTCCCCAGAAATCGGGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTaa                            <  1:598211/150‑1 (MQ=255)
                                                                                                                          tttAACCGTCCCCGGGGGGGCAG‑T‑ACCCCCCG‑CCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTaa                            <  1:571270/150‑1 (MQ=255)
                                                                                                                                    cccAGGGGGGCAG‑T‑CCCCCCCG‑CCCTTTTTTTCCCCCAGATATCGGGATGCAAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCg                 <  2:43834/151‑1 (MQ=255)
                                                                                                                                    cccAGGCGGGCAG‑A‑ACCCTCCG‑CCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGc                   <  1:63107/149‑1 (MQ=255)
                                                                                                                                        ggCGGGCAG‑T‑AACCTCCG‑CCCGATTCTTTGCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACAt             >  2:461709/1‑151 (MQ=255)
                                                                                                                                          cGGGCAG‑T‑AACCTCCG‑CCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACAttt           >  1:448796/1‑151 (MQ=255)
                                                                                                                                          cGGGCAG‑T‑AACCTCCG‑CCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACAttt           >  1:710564/1‑151 (MQ=255)
                                                                                                                                           gggCAT‑T‑ACCCCCCG‑CCCGTTTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACAttt           <  1:140221/150‑1 (MQ=255)
                                                                                                                                           gggCAG‑T‑ACCCTCCG‑CCCATTTCTTCCCCCATATACCGAGAGGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACAttt           <  2:542296/150‑1 (MQ=255)
                                                                                                                                           gggCAG‑T‑ACCCCCCG‑CCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACAttt           <  1:562385/150‑1 (MQ=255)
                                                                                                                                           gggCAG‑T‑ACCCCCCG‑CCCGATTCTTCCCCCAGATAGCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACAtttt          <  1:597183/151‑1 (MQ=255)
                                                                                                                                           gggAAG‑T‑ACCCCCCC‑CCCGATTCTTCCCCCAGTTATCGAGAGCCCGGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACAttt           <  1:272505/150‑1 (MQ=255)
                                                                                                                                           ggccAG‑C‑AACCCCCG‑CCCGTTTTTTCCCCCAGATTTCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACAtttt          <  2:671910/148‑1 (MQ=255)
                                                                                                                                                  t‑AACCTCCG‑CCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCg    >  1:242642/1‑151 (MQ=255)
                                                                                                                                                     aCCTCCG‑CCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCt  >  1:403062/1‑151 (MQ=255)
                                                                                                                                                        |                                                                                                                                                    
TCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCG‑CG‑TTCCGCCGCCCGTGC‑TTTAACCGTACCCAGGCGGGCAG‑T‑AACCTCCG‑CCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCT  >  NC_000913/2228128‑2228422

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