Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913987,4940TG47.3% 163.9 / 26.2 93intergenic (‑512/+91)ompF/asnSouter membrane porin 1a (Ia;b;F)/asparaginyl tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base T (1/48);  new base G (44/0);  total (45/48)
Fisher's exact test for biased strand distribution p-value = 6.29e-26
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACG  >  NC_000913/987346‑987636
                                                                                                                                                    |                                                                                                                                               
tcttTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGtt                                                                                                                                                <  1:1269596/150‑1 (MQ=255)
tcttTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGt                                                                                                                                                 <  3:885824/149‑1 (MQ=255)
    tCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCaa                                                                                                                                            <  1:1148497/150‑1 (MQ=255)
      aGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGaaa                                                                                                                                           >  3:1281789/1‑149 (MQ=255)
      aGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGaa                                                                                                                                            >  2:147281/1‑148 (MQ=255)
         cATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGAt                                                                                                                                      <  3:1120606/151‑1 (MQ=255)
         cATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGAt                                                                                                                                      <  2:370390/151‑1 (MQ=255)
          aTCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGaaaaaaa                                                                                                                                       >  1:388166/1‑145 (MQ=255)
             cAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCaaaaaaaatttt                                                                                                                                  >  4:13050/1‑142 (MQ=255)
              aGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCtt                                                                                                                                 <  1:74917/151‑1 (MQ=255)
                 aTGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGAc                                                                                                                              <  3:710124/151‑1 (MQ=255)
                   gCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACtt                                                                                                                            <  4:517302/151‑1 (MQ=255)
                   gCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACtt                                                                                                                            <  1:1067014/151‑1 (MQ=255)
                   gCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACtt                                                                                                                            <  1:1018370/151‑1 (MQ=255)
                   gCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACtt                                                                                                                            <  3:382255/151‑1 (MQ=255)
                    ccAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGTAAAAAAATTGTTTGTCtt                                                                                                                            >  2:527564/1‑150 (MQ=255)
                    ccAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAATAGaatttttttcttt                                                                                                                           >  1:249543/1‑139 (MQ=255)
                    ccAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAAAAGAATTCtttgctt                                                                                                                            >  3:676305/1‑145 (MQ=255)
                     cAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTaa                                                                                                                          <  3:259984/151‑1 (MQ=255)
                      aaaaaaTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGTGGTGAAAATAGATTACttctctt                                                                                                                            >  3:1056625/1‑143 (MQ=255)
                       aaaaaTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAAAAAAtttgttgttttaa                                                                                                                          >  1:1352885/1‑138 (MQ=255)
                         aaaTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTaacaaa                                                                                                                      <  3:1357608/151‑1 (MQ=255)
                         aaaTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGGGGAAAAAAATTTATTGAATTtacaa                                                                                                                       >  4:146872/1‑150 (MQ=255)
                          aaTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTTCAAATAAATTTCTTTGCTTtacaaa                                                                                                                      >  3:231436/1‑150 (MQ=255)
                          aaTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCGAATAAAATTCTTTACTTtaaaaa                                                                                                                      >  1:823177/1‑150 (MQ=255)
                             aGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAg                                                                                                                   <  2:100153/150‑1 (MQ=255)
                             aGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAATAAATTGCTTGACTTaaaaaa                                                                                                                      >  3:1238242/1‑147 (MQ=255)
                              gTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGAt                                                                                                                 <  1:4535/151‑1 (MQ=255)
                              gTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAATAGAATGGttttatttacaaaaa                                                                                                                    >  2:342350/1‑135 (MQ=255)
                               tATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAATAAATTGCTTGGCTTaacaaaaa                                                                                                                    >  3:741615/1‑145 (MQ=255)
                               tATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAAAAAATTTATTTACTTAACAAACAGa                                                                                                                  >  4:243965/1‑149 (MQ=255)
                                  gAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAATAAAATGATTTAATTaaaaaaaaaaaaaa                                                                                                              >  1:910632/1‑142 (MQ=255)
                                    aaaTTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTg                                                                                                           <  3:1220191/151‑1 (MQ=255)
                                       ttGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGGGCGAATAAATTGTTTTACTTAAAGAACAAATaa                                                                                                               >  2:27166/1‑144 (MQ=255)
                                        tgtgCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACa                                                                                                        <  1:1336784/150‑1 (MQ=255)
                                        tgtgCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAAAAAATTTCTTTACTTTAAAAACAAATAAATTTAc                                                                                                         >  4:27101/1‑149 (MQ=255)
                                          tgCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAAAAAATTGGTTTACTTaaaaaaaacaaaaatttaaa                                                                                                        >  1:1015167/1‑134 (MQ=255)
                                           gCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGaa                                                                                                     <  3:423544/150‑1 (MQ=255)
                                           gCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAAAAAAATTCTTTTCTTAAAAAAAAAAAAAATTTACAGaa                                                                                                     >  3:708519/1‑150 (MQ=255)
                                                tAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGAt                                                                                                <  3:404074/150‑1 (MQ=255)
                                                     aaGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATaaaa                                                                                            <  4:206090/149‑1 (MQ=255)
                                                          aaaaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGc                                                                                     <  3:192866/151‑1 (MQ=255)
                                                             aCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAg                                                                                  <  3:152064/151‑1 (MQ=255)
                                                             aCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAg                                                                                  <  1:1023013/151‑1 (MQ=255)
                                                             aCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGCAATAAAATGGTTTACTTTAAAAAAAGATAAATTTAAAGaaaaaaaaaaaaaaaaaac                                                                                    >  3:334864/1‑144 (MQ=255)
                                                             aCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAAAAAAATGCTTGTCTTTAAAAAAAAAAAACTTGACAAAAAAGAaaaaaaaaaa                                                                                       >  1:1038993/1‑144 (MQ=255)
                                                             aCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAAAAGATTGCTTTTCTTTAAAAACAAAAAAATTTAAAGAAAAAATaaaaaaaaa                                                                                       >  3:1082963/1‑144 (MQ=255)
                                                              cTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGTAAAAAAATTGCTTTTCTTTAAAAAAAAAAAAATTTTCAGaaaaaaaaaaaaaaaaacc                                                                                    >  4:148223/1‑143 (MQ=255)
                                                                aaTCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCgg                                                                               <  1:281137/151‑1 (MQ=255)
                                                                 aTCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCggg                                                                              <  1:604005/151‑1 (MQ=255)
                                                                 aTCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCggg                                                                              <  1:1337559/151‑1 (MQ=255)
                                                                    cGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAAAAAATAGGTTGACTTAACAAAAAAATAACTTTACAGaaaaaaaaaaaaaaaaaaaaacggg                                                                              >  1:837371/1‑137 (MQ=255)
                                                                     gCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAAAAAAATGGTTTGCTTTACAAAAAAAAAAAATTACAGAAAAAAAAAAAAAAAAAGCAGCGGgcg                                                                            >  2:372972/1‑149 (MQ=255)
                                                                      cAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGGAAAAAAATTCTTTACTTTACAAAAAAAAAAAATTACAGAAAAAAAAAAAAAAAAAACAACGGGGGCTg                                                                         >  4:235815/1‑151 (MQ=255)
                                                                        aTAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAAAAAAATTTTTTACTTTACAAAAAAAAAACATTATAGAAAAAAAAAAAAAAAAACCaaaggggg                                                                            >  4:490244/1‑139 (MQ=255)
                                                                         tAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAATAAATTGCTTTTCTTTACAAAAAAATAAATTTACAGAAAAAAAAAAAAAAAAACCAACggggg                                                                            >  3:1284438/1‑143 (MQ=255)
                                                                          aaGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGctc                                                                     <  1:1311705/151‑1 (MQ=255)
                                                                          aaGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGTAAAAAAAATTTTTTACTTAAAAAAAAAAAAAATTGATAGAAAAAAAAAAAAAAAAACCAGCGGGGGCCgggt                                                                      >  2:32701/1‑148 (MQ=255)
                                                                            gAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAAAAAAATGTTTTACTTTACAAAAAAATAAATTGACAGAAAAGATAAAAAAAAAAACAGCggggggcgggtcct                                                                   >  1:59471/1‑140 (MQ=255)
                                                                             aaTGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTg                                                                 <  1:734843/151‑1 (MQ=255)
                                                                             aaTGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTg                                                                 <  1:88857/151‑1 (MQ=255)
                                                                               tGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAAAAGAATGCTTGTCTTAAAAAAAAAAAAAAATTAAAGAAAAAAAAAAAAAAAAGACAGCGGGCGCGGGCGCtt                                                                  >  1:716056/1‑148 (MQ=255)
                                                                               tGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAAAAAAATGCTTTACTTAAAAAAAAAAAAAATTTACAGAAAAAAAAAAAAAAAAAACAGCGGGCGCCGGGTCtt                                                                  >  2:343405/1‑148 (MQ=255)
                                                                                       aGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTa                                                        <  3:40410/150‑1 (MQ=255)
                                                                                        gTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAATAGATTTCTTGTCTTTACAAACAGATAAATTTACAGAAAAAATAAAAAAAAAGACCGCGGGCGCTGGCTCTTTAAGAAGCTaa                                                       >  2:125491/1‑150 (MQ=255)
                                                                                          ttttGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGGAATAAATTGTTTGACTTTACAAAAAAAAAAATTTACAGAAAAAATAAAAAAAAAGACAGCGGGGGCGGGCTCTTGAAGAAGCCaaa                                                      >  2:288011/1‑149 (MQ=255)
                                                                                            ttGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAAAAAAATTGTTGACATAAAAAAACAAAAAATTTTCAGAAAAAAAAAAAAAAAAACAAAAGGGCGCGGGCTCCTTAAGAAGCCaaa                                                      >  4:158399/1‑147 (MQ=255)
                                                                                             tGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAg                                                  <  3:1337688/150‑1 (MQ=255)
                                                                                                gcgTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTGGGAAATAAATTTTTTTACTTTAAAAAAAAATAAATTTACAGAAAAAAAAAAAAAAAAAACAGCGGGGGCTGGCTCTTTTAAAAGCtaaaataaaa                                                >  1:381494/1‑149 (MQ=255)
                                                                                                 cgTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCt                                             <  1:10167/151‑1 (MQ=255)
                                                                                                 cgTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCt                                             <  3:422269/151‑1 (MQ=255)
                                                                                                 cgTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGTAAAAAAAATTCTTGAATTAAAAAACAAAAAAATTGACAGAAAAAATAAAAAAAAAGACAGCGGGCGCCGGCGCTTTTAGAAGCTAAATTAGaa                                                >  1:1128646/1‑148 (MQ=255)
                                                                                                        gTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGtt                                       <  3:1362591/150‑1 (MQ=255)
                                                                                                        gTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTa                                      <  3:70174/151‑1 (MQ=255)
                                                                                                        gTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTCGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTa                                      <  2:146874/151‑1 (MQ=255)
                                                                                                            cAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTAcgcg                                  <  3:548207/151‑1 (MQ=255)
                                                                                                               tGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAAAAAATTTTTTTTATTTAAAAAAAGATAAATTGACAAAAAAAAAAAAAAAAAAAACAGCGGGCGCCGGGTCTTGAAGAAGCGAAAATAAAAAATGGGGTTa                                      >  1:301568/1‑144 (MQ=255)
                                                                                                                gTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGt                              <  4:340562/151‑1 (MQ=255)
                                                                                                                gTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGt                              <  3:736685/151‑1 (MQ=255)
                                                                                                                     atTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTg                         <  3:82244/151‑1 (MQ=255)
                                                                                                                     atTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAAAAAAATTTTTTACTTTAAAAAAAAAAAACTTTACAGAAAAAAAAAAAAAAAAACAAGCGGGGGCCGGGGCCTTTAGACGCCAAAATAAAAAATGGCGGTACGCGGAGGACggg                         >  1:1176120/1‑149 (MQ=255)
                                                                                                                     atTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAAAAAATTTTTTTGCTTAACAAAAAGATAAAATTACAGAAAAGATAAAAAAAAAAACAGCGGGGGGTGGGCCTTTTAAAAGCGAAAATAAAAACAGGCGGTACGCGGGGGACGTg                         >  1:187763/1‑151 (MQ=255)
                                                                                                                       tCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTgg                        <  3:10062/150‑1 (MQ=255)
                                                                                                                        ccGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTggg                       <  1:1272033/150‑1 (MQ=255)
                                                                                                                         cGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGaa                     <  1:83342/151‑1 (MQ=255)
                                                                                                                          gTAACAGGATGATCGTTATGGCGGGGGTGGGAAAAAATTTTTTTAATTATAAAAAAGATAAATTTTCAGAAAAAAAAAAAAAAAAAGCAGCGGGGGCCGGCTCTTTAAGACGCCAAAATAAAAAATGGGGGTACgggggggaagggggg                      >  4:377799/1‑138 (MQ=255)
                                                                                                                            aaCAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACg                   <  3:716159/150‑1 (MQ=255)
                                                                                                                              cAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGaa                <  3:736770/151‑1 (MQ=255)
                                                                                                                              cAGGATGATCGTTATGGCGGGGGTGGCAAAAAATTTCTTTGCTTTACAAACAAATAAATTTACAAAAAAAAAAAAAAAAAAAGCAGCGGGGGCTGGGTCTTGTAGAAGCTAAAATTAAAACTGGCGGTACGCCGAGTAAGTggggaagg                  >  3:1369920/1‑144 (MQ=255)
                                                                                                                              cAGGATGATCGTTATGGCGGGGGTGCAAATATATTTTTTTATTTTACAAAAAGAAAACTTGTCAAAAAAGAAAAAAAAAAAGGCAACGGGGGCCGGGTCTTTGAAAAGCCAAAATTAAAACAGGCGGTACGCGGAggaaggggggaagg                  >  1:393638/1‑136 (MQ=255)
                                                                                                                                 gatgatCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatca             <  3:551709/151‑1 (MQ=255)
                                                                                                                                     atCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatcacatc         <  3:64623/151‑1 (MQ=255)
                                                                                                                                        gTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCacg      <  2:91428/151‑1 (MQ=255)
                                                                                                                                             ggCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCacgtacg  <  1:953567/150‑1 (MQ=255)
                                                                                                                                             ggCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCacgtacg  <  3:249562/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
TCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACG  >  NC_000913/987346‑987636

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