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.8%	114.2 / 33.2	71	P125S (CCA→TCA)	emrY	putative multidrug efflux system
Reads supporting (aligned to +/- strand):  ref base G (38/4);  new base A (0/29);  total (38/33)
Fisher's exact test for biased strand distribution p-value = 5.98e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCC‑TTAATAATAAA‑CTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACC  >  NC_000913/2481654‑2481953                                                                                                                                                         |                                                                                                                                                      atatCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCC‑TTAATAATAAA‑CTCTGTGACAGTgg                                                                                                                                                       >  1:318436/1‑151 (MQ=255)     ccAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCC‑TTAATAATAAA‑CTCTGTGACAGTGGAAtt                                                                                                                                                   >  1:337745/1‑151 (MQ=255)     ccAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCC‑TTAATAATAAA‑CTCTGTGACAGTGGAAtt                                                                                                                                                   <  2:563265/151‑1 (MQ=255)         ccccAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCC‑TTAATAATAAA‑CTCTGTGACAGTGGAATTAACg                                                                                                                                               >  2:1156861/1‑151 (MQ=255)          cccAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCC‑TTAATAATAAA‑CTCTGTGACAGTGGAATTAACgg                                                                                                                                              <  2:949995/151‑1 (MQ=255)           ccAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCC‑TTAATAATAAA‑CTCTGTGACAGTGGAATTAACGGc                                                                                                                                             >  2:705743/1‑151 (MQ=255)             aGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCC‑TTAATAATAAA‑CTCTGTGACAGTGGAATTAACGGccc                                                                                                                                           <  1:1083250/151‑1 (MQ=255)              gCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCC‑TTAATAATAAA‑CTCTGTGACAGTGGAATTAACGGccca                                                                                                                                          <  2:32441/151‑2 (MQ=255)                tAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCC‑TTAATAATAAA‑CTCTGTGACAGTGGAATTAACGGCCCCGc                                                                                                                                        >  1:252627/1‑151 (MQ=255)                    gTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCC‑TTAATAATAAA‑CTCTGTGACAGTGGAATTAACGGCCCCGCCAtt                                                                                                                                    >  1:33163/1‑151 (MQ=255)                      tATCACAAATATAACCGCCTAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCC‑TTAATAATAAA‑CTCTGTGACAGTGGAATTAACGGCCCCGCCATTa                                                                                                                                   >  1:482169/1‑150 (MQ=255)                         cacaAAAATACCCGCCAAATACCGGCCCAAATATGGGAGCGTTAATCACGGTAATTGACCATAATGCCAGAGCAAATTTTTTTTTTTCGGGCGAATACTTCC‑TAAAAAATAAC‑CTCTTAGACGATGAAATTAACGGCCCCGCCATTAAcc                                                                                                                                <  2:68492/150‑1 (MQ=255)                                ataACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCC‑TTAATAATAAA‑CTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACg                                                                                                                        >  2:387940/1‑151 (MQ=255)                                 taACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCC‑TTAATAATAAA‑CTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACg                                                                                                                        >  1:295122/1‑150 (MQ=255)                                    ccGCCCATTATCGGCCCACATATCGGACCGAAAGCCACGTTTATTTCCCATAATGCCAGAGCAAATTTTTTTTTTTGGGGGGGAAAATTCC‑TTAAAAATAAA‑CTCTGGAACAGTGAAATAAGCGGCCCCGCCATTAACCCCTGAACGAct                                                                                                                     <  1:18332/150‑1 (MQ=255)                                    ccGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCC‑TTAATAATAAA‑CTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGActc                                                                                                                    >  1:481719/1‑151 (MQ=255)                                       cccAATATCGCCCCACATATCGGAGCAATAATCACGGTCATTGACCATAATGCCAGAGAAAATTTTCTTTTTTCTGGCGGATAATTCC‑TAAATAAAAAC‑CTCTGTGACAGTGAAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTa                                                                                                                  <  1:184235/150‑1 (MQ=255)                                        ccAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCC‑TTAATAATAAA‑CTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTaa                                                                                                                 >  1:121983/1‑150 (MQ=255)                                                 gCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCC‑TTAATAATAAA‑CTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATc                                                                                                         >  1:22126/1‑149 (MQ=255)                                                  cccACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCC‑TTAATAATAAA‑CTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGc                                                                                                      >  1:230869/1‑151 (MQ=255)                                                   ccACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCC‑TTAATAATAAA‑CTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGc                                                                                                      >  2:277787/1‑150 (MQ=255)                                                    caaatatCGGAGCGATAATCACGTTCATTAACAATAATGCAAGAAAAAATTTTTTTTTTTCGGGCAAATAATTCC‑TAAAAAAAAAT‑CTCTGTGACAGTGAAATAAACGACCCCGCCATTAACCCCTGAACGACTCTAAAAAatat                                                                                                          <  2:1002582/142‑1 (MQ=255)                                                     acaTATCGGGGCAATAATCCCGGAAATTACCCTGAAAGCCAGAACAAATGTTTTTTTTCCTGGCGAAT‑ATTCC‑ATAAAAAAAAACCTCTGTGACAGTAAAATAACCGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCa                                                                                                        <  2:642520/146‑1 (MQ=255)                                                     acaTATCGGAGCAATAACCATGGTAATTGACAATAAGGCCAAAGAAAATGTTCTTTTTTCGGGGGAATAATCCC‑TAAAAAAAAAA‑CCCTGTGACAGGGAAATAACCGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGcac                                                                                                    <  2:348475/150‑1 (MQ=255)                                                         agcGGAGCGATAATCACGGTCATTGCCCAAAATGCCAGAGAAATTTTTTTTTTTTCTGGCGAATATTTCC‑TTAATAAAAAC‑CTCTGTGACAGTGAAATAAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCa                                                                                                <  1:165321/148‑1 (MQ=255)                                                            ggggcaaTAATCACGTTCATTACCCAAAATGCAAGAAAAATTTTTTTTTTTTCTGGCGGATATTTCC‑TAAAAAAAAAC‑CTCTGTAACAGTAAAATAACCGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAg                                                                                              <  1:166493/143‑1 (MQ=255)                                                            ggAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCC‑TTAATAATAAA‑CTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATc                                                                                                 >  2:942826/1‑146 (MQ=255)                                                              aGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCC‑TTAATAATAAA‑CTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTg                                                                                          >  2:1078320/1‑151 (MQ=255)                                                                  aTAATCACGGTCATTGCCCAAAATGCAAGAGAAAATGTTTTTTTTTCTGGCGGATATTTCC‑TAAAAAAAAAA‑CTCTGGGACAGTGAAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGt                                                                                        <  1:1028673/149‑1 (MQ=255)                                                                    aaTCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCC‑TTAATAATAAA‑CTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGAt                                                                                    >  1:829248/1‑151 (MQ=255)                                                                        aCGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCC‑TTAATAATAAA‑CTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGc                                                                                >  1:412417/1‑151 (MQ=255)                                                                          gttCATTGACAAAAATGCCAGAGCAAATGTTTTTTTTTCTGGCGAATAATTCC‑TTAAAAAAAAA‑CTCTGTGAAAGTGAAATAACCGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAg                                                                                  <  2:863946/145‑1 (MQ=255)                                                                             cATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCC‑TTAATAATAAA‑CTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTacac                                                                           >  1:1027368/1‑151 (MQ=255)                                                                             cATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCC‑TTAATAATAAA‑CTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTacac                                                                           >  2:512301/1‑151 (MQ=255)                                                                               ttGACCAAAATCCCAGAG‑AATTTTTTTTTTTTTCGGGGGATTATTCCTTTAATAATAAC‑CTCTGGGACAGTGAAATAACCGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTaca                                                                            <  1:124765/148‑1 (MQ=255)                                                                                 gACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCC‑TTAATAATAAA‑CTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACAtt                                                                        >  2:638849/1‑150 (MQ=255)                                                                                     aaaaggcaagagAAAATGTTTTTTTTTCTGGCGGATAATTCC‑TTAAAAAAAAC‑CTCTGTGACAGTGAAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGa                                                                   <  1:562571/143‑1 (MQ=255)                                                                                     aTAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCC‑TTAATAATAAA‑CTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATg                                                                    >  2:727302/1‑150 (MQ=255)                                                                                          actagagAAAATGTTCTTTTTTCGGGCGAATAATTCT‑TAAAAAAAAAA‑CTCTGAGAAAATGAAATAACCGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGa                                                                   <  2:476010/143‑1 (MQ=255)                                                                                             agagCAAATGTTCTTTTTTCTGGCGGATAATTCC‑TTAATAATAAA‑CTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCg                                                           >  1:484610/1‑151 (MQ=255)                                                                                               agCAAATGTTCTTTTTTCTGGCGGATAATTCC‑TTAATAATAAA‑CTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGa                                                          >  1:1159010/1‑150 (MQ=255)                                                                                                gCAAATGTTCTTTTTTCTGGCGGATAATTCC‑TTAATAATAAA‑CTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGaa                                                         >  2:605017/1‑150 (MQ=255)                                                                                                  aattttttttttttCTGGCGGATAATTCC‑TAAAAAAAAAC‑CTCTGGGACAGTGAAATAAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGaa                                                         <  1:698748/140‑1 (MQ=255)                                                                                                      gTTCTTTTTTCTGGCGGATAATTCC‑TTAATAATAAA‑CTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGaaaaaaa                                                    >  2:839888/1‑149 (MQ=255)                                                                                                              ttCTGGCGGATAATTCC‑TTAATAATAAA‑CTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGGTAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGa                                           >  2:814247/1‑150 (MQ=255)                                                                                                              ttCTGGCGAATAATTCC‑TTAAAAATAAA‑CTCGGTGACAGGGAAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGa                                           <  1:352897/150‑1 (MQ=255)                                                                                                                cgggCGGATAATTCC‑TAAAAAATAAA‑CTCGGTGACAGAGAAATAAACGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGaaa                                         <  2:215180/148‑1 (MQ=255)                                                                                                                 gggCGGATAATTCC‑TAAAAAAAAAC‑CTCTGTCACGGTGAATTAAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAg                                        <  1:943272/149‑1 (MQ=255)                                                                                                                 ggacgaatatTTCC‑TAAAAAAAAAC‑CTCTGTGACAGGGAAATAACCGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGaaa                                         <  2:295050/143‑1 (MQ=255)                                                                                                                  ggCGGATAATTCC‑TAAAAAAAAAA‑CTCTGAAACAGGAAAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTa                                      <  2:630560/151‑1 (MQ=255)                                                                                                                    cGGATAATTCC‑TTAATAATAAA‑CTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTaa                                     >  1:580841/1‑150 (MQ=255)                                                                                                                     ggATAATTCC‑TTAATAATAAA‑CTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTaaa                                    >  2:541382/1‑150 (MQ=255)                                                                                                                       aTAATTCC‑TTAATAATAAA‑CTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAaat                                 >  2:364815/1‑151 (MQ=255)                                                                                                                         aaTTCC‑TTAAAAAAAAA‑CTCTGTGACAGTGAAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTaaa                                    <  1:808700/146‑1 (MQ=255)                                                                                                                          aTTCC‑TTAATAATAAA‑CTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAaataat                              >  1:745732/1‑151 (MQ=255)                                                                                                                            tCC‑TTAATAATAAA‑CTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCt                            >  2:74998/1‑151 (MQ=255)                                                                                                                             cc‑TTAATAATAAA‑CTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCtt                           >  2:384511/1‑151 (MQ=255)                                                                                                                             cc‑TTAATAATAAA‑CTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCtt                           >  1:684379/1‑151 (MQ=255)                                                                                                                                ttAATAATAAA‑CTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTa                          >  2:601941/1‑150 (MQ=255)                                                                                                                                  aaaaataaC‑CTCTGTAACAGAGAAATAACCGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCtt                           <  2:642675/144‑1 (MQ=255)                                                                                                                                  aaaaaaaaa‑CTCTGTAACAGTGAAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTa                          <  1:380770/142‑1 (MQ=255)                                                                                                                                   ataataaA‑CTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTc                      >  2:952295/1‑151 (MQ=255)                                                                                                                                   aaaaaaaacctctGTGACGGTAAAATAACCGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTaa                         <  1:859258/140‑1 (MQ=255)                                                                                                                                   aaaaaaaacctctGTGACAGTGAAATAAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTaa                         <  1:113937/140‑1 (MQ=255)                                                                                                                                          aattctGTGAAAGGGAAATTACCGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGcc                   <  2:1027086/145‑1 (MQ=255)                                                                                                                                            ctttgtgAAAGGGAAATAACCGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTAt                <  1:1124859/146‑1 (MQ=255)                                                                                                                                               tgtgACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCttt           >  1:811295/1‑151 (MQ=255)                                                                                                                                                gtgACAGTGAATTTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCttt           <  1:545090/150‑1 (MQ=255)                                                                                                                                                   aCAGTGAAATAAACGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGc       <  2:1148171/151‑1 (MQ=255)                                                                                                                                                   aCAGTAAAATAAACGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTtgtg        <  2:401853/150‑1 (MQ=255)                                                                                                                                                     aGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGcc      >  2:1103336/1‑150 (MQ=255)                                                                                                                                                       tGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAAc   >  1:533286/1‑151 (MQ=255)                                                                                                                                                        ggAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAAcc  >  1:780870/1‑151 (MQ=255)                                                                                                                                                         |                                                                                                                                                      ATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCC‑TTAATAATAAA‑CTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACC  >  NC_000913/2481654‑2481953

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