Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	6,748	G→A	S404S (AGC→AGT)	yaaJ ←	putative transporter YaaJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	6,748	0	G	A	98.4%	164.7 / NA	63	S404S (AGC→AGT)	yaaJ	putative transporter YaaJ
Reads supporting (aligned to +/- strand):  ref base G (1/0);  major base A (33/27);  minor base C (1/0);  total (36/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

GAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACA‑CCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACAAGGGTAACAAACTCAGCACCCCAGGAACCCATT‑AGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTTGT  >  NC_000913/6503‑6987                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  gagaatgagTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCCAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACtt                                                                                                                                                                                                                                                >  1:75257/1‑249 (MQ=255)    aatgagTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCa                                                                                                                                                                                                                                            >  1:492718/1‑250 (MQ=255)            gATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTgccg                                                                                                                                                                                                                                    <  1:607216/250‑1 (MQ=255)             aTCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCCGGCCATTATGATATCTGCCAGTTGCCCCATCCGCGGAAGACTTAGCAAGGTgccgcc                                                                                                                                                                                                                                  >  2:103169/1‑251 (MQ=255)             aTCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTgccgcc                                                                                                                                                                                                                                  >  1:457917/1‑251 (MQ=255)               cgataggtgagatgccTCCTGAGAAAAATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTAAGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGGTAAATTGGAAATAACCATGCAGGGCATTATGATATCTGCCAGTTGCCAAATCAGCGGTAGACTAAGCAAGGTGgcgccg                                                                                                                                                                                                                                 >  2:431859/16‑250 (MQ=255)                         gATTACTCCTGAGAACCATAATCCCACGCGTCCGGAGAAAGCTGGCGGCAGATATCCGGCTAACGCAACGGATCAAACACCGGGCGCGAGCCGAGTTTACGCTGGCGTAGATAATACCTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGACATGAAAGCCATTAAGATATCTGCCAGATGCCCCATCAGCGGAAAACGTCGAACGGTGCCGCCGACGACCGTTGc                                                                                                                                                                                                                      >  2:375006/1‑251 (MQ=255)                           ttACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCAAATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCg                                                                                                                                                                                                                     >  1:67441/1‑250 (MQ=255)                             aCTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATCCAGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTTGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGCAGAATTAGCAAGGTGCCGCCGATGACCGTTGCGa                                                                                                                                                                                                                    >  1:46215/1‑249 (MQ=255)                                   gCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCa                                                                                                                                                                                                             <  1:117745/250‑1 (MQ=255)                                       aaCATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCTGAATACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATcc                                                                                                                                                                                                        <  1:301475/251‑1 (MQ=255)                                        aCATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCg                                                                                                                                                                                                       <  2:447593/251‑1 (MQ=255)                                          atcatcCCACGCGTCCGGAGAAAGCTGGCGGCCGATATCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTCTGACCCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTACCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGc                                                                                                                                                                                                      >  1:89815/1‑250 (MQ=255)                                                 cACGCGTCCGGCGAAAGCTGGCGGCCGATATCCGGATAACGCAACGGATCAAACACCGGGCGCACGCAGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGGAATAGCCATGCAGGCCATTATGATATTCGCCAGTTGCCCCATCAGCGGAAGACTTTGCAAAGTGCCGCCGATGACCGCTGCGAAGTGGCAGAACCGCacacgcc                                                                                                                                                                                               >  2:55429/1‑244 (MQ=255)                                                        ccGGAGAAAGCTGGCGGCCGATATCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACA‑CCAGATCGCt                                                                                                                                                                                      <  1:247606/251‑1 (MQ=255)                                                         cGGAGAAAGCTGGCGGCCGATATCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGGGCAGATCCGCAAACA‑CCAGATCGCtt                                                                                                                                                                                     >  1:342519/1‑251 (MQ=255)                                                              aaaGCTGGCGGCCGATATCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACA‑CCAGATCGCTTTAggg                                                                                                                                                                                >  1:236954/1‑251 (MQ=255)                                                              aaaGCTGGCGGCCGATATCCGGATAACGCAACGGATAAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTCCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACA‑CCAGATCGCTTTAggg                                                                                                                                                                                <  2:95768/251‑1 (MQ=255)                                                                 gCTGGCGGCCGATATCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGACATGCAGGCCATTTTGATATCTGCCAGTTGCCACCTCAGCGGGAGAATTAGCAAGGTGCCGCCGATGACCGTTTCGAAGGGGCAGATCCGCAAACA‑CCCGGTCGCTTTAGGgttt                                                                                                                                                                             >  2:468409/1‑250 (MQ=255)                                                                  cTGGCGGCCGATATCCGGATAACGCAACGGATGAACCACCGGGCGAACGCCGCGTTTCCGCTGGCGTCGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGGCCAATTGGATATAGCCCGGGAGGCCATTTTGAGATCTGCCAGTTGCCACGTTAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGGGAAGGTGCAGATCCGGAAACA‑CCAGCTCGCTTTAGGgtagg                                                                                                                                                                            >  2:45996/1‑248 (MQ=255)                                                                  cTGGCGGCCGATATCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACA‑CCAGATCGCTTTAGGgt                                                                                                                                                                               <  1:45996/248‑1 (MQ=255)                                                                              atCCGGATAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGGGCCGCCGATGACCGTTGCGCAGGGGCAGAACCGCAAACA‑CCAGATCGATTTAGGGGTGTTCAGGGGTaaa                                                                                                                                                                 >  2:484517/1‑250 (MQ=255)                                                                                    aTAACGCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACA‑CCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAagaaga                                                                                                                                                          <  1:103169/251‑1 (MQ=255)                                                                                      aaCGCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACA‑CCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAagag                                                                                                                                                         <  1:87428/250‑1 (MQ=255)                                                                                         gCAACGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACA‑CCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTg                                                                                                                                                     <  1:320958/251‑1 (MQ=255)                                                                                             cggATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAACATGGCGGTCAAATTGGTAATAGCCATGCAAGCCATTTTAATATCTGCCAGTTGCCACAACATCGGAAAACTTAGCAAGGGGACGCCGATGAACGTTGCGAAGGTGCAGACACGCAAACC‑CCCGCTCGCATTAGGGGGGTTCAGGAGTGAAAAGAAGAGATTGtttt                                                                                                                                                 >  2:446226/2‑251 (MQ=255)                                                                                             cggATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACA‑CCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGtttt                                                                                                                                                 <  2:342519/250‑1 (MQ=255)                                                                                              ggATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACA‑CCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTc                                                                                                                                                <  2:470539/251‑1 (MQ=255)                                                                                              ggATCAAACACCGGGCGCACGCCGAGTTTACGCGGGCGTAGATAACCACTGGCAATGGTACGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGCCCATTATGATATCTGCCAGGTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTAGCGAAGTTGCAGACTCGAAAACA‑CCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAATAGATGGTTTTc                                                                                                                                                <  2:134444/251‑1 (MQ=255)                                                                                                  caaacaaagggCGCACGCCGAGTTTACGATGGAGGAGATAATCACTGGCAATGGCCTGAACCACAGGCGAGAGCAGTAAAATGGCGGAAAAATTGGGAATAGCCATGAAGGCCATTATGATATTTGCCAGTTGCCACATCAGCGGAAGACCTCGCAAGGTGCCGCCGATGACCGTTGCGAAGGAGCAGATCCGCCAACAACCA‑ATAGCTTTAGGGTTGTTCCGGCGTAAAAAGAAGAGAtttgtttaggcg                                                                                                                                            >  2:402158/9‑241 (MQ=255)                                                                                                   aaaCACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACA‑CCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGcgt                                                                                                                                           <  1:345912/251‑3 (MQ=255)                                                                                                   aaaCACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACA‑CCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGcgt                                                                                                                                           <  1:91775/251‑3 (MQ=255)                                                                                                          gggCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACA‑CCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTa                                                                                                                                    <  2:236954/251‑1 (MQ=255)                                                                                                           ggCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACA‑CCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATgg                                                                                                                                     >  2:155584/1‑248 (MQ=255)                                                                                                            gcgcACGCCGAGTTTACGCTGGCGCAGATAATCACTGGCAATGGTATGAACCACAGGCGAGATCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGTCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTTCAGATCCGCAAACA‑CCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGt                                                                                                                                  <  1:123991/251‑1 (MQ=255)                                                                                                                   cctaGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGAGATCTGCCAGTTGCCACATCATCGGAAGACTTAGCAAGGTTCCGCCGATGACCGTTGCGAAGGTGCAGAGCCGCAATCA‑CCAGATCTCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAAc                                                                                                                           <  1:402158/248‑1 (MQ=255)                                                                                                                     gAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGTAGATCCGGAAACA‑CCAGATCGCTTTAGGGTTGTTCAGGGGGAAAAAGAAGAGATTGTTTTCGGCGTAAAAGTAGTTGGAAACg                                                                                                                          >  2:290581/1‑250 (MQ=255)                                                                                                                            cGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAACTTTGTAATAGCCATGCAGGCCATTATGATATCTGCCCGTTGCCACATCAGCGGAAGAATTAGCAAGGGGCCGCCGAAGAACGTTGCGAAGGTGCAAATCCGCAAACA‑CCAGATCGCGTTCGGGATGTTCAGGCGTAAAAAGAAGAGAATGATTTCGGCGTTAATGTAGTTGGCAAAGATGGAGc                                                                                                                   >  2:190573/1‑250 (MQ=255)                                                                                                                               gggCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGATAGCAGTAAAATGGCGGTCAAATTGGTAAGAGCCATGCAGGCCATTATGATATCTGCCATGTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACA‑CCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGaa                                                                                                               <  2:205916/250‑1 (MQ=255)                                                                                                                                  cGTAGCTAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGGAATAGCCATGCACGCCATTATGATATCTGCCCGTTTCCCCATCAGCGGAATACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGGGCAGATCCGCAAACA‑CCAGATCGGTTTAGGGATGTTCAGGCGGAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGGTGGCAACGATGGAGCTGAAgg                                                                                                             >  1:42441/1‑250 (MQ=255)                                                                                                                                      gATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCCGTTGCCCCATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGAGAAGGGGCAGATCCGCAAACA‑CCCGATCGCTTTAGGGGTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGGTGGCAACGAAGGAGCTGAAGGCAAAc                                                                                                        >  2:345721/1‑251 (MQ=255)                                                                                                                                              cTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGGAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACA‑CCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGGTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCGGAAGGCAAACAGAATAAc                                                                                                >  1:337547/1‑251 (MQ=255)                                                                                                                                                ggCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACA‑CCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAAcc                                                                                               <  2:457917/250‑1 (MQ=255)                                                                                                                                                                  aGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACA‑CCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTc                                                                              >  1:230600/1‑249 (MQ=255)                                                                                                                                                                      gagaGCAGTAAAATGGCGGTCAAATTGGTAAGAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACA‑CCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCAcc                                                                        <  2:75257/251‑1 (MQ=255)                                                                                                                                                                               aaaaTGGCGGTCAAATTGGTAATAGCCCTGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACA‑CCCGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCa                                                                                                                             >  2:257558/1‑189 (MQ=255)                                                                                                                                                                               aaaaTGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACA‑CCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCa                                                                                                                             <  1:257558/189‑1 (MQ=255)                                                                                                                                                                                aaaTGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACA‑CCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAAccc                                                              >  2:542506/1‑251 (MQ=255)                                                                                                                                                                                 aaTGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACA‑CCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGATGGTAACAAACTCAGCACCCCAGGAACCCa                                                             <  2:308851/251‑1 (MQ=255)                                                                                                                                                                                    ggCGGTAAAATTGGTAATAGCCATGCAGGCCCTTATGATATCGGGCAGTTAACACATCAGCGGAAGTATTCGCCTGATGACGCCGATGAACGTTGCGAAGGTGCAGATCCGCAAACA‑CAAGATAGCTTTAGGGTTGTTCTGGCGTACACAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAATAGAATAACCACGAGGCTGACAAACTCAGCACCCCAGGAACCCAtt                                                           >  2:59482/1‑250 (MQ=255)                                                                                                                                                                                     gCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGACGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACA‑CCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATT‑a                                                         >  1:272819/1‑250 (MQ=255)                                                                                                                                                                                       ggTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACA‑CCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCACAGGAACCCATT‑aac                                                       >  1:51876/1‑248 (MQ=255)                                                                                                                                                                                       ggTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACA‑CCAGATCGCTTTAGGGGTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAAc                                                                                    >  1:426211/1‑222 (MQ=255)                                                                                                                                                                                                      aGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACA‑CCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGTGTAACAAACTCAGCACCCCAGGAACCCATT‑AACACCCGCATCGCCTTc                                        <  1:550022/250‑1 (MQ=255)                                                                                                                                                                                                      aGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACA‑CCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCAACCCAGGAACCCATT‑AACACCCGCATCGCCTTCt                                       >  1:19329/1‑251 (MQ=255)                                                                                                                                                                                                      aGCCATGCAGGCCATTATGATAGCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACA‑CCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATT‑AACACCCGCAACGCCTTc                                        <  2:2197/250‑1 (MQ=255)                                                                                                                                                                                                             cAGGCCATTACGATCTCTGCCGGTTGCCACATCAGCGGAATACTTAGCAAGGTGACGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACA‑CCAGCTCGCTTTAGGGATGTTCCGGCGTACAAAGCAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGAAAACAGAATAACCACGAGGGGAACAAACTCAGCACCCCAGGAACCCATT‑AACACCGGCATCGCCTTCTGGATa                                  >  2:304487/1‑249 (MQ=255)                                                                                                                                                                                                                           atCTGCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCAGCAAACA‑CCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGAAAAAAACTAAGCACCCCAGGA‑CCCATTTAACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTcc                   >  1:572259/1‑250 (MQ=255)                                                                                                                                                                                                                               gCCAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGAGCCGCAAACA‑CCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATT‑AACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCgg              <  1:399113/251‑1 (MQ=255)                                                                                                                                                                                                                                 cAGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACA‑CCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATT‑AACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCa            <  1:155584/251‑1 (MQ=255)                                                                                                                                                                                                                                  aGTTGCCACATCAGCGGACGCCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACA‑CCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATT‑AACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCAt           >  1:472022/1‑251 (MQ=255)                                                                                                                                                                                                                                  aGTTGCCACATCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACA‑CCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATT‑AACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCAt           >  2:394866/1‑251 (MQ=255)                                                                                                                                                                                                                                           aTCAGCGGAAGACTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACA‑CCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCGTAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACGAGGGTAACAAACTCAGCACCCCAGGAACCCATT‑AACACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTtgt  >  1:136287/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  GAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACA‑CCAGATCGCTTTAGGGTTGTTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACAAGGGTAACAAACTCAGCACCCCAGGAACCCATT‑AGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTTGT  >  NC_000913/6503‑6987

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