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

ATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACA‑GTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAAC  >  NC_000913/2481656‑2481952                                                                                                                                                        |                                                                                                                                                     atCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACA‑GTGGa                                                                                                                                                     <  2:531267/150‑1 (MQ=255)    cAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACA‑GTGGAAtt                                                                                                                                                  >  2:280875/1‑150 (MQ=255)    cAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACA‑GTGGAATTa                                                                                                                                                 <  2:523186/151‑1 (MQ=255)     aaCCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACA‑GTGGAAtt                                                                                                                                                  >  1:683556/1‑149 (MQ=255)     aaCCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACA‑GTGGAATTa                                                                                                                                                 >  1:899302/1‑150 (MQ=255)      aCCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACA‑GTGGAATTAAc                                                                                                                                               <  2:252493/151‑1 (MQ=255)      aCCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACA‑GTGGAATTAAc                                                                                                                                               >  1:403967/1‑151 (MQ=255)       ccccAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACA‑GTGGAATTAACg                                                                                                                                              <  2:407062/151‑1 (MQ=255)           aGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACA‑GTGGAATTAACGGccc                                                                                                                                          <  1:66004/151‑1 (MQ=255)                                 cccccccAATACCGCCCCACATAGGGGAGCAATAATCCCGGTCATTGCCAATAATGCAAGAAAAATTGTTTTTTTTTCGGGCGAATATTTCCTA‑AAAAAAAAC‑CTCTGGGAAG‑GTGGAATAACCGCCCCCGCCATTAACCCCTGAACGAct                                                                                                                    <  1:407031/147‑1 (MQ=255)                                  ccGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACA‑GTGGAATTAACGGCCCCGCCATTAACCCCTGAACGActc                                                                                                                   >  1:73043/1‑151 (MQ=255)                                  ccGCCAAAAACCGCCCCACATAGCGGAGCGATAATCGCGGTCTTTGCCAAAAATGCCAGAGAAATTGTTCTTTTTTCGGGGGAATAATTCCTA‑AAAAACAAC‑CTCTGTGACG‑GTGAAATTAGCGCCCCCGCCATTAACCCCTGAACGAct                                                                                                                    <  1:50483/150‑1 (MQ=255)                                   ccccaaaTACCGCCCCACAAACCGGAGCAATAATCCCGGTCATTGACCATAATGCCAGAGAAATTTTTTTTTTTTCGGGCGAATAATTCTTA‑AATAAAAAC‑CTCTGGGACA‑GTGGAATTACCGGCCCCGCCATTAACCCCTGAACGActct                                                                                                                  <  2:858841/146‑1 (MQ=255)                                   cGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACA‑GTGGAATTAACGGCCCCGCCATTAACCCCTGAACGActct                                                                                                                  >  1:119237/1‑151 (MQ=255)                                       cAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACA‑GTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTaaa                                                                                                               >  2:541989/1‑150 (MQ=255)                                         atatCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACA‑GTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTaaaa                                                                                                              >  1:1038045/1‑149 (MQ=255)                                          taccgccccACATATGGGGGCGATAATCACGTTCATTAACAATAAGGCCAGAGCAAGTTTTTTTTTTTCGGGAGAATAATCCCTA‑AAAAATAAA‑CTCCGTAACA‑GAAAAATTACCGGCCCCGCCATTAACCCCTGAACGACTCTaaaaaa                                                                                                            <  1:880911/144‑1 (MQ=255)                                           atCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACA‑GTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAat                                                                                                           >  1:511042/1‑150 (MQ=255)                                                cccACATAGCGGAGCAATAACCACGGTCATTGACAAAAATGCAAGAGAAATTGTTTTTTTTTCGGGCGAATAATTCCTT‑AAAAAAAAA‑CTCTGTAAAA‑AAAAATTTAACGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCa                                                                                                       <  1:179410/149‑1 (MQ=255)                                                        tCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACA‑GTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCa                                                                                               >  2:193597/1‑149 (MQ=255)                                                          ggAGCAATAATCACGGTCATTGCCCAAAATGCAAGAGCAAATTTTTTTTTTTCGGGCGAATAATTCTTA‑AAAAAAAAA‑CTCTGGAAAA‑GTGAAATAAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAg                                                                                             <  2:601783/149‑1 (MQ=255)                                                             gtaaTAATCCCGGTCATTGACTATAATGCCAGAGAAATTTTTTTTTTTTCGGGGGGATATTTCTTT‑AAAAAAAAC‑CTTGGTGACG‑GTGAAATAAACGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTg                                                                                         <  2:680544/147‑1 (MQ=255)                                                             gCGATAATCGCGTTCATAGACCATAATCCCAGAAAAATTTTTTTTTTTTCGGGGGAATAA‑TCCTT‑AAAAAAAAACCTCTGGGACA‑GAGAAATAACCGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGAtt                                                                                          <  1:775074/149‑1 (MQ=255)                                                             gCGATAATCCCGGTCATTAACAATAATGCAAGAGAAATTTTTTTTTTTTTGGGCGAATATTTCTTA‑AAAAAAAAC‑CTCTGTGACA‑GTGAAATTAACGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTg                                                                                         <  2:94917/150‑1 (MQ=255)                                                             gCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACA‑GTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTgg                                                                                        >  1:144780/1‑151 (MQ=255)                                                                    tCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACA‑GTGGAATTAACGGCCCCGCAATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGAt                                                                                   >  1:982446/1‑149 (MQ=255)                                                                      aCGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACA‑GTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGc                                                                               >  2:395638/1‑151 (MQ=255)                                                                      aCGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACA‑GTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGc                                                                               >  1:1017680/1‑151 (MQ=255)                                                                             ttGACAAAAATCCCAAAAAAAATTTTTTTTTTTCGGGGGAATAATTCCTA‑AAAAATAAC‑CTCTTTGAAA‑GAGAAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTAcaca                                                                         <  2:786247/150‑1 (MQ=255)                                                                               gACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACA‑GTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACAtt                                                                       >  2:770566/1‑150 (MQ=255)                                                                                  caaaatgcagaagCAATTTTTTTTTTTTGGGGCGAATAATCACTA‑AAAAACAAC‑CTCTGTGAGA‑GGGAAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATg                                                                   <  2:315521/140‑1 (MQ=255)                                                                                  cATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACA‑GTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATg                                                                   >  2:891574/1‑151 (MQ=255)                                                                                   aaaataaaagagCAATTTTTTTTTTTCCGGGGGGATAATTCATA‑AATAAAAAA‑CTCGGTGACA‑GTAAAATTAGCGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGt                                                                                       <  1:214244/122‑1 (MQ=255)                                                                                   aTAATGCCAGAGAAATTTTTTTTTTTTCGGGCGAATATTTCCTA‑AAAAAAAAC‑CCCTGTGACA‑GTGAAATAAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATg                                                                   <  2:49689/150‑1 (MQ=255)                                                                                   aTAATCCCAGAGAAATTTTTTTTTTTTCTGGCGGATAATTCCTA‑AAAAAAAAC‑CCCTGTGACA‑GTGAAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACAtt                                                                       <  2:263090/146‑1 (MQ=255)                                                                                     aaTGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACA‑GTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGa                                                                  >  1:851085/1‑149 (MQ=255)                                                                                                aattttttttttttCGGGGGAATAATTCCTA‑AAAAAAAAC‑CTTGGTAACA‑GTGAAATAACCGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGaaaaa                                                     <  2:335117/143‑1 (MQ=255)                                                                                                aattttttttttttCGGGCGAATAATTCATT‑AATAAAAAC‑CTCTGGGAAA‑GTGAAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGaaa                                                       <  2:44526/141‑1 (MQ=255)                                                                                                aattttttttttttCGGGAGGATAATCCTTA‑AAAAATAAA‑TTCTGTGAAA‑GTGAAATAAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGa                                                         <  1:268157/139‑1 (MQ=255)                                                                                                aaaTTTTTTTTTTTCTGGGGAAAAATTCATT‑AAAAAAACC‑CTC‑GTGAAAGGGGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGaaaa                                                      <  2:587583/150‑1 (MQ=255)                                                                                                 aaTTTTTTTTTTTCGGGGGAATAATTCCTA‑AAAAAAAAC‑CTCTGGGACA‑GTGAAATAACCGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGaaaaaa                                                    <  2:1029351/151‑1 (MQ=255)                                                                                                 aaTGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACA‑GTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGaaaaa                                                     >  1:708860/1‑150 (MQ=255)                                                                                                   tgttgtTTTTTCTGGGGGATTATTCCTT‑AAAACCAAC‑CTCTGCGAGA‑GTGAATATAGCGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGa                                                         <  2:792131/144‑1 (MQ=255)                                                                                                        ttttttCTGGCGGAAAGTTCCTT‑AATAATAAA‑CTCTCTGACA‑GTGGAATTAACGGCCCCGCCATTACCTCCTGAACGACTCTAAAAAATATCAGCACATCAGGATTGGTCGATAGGCAACACATTAATGAAGACAGCGtaaaaattgcga                                              >  1:176396/1‑144 (MQ=255)                                                                                                          ttttgtgaaaaaaaaTCCCTA‑AAAAATAAA‑CTATGTAACA‑GTGAATTAAACTGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACt                                            <  1:524897/137‑1 (MQ=255)                                                                                                          ttttCGGGCGAATAATTCCTA‑AAAAAAAAC‑CTCTGGAACA‑GGGGAATAACCGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACt                                            <  1:214250/150‑1 (MQ=255)                                                                                                           tttCTGGCGAATATTTCCTA‑AAAAAAAAC‑CTCTGAGAAA‑GGAAAATTACCGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTg                                           <  1:414879/150‑1 (MQ=255)                                                                                                                 ggcgaataTTCCTTAAAAAAAAAA‑CTCTGTGACA‑GGGAAATTACCGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAg                                       <  1:1002127/142‑1 (MQ=255)                                                                                                                 ggaaaaaaTTTCTTAAATAATAAA‑CTCGTGGAAA‑ATAAAATTACGGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGt                                      <  1:502338/144‑1 (MQ=255)                                                                                                                   ggATAATTCCTT‑AATAATAAA‑CTCTGTGACA‑GTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTaaa                                   >  1:298939/1‑150 (MQ=255)                                                                                                                       aaTTCCTT‑AATAAAAAC‑CTCTGTGACA‑GTGAAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAaataa                              <  2:458880/151‑1 (MQ=255)                                                                                                                         ttCCTT‑AATAATAAA‑CTCTGTGACA‑GTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATc                            >  1:461102/1‑151 (MQ=255)                                                                                                                         ttCCTA‑AAAAAAAAA‑CTGTGTGACA‑GAGAAATTACCGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAaataa                              <  1:867394/149‑1 (MQ=255)                                                                                                                                aaaaataaC‑CTCTGTGACA‑GTGAAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAaat                                <  2:32980/138‑1 (MQ=255)                                                                                                                                 ataataaA‑CTCTGTGACA‑GTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTc                     >  1:364127/1‑151 (MQ=255)                                                                                                                                 ataaaaaC‑CTCTGGGACA‑GGGAAATTACCGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTc                     <  2:133149/151‑1 (MQ=255)                                                                                                                                   aataaA‑CTCTGTGACA‑GTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGc                   >  1:1045952/1‑151 (MQ=255)                                                                                                                                   aaacaacccctgcaaCA‑ATGAAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCg                    <  2:526232/137‑1 (MQ=255)                                                                                                                                   aaaaaa‑CGCTGTAACA‑GAAAAATTAGCGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTc                     <  1:76681/146‑1 (MQ=255)                                                                                                                                    ataaA‑CTCTGTGACA‑GTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCATAATTCGcc                  >  1:948519/1‑151 (MQ=255)                                                                                                                                      aaa‑CTCTGTGACA‑GTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCt                 >  2:515352/1‑150 (MQ=255)                                                                                                                                        a‑CTCTGTGACA‑GTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTAtt              >  1:978982/1‑151 (MQ=255)                                                                                                                                            ctGTGACA‑GTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCtt           >  1:162649/1‑151 (MQ=255)                                                                                                                                              gtgACA‑GTGAAATAACCGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCt            <  2:362967/148‑1 (MQ=255)                                                                                                                                               tgACA‑GTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTtg         >  1:385593/1‑150 (MQ=255)                                                                                                                                                 aCA‑GTGAAATAAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGc      <  2:852206/151‑1 (MQ=255)                                                                                                                                                   a‑GTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTAAACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGcc     >  2:388221/1‑150 (MQ=255)                                                                                                                                                      tGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAAc  >  2:938523/1‑151 (MQ=255)                                                                                                                                                        |                                                                                                                                                     ATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTT‑AATAATAAA‑CTCTGTGACA‑GTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAAC  >  NC_000913/2481656‑2481952

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