Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,856,243	A→G	T636A (ACG→GCG)	fhlA →	DNA‑binding transcriptional activator FhlA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,856,243	0	A	G	98.0%	146.7 / NA	49	T636A (ACG→GCG)	fhlA	DNA‑binding transcriptional activator FhlA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (18/30);  minor base T (1/0);  total (19/30)
Fisher's exact test for biased strand distribution p-value = 3.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GAGCGTCCGGAAGATATTCCGCTGCTGGCGAAAGCCTTTACCTTCAAAATTGCCCGTCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGCGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTC  >  NC_000913/2855997‑2856485                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   gAGCGTCCGGAAGATATTCCGCTGCTGGCGAAAGCCTTTACCTTCAAAATTGCCCGTCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGCGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAGCGGt                                                                                                                                                                                                                                                <  1:236339/251‑1 (MQ=255)    cGTCCGGAAGATATTCCGCTGCTGGCGAAAGCCTTTACCTTCAAAATTGCCCGTCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGCGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAGCGGTTGt                                                                                                                                                                                                                                             <  1:288266/251‑1 (MQ=255)     gTCCGGAAGATATTCCGCTGCTGGCGAAAGCCTTTACCTTCAAAATTGCCCGTCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGCGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAGCGGTTGTc                                                                                                                                                                                                                                            <  2:49706/251‑1 (MQ=255)           aaGATATTCCGCTGCTGGCGAAAGCCGTTACCTTCAAAATTGCCCGTCGTCTGGGGCGCAATAGCGACAGCATTCCTGCCGAGACGCTGCGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAGCGGTTGTCGCCCTg                                                                                                                                                                                                                                      <  2:427058/251‑1 (MQ=255)                atTCCGCTGCTGGCGAAAGCCTTTACCTTCAAAATTGCCCGTCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGCGCACCTTGAGCAAAATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAGCGGTTGTCGCCCTGGAggg                                                                                                                                                                                                                                 >  1:296833/1‑251 (MQ=255)                   ccGCTGCTGGCGAAAGCCTTTACCTTCAAAATTGCCCGTCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGCGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAGCGGTTGTCGCCCTGGAGGGCg                                                                                                                                                                                                                               <  2:345209/250‑1 (MQ=255)                                                 aTTGCCCGTCGTCTGGGGCGCAATAGCGACAGCATTCCTGCCGAGACGCTGCGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGATAACACGCGGTAACGGGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAGCGCCGACTGCAGCATCGGTTGTCGCCCTGGAGGGCGTAGATGAATATCAGTTGATTGGGCGCGTgtt                                                                                                                                                                                                >  2:12126/1‑249 (MQ=255)                                                     cccGTCGTCTGGGGAGCAATACCGACAGCATTCCTGCCGAGACGCTGCGCACCTTGAGCAACATGGAGGGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCACTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAGCGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTgaaa                                                                                                                                                                                            <  1:104388/251‑1 (MQ=255)                                                           gtcTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGCGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGACCTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAGCGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAAcc                                                                                                                                                                                      >  1:431650/1‑251 (MQ=255)                                                               ggggCGCAATATCGACAGCATTCCTGCCGAGACGCTGCGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAGCGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAAc                                                                                                                                                                                   <  1:153003/250‑1 (MQ=255)                                                               ggggCGCAATATCGACAGCATTCCTGCCGAGACGCTGCGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAGCGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACg                                                                                                                                                                                  <  1:97944/251‑1 (MQ=255)                                                               ggggCGCAATATCGACAGCATTCCTGCCGAGACGCTGCGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAGCGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACg                                                                                                                                                                                  >  1:310218/1‑251 (MQ=255)                                                                gggCGCAATATCGACAGCATTCCTGCCGAGACGCTGCGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAGCGGTTGTCGCCCTGGAGGGCGCAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACgg                                                                                                                                                                                 >  2:365892/1‑251 (MQ=255)                                                                gggCGCAATATCGACAGCATTCCTGCCGAGACGCTGCGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAGCGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACgg                                                                                                                                                                                 <  2:35842/251‑1 (MQ=255)                                                                             aCAGCATTCCTGCAGAGACGCTCCGCCCATTGCGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAGCGGTTGTCGCCCTGGAGGGCGAAGATGAATATCCGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCggg                                                                                                                                                                    >  2:324442/1‑251 (MQ=255)                                                                               aGCATTCCTGCCGAGACGCTGCGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAGCGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGcc                                                                                                                                                                  <  1:305343/251‑1 (MQ=255)                                                                                     ccTGCCGAGACGCTGCGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAGCGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAgg                                                                                                                                                            <  1:220705/251‑1 (MQ=255)                                                                                          cGAGACGCTGCGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAGCGGTTGTCGCCCCGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTg                                                                                                                                                       <  2:334371/251‑1 (MQ=255)                                                                                                   gcgcACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATGGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAGCGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGt                                                                                                                                               <  1:324442/250‑1 (MQ=255)                                                                                                      cACCTTGAGCAACATGGAGTGGCCGGGTAACGTAAGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTCACGAGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAGCGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTgg                                                                                                                                           >  2:435404/1‑251 (MQ=255)                                                                                                           tGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCAGTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAGCGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCt                                                                                                                                       <  2:132430/250‑1 (MQ=255)                                                                                                               caccATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGCGCGCGCGGTATTGCTCACACACGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGACCCTGAAACGCCGCCTGCCGCAGCGGTTGTCGCCCTGGCGGGCGCAGTTGAATATCAGTTGATTGTGCGCGTGCTGCACGAAACCAACGGCGTGGTTGCCGGGCCTAACGTCGCTGCGCAACGTCTGGGGCTGaaa                                                                                                                                   >  2:406510/4‑250 (MQ=255)                                                                                                                 acagggAGTGGCCGGGTAACGTTCGCGAACTGGAAAACGTCATGGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAGCGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGTCGCTGCGCACCGTCTGGGGCt                                                                                                                                       <  1:335227/240‑1 (MQ=255)                                                                                                                 aCATGGAGTGGCCGGGTAACGTAAGCGAACTGGAAACCGTCATTGAGCGCGCGGTATTGCTAACAAGCGGTACCGTGCTGCAGCTGTCTTTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAGCGGTTGTCGCCCTGGAGGGCGCCGATGAATCTCCGGTGATTGTGCGCGGGCTGAAAGGAACCAACGGCGGGGGTGCCGGGCCTAAAGGAGATGCGCCAAGTCTGGGGCtaggtcgg                                                                                                                                >  2:335227/1‑244 (MQ=255)                                                                                                                      gAGTGGCCGGGTAACGTACGCGAACTGGAAATCGTCATTGGACGAGCGGTATTGCTAACACGCGGTAACGTGCTGCAGATGTCATTGCCAGATATGGTTGTACCGGAACCTGAAGCGCCGCCTGCCGCAGCGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGAc                                                                                                                           <  2:102815/251‑1 (MQ=255)                                                                                                                       ggTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAGCGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACTTCTGGGGCTGAAACGCACGAcc                                                                                                                          <  1:168532/250‑1 (MQ=255)                                                                                                                                          cgAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAGCGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGc                                                                                                       <  1:407873/251‑1 (MQ=255)                                                                                                                                           gAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAGCGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCg                                                                                                      <  1:141270/251‑1 (MQ=255)                                                                                                                                                  aaaaCGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAGCGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGACAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGACTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTggg                                                                                                >  1:244638/1‑250 (MQ=255)                                                                                                                                                   aaaCGTCATGGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAGCGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGaa                                                                                              <  2:145018/251‑1 (MQ=255)                                                                                                                                                      cGTCATTGAGCGCGCGGTATTGCTAACCCGCGGTAACGTGCTGCAGCTGTCCTTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAGCGGTTGTCGCCCTGGAGGGCGACGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGTTGAAGCGGCTGGGAATTg                                                                                           >  2:154513/1‑251 (MQ=255)                                                                                                                                                      cGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAGCGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCaa                                                                                                                                                  <  1:429046/196‑1 (MQ=255)                                                                                                                                                      cGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAGCGGTTGTCGCCCTGGAGGGCGAAGATGAATATAAGTTGATTGTGCGCGTGCTGCAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCaa                                                                                                                                                  >  2:429046/1‑196 (MQ=255)                                                                                                                                                          agggAGCGTGCGGTATTGCTAACACTCGTTAACGTGCTGCAGCTGTAATTGCCAGATATTGTTTTACCGGAACCTGCAACGCTGCCTGCCGCAGCGGTTGTCGCACTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGGCTGGGGCTGAAACGCACGGCCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATaa                                                                                       <  2:343220/248‑1 (MQ=255)                                                                                                                                                           ttGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGATGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAGCGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATaa                                                                                       >  1:137440/1‑250 (MQ=255)                                                                                                                                                                  cgcgGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCGGAAACGCCGCCTGCCGCAGCGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCAt                                                                               <  1:181073/251‑1 (MQ=255)                                                                                                                                                                          tgctaacgcgcgctcccGTGCTGAGGCTGTCATTACCAGATATTGTTTTACCGGACCCTGAACCGCCGCCTGCCGCCGCAGTTGTCGACCTGGAGGGCGCAGATGACTATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGAGCCTCAAGGCGCTGCGCCACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGGTGATGAGGCTGGGGATTGATATATCGGAATTGATTTa                                                                        >  2:50988/17‑250 (MQ=255)                                                                                                                                                                          gggtAACACGCGGTTACGTGCTGCCGCTGTCAATGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAGCGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATGGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCTCGGAGGAAGCGGCTGGGAAGTGATAAATCGGCATTGATTTa                                                                        <  2:334655/247‑1 (MQ=255)                                                                                                                                                                               acacGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAGCGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTgaaagaaa                                                                                                                                                                                        >  2:316416/1‑133 (MQ=255)                                                                                                                                                                               acacGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAGCGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTgaaagaaa                                                                                                                                                                                        <  1:316416/133‑1 (MQ=255)                                                                                                                                                                                cacGCGGTAACGTGCTGCAGCGGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAGCGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCaa                                                                 <  1:80919/251‑1 (MQ=255)                                                                                                                                                                                       tAACGAGCTGCAGCTGTCATTGCCAGATATTGATTTACCGGAACCTGAAACGCCGCCTGCCGCAGCGGTTGTCGCCCTGGAGGGCGCAGATGAATATCCGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCGGCTGTCACGGATGACGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCg                                                          >  1:97092/1‑251 (MQ=255)                                                                                                                                                                                        aaCGTGCGGCAGCTGTCAGTGCCAGCTTTCGTTTTCCCGGACCCGGAAACGCCGCCGGCCGCAGCGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGGTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGcca                                                          <  2:11686/250‑2 (MQ=255)                                                                                                                                                                                              ctgcAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAGCGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCCCGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGAc                                                       >  1:210469/1‑247 (MQ=255)                                                                                                                                                                                                               aGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAGCGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATAc                                  <  2:310218/251‑1 (MQ=255)                                                                                                                                                                                                                 atatTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAGCGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTa                                <  1:202801/251‑1 (MQ=255)                                                                                                                                                                                                                      gTTTTACCGGAACCTGAAACGCCGCCTGCCGCAGCGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCAt                           <  1:105410/251‑1 (MQ=255)                                                                                                                                                                                                                       ttttACCGGAACCTGAAACGCCGCCTGCCGCAGCGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATg                          >  1:93658/1‑251 (MQ=255)                                                                                                                                                                                                                                               cTGCCGCAGCGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTc  >  1:229401/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   GAGCGTCCGGAAGATATTCCGCTGCTGGCGAAAGCCTTTACCTTCAAAATTGCCCGTCGTCTGGGGCGCAATATCGACAGCATTCCTGCCGAGACGCTGCGCACCTTGAGCAACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTC  >  NC_000913/2855997‑2856485

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