Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,133,175	G→A	T13I (ACC→ATC)	yghR ←	putative ATP‑binding protein YghR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,133,175	0	G	A	96.7%	193.5 / NA	61	T13I (ACC→ATC)	yghR	putative ATP‑binding protein YghR
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (23/36);  minor base C (2/0);  total (25/36)
Fisher's exact test for biased strand distribution p-value = 1.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGGTTTTCATGCTTTGGGTTTTCGAGGATTTGGCCGCCAGTCGCCGTTCCAGCCAGACGCCAACCAACGGCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAAC‑GGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGT  >  NC_000913/3132932‑3133418                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     ttGGTTTTCAGTCTGTGGTTTTGCGAGGATGTGGCCTCCAGGCGCCGTGCCAGCCAGACGCCAACCAACGGCAATCGTGTGATGTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTGTGCCAGTGTTGGTGCAGCGATTTCACCAGGGCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGTCCTCCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTAc                                                                                                                                                                                                                                                <  1:421339/250‑1 (MQ=255)        tCATGCTTTGGGTTTTCGAGGATTTGGCCGCCAGTCGCCGTTCCAGCCAGACGCCAACCAACGGCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAAc                                                                                                                                                                                                                                         <  1:356880/250‑1 (MQ=255)          aTGCTTTGGGTTTTCGGTGATTTGGCCGCCAGTCGCCGTTCCAGCCAGACGCCAACCAACTGCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCGCCAGATAGCGCCGCTCGGTTTGCCAGGGTTGTTGCAGCGATTTCACCAGGGCGGGGGGCAGTGTGGAGTTACCTGGGCCATCACACCCGACCACCTCAATCAGCCACGGAATGTAGTTGGGCGGCGGCGGGATTTTACTATTAACagt                                                                                                                                                                                                                                      <  2:221962/251‑1 (MQ=255)              tttGGGTTTTCGAGGATTTGGCCGCCAGTCGCCGTTCCAGCCAGACGCCAACCAACGGCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACACCCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACagtttgc                                                                                                                                                                                                                                  >  1:268596/1‑250 (MQ=255)               ttGGGTTTTCGAGGATTTGGCCGCCAGTCGCCGTTCCAGCCAGACGCCAACCAACGGCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACagtttg                                                                                                                                                                                                                                   >  1:238405/1‑249 (MQ=255)                 gggTTTTCGAGGATTTGGCCGCCAGTCGCCGTTCCAGCCAGACGCCAACCAACGGCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACagtttgagtt                                                                                                                                                                                                                               >  1:110694/1‑251 (MQ=255)                   gTTTTCGAGGATTTGGCCGCCAGTCGCCGTTCCAGCCAGACGCCAACCAACGGCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACagtttgagtttg                                                                                                                                                                                                                             >  1:347168/1‑251 (MQ=255)                    ttttCGCGGACTTGGCCGCCAGTCGCCGTTCCAGCCAGACGCCAACCAACGGAAATCGTTTGATTTTGTCGCAGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGCTTGTACTATTAAAAGTTTGAGTTTGt                                                                                                                                                                                                                            >  2:432997/1‑251 (MQ=255)                        cGAGGATTTGGCCGCCAGTCGCCGTTCCAGCCAGACGCCAACCAACGGCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTg                                                                                                                                                                                                                        <  1:307110/251‑1 (MQ=255)                                tGGCCGCCAGTCGCCGTTCCAGCCAGACGCCAACCAACGGCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGGGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTGGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCAt                                                                                                                                                                                                                 <  1:286439/250‑1 (MQ=255)                                 ggCCGCCAGTCGCCGTTCCAGCCAGACGCCAACCAACGGCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTc                                                                                                                                                                                                               <  1:437459/251‑1 (MQ=255)                                  gccgccAGTCGCCGTTCCAGCCAGTCGCCAACCAACGGCAATCGTTTGATTTTGGCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCTGAATGTAGTTGGGCTGCGGCGCGAGTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCg                                                                                                                                                                                                              <  2:157016/251‑1 (MQ=255)                                   ccgccAGTCGCCGTTCCAGCCAGCCGCCAACCAACGGCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGt                                                                                                                                                                                                             >  1:122283/1‑251 (MQ=255)                                     gccAGTCGCCGTTCCAGCCAGACGCCAACCAACGGCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGtt                                                                                                                                                                                                            >  2:289160/1‑250 (MQ=255)                                       cAGTCGCCGTTCCAGCCAGACGCCAACCAACGGCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTaa                                                                                                                                                                                                          >  1:197858/1‑250 (MQ=255)                                        aGTCGCCGTTCCAGCCAGACGCCAACCAACGGCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAat                                                                                                                                                                                                        <  1:27332/251‑1 (MQ=255)                                            gCCGTTCCCGCCCGACGCCAACCAACGGCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAatcatc                                                                                                                                                                                                    >  2:378791/1‑251 (MQ=255)                                            gCCGTTCCAGCCAGACGCCAACCAACGGCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTTGGCTGCGGCGCGATTGTa                                                                                                                                                                                                                                                 <  2:427017/206‑1 (MQ=255)                                            gCCGTTCCAGCCAGACGCCAACCAACGGCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTa                                                                                                                                                                                                                                                 >  1:427017/1‑206 (MQ=255)                                               gTTCCAGCCAGACGCCAACCAACGGCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATcccc                                                                                                                                                                                                 >  2:86820/1‑251 (MQ=255)                                                      ccagACGCCAACCAACGGCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCAt                                                                                                                                                                                          >  2:5927/1‑251 (MQ=255)                                                            gCCAACCAACGGCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGt                                                                                                                                                                                    <  2:43667/251‑1 (MQ=255)                                                              caaccaacGGCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGa                                                                                                                                                                                  <  1:49639/251‑1 (MQ=255)                                                               aaccaacGGCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGTTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGAt                                                                                                                                                                                 <  1:245804/251‑1 (MQ=255)                                                                accaacGGCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATc                                                                                                                                                                                >  1:222912/1‑251 (MQ=255)                                                                   aacGGCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAgg                                                                                                                                                                             <  1:129799/251‑1 (MQ=255)                                                                     cGGCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGGTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTATTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTAGCCATTTTGGTGATCAGGaa                                                                                                                                                                           <  2:162828/251‑1 (MQ=255)                                                                      ggCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGaaa                                                                                                                                                                          <  2:347168/251‑1 (MQ=255)                                                                      ggCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGaa                                                                                                                                                                           >  1:370385/1‑250 (MQ=255)                                                                        cAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATCGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCAATCAAACCCGACCACCGCAATCAGCCCCGGACTGTAGTTGGGCTGCGGCGCGCTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATACCCTATCCATTTTGGTGATCAGGaacgc                                                                                                                                                                        >  2:174335/1‑248 (MQ=255)                                                                        cAATCGTTTGAGTTTGGCGCCGTCTTCGCCTGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAgc                                                                                                                                                                        <  1:356151/251‑1 (MQ=255)                                                                                   ttttGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAAt                                                                                                                                                             >  2:389603/1‑251 (MQ=255)                                                                                        tCGCAGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTgagc                                                                                                                                                        >  2:322303/1‑248 (MQ=255)                                                                                         cGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCATCGATTTCACCAGGTCGGTGGTCAGTGTTGATTTACCGTTGCCATCACACCCTACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTgcgc                                                                                                                                                        <  1:420793/250‑1 (MQ=255)                                                                                           ccGTCTTCGACGGAGAGCAACCCCAGAAAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTCCTATTAACAGTTTGAGTTTGTAGTGTATCCATTCGTTAAATCATCCCCTATCCATTTTGATGATCAGGAAAGCGCGGCGTGAATTGCGCGt                                                                                                                                                      >  2:305539/1‑250 (MQ=255)                                                                                                 tCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCCCCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGa                                                                                                                                                 <  2:110694/249‑1 (MQ=255)                                                                                                 tCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGAc                                                                                                                                                <  1:86820/250‑1 (MQ=255)                                                                                                    ccGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTCGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCtt                                                                                                                                            >  2:453873/1‑251 (MQ=255)                                                                                                            cagcCCCAGATAGAGCCGCTCGGGTTGCCAGTGTTGGTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCAGTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAATCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAAtt                                                                                                                                    <  2:184814/251‑1 (MQ=255)                                                                                                                   aGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGt                                                                                                                                                      <  1:174149/226‑1 (MQ=255)                                                                                                                   aGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGt                                                                                                                                                      >  2:174149/1‑226 (MQ=255)                                                                                                                                              tgttgCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAgg                                                                                                   <  1:391106/250‑1 (MQ=255)                                                                                                                                               gttgCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGa                                                                                                  <  1:351520/250‑1 (MQ=255)                                                                                                                                               gttgCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCAGTTCAAGGa                                                                                                  >  1:213649/1‑250 (MQ=255)                                                                                                                                                       gATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCCTCCAGTTCACAGATCTTTgcgc                                                                                         >  2:256567/1‑251 (MQ=255)                                                                                                                                                        aTTTCACCAGGTCGGAGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCCCCGCAATCAGCCCCGGATTGTAGTTGGACTGCGTCGCAACTGTACTAGTCACAGTTTGAGTTTGTAGTGTATCCATTCGTTACATAATCCCCTATACATTTTGGTGAGCAGGCAAGCGCGGTGCGCATTGTACGTAGTGACGCTTGCAGAATTTCATCGGCCGGATGTCTccccacccgttcaatgaatttttcgcc                                                                                        >  2:465576/1‑227 (MQ=255)                                                                                                                                                              ccAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGaa                                                                                  <  1:267123/251‑1 (MQ=255)                                                                                                                                                                                gggATTTACCGGTGCCCTCCCACCCGACCACCGCAATCAGCCTCGGAATGTAGTGGGGCTGCTGCTCGATTGTACTATTAACAGTTTGAGGTTGTAGTGCATCCATTCGTTAAATCATCCCCTAACCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCATAATTTCATCGGCGTGGTTTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCATTTGCGGCGTAAc                                                                <  2:405739/250‑1 (MQ=255)                                                                                                                                                                                 ggATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAAc                                                                <  1:174335/250‑1 (MQ=255)                                                                                                                                                                                 ggATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAAc                                                                <  2:370385/250‑1 (MQ=255)                                                                                                                                                                                   atgtACCGGTGCCATCCCACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCAGTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTTACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTTCGCCATTGAATGTCAGTTGCGGCGTAAC‑GGc                                                            <  2:312666/248‑1 (MQ=255)                                                                                                                                                                                        ccGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAAC‑GGCGAttt                                                       <  1:256567/251‑1 (MQ=255)                                                                                                                                                                                                  cacaCCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAAC‑GGCGATTTTTTCCTGTa                                              <  1:377387/250‑1 (MQ=255)                                                                                                                                                                                                          accaccGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCCGTTCAAGGATCTTTTCGCCCTTGTATTTCATGTGTGGCGTAAC‑GGCGGTTTTTTTCTTTTACGCTGccc                                     >  2:8001/1‑250 (MQ=255)                                                                                                                                                                                                                gCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCCGCCATCCCGGTCAAAGATCTTTGCGCCCTTGTATGTAAGGTGCGG‑GTGACAGGCGATTTTTTTCTTGAACGCTgaccggtggg                               >  2:247563/1‑242 (MQ=255)                                                                                                                                                                                                                         cccGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAAC‑GGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGtt                      <  1:5927/251‑1 (MQ=255)                                                                                                                                                                                                                            ggAATGTAGTGGGTCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGAGCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAAC‑GGCGATTTTTTCCTGTAACGCTGCCAGTGGGTGGTCAGGTTTAc                   <  2:268596/251‑1 (MQ=255)                                                                                                                                                                                                                               aTGTCGTTGGGCTGCGGCGCGAGTGTACTAGTCACAGGTTGAGTGTGTAGTGCATCCATTCGTTAAAGCATCCCCTGTCCCTTTTGGGGATCAGGCAATCGCGGCGTTAATTGCGCGTAGTTACGCTTGCAGAATTTCATCGGTGGGGTGGCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAAC‑GGCGATTTTTTCCTGTAACGCTGCCAGTTGGTCGGCATGTTTAcgag                <  2:248943/251‑3 (MQ=255)                                                                                                                                                                                                                                    gTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAAC‑GGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAAc           <  1:432997/251‑1 (MQ=255)                                                                                                                                                                                                                                            gcggcgCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAAC‑GGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGtt     <  2:238405/249‑1 (MQ=255)                                                                                                                                                                                                                                             cggcgCGATTGTACTATTAACAGTGTGAGTTTGTAGTGCATCCAGTCGTTAAATCATCCCCTAGCCATTTTGGGGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAAC‑GGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGt  <  2:122283/251‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     TTGGTTTTCATGCTTTGGGTTTTCGAGGATTTGGCCGCCAGTCGCCGTTCCAGCCAGACGCCAACCAACGGCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAAC‑GGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGT  >  NC_000913/3132932‑3133418

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