Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,281,093	G→A	A372T (GCG→ACG)	yjcE →	putative transporter YjcE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,281,093	0	G	A	98.2%	167.5 / ‑4.4	56	A372T (GCG→ACG)	yjcE	putative transporter YjcE
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (28/27);  total (28/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.11e-01

AGGTCTGCAGCTGCCGGGTATTCTGGAGACGTCGCTGATGGCAGCAGAAATCGACCCTAACGTCGAAATCTGGATGCTGTTTACCAATATTATTCTGATATATGCGGCGCTGATGCTGGTCCGTTTCGGCTGGCTGTGGACGATGAAAAAGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCGTCATTCTCTTCTCGCTGTTTGTCGGCGTGGTGATGTTGCCTATTCTGCTACAACACATTGAAGTCGCGGATCATTCGCAACAATTGAAAGAGGAACGTATTGCGCGAGCGGCAACGGCAGAAGTGGCGATTGTGGCGATCCAGAAAATGG  >  NC_000913/4280849‑4281339                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       aGGTCTGCAGCGGCCGGGTATTCTGGAGAAGTCGCTGATGGCAGCAGAAATCGACCCTAACGTCGAAATCTGGATGCTGTTTACCAATATTATTCTGATATATGCGGCGCTGATGCTGGTCCGTTTCGGCTGGCTGTGGACGATGAAAAAGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGACGATc                                                                                                                                                                                                                                                   <  2:254147/250‑1 (MQ=255)   gTCTGCAGCTGCCGGGTATTCTGGAGACGTCGCTGATGGCAGCAGAAATCGACCCTAACGTCGAAATCTGGATGCTGTTTACCAATATTATTCTGATATATGCGGCGCGGATGCTGGTCCGTTTCGGCTGGCTGTGGACGATGAAAAAGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGACGATCAct                                                                                                                                                                                                                                                <  2:371883/251‑1 (MQ=255)       gcatcagccgtgaaTTCTGGAGACGACGCTGATGGCAGCAGAAATCGACCCTAACGGCGAAATCTGGAGGCTGTTTACCAATATTATTCTGATATATGCGGCGCTGATGCTGGTCCGTGTCGGCGGGCTGTGGACGATGAAAAAGGTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGCGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGACGATCCctc                                                                                                                                                                                                                                               <  2:90305/235‑1 (MQ=255)               cGGGTATTCTGGAGACGTCGCTGATGGCAGCAGAAATCGACCCTAACGTCGAAATCTGGATGCTGTTTACCAATATTATTCTGATATATGCGGCGCTGATGCTGGTCCGTTTCGGCTGGCTGTGGACGATGAAAAAGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGACGATCACTCTGGCCGGtgtg                                                                                                                                                                                                                                    <  1:245494/251‑1 (MQ=255)                          gagaCGTCGCTGCTGGCAGCAGAAATCGACCCTAACGTCGAAATCTGGATGCTGTTTACCAATATTATTCTGATATATGCGGCGCTGATGCTGGTCCGTTTCGGCTGGCTGTGGACGATGAAAAAGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGACGGGGTGCGTGGGACGATCACTCTGGCCGGTGTGCTCTCTATccc                                                                                                                                                                                                                         >  1:288889/1‑251 (MQ=255)                                      aTGGCAGCAGAAATCGACCCTAACGTCGAAATCTGGATGCTGTTTACCAATATTATTCTGATATATGCGGCGCTGATGCTGGTCCGTTTCGGCTGGCTGTGGACGATGAAAAAGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGACGATCACTCTGGCCGGTGAGCTCTCTATCCCGCTGCTCTTGcc                                                                                                                                                                                                             >  1:243112/1‑251 (MQ=255)                                      aTGGCAGCAGAAATCGACCCTAACGTCGAAATCTGGATGCTGTTTACCAATATTATTCTGATATATGCGGCGCTGATGCTGGTCCGTTTCGGCTGGCTGTGGACGATGAAAAAGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGACGATAACTCTGGCAGGTGTGCTCTCTATCCCGCTGCTCTTGcc                                                                                                                                                                                                             >  1:121131/1‑251 (MQ=255)                                       tGGCAGCAGAAATCGACCCTAACGTCGAAATCTGGATGCTGTTTACCAATATTATTCTGATATATGCGGCGCTGATGCTGGTCCGTTTCGGCTGGCTGTGGACGATGCAAAAGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGACGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCg                                                                                                                                                                                                            >  1:71411/1‑251 (MQ=255)                                        ggCAGCAGAAAGCGACCCTAACGTCGAAATCTGGACGCTGTTTACCAATATTATTCTGATATATGCGGCTCTGATGCTGGTCCGTTTCGGCTGGCTGTGGACGATGAAAAAGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGACGATCACTCTGGCCGTTGTGCTCTCTATCCCGCTGCTCTTGCCg                                                                                                                                                                                                            <  2:383400/250‑1 (MQ=255)                                             cagAAATCGACCCTAACGTCGAAATCTGGATGCTGTTTACCAATATTATTCTGATATATGCGGCGCTGATGCTGGTCCGTTTCGGCTGGCTGTGGACGATGAAAAAGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGACGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGt                                                                                                                                                                                                      >  1:390163/1‑251 (MQ=255)                                                     gACCCTAACGTCGAAATCTGGTTGCTGTTTACCAATATTATTCTGATAAATGCGGCGCTGATGCTGGTCCGTTTCGGCTGGCTGTGGACGATGAAAAAGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGACGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCtt                                                                                                                                                                                              <  2:375993/251‑1 (MQ=255)                                                      aCCCTAACGTCGAAATCTGGATGCTGTTTACTAATATTATTCTGATATATGCGGCGCTGATGCTGGTCCGTTTCGGCTGGCTGTGGACGATGAAAAAGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGACGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTc                                                                                                                                                                                             <  1:308107/251‑1 (MQ=255)                                                      aCCCTAACGTCGAAATCTGGATGCTGTTTACCAATATTATTCTGATATATGCGGCGCTGATGCTGGTCCGTTTCGTCTTGCTGTGGACGATGAAAAAGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGACGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCt                                                                                                                                                                                               <  2:254963/249‑1 (MQ=255)                                                         cTAACGTCGAAATCTGGATGCTGTTTACCAATATTATTCTGATATATGCGGCGCTGATGCTGGTCCGTTTCGGCTGGCTGTGGACGATGAAAAAGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGACGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGGAACGTCTTCCCg                                                                                                                                                                                          >  2:439791/1‑251 (MQ=255)                                                           aaCGTCGAAATCTGGATGCTGTTTACCAATATTATTCTGATATATGCGGCGCTGATGCTGGTCCGTTTCGGCTGGCTGTGGACGATGAAAAAGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGACGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGgc                                                                                                                                                                                        >  1:258636/1‑251 (MQ=255)                                                              gTCGAATTCTGGAGGCGGTGTACCAATATTATTCTGATATATGCGGGGCTGATGCTGGTCCGTTTCGGCTGGTGGTGGACGATGAAAAAGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGTTGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGACGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGgcgcg                                                                                                                                                                                     <  2:501078/251‑1 (MQ=255)                                                                 aaaaTCTGGAGTCTGTTTACCAATATTATTCCGAGATATGTGGCGCGGATGCTGGGCCGTTTCGGCTTGCTGTGGACGATGAAAAAGTTCAGCAACCGCTTCCGTAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCtc                                                                                                                                                                                  <  2:109648/250‑2 (MQ=255)                                                                        ggATGCTGTTTACCAATATTATTCTGATATATGCGGCGCTGATGCTGGTCCGTTTCGGCTGGCTGTGGACGATGAAAAAGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGACGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTg                                                                                                                                                                           >  1:74239/1‑251 (MQ=255)                                                                          aTGCTGTTTACCAATATTATTCTGATATATGCGGCGCTGATGCTGGTCCGTTTCGGCTGGCTGTGGACGATGAAAAAGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGACGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGgt                                                                                                                                                                         >  1:351950/1‑251 (MQ=255)                                                                               gTTTACCAATATTATTCTGATATATGCGGCGCTGATGCTGGTCCGTTTCGGCTGGCTGTGGACGATGAAAAAGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGACGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTcc                                                                                                                                                                    <  2:209802/251‑1 (MQ=255)                                                                                       atatTATTCTGATATATGCGGCGCTGAGGCTGGTCCGTTTCGGCTGGCTGTGGACGATGAAAAAGTTCAGCACCCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGACGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCt                                                                                                                                                            <  1:112995/251‑1 (MQ=255)                                                                                               cTGATATATGCGGCGCTGATGCTGGTCCGTTTCGGCTGGCTGTGGACGATGAAAAAGTTCAGCAACCGCTTCCTGAAGACGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGACGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCGTCAt                                                                                                                                                    >  1:34910/1‑251 (MQ=255)                                                                                                    atatGCGGCGCTGATGCTGGTCCGTTTCGGCTGGCTGTGGACGATGAAAAAGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGACGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCGTCAttctct                                                                                                                                               <  2:351950/251‑1 (MQ=255)                                                                                                       tGCGGCGCTGATGCTGGTCCGTTTCGGCTGGCTGTGGACGATGAAAAAGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGACGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCGTCAttctcttct                                                                                                                                            >  1:363178/1‑251 (MQ=255)                                                                                                            cgcTGATGCTGGTCCGTTTCGGCTGGCTGTGGACGATGAAAAAGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGACGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCGTCATTCTCTTCTCGCTg                                                                                                                                       <  1:17936/251‑1 (MQ=255)                                                                                                            cgcTGATGCTGGCCAGTTTCGGCCGGCTGTGAACGATGCATTTGTTCAGCACGCACTTCATGAAGACGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGCAATACTGATTGCGTCTTTCGCCGGTGTGCGTGGGACGATCACTCTGGCCGGTGAGCTCTCTATCCAGATGCTCTTGCCGGGTGGTAACGTCTTCCCGGCGCGCTCTGAGCTGGTGTTCCTGGCGGCTGGCGTAATTTTCTTCTTGCTg                                                                                                                                       >  2:45582/1‑251 (MQ=255)                                                                                                               ggATGCTGGTCCGTTTCGGATGTCTGTGGAAGACGAAAAAGTTCAGCACCCGCTTCCTGAAGAAGAAGCCGATGGAGTTGGGTTCGTGGACCACCCGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGACGATCTCTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCGTCATTCTCTTCTCGCTGttt                                                                                                                                    <  2:124352/250‑1 (MQ=255)                                                                                                                   gCTGGTCCGTTTCGGCTGGCTGTGGTCGATGAAAAAGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGACGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTGGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCGTCATTCTCTTCTCGCTGTTTGTc                                                                                                                                 <  2:177831/250‑1 (MQ=255)                                                                                                                    cTGGTCCGTTTCGGCTGGCTGTGGACGATGAAAAAGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGACGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCGTCATTCTCTTCTCGCTGTTTGTCgg                                                                                                                               >  1:263503/1‑251 (MQ=255)                                                                                                                     ggggccGTTTCGTCTGGCGGGGGACGATGAAAAAGTTAAGTAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGGGGACCACACGAGAAACCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGACGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCGTCATTCTCTTCTCGCTGTTTGTCGGc                                                                                                                              <  2:6582/247‑1 (MQ=255)                                                                                                                                   tggctgTGGACGATGAAAAAGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGACGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTCTGAGCTGGTGTTCCTGGCGGCTGGCGTCATTCTCTTCTCGCTGTTTGTCGGCGTGGTGATGTTGcc                                                                                                                >  1:320835/1‑251 (MQ=255)                                                                                                                                    ggctgTGGACGATGAAAAAGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGACGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCGTCATTCTCTTCTCGCTGTTTGTCGGCGTGGTGATGTTGCCt                                                                                                               >  1:419845/1‑251 (MQ=255)                                                                                                                                        gtgGACGAAGAAAAAGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGTAGTTTGGTTCGTTGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGACGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCGTCATTCTCTTCTCGCTGTTTGTCGGCGTGGTGATGTTGCCTAt                                                                                                             <  1:355591/249‑1 (MQ=255)                                                                                                                                                     aaGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGGGTTTGGTTCGCGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGACGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCGTCATTCTCTTCTCGCTGTTTGTCGGCGTGGTGATGTTGCCTATTCTGCTACAACACAt                                                                                              >  2:207217/1‑251 (MQ=255)                                                                                                                                                      aGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGACGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCGTCATTCTCTTCTCGCTGTTTGTCGGCGTGGTGATGTTGCCTATTCTGCTACAACACAtt                                                                                             <  1:196940/251‑1 (MQ=255)                                                                                                                                                       gTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGACGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCTCTATGAGCTGTGGTGCCTGGCGGCTGGCGTCATTCTCTTCTCGCTGTTTGTCGGCGTGGTGATGTTGCCTATTCTGCTACAACACATTg                                                                                            <  1:328833/251‑1 (MQ=255)                                                                                                                                                          cAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGACGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCGTCATTCTCTTCTCGCTGTTTGTCGGCGTGGTGATGTTGCCTATTCTGCTACAACACATTGa                                                                                           >  2:153545/1‑249 (MQ=255)                                                                                                                                                            gcatccgtttCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGACGATCACTCTGGCCGGTGTGCTCGCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCGTCATTCTCTTCTCGCTGTTTGTCGGCGTGGTGATGTTGCCTATTCTGCTACAACACATTGAAGTc                                                                                       <  2:337533/243‑1 (MQ=255)                                                                                                                                                                                          aGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGACGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCGTCATTCTCTTCTCGCTGTTTGTCGGCGTGGTGATGTTGCCTATTCTGCTACAACACATTGAAGTCGCGGATCATTCGCAACAATTGAAAGAGGaa                                                         <  1:153545/251‑1 (MQ=255)                                                                                                                                                                                          aGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGACGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCGTCATTCTCTTCTCGCTGTTTGTCGGCGTGGTGATGTTGCCTATTCTGCTACAACACATTGAAGTCGCGGATCATTCGCAACAATTGAAAGAGGa                                                          >  1:177298/1‑250 (MQ=255)                                                                                                                                                                                            tttGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGACGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTc                                                                                                                                                                     <  2:165162/141‑1 (MQ=255)                                                                                                                                                                                            tttGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGACGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTc                                                                                                                                                                     >  1:165162/1‑141 (MQ=255)                                                                                                                                                                                              tGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGACGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCGTCATTCTCTTCTCGCTGTTTGGCGGCGTGGTGATGTTGCATATTCTGCTACAACACATTGCAGTCGCGGATCATTCGCAACAATTggcagatgaacgg                                                      >  2:194381/1‑238 (MQ=255)                                                                                                                                                                                                   cGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGACGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCACTTGCCGGCTGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCTGCTGGCGTCATTCTCTTCTCGCTGTTTGTCGGCGTGGTGATGGTGCCTATTCTGCTACAAACCCTTGAAGTCGCGGATCATTCGCAACCATTGAACGAggcaagttttg                                                  >  2:295174/1‑240 (MQ=255)                                                                                                                                                                                                                agaAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGACGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCGTCATTCTCTTCTCGCTGTTTGTCGGCGTGGTGATGTTGCCTTTTCTGCTACCACACATTGAAGTGGCGGATCATTCGCAACAATTTAAAGAGGAACGTATTGCGCGTGCGGAAACgg                                   >  2:391765/1‑251 (MQ=255)                                                                                                                                                                                                                  aaaTCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGACGATCACTCTTGCCGGCATGCTCTCAATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCGTCATTCTCTTCTCGATGTTTGTCGGCGTGGTTATGTTGCCTATTCTGCTACAACACATTGAAGACGCGGATCATTCGCAACAATTGCACGAGGAACGGATTGCGCGAGCGGCCACGGcc                                 >  1:246376/1‑250 (MQ=255)                                                                                                                                                                                                                     tCCTGATCGCGTCTTTCGCCGGGGAGCGTGGGACGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCGTCATTCTCTTCTCGCTGTTTGTCGGCGTGGTGATGTTGCCTATTCTGCTACAACACATTGAAGTCGCGGATCATTCGCAACAATTGCAAGAGGCACGTATTGCGCGAGCGGCAACGGCAGa                               >  1:368368/1‑250 (MQ=255)                                                                                                                                                                                                                      ccTGATCGCGTCTTTCGCCGGGGTGCGTGGGACGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCGTCATTCTCTTCTCGCTGTTTGTCGGCGTGGTGATGTTGCCTATTCTGCTACAACACATTGAAGTCGCGGATCATTCGCAACAATTGAAAGAGGAACGTATTGCGCGAGCGGCAACGGCAGa                               >  1:444073/1‑249 (MQ=255)                                                                                                                                                                                                                        tGATCGCGTCTTTCGCCGGGGTGCGTGGGACGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCGTCATTCTCTTCTCGCTGTTTGTCGGCGTGGTGATGTTGCCTATTCTGCTACAACACATTGAAGTCGCGGATCCTTCGCAACAATTTGAAGAGGAACGTATTGCGCGAGCGGCCACGGCAGAAGTg                           >  2:206525/1‑251 (MQ=255)                                                                                                                                                                                                                           tCGCGTCTTTCGCCGGGGTGCGTGGGACGTTCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCGTCATTCTCTTCTCGCTGTTTGTCGGCGTGGTGATGTTGCCTATTCTGCTACAACACATTGAAGTCGCGGATCATTCGCAACAATTGAAAGAGGAACGTATTGCGCGAGCGGCAACGGCAGAAGTGGc                         <  2:363178/250‑1 (MQ=255)                                                                                                                                                                                                                              cgTCTTTCGCCGTGGTGCGTGGTACGATCACTCTGGGCGGTGTGCTCTCTATCCCGCTGCTCTTGTCGGATGGTAACGCCTTCCTGGCGCGCTATGAGCTGGTGTTCCCGGCGGCTGGCGTCATTCGCTTCTCGCTGTTTGGCGGCGTGGTGATGTTGCCTATTCTGCTACAACACATTGAAGTCGCGGATCATTCGCAACAATTGAAAGAGGAACGTATTGCGCGAGCGGCAACGGCAGAAGTGGCGAtt                     <  2:305884/251‑1 (MQ=255)                                                                                                                                                                                                                                   ttCGCCGGGGTGCGTGGGACGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCGTCATTCTCTTCTCGCTGTTTGTCGGCGTGGTGATGTTGCCTATTCTGCTACAACACATTGAAGTCGCGGATCATTCGCAACAATTGAAAGAGGAACGTATTGCGCGAGCGGCAACGGCAGAAGTGGCGATTGTgg                 <  1:439791/250‑1 (MQ=255)                                                                                                                                                                                                                                     cGCCGGGGTGCGTGGGACGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCGTCATTCTCTTCTCGCTGTTTGTCGGCGTGGTGATGTTGCCTATTCTGCTACAACACATTGAAGTCGCGGATCATTCGCAACAATTGAAAGAGGAACGTATTGCGCGAGCGGCAACGGCAGAAGTGGCGATTGTGGCGa              <  2:177298/251‑1 (MQ=255)                                                                                                                                                                                                                                           ggTGCGTGGGACGATCACTCTGGCCGGTGTGCTCTCTATCACGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCGTCATTCTCTTCTCGCTGTTTGTCGGCGTGGTGATGTTGCCTATTCTGCTACAACACATTGAAGTCGCGGATCATTCGCAACAATTGAAAGAGGAACGTATTGCGCGAGCGGCCACGGCAGAAGCGGCGATTGTGGCGGTCCAgg        >  2:423455/1‑250 (MQ=255)                                                                                                                                                                                                                                            gTGCGTGGGACGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCGTCATTCTCTTCTCGCTGTTTGTCGGCGTGGTGATGTTGCCTATTCTGCTACAACACATTGAAGTCGCGGATCATTCGCAACAATTGAAAGAGGAACGTATTGCGCGAGCGGCAACGGCAGAAGTGGCGATTGTGGCGATCCAGaa       <  2:155411/251‑1 (MQ=255)                                                                                                                                                                                                                                                 tGGGACGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCGTCATTCTCTTCTCGCTGTTTGTCGGCGTGGTGATGTTGCCTATTCTGCTACAACACATTGAAGTCGCGGATCATTCGCAACAATTGAAAGAGGAACGTATTGCGCGAGCGGCAACGGCAGAAGTGGCGATTGTGGCGATCCAGAAAATgg  >  1:134098/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       AGGTCTGCAGCTGCCGGGTATTCTGGAGACGTCGCTGATGGCAGCAGAAATCGACCCTAACGTCGAAATCTGGATGCTGTTTACCAATATTATTCTGATATATGCGGCGCTGATGCTGGTCCGTTTCGGCTGGCTGTGGACGATGAAAAAGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCGTCATTCTCTTCTCGCTGTTTGTCGGCGTGGTGATGTTGCCTATTCTGCTACAACACATTGAAGTCGCGGATCATTCGCAACAATTGAAAGAGGAACGTATTGCGCGAGCGGCAACGGCAGAAGTGGCGATTGTGGCGATCCAGAAAATGG  >  NC_000913/4280849‑4281339

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