Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	501,564	0	A	C	45.3%	97.9 / 22.8	64	*559E (TAA→GAA)	ybaL	inner membrane putative NAD(P)‑binding transporter
Reads supporting (aligned to +/- strand):  ref base A (31/4);  new base C (0/29);  total (31/33)
Fisher's exact test for biased strand distribution p-value = 2.98e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTC  >  NC_000913/501415‑501711                                                                                                                                                      |                                                                                                                                                    tACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTa                                                                                                                                                     <  1:35131/150‑1 (MQ=255) tACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTa                                                                                                                                                     <  2:608004/150‑1 (MQ=255)    tGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTAcccc                                                                                                                                                 >  1:84990/1‑151 (MQ=255)       gATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTc                                                                                                                                              <  1:19520/151‑1 (MQ=255)            ttAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTcaccac                                                                                                                                         >  2:399372/1‑151 (MQ=255)                    tcgtcgGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCAccc                                                                                                                                  <  2:330018/150‑1 (MQ=255)                    tcgtcgGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACccg                                                                                                                                 >  2:147551/1‑151 (MQ=255)                        cgtttCGTAGCCCAGATAAGGCGTTCCCGCCGCATCTGGCATTTGGCTCTCGATGCTTGATGCGACGCTGGCGCGTCTTATCATGCCTAAATATTTTTCATATTTTACACCCGGCAACCCCCGTTTCCCCCGTCCCCCCCTCTCccgccg                                                                                                                              <  2:610874/147‑1 (MQ=37)                         ggTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCgg                                                                                                                             >  2:698932/1‑150 (MQ=255)                         ggTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGt                                                                                                                            >  1:838469/1‑151 (MQ=255)                                 ggCCAGATAGGGCTTTCCCCCCGCATCTGGCTTTTGTCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATTTTTTTCATTTTTTACACCCGGCACCCCCCGTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGttt                                                                                                                     <  2:747183/150‑1 (MQ=25)                                   ccaaataggcctttCCCGCCGCATCTGGCATTTGGCTCTCGATGCCTGAGGCGCCGCGGGCGCGTTTTTTCATGCCTCCTTTTTTTTCATTTTTTCCACCCGGCACCCCCCTTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTcc                                                                                                                   <  2:365383/138‑1 (MQ=255)                                   ccAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTcc                                                                                                                   >  1:229529/1‑150 (MQ=255)                                    cAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTAAATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTcc                                                                                                                   >  2:768512/1‑149 (MQ=255)                                     aGATAAGGCGTTCCCGCCGCATCTGGCTTTTGGCTCTCGATGCCTGATGCGCCGCTGGCGCGTCTTATCATGCCTACATATTTTTCATTTTTTACACCCGGCAACCACCGTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCcagc                                                                                                                <  1:906009/151‑1 (MQ=38)                                             cGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTcc                                                                                                          >  2:599087/1‑149 (MQ=255)                                             cGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTcc                                                                                                          >  1:570879/1‑149 (MQ=255)                                              gTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGc                                                                                                       >  1:1075728/1‑151 (MQ=255)                                               ttCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGc                                                                                                       >  2:953234/1‑150 (MQ=255)                                                 cccgccgcATCTGGCATTTGGCTCTCGATGCCTGATGCGCCGCTGGCGCGTCTTATCATGCCTACATTTTTTTCATATTTTACACCCGGCAACCCCCGTTTCCCCCGTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATg                                                                                                    <  2:331332/149‑1 (MQ=38)                                                 cccgccgcATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTTTTATCATGCCTACATTTTTTTCATTTTTTACACCCGGCACCCCCCTTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATg                                                                                                    <  1:941741/149‑1 (MQ=37)                                                     ccgcATCTGGCTTTTGCCTCTCGATGCCTGATGCGACGCTGGCGTGTCTTATCATGCCTACATATTTTTAATTTTTTACACCGGGCAACCCCCTTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTAc                                                                                                <  1:74957/151‑1 (MQ=37)                                                        cATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACggg                                                                                             >  2:136785/1‑151 (MQ=255)                                                        cATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACggg                                                                                             >  2:987228/1‑151 (MQ=255)                                                                 tttGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATTTTTTTCATTTTTTACACCCGGAACCCCCCTTTTCCCCCGTCACCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCa                                                                                     <  1:891866/150‑1 (MQ=37)                                                                     gCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTc                                                                                >  2:688373/1‑151 (MQ=255)                                                                       tctggctGCTTGATGCGAGGCGGGCGCGTCTTTTCATCCCTACATTTTTTTCATATTTTACATCCGGAACCCACCGTTTCCCCCGTCCCCACCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTc                                                                                <  2:97189/143‑1 (MQ=255)                                                                        ctcGATGCCTGATGCGCCGCGGGCGCGTCTTATCATCCCTACATATTTTTCATATTTTACACCCGGCAACCCCCGTTTCCCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTc                                                                                <  1:530981/148‑1 (MQ=255)                                                                          cGATGCCTGGTGCGCCGCGGGCGCGTTTTATCATCCCTACTTTTTTTTCATTTTTTACACCCGGCACCCCCCTTTTCCCCCTCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTc                                                                                <  2:20321/146‑1 (MQ=255)                                                                          cGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGccc                                                                            >  1:632868/1‑150 (MQ=255)                                                                           gATGCCTGAGGCGACGCTGGCGCTTTTTATCATGCCTACATATTTTCCATTTTTTACACCCGGCACCCCCCTTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCAt                                                                          <  1:1036479/151‑1 (MQ=255)                                                                               ccTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACt                                                                      >  2:461971/1‑151 (MQ=255)                                                                                 tGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTAc                                                                    >  1:325765/1‑151 (MQ=255)                                                                                         cgcgGGCGCGTTTTATCATGCCTACATATTTTTCATATTTTACACCCGGCAACCCCCTTTTCCCCCGTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGAtt                                                              <  1:837088/149‑1 (MQ=255)                                                                                         cGCTGGCGCGTTTTATCATGCCTACATTTTTTTCATATTTTACACCCGGCACCCCCCTTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTc                                                             <  1:21485/150‑1 (MQ=255)                                                                                             ggCGCGTCTTATCATGCCTACATTTTTTTCATATTTTACACCCGGCACCCCCCTTTTCCCCCGTCACCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGcac                                                         <  1:876536/150‑1 (MQ=255)                                                                                               cgcgTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGcaccac                                                      >  1:365274/1‑151 (MQ=255)                                                                                                     ttATCATCCCTACATTTTTTTCATTTTTTACATCCGCCACCCCCCGTTTCCCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCg                                                 <  2:545818/150‑1 (MQ=255)                                                                                                         cATGCCTACTTATTTTTCTTTTTTTCCACCCGCCACCCTCTTTTTCCCCCCTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGa                                             <  2:699676/150‑1 (MQ=255)                                                                                                          aTGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGat                                            >  1:759885/1‑150 (MQ=255)                                                                                                                aCATTTTTTTCATTTTTTACACCCGGCAACCCCCTTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGc                                      <  1:401217/150‑1 (MQ=255)                                                                                                                  atatTTTTCATATTTTACATCCGGCACCCCCCGTTTACCCCGCCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCAc                                   <  1:61500/151‑1 (MQ=255)                                                                                                                   tttttttatttttttCCACCCGGCACCCCCCGTTTACCCCCCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGcc                                     <  2:388178/137‑1 (MQ=255)                                                                                                                   tatTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACt                                  >  1:999624/1‑151 (MQ=255)                                                                                                                   tatTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACt                                  >  2:49619/1‑151 (MQ=255)                                                                                                                     tttttCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTc                                >  2:732360/1‑151 (MQ=255)                                                                                                                      ttttAATATTTTACACCCGCCAACCCCCTTTTCCCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTc                                <  2:873628/150‑1 (MQ=255)                                                                                                                       tttCATATTTTACACCCGGCACCCCCCGTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCa                               <  1:222217/150‑1 (MQ=255)                                                                                                                          cATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCg                            >  2:627910/1‑150 (MQ=255)                                                                                                                            ttttttccccccagCAACCCCCTTTTCCCCCGCCCCCCCCACCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCg                            <  2:749261/135‑1 (MQ=255)                                                                                                                            ttttttACACCCGGCACCCCCCGTTTCCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATc                             <  1:910972/145‑1 (MQ=255)                                                                                                                              ttttACACCCGCCACCCCCCTTTCCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATc                             <  2:767987/145‑1 (MQ=255)                                                                                                                                ttACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATaa                      >  1:1040817/1‑150 (MQ=255)                                                                                                                                   cacccGGCACCCCCCGTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATgg                   <  2:990925/147‑1 (MQ=255)                                                                                                                                    aTCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATgg                   >  2:389723/1‑149 (MQ=255)                                                                                                                                    aTCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGgc                 >  1:227382/1‑151 (MQ=255)                                                                                                                                             cccccttttCCCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaa         <  2:476321/144‑1 (MQ=255)                                                                                                                                              ccccGTTTCCCCCGTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaata       <  1:970544/149‑1 (MQ=255)                                                                                                                                              ccccGTTTCCCCCGCCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaata       <  2:420476/149‑1 (MQ=255)                                                                                                                                              caccGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaat        >  1:378835/1‑150 (MQ=255)                                                                                                                                               accGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaata       >  2:685561/1‑150 (MQ=255)                                                                                                                                               accGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaata       >  1:77160/1‑150 (MQ=255)                                                                                                                                                 cttttCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaataat     <  2:998444/148‑1 (MQ=255)                                                                                                                                                  gTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAAtct   >  2:183038/1‑151 (MQ=255)                                                                                                                                                   tttCCCCCGTCCCCCCCCCTCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAAtc    <  2:957873/149‑1 (MQ=255)                                                                                                                                                   tttACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAAtctc  >  1:657433/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                    TACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTC  >  NC_000913/501415‑501711

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