Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,483,254	0	G	A	44.9%	76.2 / 34.6	69	S29L (TCA→TTA)	emrK	multidrug resistance efflux pump membrane fusion protein
Reads supporting (aligned to +/- strand):  ref base G (25/13);  new base A (0/31);  total (25/44)
Fisher's exact test for biased strand distribution p-value = 1.33e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCAT  >  NC_000913/2483108‑2483403                                                                                                                                                   |                                                                                                                                                      gTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGaaaa                                                                                                                                                   <  1:640824/151‑1 (MQ=255)   cGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACg                                                                                                                                                 <  2:570907/151‑1 (MQ=255)     aaCGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCa                                                                                                                                               <  2:475730/151‑1 (MQ=255)      aCGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCaa                                                                                                                                              <  1:523793/151‑1 (MQ=255)        gTAGTTCGTATTTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCaat                                                                                                                                             >  2:76847/1‑150 (MQ=255)          aGTTCGTATCTTTATGATTAACGACAGTGACACTCCCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCaataa                                                                                                                                           <  2:747123/150‑1 (MQ=255)          aGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCaataa                                                                                                                                           <  2:388774/150‑1 (MQ=255)             tCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCaataaata                                                                                                                                        <  1:1049900/150‑1 (MQ=255)                   cTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACACAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTAc                                                                                                                                 >  1:232595/1‑151 (MQ=255)                     ttATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTAcc                                                                                                                                <  1:643077/150‑1 (MQ=255)                        tGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTTTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGcc                                                                                                                             >  1:1039871/1‑150 (MQ=255)                                        cccTACCTGAGACTTGTGCAAAAATGAGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCAAATCCCATTGCCCAATAGACATAGGCACCTAAAAACGCAATAAAAAAAACTACCGCCAATAAAGAAAAGTAttt                                                                                                            <  2:483741/149‑1 (MQ=255)                                        cacTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTATAATTCCATTGACCAATAGGCAAAGACACCTAAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTAtt                                                                                                             <  2:791983/150‑1 (MQ=255)                                        cacTACCTGAGACTTGTACAAAAATTGGATCTGCTTCCCCCGTGACATAGGCGTCATTTGTACAAATCATTTCTTCTAATTCCATTGACAAATAGGAAAAGCCACCTAAAAACGCAAAAAAAAAAACTACCGCCAATAAAGAAAAGTAtt                                                                                                             <  2:1057779/150‑1 (MQ=255)                                         acTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTAtttt                                                                                                           >  1:391454/1‑151 (MQ=255)                                         acTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTAtttt                                                                                                           >  2:1041223/1‑151 (MQ=255)                                            acttGAGATTTGTACAAAAATTGAATCTGCATTCCCCTGGACATGGGGGTCATCTGTTCTATTCATGTTTTCTAATCCCATTGACCAATAGAAATAGGCACTTAAATACTCAATAAAAAAACTTACCGCCAATAAAGAAAAGTATTttctt                                                                                                        <  2:574475/148‑1 (MQ=255)                                             ccTGAGCTTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATCCCATTGCCCAATAGGCATAGGCACCTAAAACCGCAATAAAAAAAACTACCGCCAATAAAGAAAAGTATTttct                                                                                                         <  1:827248/149‑1 (MQ=255)                                             ccTGAGATTTGTGAAAAAATTGGTTCTGCATTCCCCGTGACATAGGCGTCATTTGTACTAATAATGTCTTCTAATCCCATTGACAAATAGGAATAGGCACCTAAAAACGAAAAAAAAAAAACTACCGCCAATAAAGAAAAGTATTttcttc                                                                                                       <  1:540776/151‑1 (MQ=255)                                             ccTGAGACTTGTGCAAAAATTGCATCTGTATCCCCCGAGACATAGGCGTCATCTGTACTAATCATTTCTTCAAAGTCCATTAAACAAAAGACATAGGCACTTGAAAACACAATAAATAAAACTACCGCCAATAAAGAAAAGTAttt                                                                                                            <  2:380112/146‑1 (MQ=255)                                               tGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTg                                                                                                     >  1:601874/1‑151 (MQ=255)                                               tGAGACTTGTGCAGAAATTGGATCTACATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTAAAAACGCAATAAAAAAAACTACCGCCAATAAAGAAAAGTATTttctt                                                                                                        <  1:868607/148‑1 (MQ=255)                                                 agaCTTGTGCAAAAATTGGATCTGCATTCCCCGTAACATAGGCTTCATCTTTACTAATCATGTCTTCAAATCCAATTGCCAAATAGGCAAAGCCACTTAAAAACCCAATAAATAAACCTACCGCCAATAAAGAAAAGTATTTTCTTCTGtt                                                                                                   <  2:300601/151‑1 (MQ=255)                                                    cTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGtt                                                                                                   >  1:544512/1‑148 (MQ=255)                                                     ttGTGCAAAAATGGGATCTGCATCCCCCGTGAAATAGGCTTCATCTGTACTAATCATGTCTTCTAATCCCATTGACCAATAGGCATAGGCCCCTAAAAACGAAAAAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAg                                                                                                 <  1:893678/149‑1 (MQ=255)                                                     ttGTGAAAAAATTGTATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTATAATTCCATTGACCAATAGACATAGGCCCCTAAAACCGAAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGtt                                                                                                   <  2:431887/147‑1 (MQ=255)                                                       ttgCAAAAATTGTATCTCCATTCCCCTTAACATAGCCTCCATCTTTACTAATAATTTCTTAAAATACAATTGACCAATAGACATAGGCACCTAAAACCGAAATAAAAAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAAt                                                                                              <  2:370166/149‑1 (MQ=255)                                                        tgCAGAAATTGGATCTGCTTCCCCCGTGACATAGGCTTCATCTGTACTAATCATTTCTTCTAATCCCATTAACCAATAGGCATAGCCCCCTAAAAACGCAATAAAAAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATg                                                                                             <  1:251912/150‑1 (MQ=255)                                                        tgCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGt                                                                                            >  2:939152/1‑151 (MQ=255)                                                         gCAGAATTGGGTTTTGCATCCCCCGTGACATAGGCGTTATCTGTTCTAAGCATGTTTTCTATTTCCATTGACCAAGAGGCAAAGGCACTAAAAAACGAAATAAAAAAACCTCCCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGaa                                                                                               <  1:723408/147‑1 (MQ=255)                                                              atttGGATCTGCATCCCCCGTAACATAGGCGTCATCTTTACTAATCATGTCTTCTAATTCCATTGACCAATAGGAAAAGGCCCCGAAAACCCCAATAAAAAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTa                                                                                      <  1:570273/149‑1 (MQ=255)                                                                ttGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTAtt                                                                                    >  1:229346/1‑151 (MQ=255)                                                                 tGGATCTGCATTCCCCGTAACATAGGCGTCATCTGTACTAATCATTTTTTCTAATCCCATTAACCAATAGGCATAGCCACCTAAAAACTCAATAAAAAAAACTTCCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTAtt                                                                                    <  2:990873/150‑1 (MQ=255)                                                                    aTCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGa                                                                                 >  1:185187/1‑150 (MQ=255)                                                                      atgctgtcaccGTGACATACGCCTTAACTGTATTAATCATGTCTTCTACTTCCATTGAACAAGAGGGCTAGGCAGATGAAAACGCAATAAATAAAACTAACGCAAATAAAAACAAGTATTTTCTTGTGTAAGAATGATTTGTATTTGaaat                                                                              >  2:468182/10‑149 (MQ=255)                                                                        gCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCCCCTAAAAACGCAAAAAAAAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGa                                                                                 <  1:408318/146‑1 (MQ=255)                                                                        gCATCCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCAAATCCCATTGACAAATAGGCATAGGCCCCTAAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTa                                                                             <  1:1055296/150‑1 (MQ=255)                                                                          aTTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATGGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATc                                                                          >  1:180853/1‑151 (MQ=255)                                                                           ttCCCCGTAAAATAGGCCTTATTTGTACAAATCATGTCTTCTAATCCCATAGACCAATAGACATAGGCACCTAAAAACGAAAAAAATAAAATTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAAt                                                                               <  1:882790/145‑1 (MQ=255)                                                                           tcccccGTGACATAGGCGTCATCTGTACTAATCATTTCTTCTAATTCCATTGCCCAATAGGCATAGCCCCCTAAAAACGCAATAAAAAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTaa                                                                            <  1:334870/146‑1 (MQ=255)                                                                            tCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGAATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTg                                                                        >  1:88956/1‑151 (MQ=255)                                                                             ccccGTGACATAGGGGTCATCTGTCCAAATCATGTCTTATAATTCCATTGACCAAAAGGTATAGGCACCTAAAAACCAAAAAAATAAACTTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGt                                                                       <  2:815316/151‑1 (MQ=255)                                                                                 gTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTcc                                                                    >  1:869313/1‑150 (MQ=255)                                                                                    aCATAGCCGTCATCTATACAAATAATGTCTTCTAATTCAATTAACCAATAGGAATAGCCACCTAAAAACGCAATAAAAAAAACTTCCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCAc                                                                  <  2:453733/149‑1 (MQ=255)                                                                                      aTAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACtatt                                                              >  2:598933/1‑151 (MQ=255)                                                                                               aTCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGAAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCa                                                      <  1:599703/150‑1 (MQ=255)                                                                                                tCTGTACTAATCATGTCTTCTAATTCCATTGACAAATAGGCATAGGCACCTGAAAACGAAATAAAAAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTAtctctc                                                       <  1:318413/148‑1 (MQ=255)                                                                                                 cTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCAttt                                                   >  1:968801/1‑151 (MQ=255)                                                                                                      cTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCAt                                              >  2:883749/1‑151 (MQ=255)                                                                                                        aataatGTCTTCTAATCCAATTGACCAATAGGCAAAGGCACCTGAAAACGAAATAAAAAAAATTTCCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAg                                            <  2:81001/151‑1 (MQ=255)                                                                                                        aaTCTTGTCTTCTAATCCCATTGACCAATAGGCATAGGCACCTGAAAACACAATAAAAAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATa                                             <  1:940197/150‑1 (MQ=255)                                                                                                        aaTCATGTTTTCTAATCCCATTGCCCAAAAGGCATAGGCACCTAAAAACGAAAAAAAAAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCa                                               <  1:66326/148‑1 (MQ=255)                                                                                                        aaTCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCAt                                              >  1:726926/1‑149 (MQ=255)                                                                                                        aaTCATGTCTTCTAATTCCATTGACAAATAGACATAGGCACCTAAAAACGCAAAAAAAAAAACTTCCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCAttt                                                   <  2:828931/144‑1 (MQ=255)                                                                                                             tttcttctAATTCAATTACCCAATAGGCATAGACCCTTAAAACCACAATAAATAAAATTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCAt                                              <  2:32833/142‑1 (MQ=255)                                                                                                                  tctAATTCCATTGACCAATAGGCATAGGCACCTAAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAAt                                     <  1:703695/148‑1 (MQ=255)                                                                                                                     aaTTCCATTGACCAATAGGCATAGGCACCTAAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTa                                  <  1:377401/148‑1 (MQ=255)                                                                                                                       ttCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTAtgaat                             >  1:350096/1‑151 (MQ=255)                                                                                                                           aTTGACCAAAAGGCATAGCAACCTAAAAACCAAAAAAAAAAAATTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTAtgaattgaa                         <  2:423991/151‑1 (MQ=255)                                                                                                                            ttGACCAATAGGCATAGGCACCTGAAAACGCAATAACTAAAACTACAGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCACTTCTCATAGATGAAATTTAtgaattgaa                         >  2:1018595/1‑150 (MQ=255)                                                                                                                            ttGACCAATAGGCATAGGCACCTGAAAAAGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTAtgaattgaat                        >  1:1045499/1‑151 (MQ=255)                                                                                                                              aacaaataggaataggCACCTAAAAACGAAAAAAAAAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATtgaatga                      <  2:331407/147‑1 (MQ=255)                                                                                                                                 cAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGaaa                    >  1:842515/1‑150 (MQ=255)                                                                                                                                  aaTAGAAATAGGAACTTAAAAACGCAATAAATAAAACTTCCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGaaaa                   <  2:423656/150‑1 (MQ=255)                                                                                                                                    taggcataggcaCCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAag                  >  1:741649/1‑149 (MQ=255)                                                                                                                                           agcccccTAAAACCACAATAAAAAAAACTCCCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGATCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGcattca          <  1:76856/145‑1 (MQ=255)                                                                                                                                            ggcaCCTAAAAACGCAATAAAAAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGcattcat         <  1:413087/150‑1 (MQ=255)                                                                                                                                              caCCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCAttct      >  1:748249/1‑151 (MQ=255)                                                                                                                                                  tGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTACTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCat  >  2:186623/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                                      GTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCAT  >  NC_000913/2483108‑2483403

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