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

TCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCT  >  NC_000913/2228128‑2228422
                                                                                                                                                        |                                                                                                                                                    
tctTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTACCGTTTCCCCCCCCCGTGTTTTAACCTTACCCCGCGGGGAGTAAACTc                                                                                                                                                     >  2:597123/1‑149 (MQ=255)
    tttCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGCGTTCCGCCCCCCCCGCTTTTACCCCACCCACGCGGGCAGGAAAcccccccc                                                                                                                                               >  2:685731/1‑150 (MQ=255)
    tttCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCGCGTTCCCCCCCCCCTGCTTTAAACGTACCCCAGCGGGCAGTAACCCCCGccc                                                                                                                                               >  1:111230/1‑150 (MQ=255)
     ttCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCg                                                                                                                                                   <  2:334183/146‑1 (MQ=255)
       cTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTTGCGGGTTTTCTCCCGGTCCCCCGCCCGTGCTTTTACCCGACCCAGGCGGGCGGTAACCTcccccc                                                                                                                                               >  2:299607/1‑147 (MQ=255)
        tATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCCCGGTTCCCCGCCCGTGCTTTTACCCGAACCCGGCGGGGAGGAAACCCCGccccctt                                                                                                                                           >  1:330131/1‑146 (MQ=255)
           ttGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGCGTTCCCCCGCCCCCGCTTTTACCGCACCCCGGCGGGGCGTAAACTCCCCCCCATTCttt                                                                                                                                       >  2:567680/1‑151 (MQ=255)
              ctgcAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTTCTCGCCTTTCCCCCCCCCCGCTTTTAACGCACCCCGGCGGGGAGGAAACCCCCCCCCATTCTTTtcc                                                                                                                                    >  2:122102/1‑151 (MQ=255)
              ctgcAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCACCCAGCAACTGGCGGGCTTGCCCGCCGTTCCCCCCCCCTTCTTTAACCCTACCCACACGGGGAGTTATCCTCCCCCGGGTTCTTTcc                                                                                                                                     >  2:59692/1‑151 (MQ=255)
                    gACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCGCGGTCCGCCCCCCCGGCTTTAAACGTACCCCGGCGGGGAGTAAACCCCGCCCCGTTTTTTTCCCCGat                                                                                                                               >  1:620972/1‑150 (MQ=255)
                     aCCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCCCCGCCCGTGCTTTAAACGTAACCACGCGGGCAGGTAACTCCGCCCGATTCTTTtccccaat                                                                                                                               >  1:43837/1‑145 (MQ=255)
                     aCCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGCGTTTCGCCGCCCGGGCTTTAACCGTACCCAGGCGGGGAGGAAACCCCGCCCCCTTCTTTTCCCCGata                                                                                                                              >  1:121768/1‑150 (MQ=255)
                      ccGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGatat                                                                                                                             <  2:32857/150‑1 (MQ=255)
                      ccGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGCGGTCCCCCCCCCCTGCTTTTACCCTACCCCCGCGGGCCGTAAACCCCCCCCCCTTCTTTtccc                                                                                                                                   >  2:955780/1‑144 (MQ=255)
                       cGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCGCGTTTCCCCGCCCCCGCCTTTACCCGACCCCAGCGGGCCGTAAACCCCCCCCCCATCTTTTcccccgt                                                                                                                                >  2:668717/1‑144 (MQ=255)
                              aaCGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCGCGTTCCGCCGCCCGGGCTTTAAACGTACCCAGGCGGGCAGTAACCTCCCCCCGGTTTTTTCCCCCGATATCGAGATGcc                                                                                                                    >  1:1001181/1‑151 (MQ=255)
                               aCGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTTCTCCCGGTCCCCCCCCCCTGCTTTTACCCTACCCCAGCGGGGAGGAAACCCCGCCCCCATTTTTTCCCCAATATCCAGATGcc                                                                                                                    >  2:962062/1‑150 (MQ=255)
                                cgccgcCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAg                                                                                                                  <  2:73533/151‑1 (MQ=255)
                                 gccgccTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGc                                                                                                                 <  2:182746/151‑1 (MQ=255)
                                    gccTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCCCCCGTTTCGCCCCCCCTGCTTTTACCGTACCCAAGCGGGCAGTAAACTCCCCCCCCTTTTTTTCCCCGATATCGAGATGCCAGccc                                                                                                               >  1:268025/1‑150 (MQ=255)
                                            ttGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCCCCCGCGCTTTTAACGTACCCAGGCGGGCAGGAAACTCCCCCCCATTCTTTCCCCCGATATCGAGAAGCCAACccgcccgtaa                                                                                                       >  1:815432/1‑145 (MQ=255)
                                            ttGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCCCGTTCCCCCCCCCCCGCTTTTAACCGACCCCCGCGGGCAGGAAACCCCCCCCCCTTCTTTTCCCCGATATCCAGAAGCCCGcccccccgt                                                                                                         >  2:192409/1‑145 (MQ=255)
                                             tGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTTCTCCCGTTTCGCCGCCCCCGCTTTTACCCTACCCCGGGGGGGCGGAACCCCCCCCCGGTTTTTTTCCCCCATATCGAGATGCCCGCCCCCCCTTaa                                                                                                       >  2:240944/1‑149 (MQ=255)
                                                  cTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTTCTCCCGGTCCGCCGCCCCTGCTTTAAACCGACCCCGGCGGGCCGTAAACCCCGCCCCGATCTTTTCCCCGATATCGAGATGCCAACCCCCCCTTTAGCCaa                                                                                                  >  1:823239/1‑149 (MQ=255)
                                                         tACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTGTCCCGGTCCGCCGGCCGTGGTTTTACCGGAACCCGGCGGGGAGGAAACCCCCCCCCAGTCTTTTCCCCGATAATCAGATACCCGCCCCCCGATTAACCaaaaagaac                                                                                           >  2:118309/1‑144 (MQ=255)
                                                            aGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCCGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCCAAAGTACCgct                                                                                        >  2:379400/1‑147 (MQ=255)
                                                            aGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCCGATATCGAGATGGCAGCCCGCCGGTTAGCCCAAAGTACCGTTgg                                                                                      >  2:513663/1‑151 (MQ=255)
                                                              gCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGgtg                                                                                    >  1:719256/1‑151 (MQ=255)
                                                                 ggCTACGCTGCCAGCCAGCAATGGGGGGGTTTGCTCGCGTCCCCCCGCCCGTGCTTTACCCGTCCCCGGGGGGGAAGTACCCCCCGCCCGTTTCTTCCCCAAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTAc                                                                                 <  2:849786/151‑1 (MQ=255)
                                                                       gCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTCCCCGGGCGGGCAGTACCCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACgg                                                                            <  1:561624/150‑1 (MQ=255)
                                                                            cagccagcAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTTCCCAGATATCGAGATGCCAGCCCGCCGTTAAACCAAGAGTACCGTTGGGGTACCACGGGCCgg                                                                       >  2:674957/1‑150 (MQ=255)
                                                                                 agcAATTGGCGGGGTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCg                                                                        >  1:546549/1‑144 (MQ=255)
                                                                                         gCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGcc                                                         >  1:712419/1‑151 (MQ=255)
                                                                                           ggttttGCCGGCGTCCCGCCCCCGGTGTTTAACCCGTCCCCCGGCGGACAGAACCCCCCCCCCAATTTTTTCCCCAGTTATCGAGACGCCGCCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCgg                                                           <  2:18026/144‑1 (MQ=255)
                                                                                                       cgTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCag                                              >  2:775745/1‑149 (MQ=255)
                                                                                                        gTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGc                                           >  1:729540/1‑151 (MQ=255)
                                                                                                          ttCCGCCGCCCGTGCTTAACCCGTCCCCAGGCGGGAAGTACCCCCCCCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCag                                              <  2:79412/147‑1 (MQ=255)
                                                                                                            ccgccgCCCGTGCTTTACCCGTACCCAGGCGGGCAGTACCCCCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCa                                          <  1:514250/149‑1 (MQ=255)
                                                                                                               ccgccCGTGCTTTAACCGTCCCCAGGCGGGCAGTACCCCCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACgg                                     <  1:269819/151‑1 (MQ=255)
                                                                                                               ccgccCATGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACgg                                     >  1:403388/1‑151 (MQ=255)
                                                                                                                cgccCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGc                                    >  1:89913/1‑151 (MQ=255)
                                                                                                                 gccCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGcc                                   >  1:333951/1‑151 (MQ=255)
                                                                                                                 gccCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGcc                                   >  1:435750/1‑151 (MQ=255)
                                                                                                                      tGCTTTACCCGTCCCCAGGGGGGCAGTACCCCCCCCCCGATTCTTCCCCCAGATAGCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCa                                <  1:675589/149‑1 (MQ=255)
                                                                                                                          tttACCCGTACCCTGCCGGGATTTACCCCCCGCCCATTTTTTTCCCCAGAAATCGGGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTaa                            <  1:598211/150‑1 (MQ=255)
                                                                                                                          tttAACCGTCCCCGGGGGGGCAGTACCCCCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTaa                            <  1:571270/150‑1 (MQ=255)
                                                                                                                                    cccAGGGGGGCAGTCCCCCCCGCCCTTTTTTTCCCCCAGATATCGGGATGCAAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCg                 <  2:43834/151‑1 (MQ=255)
                                                                                                                                    cccAGGCGGGCAGAACCCTCCGCCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGc                   <  1:63107/149‑1 (MQ=255)
                                                                                                                                        ggCGGGCAGTAACCTCCGCCCGATTCTTTGCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACAt             >  2:461709/1‑151 (MQ=255)
                                                                                                                                          cGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACAttt           >  1:448796/1‑151 (MQ=255)
                                                                                                                                          cGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACAttt           >  1:710564/1‑151 (MQ=255)
                                                                                                                                           gggCATTACCCCCCGCCCGTTTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACAttt           <  1:140221/150‑1 (MQ=255)
                                                                                                                                           gggCAGTACCCTCCGCCCATTTCTTCCCCCATATACCGAGAGGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACAttt           <  2:542296/150‑1 (MQ=255)
                                                                                                                                           gggCAGTACCCCCCGCCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACAttt           <  1:562385/150‑1 (MQ=255)
                                                                                                                                           gggCAGTACCCCCCGCCCGATTCTTCCCCCAGATAGCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACAtttt          <  1:597183/151‑1 (MQ=255)
                                                                                                                                           gggAAGTACCCCCCCCCCGATTCTTCCCCCAGTTATCGAGAGCCCGGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACAttt           <  1:272505/150‑1 (MQ=255)
                                                                                                                                           ggccAGCAACCCCCGCCCGTTTTTTCCCCCAGATTTCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACAtttt          <  2:671910/148‑1 (MQ=255)
                                                                                                                                                  tAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCg    >  1:242642/1‑151 (MQ=255)
                                                                                                                                                     aCCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCt  >  1:403062/1‑151 (MQ=255)
                                                                                                                                                        |                                                                                                                                                    
TCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCT  >  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: