Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,481,804	0	G	A	40.5%	109.7 / 35.5	75	P125S (CCA→TCA)	emrY	putative multidrug efflux system
Reads supporting (aligned to +/- strand):  ref base G (38/6);  new base A (0/30);  total (38/37)
Fisher's exact test for biased strand distribution p-value = 4.15e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATG‑CCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTG  >  NC_000913/2481657‑2481945                                                                                                                                                        |                                                                                                                                               tCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATG‑CCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAAt                                                                                                                                             >  1:482832/1‑151 (MQ=255)  ccAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATG‑CCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAAtt                                                                                                                                            <  1:471297/151‑1 (MQ=255)     aCCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATG‑CCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTAA‑c                                                                                                                                        >  2:62622/1‑151 (MQ=255)     aCCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATG‑CCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTAA‑c                                                                                                                                        <  2:255297/151‑1 (MQ=255)      ccccAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATG‑CCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTAA‑Cg                                                                                                                                       >  1:614509/1‑151 (MQ=255)        ccAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATG‑CCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTAA‑CGGc                                                                                                                                     <  1:784350/151‑1 (MQ=255)            cTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATG‑CCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTAA‑CGGcccc                                                                                                                                  >  2:709978/1‑150 (MQ=255)            cTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATG‑CCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTAA‑CGGCCCCg                                                                                                                                 >  2:418768/1‑151 (MQ=255)              aaaGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATG‑CCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTAA‑CGGCCCCGcc                                                                                                                               >  1:968187/1‑151 (MQ=255)              aaaGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATG‑CCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTAA‑CGGCCCCGcc                                                                                                                               >  1:899526/1‑151 (MQ=255)                  ttATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATG‑CCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTAA‑CGGCCCCGCCAtt                                                                                                                            >  2:498767/1‑150 (MQ=255)                  ttATAACAAATAAAACCGCCCAATATCGGCCCATATATGGGGGCAATAATCACGGTCATTGACCATAATG‑CCAGAGCATTTTTTTTTTTTTCGGGCGGATACTTCCTT‑AATAATAAA‑CTCTGTGAGAGT‑GGACTTAA‑CGGCCCCGCCAtt                                                                                                                            <  1:772869/150‑1 (MQ=255)                    aTCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATG‑CCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTAA‑CGGCCCCGCCATTAAc                                                                                                                         >  1:618535/1‑151 (MQ=255)                     tCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATG‑CCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTAA‑CGGCCCCGCCATTAAcc                                                                                                                        >  2:346049/1‑151 (MQ=255)                       acaAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATG‑CCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTAA‑CGGCCCCGCCATTAAcccc                                                                                                                      >  1:325194/1‑151 (MQ=255)                            tataACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATG‑CCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTAA‑CGGCCCCGCCATTAACCCCTGAAc                                                                                                                 >  2:266779/1‑151 (MQ=255)                               acccGCCAAATATCGCCCCACATATCGGAGCAATAACCACGGTCTTGAACCAAAATG‑CCAGAGCAAATTTTTTTTTTTCTGGCGAATATTTCTTT‑AATAAAAAC‑CTCTGTGACGGT‑GAAATAAA‑CGCCCCCGCCATTAACCCCTGAACGa                                                                                                               <  2:874197/148‑1 (MQ=255)                                aCCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATG‑CCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGAct                                                                                                             >  2:757335/1‑151 (MQ=255)                                      cAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATG‑CCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTaaa                                                                                                        >  1:195751/1‑150 (MQ=255)                                       aatgccgccccAAATAGGGGAGCGATAATCAGGGTAATTGACAATAAGC‑CAAAAGCAATTTTTTTTTTTGCGGGCGAATATTTCCTA‑AATAATAAC‑CTCGGTGAAAGG‑GAAAATAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTaaa                                                                                                        <  1:566403/141‑1 (MQ=255)                                          atCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATG‑CCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAat                                                                                                    >  2:258339/1‑150 (MQ=255)                                            cGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATG‑CCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAatat                                                                                                  >  2:107751/1‑150 (MQ=255)                                               cccACATATCGGAGCGATAATCACGGTCATTGACCATAATG‑CCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGc                                                                                              >  2:736121/1‑151 (MQ=255)                                                 cacaTATCGGAGCGATAATCACGGTCATTGACCATAATG‑CCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGcac                                                                                            >  1:722099/1‑151 (MQ=255)                                                      gtCGGAGCAATAACCACGGTCTTTGACAAAAAGG‑CAGGAGAAATTGTTTTTTTTTCGGGCGAATATTTCTTT‑AAAAAAAAA‑CTCTGAAAAAGA‑GAAATAAC‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCa                                                                                        <  1:219547/149‑1 (MQ=255)                                                         ggggtaaaaaTCACGTTATTTGCCCAAAATA‑CAAGAAAAATTTTTTTTTTTTCGGGCGTATATTTCTTA‑AAAAAAAAC‑CCCGGTGACAGT‑AAAATTAC‑CGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCaa                                                                                       <  1:125706/140‑1 (MQ=255)                                                         ggggcaaTAATCAGGTTAATTGCCCA‑AATGCCCAGAGAAATTTTTTTTTTTTCGGGCGGAT‑ATTTCTT‑AA‑AAAAAA‑CGCCCGGAAAATGGGAAATAAGGGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATc                                                                                         <  2:210046/140‑1 (MQ=255)                                                            gCGATAATCACGGTCATTGACCATAATG‑CCAGAGCAAATGTTCTTTTTTCTGGCGGGTAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTgg                                                                                 >  1:183945/1‑151 (MQ=255)                                                             caaTAATCCGGTTCATTGACCATAAGG‑CAAGAAAAAATTTTTTTTTTTCGGGGGAATAATTCTTA‑AAAAAAAAC‑CTCTGGAACGGT‑GAAATAAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTg                                                                                  <  2:159408/147‑1 (MQ=255)                                                              gATAATCACGGTCATTGACCATAATG‑CCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTc                                                                               >  2:487366/1‑151 (MQ=255)                                                              gATAATCACGGTCATTGACCATAATG‑CCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTc                                                                               >  1:525983/1‑151 (MQ=255)                                                               aTAATCACGGTCATTAACCATAAGG‑CAAAAGAAAATGTTTTTTTTTCTGGCGAATATTTCCTA‑AATAAAAAA‑CTCTGGGACAGG‑GAAATTAC‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCg                                                                              <  2:1056404/151‑1 (MQ=255)                                                               aTAATCACGGTAATTGACAATAATG‑CCAGAACAAATTTTTTTTTTTCGGGCGGATATTTCCTT‑AAAAATAAA‑CTCTGTGACAGT‑GAAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGt                                                                                <  1:484665/149‑1 (MQ=255)                                                               aTAACCACGGTCATTGACCATAATG‑CCAGAGAAATTGTTTTTTTTTCGGGGGAATAATTCCTA‑AAAAAAAAA‑CTCTGTGACAGT‑GAAATTAC‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCg                                                                              <  2:1089962/151‑1 (MQ=255)                                                                 aaTCACGTTTATTGACCATAATG‑CAAGAGAAAATGTTTTTTTTTCGGGCGAATAATCCTTA‑AAAAATAAC‑CTTTGGGACAGT‑GAAATAAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTg                                                                                  <  2:409011/145‑1 (MQ=255)                                                                 aaTCACGGTCATTGACCATAATG‑CCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTAA‑CGGCCCAGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCg                                                                              >  1:648823/1‑149 (MQ=255)                                                                      cGGTCATTGACCATAATG‑CCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCt                                                                       >  2:67784/1‑151 (MQ=255)                                                                      cGGTCATTGACCATAATG‑CCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCt                                                                       >  1:382342/1‑151 (MQ=255)                                                                       gatcatgaaCCATAATG‑CAAGAGAAATTTTTTTTTTTCCGGGCGAATAATCCTTA‑AATAAAAAA‑CTCTGAGCGGGT‑GTAATTAC‑GGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGc                                                                        <  1:529260/141‑1 (MQ=255)                                                                              taacaataagccaagaaCAAATTTTTTTTTTTCGGGCGAATAATCCCTT‑AAAAAAAAC‑CTCGTTGAAAGT‑GAAATTAC‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACAt                                                                 <  1:719223/137‑1 (MQ=255)                                                                               aCCATAATG‑CAAGAGAAATTTTTTTTTTTTCGGGCGGATATTTCTTA‑AAAAAAAAA‑CTTTGTGAAGGT‑GAAATTAA‑CGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACAt                                                                 <  1:886793/148‑1 (MQ=255)                                                                                  aTAATG‑CCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAAATTGGTCGATAGGCTACACATTAATg                                                            >  2:217073/1‑150 (MQ=255)                                                                                    aagg‑CCGGAGCTATTGTTTTTTTTTC‑GGCTGTAAAATCCTTCAATAAAAAA‑CTTGGCGAAGGT‑GAAATACG‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATg                                                            <  2:668880/145‑1 (MQ=255)                                                                                    aaTG‑CCAGAGAAATTTTTTTTTTTTCTGGGGGATAATTCCTT‑AATAATAAC‑CTCTGTGACAGT‑GAAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAAt                                                             <  2:761455/147‑1 (MQ=255)                                                                                    aaTG‑CAAGAGAATTTTTTTTTTTTTCGGGGGAATAATTCATA‑AATAAAAAA‑CTCTGTAAAATA‑AAAATTAG‑CGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGa                                                           <  1:448126/149‑1 (MQ=255)                                                                                      tG‑CCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGAc                                                       >  1:298694/1‑151 (MQ=255)                                                                                      tG‑CCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGAc                                                       >  2:729975/1‑151 (MQ=255)                                                                                         ccAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGAc                                                       >  2:948820/1‑149 (MQ=255)                                                                                             agaaatttttCTTTTTTCTGGCGAATAATTCCTA‑AAAAAAAAA‑CTCTGTGACAGT‑GAAATTAC‑CGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGaa                                                 <  2:510453/143‑1 (MQ=255)                                                                                             agCAATTTTTTTTTTTTCGGGGGGATATTTCTAT‑AAAAATAAT‑CTCTGAGACAGG‑GAAATTAG‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCg                                                   <  1:688908/149‑1 (MQ=255)                                                                                             agCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGa                                                  >  1:727187/1‑150 (MQ=255)                                                                                             aaaaattttttttttttCGGGCAGA‑AATTTCTA‑AAAAAAAAACCTCTGTAAAGGT‑AAAATAAC‑CGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACATCGaa                                                 <  1:785546/140‑1 (MQ=255)                                                                                                  aTGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGaaaaaa                                             >  2:930694/1‑150 (MQ=255)                                                                                                  aTGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGaaaaaa                                             >  1:190827/1‑150 (MQ=255)                                                                                                   tttttttttttCTGGCGAATATTTCCTA‑AAAAAAAAC‑CTCGGTGACAGT‑GAAATTAA‑CGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGaaaaaa                                             <  2:602875/144‑1 (MQ=255)                                                                                                   tttttttttttCTGGCGAATAATTCCTT‑AAAAATAAC‑CTGTGTGAAAGT‑AAAATAAA‑CGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGaa                                                 <  1:1061971/140‑1 (MQ=255)                                                                                                        ttttttCGGGGGAATATTTCCTA‑AAAAAAAAC‑CTCTGGGACAGT‑AAAATAAA‑CGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTg                                        <  2:685737/149‑1 (MQ=255)                                                                                                          ttttGTGGGAGATAATTCTTA‑AATAAAAAC‑CTCTGAAAAAGA‑GAAAATAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACt                                     <  1:308254/150‑1 (MQ=255)                                                                                                            ttCTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGa                                   >  2:664814/1‑150 (MQ=255)                                                                                                              cgggCGGAAATTTCTTT‑AAAAATAAC‑CTCTGGAAGAGG‑AAAATTAA‑CGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGaaa                                 <  1:626014/148‑1 (MQ=255)                                                                                                               tGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAg                                >  1:846781/1‑150 (MQ=255)                                                                                                                ggggaatatTTCCTT‑AAAAAAAAA‑CTCGGGGACAGT‑GAAATTAA‑CGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTa                              <  1:859256/146‑1 (MQ=255)                                                                                                                ggagaatatTTCTTA‑AAAAAAAAA‑CTCTGTGACAGT‑GAAATAAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGt                               <  2:582806/145‑1 (MQ=255)                                                                                                                   ggATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTaaa                            >  2:518202/1‑150 (MQ=255)                                                                                                                     atatTTCCTT‑AATAATAAA‑CTCTGTACAAGA‑GGAATAAA‑CCGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTaaaa                           <  1:801767/149‑1 (MQ=255)                                                                                                                     aTAATTCTTT‑AAAAAAAAC‑CTCGGTGAAAGT‑GAAATAAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTaaaaa                          <  2:1038658/150‑1 (MQ=255)                                                                                                                     aTAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAaat                         >  2:182510/1‑151 (MQ=255)                                                                                                                            cTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGGCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCtt                   >  2:187939/1‑150 (MQ=255)                                                                                                                                aaaaaaaaa‑CTTTGTGAAAGT‑GAAATAAC‑CGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTaa                 <  1:498684/143‑1 (MQ=255)                                                                                                                                    ataaC‑CTCTGTGAAAGT‑GAAATAAC‑CGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGcc           <  1:876476/151‑1 (MQ=255)                                                                                                                                      aaa‑CTCTGTGACAGT‑GGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCt          >  1:628296/1‑150 (MQ=255)                                                                                                                                        a‑CTCTGTGACAGT‑GGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTAtt       >  1:506314/1‑151 (MQ=255)                                                                                                                                              gtgACAGT‑GAAATAAC‑CGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTc      <  1:658619/147‑1 (MQ=255)                                                                                                                                              gagaaaGT‑GAAATTAC‑CGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTtg  <  1:390117/146‑1 (MQ=255)                                                                                                                                                 aCAGC‑GAAATAAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCGGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTtg  <  2:686364/148‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                               TCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATG‑CCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTG  >  NC_000913/2481657‑2481945

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