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

ATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTT‑CTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTA‑ACG‑GCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGC  >  NC_000913/2481656‑2481948                                                                                                                                                         |                                                                                                                                                    atCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTT‑CTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGa                                                                                                                                                    >  1:1047131/1‑150 (MQ=255)    cAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTT‑CTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTa                                                                                                                                                <  1:700648/151‑1 (MQ=255)     aaCCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTT‑CTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTa                                                                                                                                                >  1:953581/1‑150 (MQ=255)      aCCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTT‑CTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTA‑Ac                                                                                                                                             >  1:648850/1‑151 (MQ=255)        cccAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGATCTATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTT‑CTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTTACAGT‑GGAATTA‑ACg‑g                                                                                                                                          <  2:292871/151‑1 (MQ=255)        cccAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTT‑CTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTA‑ACg‑g                                                                                                                                          >  1:174644/1‑151 (MQ=255)        cccAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTA‑CTGGCGGATAATTCCTT‑AATAATAAA‑CTGTGTCACAAT‑GGAATTA‑ACg‑g                                                                                                                                          <  2:519889/151‑1 (MQ=255)          cAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTT‑CTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTA‑ACG‑Gcc                                                                                                                                        >  2:679366/1‑151 (MQ=255)             cTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTT‑CTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTA‑ACG‑GCCCCg                                                                                                                                     >  1:832343/1‑151 (MQ=255)             cTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTT‑CTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTA‑ACG‑GCCCCg                                                                                                                                     >  1:404181/1‑151 (MQ=255)              tAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTT‑CTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTA‑ACG‑GCCCCGc                                                                                                                                    >  1:397636/1‑151 (MQ=255)                   ttATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTT‑CTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTA‑ACG‑GCCACGCCAtt                                                                                                                                >  2:177394/1‑150 (MQ=255)                          aaaTATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTT‑CTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTA‑ACG‑GCCCCGCCATTAACCCCTg                                                                                                                        >  1:1287206/1‑151 (MQ=255)                              aaacccGCCCAATATCGCCCCAAATATGGGAGCGATAATCGCGCTCTTACACTATAATGCCAGAGCAAATGTTCTTTTCT‑TGGGCGAATATTCCTTT‑AAAAAAAAC‑CTCTGGTAAAGT‑GAAATAA‑ACG‑GCCCCGCCATTAACCCCTGAACg                                                                                                                    <  2:568643/147‑1 (MQ=255)                               taACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTT‑CTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTA‑ACG‑GCCCCGCCATTAACCCCTGAACg                                                                                                                    >  1:568403/1‑150 (MQ=255)                               taACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTT‑CTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTA‑ACG‑GCCCCGCCATTAACCCCTGAACg                                                                                                                    >  2:618727/1‑150 (MQ=255)                                    gCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTT‑CTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTA‑ACG‑GCCCCGCCATTAACCCCTGAACGActct                                                                                                               >  1:210320/1‑150 (MQ=255)                                      ccAATACCGGCCCACATACCGGAGCGATAATCACGGTCATGGACAATAATGCCAGAGAAATTTTTTTTTTTT‑CGGGCGAATATTTCTTA‑AAAAAAAAC‑CTCTGTGACAGG‑GAAATAA‑ACG‑GCCCCGCCATTAACCCCTGAACGACTCTaaa                                                                                                            <  1:319005/151‑1 (MQ=255)                                             cGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTT‑CTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTA‑ACG‑GCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATc                                                                                                     >  1:396677/1‑151 (MQ=255)                                               gCCCAAATATCGGAGCAAAAACCACGGTCATTAACAAAAATGCCAGAGAAAAGGTTTTTTTTT‑CTGGCGGATAATTCTTT‑AAAAATAAA‑CCCTGTGAAAAT‑GATATCA‑GCT‑CCGCCGCCATTAACCCCGGAACGACTCTAAAAAATATc                                                                                                     <  2:455450/149‑1 (MQ=255)                                                cccAAATATCGGGGTGAAAATCACGGTCTTTACCCATAAGGCCAGA‑CAAATTTTTTTTTTTCCTGGAGAATATTCCCTT‑AAAAAAAAC‑CTCTTTGAAATT‑AAAATTA‑GCG‑GCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATc                                                                                                     <  2:416153/148‑1 (MQ=255)                                                 caacaattgggAGCGATAATCAGGAATTTTTACTATATTGCCAGAGAAAAATTTCTTTTTA‑CTGGGGGATAATTCCTT‑AAAAATAAA‑CTGTGTGA‑AGTAGAAATTA‑GGG‑GGGCCGCCATTAACCCCTGAACGACTCTAAAAAATATCa                                                                                                    <  2:1101362/139‑1 (MQ=255)                                                        tCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTT‑CTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTA‑ACG‑GCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCa                                                                                            >  1:973103/1‑149 (MQ=255)                                                           gAGCGATAATCTCGGTCATTGACCAAAATGCCAAAGAAATTGTTCTTTTTT‑CGGGCGAAAAATTCCTA‑AAAAATAAA‑CTCGGTAACAGT‑GAAATTA‑ACG‑GCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAg                                                                                                   <  2:830722/139‑1 (MQ=255)                                                           gAGCGATAATCACGTTCATTGACCATAATGCCAGAGCAAATTTTTTTTTTT‑CTGGCGGATAATTCCTT‑AAAAAAAAC‑CTCTGTGACAGT‑GGAATTA‑ACG‑GCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGAtt                                                                                       <  1:122209/151‑1 (MQ=255)                                                                aTAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTT‑CTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTA‑ACG‑GCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTc                                                                                   >  2:185915/1‑150 (MQ=255)                                                                aTAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTT‑CTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTA‑ACG‑GCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCg                                                                                  >  1:1283642/1‑151 (MQ=255)                                                                     cccGGTCATTGACAAAAATGCAAAAAAAATTGTTTTTTTTT‑CGGGCGGAAAATTCTTA‑AAAAAAAAA‑CCTTGGAAAGGT‑AAAATTA‑CCC‑GCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCg                                                                                  <  2:298959/144‑1 (MQ=255)                                                                        ggTCATTGACCATAATGCCAGAGAAATTTTTTTTTTTT‑CTGGCGGATAATTCCTA‑AATAAAAAC‑CTCTGTGACAGG‑GAAATTA‑ACG‑GCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAg                                                                              <  1:1275607/147‑1 (MQ=255)                                                                        gcttaTTGACCATAATCCCAGACAAAATGTTTTTTTTT‑CGGGCAGATAATCCCTT‑AAAAAAAAA‑CGGTTACAAATG‑GAAATTA‑GCA‑CCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAgg                                                                             <  2:865664/144‑1 (MQ=255)                                                                            aTTGACCATAATACAAGAGCAAATTTTTTTTTTT‑CTGGCGGATAATTCCTA‑AAAAAAAAA‑CTCTGTAACAGT‑GAAATTA‑ACG‑GCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCt                                                                           <  2:358696/146‑1 (MQ=255)                                                                                   aTAATGCCAGAGCAAATGTTCTTTTTT‑CTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTA‑ACG‑GCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATg                                                                >  1:687911/1‑150 (MQ=255)                                                                                   aTAATGCCAGAGAAAATGTTTTTTTTT‑CTGGCGGATAATTCCTA‑AAAAATAAC‑CTCGTTAAAAGA‑GAAATAA‑ACG‑GCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATg                                                                <  2:168988/150‑1 (MQ=255)                                                                                    tAATGCCAGAGCAAATGTTCTTTTTT‑CTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTA‑ACG‑GCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGa                                                               >  2:526247/1‑150 (MQ=255)                                                                                     aatgaaagagaatttGTTTTTTTTT‑CGGGGGAATAATCCCTA‑AAAAAAAAA‑CCTGGAAAAAGT‑GAAATAA‑ACG‑CCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCGACACATTAATGa                                                               <  2:1172020/135‑1 (MQ=255)                                                                                     aaggCCAGAGAAATTTTTTTTTTTT‑CGGGCGAATAATTCTTA‑AAAAAAAAC‑CTCTGTGACAGG‑GAAATTA‑ACG‑CCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAg                                                             <  1:904846/148‑1 (MQ=255)                                                                                     aaTGCAAGAGCAAATTTTTTTTTTT‑CTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GAAATTA‑ACG‑CCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGa                                                               <  2:1052046/149‑1 (MQ=255)                                                                                      aTGCCAGAGCAAATGTTCTTTTTT‑CTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTA‑ACG‑GCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGa                                                               >  1:457773/1‑148 (MQ=255)                                                                                          cAGAGCAAATGTTCTTTTTT‑CTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTA‑ACG‑GCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGc                                                        >  2:279357/1‑151 (MQ=255)                                                                                             agCAAATTTTCTTTTTT‑CTGGCGGATAATTCCTT‑AAAAAAAAA‑CTCTGTAAAAGT‑GAAATTA‑ACG‑CCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑a                                                     <  2:263412/150‑1 (MQ=255)                                                                                             agCAAATGTTCTTTTTT‑CTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTA‑ACG‑GCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑a                                                     >  1:942036/1‑150 (MQ=255)                                                                                              gcagttGTTTTTTTTG‑GGGGGGGATAATTCCTA‑AAAAAAAAC‑GGGTGTGAGGGG‑GGAATGA‑GGG‑GCGCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGGTACGCATGAATGAAGACAGc                                                        <  2:144611/142‑1 (MQ=255)                                                                                                  aTGTTCTTTTTT‑CTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTA‑ACG‑GCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑aaaaaa                                                >  1:646669/1‑150 (MQ=255)                                                                                                   ttttttttttg‑CGGGCAAATATTTCCTT‑AAAAATAAC‑CCCGGTGAAAGA‑GAAATAA‑GCG‑CCGCCGCCATTAACCCCTGAACTACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑aaaaaaa                                               <  1:1144494/145‑1 (MQ=255)                                                                                                     ttttttttt‑CTGGGGAATAATTCCTA‑AATAAAAAG‑CGCTG‑GAC‑GA‑GGAAAAATACGTCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAg                                             <  1:985555/147‑1 (MQ=255)                                                                                                     ttCTTTTTT‑CTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTA‑ACG‑GCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAGt                                            >  1:1366507/1‑151 (MQ=255)                                                                                                           ttt‑CTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTA‑ACG‑GCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑GAAAAAAAGTGACTg                                       >  2:942999/1‑150 (MQ=255)                                                                                                           ttt‑CTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTA‑ACG‑GCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAGTGACTg                                       >  1:716708/1‑150 (MQ=255)                                                                                                            tt‑CTGGGAGATAATTCTTA‑AAAAAAAAC‑CTTTGGAAAATT‑GAAATTA‑ACG‑CCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAGTGACTGa                                      <  2:797816/150‑1 (MQ=255)                                                                                                            tt‑CTGGCGTATATTTCCTT‑AATAATAAA‑CTCTGTAACAGT‑GAAATTA‑ACG‑GCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAGTGACTGa                                      <  1:750860/150‑1 (MQ=255)                                                                                                             t‑CTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTA‑ACG‑GCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAGTGACTGaa                                     >  2:882407/1‑150 (MQ=255)                                                                                                             t‑CTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTA‑ACG‑GCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAGTGACTGaa                                     >  1:641276/1‑150 (MQ=255)                                                                                                                 ggCGAATAATTCCTA‑AAAAAAAAC‑CTCTGGGACAGT‑GAAATTA‑CCG‑CCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAGTGACTGAAAg                                   <  1:679546/149‑1 (MQ=255)                                                                                                                  ggggAAAATTTCTTAAATAAAAAC‑CTTTGTAACAGT‑GAAATTA‑ACG‑GCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAAGTGACTGaaa                                    <  1:561328/147‑1 (MQ=255)                                                                                                                  gcgaaaaaTTCCTA‑AATAATAAC‑CTCTGTGACGGT‑GAAATTA‑CCG‑GCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAGTGACTGAAAGt                                  <  1:1307885/143‑1 (MQ=255)                                                                                                                   cGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTA‑ACG‑GCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAGTGACTGAAAGTaa                                >  1:152887/1‑150 (MQ=255)                                                                                                                   cGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTA‑ACG‑GCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAGTGACTGAAAGTaa                                >  1:1208819/1‑150 (MQ=255)                                                                                                                    gcaaaaTTCATT‑AATAAAAAA‑CCTTCTGACAAA‑GAAATCA‑CCG‑GCACCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACAACAAGAGTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAGTGACTGAAAGTa                                 <  2:628444/144‑1 (MQ=255)                                                                                                                     gATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTA‑ACG‑GCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAGTGACTGAAAGTaaaa                              >  2:1341792/1‑150 (MQ=255)                                                                                                                          ttCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTA‑ACG‑GCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAGTGACTGAAAGTAAAAATAATc                        >  1:641091/1‑151 (MQ=255)                                                                                                                          ttCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTA‑ACG‑GCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGCCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAGTGACTGAAAGTAAAAATAATc                        >  2:713058/1‑151 (MQ=255)                                                                                                                             cTT‑AAAAAAAAACCTCTGGGACAGG‑GAAATTA‑CCG‑GCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAGTGACTGAAAGTAAAaataat                         <  1:740757/148‑1 (MQ=255)                                                                                                                              tt‑AATAATAAA‑CTCTGTGACAGT‑GGAATTA‑ACG‑GCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAGTGACTGAAAGTAAAAATAATCTTa                     >  2:615342/1‑150 (MQ=255)                                                                                                                              tt‑AAAAAAAAC‑CTCTGTGACAGT‑GAAATAA‑ACG‑GCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAGTGACTGAAAGTAAAaataa                          <  1:434071/145‑1 (MQ=255)                                                                                                                              ctaaataaaaaA‑CTCTAGGAAAGT‑AAAATTA‑ACG‑GCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAGTGACTGAAAGTAAAAATAATCtt                      <  2:900862/147‑1 (MQ=255)                                                                                                                                 taataaaaccccttgtgAAAGT‑GAAAAAA‑ACC‑CCACCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAGTGACTGAAAGTAAAAATAATCTTaa                    <  1:678640/137‑1 (MQ=255)                                                                                                                                 aaaaaaaaa‑CTCTGTAACAGT‑AAAATAA‑ACG‑TCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAGTGACTGAAAGTAAAAATAATCTTAAtt                  <  1:162529/145‑1 (MQ=255)                                                                                                                                 aaaaaaaaa‑CTCTGTAAAGGT‑GAAAATA‑CCG‑CCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAGTGACTGAAAGTAAAAATAATCTTAAtt                  <  2:929377/145‑1 (MQ=255)                                                                                                                                    aataaA‑CTCTGTGACAGT‑GGAATTA‑ACG‑GCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGc               >  1:542250/1‑151 (MQ=255)                                                                                                                                    aataaA‑CTCTGTGACAGT‑GGAATTA‑ACG‑GCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGc               >  1:1188125/1‑151 (MQ=255)                                                                                                                                    aataaA‑CTCTGTGACAGT‑GGAATTA‑ACG‑GCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGc               >  2:686008/1‑151 (MQ=255)                                                                                                                                    aaaaaacctctGTGAAAGG‑AAAATAA‑ACG‑GCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCg                <  2:255631/144‑1 (MQ=255)                                                                                                                                     ataaA‑CTCTGTGACAGT‑GGAATTA‑ACG‑GCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGcc              >  2:380681/1‑151 (MQ=255)                                                                                                                                      taaA‑CTCTGTGACAGT‑GGAATTA‑ACG‑GCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCt             >  2:956063/1‑151 (MQ=255)                                                                                                                                        aa‑CTCTGTGACAGT‑GGAATTA‑ACG‑GCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTAt           >  2:1304563/1‑151 (MQ=255)                                                                                                                                         a‑CTCTGTGACAGT‑GGAATTA‑ACG‑GCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTCATGAAGACAGCG‑AAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTAtt          >  2:329925/1‑151 (MQ=255)                                                                                                                                            tctGTGACAGT‑GGAATTA‑ACG‑GCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCt        >  2:774576/1‑151 (MQ=255)                                                                                                                                             cggtgACAGT‑GAAATAA‑CCG‑GCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTc         <  1:66797/147‑1 (MQ=255)                                                                                                                                               gggACAGT‑GAAATAA‑ACG‑GCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCt        <  1:279423/146‑1 (MQ=255)                                                                                                                                                 gACAGT‑GGAATTA‑ACG‑GCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTtgtg   >  1:1238518/1‑151 (MQ=255)                                                                                                                                                  aCAGT‑GAATTAA‑ACG‑GCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGc  <  1:943628/151‑1 (MQ=255)                                                                                                                                                         |                                                                                                                                                    ATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTT‑CTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACAGT‑GGAATTA‑ACG‑GCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGC  >  NC_000913/2481656‑2481948

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