Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,228,2750TC44.4% 66.5 / 22.1 73pseudogene (335/1287 nt)mdtQpseudogene; putative channel/filament proteins
Reads supporting (aligned to +/- strand):  ref base T (26/14);  new base C (9/23);  total (35/38)
Fisher's exact test for biased strand distribution p-value = 2.28e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCT  >  NC_000913/2228135‑2228422
                                                                                                                                                |                                                                                                                                                         
cTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTAAGCTGCCAGCCAGCAATTGGCCGGGTTGCTCGCGGTCCGGCCCCCGTGCTTTTACCCAGCCCAGGGGGGCCCTTAACCCCGCCCGCtt                                                                                                                                                 >  1:198247/1‑151 (MQ=255)
       ctgcAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCCGCCAGCAATTGGCGGGTTTTCTCCCGTTTCGCCCCCCCTGCTTTTACCCTACCCAGGCGGGCAGTAAACTCCGCCCGGTTCTTTcc                                                                                                                                          >  2:533869/1‑150 (MQ=255)
        tgcAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTTCGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTcccc                                                                                                                                        <  2:1445565/151‑1 (MQ=255)
              aCCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCAATTCTTTCCCCAGatat                                                                                                                                 <  1:872942/151‑1 (MQ=255)
              aCCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGCCTTCCGCCCCCCCTGCTTTTACCCTACCCCGGCCGGCAGTAACCCCCCCCCCCTTCTTTccccccctat                                                                                                                                  >  1:465954/1‑145 (MQ=255)
                  tgtTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCGGCCAGCAATTGGCGCGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCgag                                                                                                                            <  1:1154652/151‑1 (MQ=255)
                   gtTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATGGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGGAACCTCCGCCCGATTCTGTCCCCAGATCTCgag                                                                                                                            <  1:493921/150‑1 (MQ=255)
                         cgccgcCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAg                                                                                                                   <  2:1113580/151‑1 (MQ=255)
                                tGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCccgccg                                                                                                            <  2:674864/151‑1 (MQ=255)
                                      tGCCACTCCCAGTACAGGCGGGCTACGCCGCCAGCCAGCAATTGGCGGGCTTGCTCGCGTTTCCCCGCCCCGGCGTTTACCCGCCCCAGGCGCGGGGTGAACCCCCCCCCAATATTTTCCCCCATCTCGGGATGCCCGCCCCCCCttttg                                                                                                      >  1:113248/1‑147 (MQ=255)
                                           cTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTTGTCCCGCCCCGCCGCCCCGGCTTTTAACCGACCCCCGCGGGCAGTAACCTCCCCCCGGTTCTTTCCCCCGATATCCAGATGCCCAGCCCCCGTTTAACCCaaa                                                                                                 >  1:1246657/1‑151 (MQ=255)
                                           cTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTTCTTCCGGTCCGCCGCCCCTGCTTTTAACCGACCCAGGCGGGGCGTTAACTCCGCCCGGTTTTTTTCCCAGAAATCTCGAATGCCGCCCGCCCTTTAGCCacaa                                                                                                >  1:47660/1‑151 (MQ=255)
                                           cTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGGTCCGCCCCCCCCGCTTTAACCCCACCCAGGCGGGCGGGAACCTTCGCCCGCTTCTTTTCCCAGATATCGAGATTCCCGCCCCCCGTTTAGCCcaaa                                                                                                >  2:255501/1‑151 (MQ=255)
                                            tCCCAGTACAGGCGGGCTACGCTGCCAGCCCGCAATTGGCGGGGTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCCTACCCCGGCGGGCAGGAAACTCCCCCCGATTTTTTTCCCCGATATTGAGAAGCCCACCCCCCCTTAAACCaaaaa                                                                                               >  2:1448071/1‑150 (MQ=255)
                                            tCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTTCTCCCGGTTCGCCCCCCCCGCTTTTACCCGACCCCGGCGGGCAGTAACCTCCCCCCCGCTCTTTTCCCCGGTAACCAGGTGCCCACCCCCCGCTTAGCCcaaaa                                                                                               >  1:183807/1‑150 (MQ=255)
                                               cAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTTTTCCCGGTCCGCCGCCCCGGCTTTTACCGTACCCCGGGGGGCCGTAACCTCCCCCCCATTCTTTTCCCCGCAATCCGGATGCCACCCCCCCCCTTAGCCCAAAATa                                                                                             >  1:1003136/1‑150 (MQ=255)
                                               cAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTTCTCGCGGTTCGCCGCCCGGGCTTTAAACGTACCCCGGCGGGCAGTAAACTCCGCCCCATTTTTTTCCCCGATATTGGGATGCCAGCCCCCCGGTTAGCCCAAAATa                                                                                             >  2:1197065/1‑150 (MQ=255)
                                                 gTACAGGCGGGCTACGCTGCCCGCCACCAATTGGCGGGGTTTCTGGCGTTCCGCCGCCCGGGCGTTTTACGGACCCACGGGGGCGGTAATCTTCCCCCCGTTTCTTCCCCCGATATTCAGACGCCCGCCCCCCGTCTAGCCcaaaagaac                                                                                           >  1:1070517/1‑145 (MQ=255)
                                                 gTACAGGCGGGCTACGCTGCCAGCCCGCAATTGGCGGGGTTTCTCCCGGTTCCCCGCCCCCTCCTTTACCCCCCCCCGGCGGGCCGTAAACCCCCCCCCATTTCTTTCCCCCGTATCGCAGAGGCCACCCCCCGCTAAGGCAAAAGTAcc                                                                                           >  1:1497517/1‑150 (MQ=255)
                                                 gTACAGGCGGGCTACGCTGCCAGCCCGCAATTGGCGGGGTTGCTCGCGTTTCGCCCCCCCGGCTTTTAACCCACCCCAGCGGGGCGTTACCCCCCCCCCGCTCTTTTCCCCGCAATCGCGAAGGCACCCCCCCCGTTAGCCAAAAATAcc                                                                                           >  1:1155585/1‑150 (MQ=255)
                                                 gTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGGGGGGTTTCTCGGGGTTCGCCCCCCCGGCTTTTAACCGACCCCAGGGGGCAGTAAACCCCCCCCCATTCTTTTTCCCGATATTGAGGTGGCCGCCCCGCGTTAAGCCcaaaacaac                                                                                           >  1:644027/1‑145 (MQ=255)
                                                 gTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCCCGGTCCGCCCCCCGCGGTTTTAACCTACCCCGGGGGGGAGTAACCCCCGCCCCCGTTTTTTCCCCCATATCCAGAAGCCAGCCCCCCGCTTAAACCaaaacaac                                                                                            >  1:1474457/1‑146 (MQ=255)
                                                 gTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCCCGGTTCCCCCCCCGTGGTTTTACCCTAACCCGGCGGGCAGGAAACCCCCCCCCCTTCTTTTCCCCGATATCGGGATGCCCACCCCCCCGTTAAGCCAAAAGAAc                                                                                            >  2:1276368/1‑150 (MQ=255)
                                                  tccAGGCGGGTTACGCGGCCGGCAAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTTTTTTAACCGTACCCAGGCGGGAAGTAACCCCCGCCCGATTCTTTCCCCAGTTATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCg                                                                                          <  2:1013378/148‑1 (MQ=255)
                                                  tACAGGCGGGCTTCGCGCCCAGCCAGAAATTGGGGGGTTTGCTCGCGTTCCCCCGCCCGTGTTTTACCCGTACCCAGGGGGGAATTACCCCCCCCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGt                                                                                         <  2:1460401/151‑1 (MQ=255)
                                                   acggggggccTACCCTGCCAGCCAGAAATGGGGGGGTTTGCCCGCGTCCCGCCGCCCGTGTTTTAACCGTCCCCGGGGGGGCAGTACCCCCCGCCCGATTTTTCCCCCAAATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGtt                                                                                        <  1:58804/142‑1 (MQ=255)
                                                   aCAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCCCGTTCCGCCCCCCGCGCTTTAACCGTACCCACGCGGGCAGTAACCCCCGCCCGATTCTTTTCCCCGATATCCAGATGCCAGCCCGCCGTTTAGCCAAAAGTACCGtt                                                                                        >  1:653883/1‑151 (MQ=255)
                                                      ggCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCCGATATCGAGATGCCCGCCCGCCGTTAAGCCAAAAGGACCCTTggg                                                                                     >  2:102557/1‑150 (MQ=255)
                                                         gggCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGGGCTTTAACCCTACCCAGGCGGGCAGGAACCTCCGCCCGATTCTTTCCCCAAATACTGACATGCCAGCCCGCCGGTCAGGCAAAAGGACCCTTGGTGTa                                                                                  >  1:1125898/1‑151 (MQ=255)
                                                           gCTCCGCTGCCAGCCAGCAATTGGGGGTTTTGCTCGCGTTCCGCCGCCCGTGCTTTACCCGTCCCCGGGCGGGAAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTAc                                                                                 <  2:783980/150‑1 (MQ=255)
                                                           gCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTACCCGTCCCCAGGCGGGCAGTAACCCCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTAc                                                                                 <  2:591062/150‑1 (MQ=255)
                                                           gCTACGCTCCCAGCCAGCAATGGGGGGGTTTGCGCGCGTTCCGCCGCCCGTGCTTTAACCGTCCCCAGGCGGGCAGTAACCCCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTAc                                                                                 <  2:537195/150‑1 (MQ=255)
                                                           ccgacccTGCAAGCGAGCATTGGCCGGTTTGCCTCCCTTCCCCCCGCCCTTGTTTCACCCGTCCCCACGCGGAGGTTACCCTCGCCCCATTCTTTTCCCAAGATATCGATTTACCAGGCCGCCGTTAAGCCAAAAGTACCGTTGGTGTAc                                                                                 <  1:1227714/144‑1 (MQ=255)
                                                            ctccGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCTTGTTTTAACCGTCCCCAGGCGGGCAGTACCCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCAc                                                                              <  2:664242/147‑1 (MQ=255)
                                                            ctccGCTCCCATCCAGCAATTGGCGGTTTTCCTCGTGTCCCTCCGCCCGCGCTTTACCTTTCCCCGGGCGGGCAGTATCCCCCGCCCTTTTCTTTCCCGGTATATCGAGATGCCACCCCGCCGTTAAGCCAAAAGTACGGTTGGTGTACCa                                                                               <  1:794935/148‑1 (MQ=255)
                                                            cgagcgtGCCAGCCGGGAGGTGGCGGGTGTGTTGGCGTCCCCCCCCCCTTTCTTTAACCGTACCCGGGCGGGCAGTACCCCCCCCCCGATCCTTCCCCAAGGTACGGAGATGCAGGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCa                                                                               <  1:1129244/145‑1 (MQ=255)
                                                             tgcgcggcaggccagaAATGGGGGGGTTTGCTCCCGTTCGGCCGCACGTCCTTTAACCGTCCCCCGGCGGGAAGACACCCCCGCCCAATTCTTCCCCCAGATAGCGAGAGGCCACACCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCAc                                                                              <  2:502889/141‑1 (MQ=255)
                                                             tccGCTGCCACCCAGCATTTGGCGGGTTTGCTCGCGTCCCGCCGCCGGTGCTTTACCCGTCCCCGGGCGGGCAGTACCCTCCGCCCGATTTTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCAc                                                                              <  2:472356/149‑1 (MQ=255)
                                                             tACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTCCCCAGGGGGGCAGTAACCTCCGCCCGATTTTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCAc                                                                              <  2:1165835/151‑1 (MQ=255)
                                                             tACGCGCCCCGCCAGAAATTGGGGGTTTTCCTCGCTTTCCGCCGCCTGTCTTTTAACCGTACCCCGGCGGGCAGTACCCCCCGCCCGTTTCTTTCCCCAGATTTAGAGATCCCAGACCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCAc                                                                              <  2:239604/151‑1 (MQ=255)
                                                                 ctcccacccagCAGTGGGCGGTTTTGCTCGCGTCCCCCCGCCCGTCCTTTACCCTTCCCCGGGCGGGCAGTAACCCCCCCCCGATTCTTCCCCCAAATACCGAGATGCCGGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGc                                                                          <  1:843098/148‑1 (MQ=255)
                                                                  tGCCAGCCAGAAATGGGGGGGTTTGCTCGCTTCCCCCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTACCCCCCGCCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGcc                                                                         <  1:770861/151‑1 (MQ=255)
                                                                   gccggccaGCATTGGGCGGGTTTGCTCGCGTCCCGCCCCCCGTGTTTTACCCGTACCCAGGGGGGAAGTAACCCCCGCCCGATTCTTCCCCCAGATATCGAGAGGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGcc                                                                         <  1:1093945/150‑1 (MQ=255)
                                                                     cacccgaaaaTGGGGGGGTTTCCCCGCGTCCCCCCCCCCGTGTTTTACCCTTCCCCAGGCGGGCAGTACCCCCCGCCCTATTCTTCCCCCAGATACCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGt                                                                      <  2:1120942/143‑1 (MQ=255)
                                                                        ccagcATTGGGCGGGTTGCCTCGCGTCCCCCCCCCCGTGTTTAACCCGACCCCAGGCGGGCAGTACCCTCCGCCCGATTCTTCCCCAAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGt                                                                      <  1:1196494/148‑1 (MQ=255)
                                                                            aaaTGGGGGGGTTGGCTCGCGTCCCGCCGCCGGTGCTTAAACCGTCCCGAGGGGGGAAGAACCCTCCGCCCGATTCTTTCCCAAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCt                                                               <  1:228846/150‑1 (MQ=255)
                                                                                       tttGCTCGCGTTCCGCCGCCCGGTCTTTCCCCGTCCCCGGGGGGGCAGTACCCCCCGCCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGAt                                                     <  1:44211/150‑1 (MQ=255)
                                                                                            cccgcgTTCCCCCGCCCTTCCTTAACCCTTACCCGGGCGGGCAGTACCCCCCGCCCTTTTCTTTCCCCAGATATCGAGATCCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTc                                                <  2:1024617/149‑1 (MQ=255)
                                                                                            cccccgTCCCCCCGCCCGTTTTTTACCGGGCCCGGGGGGGGCACAACCCCCCCCCCGATTCTTCCCCCAAATTTCGAGATCCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTc                                                <  2:124764/147‑1 (MQ=255)
                                                                                            cTCGCGTTCCGCCGCCCGTGCTTTAACCGTCCCCAGGCGGGAAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTGAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTc                                                <  2:861096/151‑1 (MQ=255)
                                                                                             caggtacgaccggccCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCa                                               >  1:589976/13‑151 (MQ=255)
                                                                                                  tcccgcccCCCGTGTTTTACCCGTCCCCGGGGGGGCAGTACCCCCCCCCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGc                                           <  2:475679/148‑1 (MQ=255)
                                                                                                      cgccccccggtgtttACCGGTCCCCATGCGGGCAGTAGCCTCCGCCCTATTCTTCCCCCAGATATCGAGATGCCGGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAAc                                       <  1:1285159/139‑1 (MQ=255)
                                                                                                        ggccCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACgg                                     >  2:906055/2‑150 (MQ=255)
                                                                                                              tGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCAtt                              >  2:1202810/1‑151 (MQ=255)
                                                                                                               gCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTa                             >  2:510496/1‑151 (MQ=255)
                                                                                                               gCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTa                             >  2:631610/1‑151 (MQ=255)
                                                                                                               gCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTa                             >  1:556349/1‑151 (MQ=255)
                                                                                                                cTTAAACCGTCCCGAGGGGGGCGATACCCTCCCCCCTATTTTTCCCCCAGATAACGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTa                             <  1:116138/150‑1 (MQ=255)
                                                                                                                 ctctAACCGTACCCAGGCGGGCAGTATCCTCCGCCCGATTCTTTCCCCAGGCATCGAGCTGCCAGCCCGCCGTTAAGCCAAAAGCACCGTCGGCGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTCCAGACCacaccacccacttaa                            >  1:1100178/4‑140 (MQ=255)
                                                                                                                    tacccGTCCCCAGGCGGGCAGTACCCCCCCCCCGATTTTTTCCCCAAATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGTCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGcc                         <  2:1071957/148‑1 (MQ=255)
                                                                                                                     acccGTCCCCGGGGGGGCAGTACCCCCCCCCCGATTTTTCCCCAAGATACCGAGAGGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCt                        <  2:1378182/149‑1 (MQ=255)
                                                                                                                     aaCCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGc                          >  1:86701/1‑149 (MQ=255)
                                                                                                                      aCCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCtt                       >  2:42705/1‑151 (MQ=255)
                                                                                                                      aCCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCtt                       >  1:413162/1‑151 (MQ=255)
                                                                                                                        cGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCtt                       >  2:118248/1‑149 (MQ=255)
                                                                                                                         gTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCtt                       >  2:207546/1‑148 (MQ=255)
                                                                                                                          tccccAGGCGGGCAGTACCCCCCGCCCGTTTTTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGc                   <  2:954472/149‑1 (MQ=255)
                                                                                                                            ccgggggggcgaGTACCCCCCGCCCGTTTCTTCCCCAAAAAAGGGAGAGGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCg                 <  2:933349/140‑1 (MQ=255)
                                                                                                                            cccAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCg                 >  2:1357721/1‑151 (MQ=255)
                                                                                                                              cAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGa                >  2:808942/1‑150 (MQ=255)
                                                                                                                                        atAACCCCCGCCCCATTTCTTCCCCAGATAATCGGGGGGCAAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTGCCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGc      <  1:644387/150‑1 (MQ=255)
                                                                                                                                              ccTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACCTTTTTTGCCGCt  >  1:638003/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                         
CTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCT  >  NC_000913/2228135‑2228422

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