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

ATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTG  >  NC_000913/5013‑5303                                                                                                                                                   |                                                                                                                                                 atcaGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATa                                                                                                                                              <  2:13359/151‑1 (MQ=255)   caGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATaac                                                                                                                                            <  2:297081/151‑1 (MQ=255)   caGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATaa                                                                                                                                             >  1:249761/1‑150 (MQ=255)       aaaaTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATaacaac                                                                                                                                         <  2:879989/150‑1 (MQ=255)       aaaaTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAAcc                                                                                                                                        >  2:8770/1‑151 (MQ=255)         aaTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAAccgc                                                                                                                                      >  1:129079/1‑151 (MQ=255)          aTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAAccgcc                                                                                                                                     >  1:658452/1‑151 (MQ=255)               ttCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTtct                                                                                                                                >  2:346479/1‑151 (MQ=255)               ttCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTtct                                                                                                                                >  2:165619/1‑151 (MQ=255)                          aaTCAGGCGGGGTTTGTTTTAACGCAGCCCGGTTTTTTTATGAAAAAATTAGGGAAAAAAAGGACGGGGAAAAAGGAAAAATTCTCAAAAAAGGCGGTACTTTAGAGATTAGGATGGCGGAAAATAACAACCGCCGTTCTCATCGAGTaa                                                                                                                      <  2:119686/150‑1 (MQ=255)                          aaTCAGGCCGGGTTTTTTTTTAGGCCGCCGGGTTTTTTTATAAAAAAATTAGGGAAAAAAATGACAGGAAAAAAGGAAAAATTCTAAAAAAAGGCGTTAATTAAGAAATTAGGATTGCGGAAAAAAACAACCGCCGTTCTCATCg                                                                                                                           <  1:916524/145‑1 (MQ=255)                             cAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCTGAGAATAACAACCGCCGTTCTCATCGAGTAAtctc                                                                                                                  >  2:794855/1‑151 (MQ=255)                               ggcgggttttttttttATGCACCCCGGTTTTTTTATAAAAAAATTATGAAAAAAAAGGACGGGAAAAAAGGAAAATTTCTCAAAAAAGGGGGTATCTTAGAGATAAGTTTTGGGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTcc                                                                                                                 <  2:599380/138‑1 (MQ=255)                               ggCCGTGTTTGTTTTTATACACCCCGGTTTTTTTATAAAAAAATTATGGAAAAAAAGGACAGGAAAAAAGGAGAATTTCTCAATAAATGCGGTATTTAAGAGATGAGGATGGCGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCg                                                                                                                <  1:421176/151‑1 (MQ=255)                               ggCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCg                                                                                                                >  1:201536/1‑151 (MQ=255)                                  cGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGat                                                                                                             >  1:173949/1‑151 (MQ=255)                                      tttGTTTTTATGCAGCCCGGTTTTTTTAGAAAAAAATTATGGAAAAAAAGGAGAGGAAAAAAGGAAAAATTCTAAAAAATTGGGGTATTTTAGAGATTAGGTTGGGGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGatat                                                                                                           <  1:610359/149‑1 (MQ=255)                                         gCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGAcc                                                                                                      >  2:62282/1‑151 (MQ=255)                                          cTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGCAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGAATGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGAccc                                                                                                     >  1:716706/1‑151 (MQ=255)                                           ttttATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGAccc                                                                                                     >  1:1001500/1‑150 (MQ=255)                                            tttATCCACCCCGTTTTTTTAATAAAAAAATTATGAAAAAAAAGGACAGGAAAAAAGGAAAAATTCTAAAAAAAGGCGGTAACTAAGAGATTAGGATTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCa                                                                                                    <  2:96176/150‑1 (MQ=255)                                             ttATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCAt                                                                                                   >  1:844501/1‑150 (MQ=255)                                               acgccccccgttttttttAAGAAAAAATTGGGAAGAAAAAGAAGAGGAAAAAAGAAAAAATTCTAAATAAATGCATTACCTTAGAATTTAGGATTGCGGAAAAAACAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAAc                                                                                                <  2:710021/139‑1 (MQ=255)                                                 gCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACgg                                                                                              >  2:804815/1‑151 (MQ=255)                                                          cTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGAt                                                                                      >  2:133955/1‑150 (MQ=255)                                                           ttttttATGAAAAAATTAGGGAAAAAAAGGACAGGGAAAAAGAAAAAATCCTCAAAAAATGGGGTATTTGAGAGATTGGGGTGGCGGAGAAAAACAGCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATa                                                                                     <  1:990723/150‑1 (MQ=255)                                                              tttATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATaaa                                                                                   >  2:531390/1‑149 (MQ=255)                                                                tATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAgg                                                                                >  2:383136/1‑150 (MQ=255)                                                                 aTGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAgg                                                                                >  1:61793/1‑149 (MQ=255)                                                                    aagaaATTATGGAAAAAAATGACAGGAAAAAAGGAAAAATTCTAAAAAAATGCGTTAATTAAGAGTTTAGGATTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGt                                                                             <  1:472759/149‑1 (MQ=255)                                                                           tATGGAAAAAAATGACAGGAAAAAAGGAAAATTTATAAAAAAAGGCGGTAATTTAGAGATTAGGATTGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtgtg                                                                     <  2:206056/150‑1 (MQ=255)                                                                              gaaaaaaaaTGAAAGGAAAAAAGGAAAAATTCTAAAAAAATGGGTTAACTTAGAGATTAGGATTGGGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGa                                                                    <  2:701436/144‑1 (MQ=255)                                                                                aaaaaaaTGAAGGGAAAAAAGGAAAAATTCTAAATAAAGGCGGTAACTTAGAGATTAGGATGGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtgtg                                                                     <  2:239186/143‑1 (MQ=255)                                                                                aaaaaaaTGAAGGGAAAAAAGAAAAAATTCTCAAAAAAGGCGTTAACTTAGAAATTAGATTTGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAg                                                               <  2:941699/149‑1 (MQ=255)                                                                                  aaaaaTAAGAGGAAAAAAGAAAAAATTCAAAAAAAAGGCGGTATTTAAGAGATTAGGATGGGGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGAt                                                             <  2:77265/151‑1 (MQ=255)                                                                                   aaaaTGACAGGAAAAAAGGAAAAATTCAAAAAAAATGCGGTAACTTAGAGATTAGGATGGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATg                                                            <  2:249821/151‑1 (MQ=255)                                                                                   aaaaGGACAGGAAAAAAGGAAAAATTCTAAATAAAGGCGGTAATTTAGAGATGAGGATGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATg                                                            <  2:525406/151‑1 (MQ=255)                                                                                   aaaaGGACAGGAAAAAAGGAAAAATTCTAAAAAAAGGCGGTAACTAAGAGATGAGGTTTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATg                                                            <  2:853486/151‑1 (MQ=255)                                                                                        gACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATc                                                       >  2:421341/1‑151 (MQ=255)                                                                                            ggaaaaaaGGAAAAATTCTAAAAAAATGCGGTAACTAAGAGATTAGGATGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCt                                                      <  2:456257/145‑1 (MQ=255)                                                                                               aaaaaGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGt                                                 >  2:425170/1‑150 (MQ=255)                                                                                               aaaaaGGAAAAATTCTAAAAAAATGCGTTAACTTAGAAATTAGGATTGGGAAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATc                                                   <  1:795988/148‑1 (MQ=255)                                                                                               aaaaaGGAAAAATTCTAAAAAAAGGCGGTAATTTAGAGATAAGGATGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCg                                                  <  2:196373/149‑1 (MQ=255)                                                                                                 aaaGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACt                                              >  2:88247/1‑151 (MQ=255)                                                                                                  aaGGAGAAATTCTCAAAAAATGCGATAACTTAGAGATGAGGATCGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTc                                             <  2:276254/151‑1 (MQ=255)                                                                                                  aaGGAAAAATTCTCAAAAAAGGCGGAACTTTAGAAATGAGGATTGCGAAAAATAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTc                                             <  1:547432/151‑1 (MQ=255)                                                                                                   aGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCg                                            >  2:401417/1‑151 (MQ=255)                                                                                                    ggAAAAATTCTCAATAAAGGCGTTAACTTAGAGATGAGGATGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCg                                            <  1:353147/150‑1 (MQ=255)                                                                                                        aaaTTCTAAAAAAGTGGGGTAACTAAGAGATTAGTATTGGGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACtt                                       <  2:124913/151‑1 (MQ=255)                                                                                                        aaaTTCAAAAAAAATGAGATAATTAAGATATTGGGATTGGGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACtt                                       <  1:165567/151‑1 (MQ=255)                                                                                                          aTTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTc                                     >  2:895611/1‑151 (MQ=255)                                                                                                           ttCTAAAAAAATGGGGTAATTTAGAGATTAGGATGGGGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTc                                     <  2:129127/150‑1 (MQ=255)                                                                                                                 aTAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGtt                              >  2:451458/1‑151 (MQ=255)                                                                                                                   aaaTGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCt                            >  1:782685/1‑151 (MQ=255)                                                                                                                    aaTGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTg                           >  2:970374/1‑151 (MQ=255)                                                                                                                    aaTGCGGTAACTTAGAGATTAGGATTGCGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTg                           <  2:201438/151‑1 (MQ=255)                                                                                                                     aTGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTgg                          >  2:146826/1‑151 (MQ=255)                                                                                                                            aattaagagATTAGGATTGTGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCt                     <  1:794739/144‑1 (MQ=255)                                                                                                                             aCTTAGAGATTAGGATTGCGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTccc                  <  2:113080/151‑1 (MQ=255)                                                                                                                                tAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATg               >  1:422933/1‑151 (MQ=255)                                                                                                                                         ggATTGCGAAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGcaca       <  1:682700/150‑1 (MQ=255)                                                                                                                                             tGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGCCAGCACAGGCTg  >  1:491065/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                                 ATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTG  >  NC_000913/5013‑5303

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: ‑