Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	5,159	0	G	A	40.5%	85.8 / 47.0	74	intergenic (+139/‑75)	thrC/yaaX	L‑threonine synthase/DUF2502 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (43/1);  new base A (1/29);  total (44/30)
Fisher's exact test for biased strand distribution p-value = 2.82e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAAT‑T‑ATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGG  >  NC_000913/5018‑5306                                                                                                                                                |                                                                                                                                                    tAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAAT‑T‑ATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATaacaac                                                                                                                                            >  1:36265/1‑151 (MQ=255)  aaaaTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAAT‑T‑ATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAAcc                                                                                                                                           >  2:913160/1‑151 (MQ=255)         aTTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAAT‑T‑ATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTAAGAGATTAGGATTGCGGGAAATaacaacgccggtt                                                                                                                                     >  1:193403/1‑143 (MQ=255)                   tcAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAAT‑T‑ATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGt                                                                                                                           >  2:222478/1‑150 (MQ=255)                     aaTCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAAT‑T‑ATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAAt                                                                                                                        >  1:1045331/1‑151 (MQ=255)                         aGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAAT‑T‑ATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTcc                                                                                                                    >  1:600492/1‑151 (MQ=255)                            ccGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAAT‑T‑ATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCgg                                                                                                                  >  1:342734/1‑150 (MQ=255)                              gggTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAAT‑T‑ATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGat                                                                                                                >  1:839993/1‑150 (MQ=255)                                   tgtttttATGCAGCCGGGTTTTTTTAAAAAAAAAT‑T‑AGGGAAAAAAAGGAGAGGAAAAAAGAAAAAATCCTAAATAAATGGGTTAACTTAGAGATTAGGTTTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGa                                                                                                           <  1:625643/147‑1 (MQ=255)                                    gCTTTTATGCAGCCCGGCTTTTTTATGAAGAAAT‑T‑ATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGAcc                                                                                                         >  2:179453/1‑151 (MQ=255)                                      ttttATGCAGCCGGTTTTTTTTATAAAAAAAT‑T‑ATGGAAAAAAATGAAAGGAAAAAAGGAAAAATTCTAAAAAAATGCGTTAATTTAGAGATTAGGATGGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGAcc                                                                                                         <  2:986421/149‑1 (MQ=255)                                      ttttATGCAGCCCGGCTTTTTTATGAAGAAAT‑T‑ATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGAccc                                                                                                        >  1:873084/1‑150 (MQ=255)                                          aggccgCCCGTTTTTTTTA‑GAAAAAAT‑TGAGGAAAAAAAAAGACAGGAAAAAAAGAGAATTTCAAAAAAAAGTCGATACATAATAATTTTGGATGGAAGAAAAAAAAACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATaa                                                                                                    <  2:319313/145‑1 (MQ=255)                                          aTGCAGCCCGGTTTTTTTATAAAAAAAT‑A‑ATGAAAAAAAAAGACAGGAAAAAAGGAAAAATTCTAAAAAAAGGCGTTAATTAAGAAATTAGGATTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATaa                                                                                                    <  1:336563/150‑1 (MQ=255)                                             cAGCCCGGCTTTTTTATGAAGAAAT‑T‑ATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACggg                                                                                                >  1:1165977/1‑151 (MQ=255)                                                cccGGCTTTTTTATGAAGAAAT‑T‑ATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCa                                                                                              >  1:1124051/1‑150 (MQ=255)                                                cccAGCTTTTTTAAAAAAAAATAT‑ATGGAAAAAAATAAAGGGAAAAAAAGAGAAATACAAAATAAGTGCGATACCTAAGAGATTAGGTTTGCGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATaa                                                                                                    <  2:255572/145‑1 (MQ=255)                                                   ggtttttttAAAAAAAAAT‑T‑AGGGAAAAAAAGGACAGGGAAAAAGGAAAAATTCAAAAAAAATGCGTTAACTGAGATTTTAGGATTGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACggg                                                                                                <  2:437281/142‑1 (MQ=255)                                                   ggCTTTTTTATGAAGAAAT‑T‑ATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATg                                                                                           >  1:374138/1‑150 (MQ=255)                                                       tttttATGAAGAAAT‑T‑ATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGAACCATAATGGGCAATGATaa                                                                                       >  1:801386/1‑150 (MQ=255)                                                       tttttATAAAAAAAT‑T‑AGGGAAAAAAAGGACAGGAAAAAAGGAAAATTTCAAAATAAAGGGGGTACTTAAGAGATTAGGTTGGCGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATaa                                                                                       <  1:1057803/150‑1 (MQ=255)                                                        ttttAAAAAAAAAT‑T‑AGGAAAAAAAATGAAAGAAAAAAAGGAAAAATTCAAAAAAAATGCGTTAATTTAGATTTGAGGATGGCGAAAAAAAAAACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATaaa                                                                                      <  2:497594/150‑1 (MQ=255)                                                          ttATGAAGAAAT‑T‑ATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAg                                                                                    >  2:1216761/1‑150 (MQ=255)                                                           tATGAAGAAAT‑T‑ATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAgg                                                                                   >  1:877187/1‑150 (MQ=255)                                                            aaaaaaaaaT‑A‑ATGGAAAAAAAGGACAGAAAAAAAGAAAAAATTCAAAATAAAGGCGATAACTAAGAGATTAGGATTGGGGAAAAAAAAACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGa                                                                                  <  2:1053954/144‑1 (MQ=255)                                                            aTGAAGAAAT‑T‑ATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAgg                                                                                   >  1:84420/1‑149 (MQ=255)                                                               aaaaaaT‑T‑ATGGAAAAAAAGGAAAGGAAAAAAGGAAAAATTCTAAATAAATGGGTTAATTTAGAGTTTAGGATTGGGGAAAATAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTa                                                                               <  1:452361/147‑1 (MQ=255)                                                                 gaaAT‑T‑ATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAAcc                                                                            >  1:400875/1‑151 (MQ=255)                                                                           ggAAAAAAAGGACAGGAAAAAAGAAAAAATTCTAAAAAAATGCGTTAACTAAGAGATTAGGATGGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtgtg                                                                        <  1:823297/147‑1 (MQ=255)                                                                             agaAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGaaaa                                                                    >  1:205955/1‑149 (MQ=255)                                                                              gaAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGaaaa                                                                    >  1:654612/1‑148 (MQ=255)                                                                               aaaaaTGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGaaaa                                                                    >  1:417790/1‑147 (MQ=255)                                                                                aaaaaaacagaaaaaaaGGAGAAATTCACAAAAAATGCGTAAACTAAAAGAGTAGAATGGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATg                                                               <  2:342720/139‑1 (MQ=255)                                                                                aaaaTGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATg                                                               >  2:831944/1‑151 (MQ=255)                                                                                aaaaTAAAAGGAAAAAAGAAAAAATACAAAAAAAAGGCGTAAACTTAGAGAGGAGGAGGGGGAAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATg                                                               <  2:41853/151‑1 (MQ=255)                                                                                  aagaaCAGGAAAAAAGAAAAAATTCAAAAAAATTGGGAAAATTTAGATATTGGGATGGGGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCa                                                             <  2:923110/147‑1 (MQ=255)                                                                                     gACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATc                                                          >  1:760410/1‑151 (MQ=255)                                                                                     gACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATc                                                          >  1:3674/1‑151 (MQ=255)                                                                                      aCAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCt                                                         >  1:638020/1‑151 (MQ=255)                                                                                      aCAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCt                                                         >  1:569309/1‑151 (MQ=255)                                                                                        aGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTAt                                                       >  1:188427/1‑151 (MQ=255)                                                                                          ggAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCg                                                     >  2:160278/1‑151 (MQ=255)                                                                                            aaaaaaGAAAAATTCACAAAAAATGCGGTAACTTAGAGTTGAGGATGGGGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTa                                                   <  2:694656/151‑1 (MQ=255)                                                                                            aaaaaGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGt                                                    >  1:449168/1‑150 (MQ=255)                                                                                            aaaaaGGAAAAATTTAAAAAAAAGGGGGTACTTAATAGATTAGGAGTGGGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTAtag                                                     <  1:916866/149‑3 (MQ=255)                                                                                            aaaaaGGAAAAATTCAAAAAAAAGGCGGTAATTGAGAAATAAGGATGGGGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGt                                                    <  2:574686/150‑1 (MQ=255)                                                                                                  gagaAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGc                                              >  2:46004/1‑150 (MQ=255)                                                                                                   aaaatttCTAAATAAAGGCGGTAACTTAGAGATTAGGTTGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGc                                              <  1:729728/144‑1 (MQ=255)                                                                                                   aaaaaTTCTAAAAAAATGCGGTAATTAAGATTTTAGGTTTGCGAAAAAAAACACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCAc                                            <  2:6587/149‑1 (MQ=255)                                                                                                    gaAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACt                                           >  2:1207722/1‑151 (MQ=255)                                                                                                            aaaaaaaGGCTGTAACTAGGAGATTAGGATTAGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTcc                                       <  1:562864/143‑1 (MQ=255)                                                                                                             aaaaaaTGCGGTAACTTAGAGATTAGGATGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTgg                                   <  1:955027/147‑1 (MQ=255)                                                                                                             aaaaaaGGCGGTAATTTAGAGATTAGGATTGGGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCt                                     <  2:737001/145‑1 (MQ=255)                                                                                                             aaaaaaGGCGGTAACTAAGAGATTAGGTTGGCGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCt                                     <  1:913157/145‑1 (MQ=255)                                                                                                              aaaaaTGCGTTACATTAGAGTTTGGGTTTGCGAAAAAAACCACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTTGtt                                 <  2:751374/149‑1 (MQ=255)                                                                                                              aTAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGtt                                 >  1:885890/1‑151 (MQ=255)                                                                                                               tAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTc                                >  1:1166830/1‑151 (MQ=255)                                                                                                               aaaaTGCGTAAACTAAGAGTTTAGGATTGCGAAAAAAAACATCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTc                                <  2:780176/150‑1 (MQ=255)                                                                                                                 aaTGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTg                              >  2:436964/1‑151 (MQ=255)                                                                                                                   tGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGt                            >  2:408320/1‑151 (MQ=255)                                                                                                                         aaCTTAGAGATTAGGTTGGCGGAGAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGt                            <  1:1076156/145‑1 (MQ=255)                                                                                                                         aaCTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTcc                      >  2:101800/1‑151 (MQ=255)                                                                                                                          aCTTAGAGATTACGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCAGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTcc                      >  2:1123063/1‑150 (MQ=255)                                                                                                                            ttAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGCTCTGGTCGCTCCCAt                   >  1:493812/1‑151 (MQ=255)                                                                                                                                agaTTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATgg                 >  2:668786/1‑149 (MQ=255)                                                                                                                                agaTTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgc                >  2:154800/1‑150 (MQ=255)                                                                                                                                agaTGAGGATGGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGTAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgca               <  1:913045/151‑1 (MQ=255)                                                                                                                                   ttAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgcagc             >  1:356768/1‑150 (MQ=255)                                                                                                                                    tgggATGGGGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgcag              <  2:449156/146‑1 (MQ=255)                                                                                                                                    tagggtgggggAAAATACCACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgcagca            <  2:502120/141‑1 (MQ=255)                                                                                                                                       gATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAgg        >  1:969901/1‑151 (MQ=255)                                                                                                                                         ttGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCt      >  1:871552/1‑151 (MQ=255)                                                                                                                                         ttGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCt      >  1:852077/1‑151 (MQ=255)                                                                                                                                            cGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGc       <  1:214441/147‑1 (MQ=255)                                                                                                                                              gagaATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCgg  >  2:287670/1‑150 (MQ=255)                                                                                                                                                |                                                                                                                                                    TAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAAT‑T‑ATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGG  >  NC_000913/5018‑5306

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