Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,601,560	0	A	C	43.9%	114.3 / 11.5	66	V176G (GTG→GGG)	ftsX	putative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (30/7);  new base C (0/29);  total (30/36)
Fisher's exact test for biased strand distribution p-value = 2.44e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTGCAACTGTGCC  >  NC_000913/3601411‑3601703                                                                                                                                                      |                                                                                                                                                ccGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATcac                                                                                                                                                >  1:323012/1‑151 (MQ=255) ccGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATcac                                                                                                                                                >  1:1274547/1‑151 (MQ=255)    aaaCGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATcaccac                                                                                                                                             <  2:629508/151‑1 (MQ=255)      aCGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCg                                                                                                                                           >  1:1089571/1‑151 (MQ=255)            ccGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTg                                                                                                                                     >  2:1236462/1‑151 (MQ=255)             cGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGc                                                                                                                                    >  1:251231/1‑151 (MQ=255)               ccccgCCCGGCCACCCCCGCCAGCCGGCAACACCAGGTGTCACCCACCCGCCCTTCTTCAATGCCGTTAATCTGCGTGATACGATCCCGCCGTGTATCCAGTGATCCCGCCCCCTGGAAATCGATTTTGGGGTTCCCCCCCGCCACTGCCg                                                                                                                                  <  2:120052/147‑1 (MQ=255)                     ccGGTCACCGCCGCCAGACGGGAAACCCAGCTGTCATCCACCCGCCCTTCGTCAATGCCGTAAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATCCCGCCCCCTGAAAATCGAGTTTCGGGATCCCCACCGCCACTGCCGGAAGCg                                                                                                                            <  1:575694/151‑1 (MQ=255)                     ccGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCg                                                                                                                            >  2:459634/1‑151 (MQ=255)                      gggTCACCGCCGCCAGAGGGGAAACCCGGCTGTCCCCCACCCGCACTTCTTAAATGCCGTTAATCTGGGTGATACGATCCCGCGGCGTATCCAGTTATCCCGCCCCCGGGAATTCGGGTTTGGGGACCCCCACCGCCACTGCCGGAAGCgg                                                                                                                           <  2:1308842/150‑1 (MQ=255)                              gccgccAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTttctt                                                                                                                   >  1:1051338/1‑151 (MQ=255)                              gccgccAGACGGGCAAACCAGCTGTCACCCACCCGCCCTTCGTCAATGCCGTTAATCTGGGTGATACGATCCCGCAGCGTATTCAGTGATCCCGCCCCCGGAAAATGGATTTTGGGGATCCCCACCGCCACTGCCGGAAGCGGGTTttctt                                                                                                                   <  1:1329130/151‑1 (MQ=255)                                    aGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGc                                                                                                              >  1:516824/1‑150 (MQ=255)                                    aGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGc                                                                                                              >  1:178558/1‑150 (MQ=255)                                      aCGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCat                                                                                                            >  2:610650/1‑150 (MQ=255)                                          gCAAACCAGTTGTCACCCATCCGCACTTCGTCAATGCCGTTAATCTGGGTGATACGATCCCGCGGCGTATTCAGTGATTCCGCCCCCGGGAAATCGAGTTGGGGGACCCCAACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCa                                                                                                       <  2:583252/151‑1 (MQ=255)                                          gCAAACCAGCTGTCATCCACCCGCCCTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCCCCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCa                                                                                                       <  1:719220/151‑1 (MQ=255)                                          gCAAACCAGCTGTCACCCACCCGCCCTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATCCCGCCCCCTGGAAATCGAGTTTCGGGATCCCCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCa                                                                                                       <  1:81964/151‑1 (MQ=255)                                             aaCCATCTGTCACCCACCCGCCCTTCGTCAATCCCTTTATTCTGCGTTATACGACCCCGCGGCGTATTCAGTGTTCCCGCCCCCGGAAAATCGATTTCCGGGTTCCCCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAgcg                                                                                                    <  2:738843/151‑1 (MQ=255)                                               ccAGCTGTCATCCACCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGATTTTCGGGATCCCCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGca                                                                                                  <  1:1061119/151‑1 (MQ=255)                                                cAGCTGTCACCCACCCCCACTTCGTCAATCCCTTAAATCTGCGTAATACGATCCCGCGGCTTATTCGGTGTTCCCGTCCCCTGAAAATCGATTTTCGGGATCCCCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGcac                                                                                                 <  1:976075/151‑1 (MQ=255)                                                  gCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGcacc                                                                                                >  1:1178497/1‑150 (MQ=255)                                                     ttCACCCACCCGCCCTTCGTCAATGCCTTTATTCTGCGTGATACGACCCCGCGGCTTATTCAGTGATCCCGCCCCCTGAAAATCGAGTTTCGGGATCCCCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCaccacc                                                                                             <  2:323016/149‑1 (MQ=255)                                                      tCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCAccacca                                                                                            >  2:27141/1‑150 (MQ=255)                                                       catccatcCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCaa                                                                                           >  2:765823/1‑150 (MQ=255)                                                            acccGCACTTCTTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGTTTCCGTCCCCGGAAAATCGTGTTTCGGGATCCCCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAaccag                                                                                     <  1:930107/149‑1 (MQ=255)                                                                ccACTTCGTCAATGCCGTTAATCTGCGTGATCCGATCACGCGGCGTATTCAGTGTTCCCGCCCCCTGAAAATCGAGTTTGGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAaccagacca                                                                                 <  1:726050/150‑1 (MQ=255)                                                                  ccTTCGTAAATGCGTTTAATCTGCGTGATACGACCACGCAGCGTATTCAGTGATCCCGTCCCCTGAAAATCGATTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGt                                                                               <  1:50842/150‑1 (MQ=255)                                                                    ttCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGtt                                                                              >  1:256412/1‑150 (MQ=255)                                                                     tCGTCAATGCCTTAAATCTGCGTGATACGATCACGCGGCGTATTCAGTGATTCCGTCCCCTGGAAATCGATTTTCGGGTTCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGt                                                                               <  1:1325198/148‑1 (MQ=255)                                                                     tCGTCAATGCCGTTAATCTGCGTGATCCGATCACGCAGCGTATTCAGTGATCCCGCCCCCGGAAAATCGAGTTTGGGGATCCCCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTa                                                                             <  1:765863/150‑1 (MQ=255)                                                                          aaTGCCGTTAATCTGCGTGATCCGATCCCGCAGCGTATTCAGTGATTCCGCCCCCTGAAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAAc                                                                       <  1:672098/151‑1 (MQ=255)                                                                              ccGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCAc                                                                   >  1:781959/1‑151 (MQ=255)                                                                               cGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCAcc                                                                  >  1:511480/1‑151 (MQ=255)                                                                                  tAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCAccc                                                                 >  2:648067/1‑149 (MQ=255)                                                                                   aaTCGGCGTGATACGATCCCGCAGCGTATTCCGTGATCCCCCCCCCGGGAATTCGAGTTTCGGCACCCCCCCCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGt                                                              <  2:37605/151‑1 (MQ=255)                                                                                       tGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGt                                                          >  1:1085929/1‑151 (MQ=255)                                                                                       ggggTGATACGACCACGGGGCGTATTCAGTGATCCCGCCCCCGGAAAATCGAGTTTCGGGATCCCCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGt                                                          <  1:1160817/148‑1 (MQ=255)                                                                                         cGTGATACGATCACGCAGCGTATTCAGTGATTCCGCCCCCGGAAAATCGAGTTTCGGGATCCCCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCt                                                        <  1:1250609/151‑1 (MQ=255)                                                                                           tGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTc                                                      >  2:685197/1‑151 (MQ=255)                                                                                               aCGATCCCGCGGTGTATTCAGTTATCCCCCCCCCGGGAAATCGGGTTTCGGGACCCCCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACga                                                  <  2:775221/151‑1 (MQ=255)                                                                                                cGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACg                                                   >  1:541673/1‑149 (MQ=255)                                                                                                  aTCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACgaga                                                >  1:1108605/1‑150 (MQ=255)                                                                                                   ccACGCAGCTTATTCAGTGATCCCGTCCCCTGAAAATCGATTTTCGGGATCCCCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGaaa                                              <  1:610675/150‑1 (MQ=255)                                                                                                    cACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGaa                                               >  2:20560/1‑149 (MQ=255)                                                                                                    cACGCAGCGTATTCAGTGATTCCGCCCCCTGGAAATCGAGTTTCGGGATCCCCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAg                                             <  1:45269/151‑1 (MQ=255)                                                                                                    cACGCAGCGTATTAAGTGATCCCGCCCCCGGAAAATGGTGTTTCGGAATCCCCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGaaa                                              <  1:84621/150‑1 (MQ=255)                                                                                                     ccgcgGCGTATTCAGTGATTCCGCCCCCTGGAAATCGAGTTTCGGGATCCCCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGa                                            <  2:5460/150‑1 (MQ=255)                                                                                                       gCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGAt                                           >  1:946579/1‑150 (MQ=255)                                                                                                            gTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTc                                      >  1:1112037/1‑150 (MQ=255)                                                                                                                tCAGTGATTCCGCCCCCTGAAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGATCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACtt                                  <  1:1236215/150‑1 (MQ=255)                                                                                                                 cAGTGATTCCGTCCCCTGGAAATCGATTTTCGGGATCCCCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACttt                                 <  1:942188/150‑1 (MQ=255)                                                                                                                 cAGTGATTCCGCCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTtc                                <  1:459651/151‑1 (MQ=255)                                                                                                                    tGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTcc                             >  2:662232/1‑151 (MQ=255)                                                                                                                        tCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGc                         >  2:671061/1‑151 (MQ=255)                                                                                                                        cccGCCCCCTGGAAATCGATTTTGGGGACCCCCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACg                          <  2:862762/149‑1 (MQ=255)                                                                                                                         ccGCCCCCTGGAAATCGATTTTCGGGACCCCCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGcc                        <  1:359564/151‑1 (MQ=255)                                                                                                                           gTCCCCTGAAAATCGATTTTCGGGATCCCCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCtt                      <  1:126643/151‑1 (MQ=255)                                                                                                                           ccccccGGAAAATCGTGTTTCGGGATCCCCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCtt                      <  1:467115/149‑1 (MQ=255)                                                                                                                             ccccTGAAAATCGAGTTTCGGGATCCCCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCtt                      <  1:458143/149‑1 (MQ=255)                                                                                                                              cccTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCt                   >  1:1288530/1‑151 (MQ=255)                                                                                                                               ccTGAAAATCGATTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTc                  <  1:935132/151‑1 (MQ=255)                                                                                                                                cgggAAATCGATTTCGGGGACCCCCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCg                 <  2:608434/149‑1 (MQ=255)                                                                                                                                cTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCg                 >  2:605832/1‑151 (MQ=255)                                                                                                                                  ggAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGc               >  2:228260/1‑151 (MQ=255)                                                                                                                                       tCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTGc           >  1:406868/1‑150 (MQ=255)                                                                                                                                        cGAGTTTCGGGATCCCCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTGCa          <  1:165666/150‑1 (MQ=255)                                                                                                                                            tttCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTGCAACtgt     >  1:219922/1‑151 (MQ=255)                                                                                                                                                gggATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTGCAACTGTGcc  >  1:959790/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                CCGAAACGCGCCCGACCAGCCCGGTCAACGCCGCCAGACGGGCAAACCAGCTGTCATCCATCCGCACTTCGTCAATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTGCAACTGTGCC  >  NC_000913/3601411‑3601703

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