Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,509,606	C→T	intergenic (‑180/‑24)	glk ← / → yfeO	glucokinase/putative transport protein YfeO

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,509,606	0	C	T	97.8%	137.4 / ‑4.8	45	intergenic (‑180/‑24)	glk/yfeO	glucokinase/putative transport protein YfeO
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (26/18);  total (26/19)
Fisher's exact test for biased strand distribution p-value = 4.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.43e-01

ATATCACACAGAGCAAGACGTGCGTTGGTGCCGCCCACATCACCGACTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAGTCAAAATAATTCTTTCTCACACTGTAAATACCTGGGGGCATAACAACAACGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAATTGCAGCGATGCCTTTTTGCAAGC‑TGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACATGCTCCATCCGCGAGCCAGAACCATGTTGTTATTATCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTGATTGTGGTGATGAAAATCGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATC  >  NC_000913/2509362‑2509851                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      atatCACCCAGAGCAAGACGTGCGTTGGTGCCGCCCTCATCACCGACTAATGCATTCTTTGTCCTTCTTCATCTGCTCCGCTAAAGTCAAAATAATTCTTTCTCTCTCTGTAAATACCTGGGTGCCTAACAACAACGCCGGAAAGGCAGGCTTCCTGTATATATCGATCTGGGTCACTCAATTACTTTATCGTTTCAGCGCCAATTCCATCGATGCCTTTTTGCAAGC‑TGGGCAAACTAAGTATTTGAccc                                                                                                                                                                                                                                                 <  2:359057/251‑1 (MQ=255)      acacaGAGCAAGACGTGCGTTGGTGCCGCCCACATCACCGACTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAGTCAAAATAATTCTTTCTCACACTGTAAATACCTGGGGGCATAACAACAACGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAATTGCAGCGATGCCTTTTTGCAAGC‑TGGGCAAACTAAGTATTTGACCCCGCAt                                                                                                                                                                                                                                            >  1:398628/1‑251 (MQ=255)            agCAAGACGTGCGTTGGTGCCGCCCACATCACCGACTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAGTCAAAATAATTCTTTCTCACACTGTAAATACCTGGGGGCATAACAACAACGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAATTGCAGCGATGCCTTTTTGCAAGC‑TGGGCAAACTAAGTATTTGACCCCGCATAAGGaa                                                                                                                                                                                                                                      <  1:271310/251‑1 (MQ=255)              cAAGACGTGCGTTGGTGCCGCCCACATCACCGACTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAGTCAAAACAATTCTTTCTCACACTGTAAATACCTGGGGGCATAACAACAACGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAATTGCAGCGATGCCTTTTTGCAAGC‑TGGGCAAACTAAGTATTTGACCCCGCATAAGGAATa                                                                                                                                                                                                                                    <  2:175166/251‑1 (MQ=255)                         atggtaccgaccgcaTCACCGCCTAATGCCTACTTTGTCATTCTTCAACTGCTACGCTAAAGTCAAAATAATTCTTTCTCACACTGTAAATACCTGGGGGCATAACAACAACGCCGGACAGGCAGGCTCCCTGAAAATATCGATCTGAGTCACACACTTACTTTATCGATTCAACACCAATTGCAGCGATGCCTTTTTGCAAGC‑TGGGCAAAATAAGTATTTGACCCCGCATAAGGAATAGAACATGCTcc                                                                                                                                                                                                                         >  2:33677/14‑251 (MQ=255)                              gcctcccACATCACCGACTAATGCATTTTTGGTCATTCTTCCACTGCTACGCTAAAGTCAAAATAATTCTTTCTCACACTGTAAATACCTGGGGGCATAACAACAACGCCGGAAAGGCAGGCGCCCTGTAAATATCGATCTGGTGCACACAATTACTTTATCGTTTCAGCACCAATTGCAGCGATGCCTTTTTTTAATC‑TGGGCAATCTAAGTATTTTACCCCGCATAAGTAATAGAACATGCTCCATcct                                                                                                                                                                                                                    <  2:225298/247‑2 (MQ=255)                                         cACCGACTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAGTCAAAATAATTCTTTCTCACACTGTAAATACCTGGGGGCATAACAACAACGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAATTGCAGCGATGCCTTTTTGCAAGC‑TGGGCAAACTAAGTATTTGACCCCGCATAAGGAATAGAACATGCTCCATCCGCGAGCCAGAAc                                                                                                                                                                                                         >  2:169052/1‑251 (MQ=255)                                            cGACTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAGTCAAAATAATTCTTTCTCACACTGTAAATACCTGGGGGCATAACAACAACGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAATTGCAGCGATGCCTTTTTGCAAGC‑TGGGCAAACTAAGTATTTGACCCCGCATAAGGAATAGAACATGCTCCATCCGCGAGCCAGAACCAt                                                                                                                                                                                                      >  1:232270/1‑251 (MQ=255)                                                               cATTCTTCAACTGCTCCGCTAAAGTCAAAATAATTCTTTCTCACACTGTAAATACCTGGGGGCATAACAACAACGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAATTGCAGCGATGCCTTTTTGCAAGC‑TGGGCAAACTAAGTATTTGACCCCGCATAAGGAATAGAACATGCTCCATCCGCGAGCCAGAACCATGTTGTTATTATCGCTcccc                                                                                                                                                                                   >  1:356814/1‑251 (MQ=255)                                                                            cTCCGCTATAGCCAAAATAATTCTTTCTCACACTGTAAATACCTGGGGGCATAACAACAACGCCGGAAGCGCAGGCACCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAATTGCAGCGATACCTTTTTGCAAGC‑TGGGCAAACTAAGAATTTGAACCCGCATAAGGAATAGAACATGCTCCATCCTCGAGCCGGTACCATCTCGTTATTATCGCTCCCCGCCGTGGCAATTg                                                                                                                                                                      >  2:88869/1‑251 (MQ=255)                                                                                 attaaGTCAAAATAATTCTTTCTCACACTGTAAATACCTGGGGGCATAACAACAACGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAATTGCAGCGATGCCTTTTTGCAAGC‑TGGGCAAACTAAGTATTTGACCCCGCATAAGGAATAGAACATGATCCATCCGCGAGCCAGAACCATGTTGTTATTATCGCTCCCCGCCGTGGCAATTGGGAtt                                                                                                                                                                 >  2:97229/4‑251 (MQ=255)                                                                                    aaGTCAAAATAATTCTTTCTCACACTGTAAATACCTGGGGGCATAACAACAACGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAATTGCAGCGATGCCTTTTTGCAAGC‑TGGGCAAACTAAGTATTTGACCCCGCATAAGGAATAGAACATGCTCCATCCGCGAGCCAGAACCATGTTGTTATTATCGCTCCCCGCCGTGGCAATTGGGATTGCg                                                                                                                                                              >  2:21085/1‑251 (MQ=255)                                                                                     agtaaaaaTAATTCTTTCTCACACTGTAAATACCTGGGGGAATAAAAACAAAGCCGAAAAGGCAGGCTCCTTGAAAAAATCGATCTGGGTCACACAATTCTTTTATCTTTTCAGCACCAATTGCAGCGATGCCTTTTTGCAAGT‑TGGGCAAACTAAGTATTTGACCCCGCATAAGTAAAAGAACATGCCCCATCCGCGAGCCGGAAGCATGTTGTTATTTCCGCCCCCCGCCGTGGCAATTGGGATTGCGt                                                                                                                                                             <  1:147661/247‑1 (MQ=255)                                                                                                   ttcctcACACGGTAAATACCTGGGGCCATAAAAAAAACGCGGGAAAGGCAGGCTCCATGTAAATATAGATCGGGGTCACAAAATTACTTTATCGTTTAAGCACCAATTGAAGCGATGCCTTTTTGAAAGC‑TGGGAAAACTAAGTATTGGACCCCGAAAAAGGAATAGAAAAGGCTACATCCGCGAGCCAGAACCAGGTTGTTATTATCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCt                                                                                                                                                <  1:225559/247‑1 (MQ=255)                                                                                                    ttctCACACTGTAAATACCTGGGGGCAAAACAACAACGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGCCACACAATTACTTTATCGTTTCAGCACCAATTGCAGCGAGGCTTTTTTGAAAGC‑TGGGAAAACTAATTATTTGCCCCCGCATAAGGAATAGAACATTCCCCATCCGCGAGCCAGAACCATTTTGTTATTATCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCt                                                                                                                                                <  1:326915/249‑1 (MQ=255)                                                                                                            cTGTAAATACCTGGGGGCATAACAACAACGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAATTGCAGCGATGCCTTTTTGCAAGC‑TGGGCAAACTAAGTATTTGACCCCGCATAAGGAATAGAACATGCTCCATCCGCGAGCCAGAACCATGTTGTTATTATCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTGATTgtggtg                                                                                                                                      >  1:9193/1‑251 (MQ=255)                                                                                                                aaaTACCTGGGGGCATAACAACAACGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAATTGCAGCGATGCCTTTTTGCAAGC‑TGGGCAAACTAAGTATTTGACCCCGCATAAGGAATAGAACATGCTCCATCCGCGAGCCAGAACCATGTTGTTATTATCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTGATTGTGGtgatg                                                                                                                                   >  1:168478/1‑250 (MQ=255)                                                                                                                 aaTACCTGGGGGCATAACAACAACGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAATTGCAGCGATGCCTTTTTGCAAGC‑TGGGCAAACTAAGTATTTGACCCCGCATAAGGAATAGAACATGCTCCATCCGCGAGCCAGAACCATGTTGTTATTATCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTGATTGTGGtgatga                                                                                                                                  >  1:444915/1‑250 (MQ=255)                                                                                                                        gggggCATAACAGCACGGCCGTAAAGCCAGGCCCCCTGTAAAGTTCAATCCTGGTCACAGAATTACTTTACCGTTGCACCACGAATGACAGCTATCCTTTTTTGCA‑GCGTGGGCAAACTAAGTATTTGTCCCCGCATAGGAAATAAACCAGGCTCCACCCGCGGGCCAGTGCCATGTTGTTATTATCGCTTCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTGATTGTGGTGATGAAAATTGcc                                                                                                                          <  1:306280/251‑1 (MQ=255)                                                                                                                                   aacaacGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAATTGCAGCGATGCCTTTTTGCAAGC‑TGGGCAAACTAAGTATTTGACCCCGCATAAGGAATAGAACATGCTCCATCCGCGAGCAAGAACCATGTTGTTATTATCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTGATTGTGGTGATGAAAATTGCCTCGGTATTAc                                                                                                                >  2:272462/1‑250 (MQ=255)                                                                                                                                    acaacGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCCCCAATTGCAGCGATGCCTTTTTGCAAGC‑TGGGCAAACTAAGTATTTGACCCCGCATAAGGAATAGAACATGCTCCATCCGCGAGCCAGAACCATGTTGTTATTATCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTGATTGTGGTGATGAAAATTGCCTCGGTATTACAg                                                                                                              <  1:443641/251‑1 (MQ=255)                                                                                                                                                    aggcTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAATTGCAGCGATGCCTTTTTGCAAGC‑TGGGCAAACTAAGTATTTGACCCCGCATAAGGAATAGAACATGCTCCATCCGCGAGCCAGAACCATGGTGTTATTATCGCTCCCCGCCGTGGCAATGGGGATTGCGTCCAGTATTATTCTGATTGTGGTGATGAAAATTGCCTCGGTATTACAGACTTTGCTCTGGCAAc                                                                                              >  1:117863/1‑251 (MQ=255)                                                                                                                                                    aggcTCCCTGTAAATATCGATCTGGGTAACACAATTACTTTATCGTTTCAGCACCAATTGCAGCGATGCCTTTTTGCAAGC‑TGGGCAAACTAAGTATTTGACCCCGCATAAGGAATAGAACATGCTCCATCCGCGAGCCAGAACCATGTTGTTATTATCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTGATTGTGGTGATGAAAATTGCCTCGGTATTACAGAATTTGCTCTGGCAAc                                                                                              >  1:346906/1‑251 (MQ=255)                                                                                                                                                              tAAATATCGATCTGGGTCACACCATTAATTTATCGTTCAAGCACAAATTACAGCGATGCCTTTTTGCAAGC‑TGGGCAAACTCAGTATTTGACCCAGCATAAGGACTAGAACATGCTCCATCAGCGAGCCAGAACCATGTTGTTATTATCGCTCCACGCTGTCGCCATTGGGATAGCGTCCAGTCTTATTCTGATTGTGGTGATGAAAATTGCGTCGGTATTCCAGAATTTGCTCTGGCAACGACTGCCggg                                                                                    >  2:440229/1‑251 (MQ=255)                                                                                                                                                                aataccgattggggTCACACATTTACTTTACCGTTTCAGCACCAATTGCAGCTATGCCTTTTTGAAAGG‑TGGGAAAATTAATTATTTGCCCCCGGAAAAGGAATAGAACATGCTCCACCCGCGAGCCAGAACCATGTTGTTATTATCCCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTGATTGTGGTGATGAAAATTGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGaa                                                                                  <  2:346906/240‑1 (MQ=255)                                                                                                                                                                          tGGGTCACACATTTACTTTATCGTTTCAGCCCCAATTGCAGCGATGCCTTTTTGCAAGC‑TGGGCAAACTAAGTATTTGACCCCGCATAAGGAATAGAACATGCTCCATCCGCGAGCCAGAACCATGTTGTTATTATCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTGATTGTGGTGATGAAAATTGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGa                                                                         <  1:377124/250‑1 (MQ=255)                                                                                                                                                                           gggTCACACAATTACTTTATCGTTTCAGCACCAATTGCAGCGATGCCTTTTTGCAAGC‑TGGGCAAACTAAGTATTTGACCCCGCATAAGGAATAGAACATGCTCCATCCGCGAGCAAGAACCATGTTGTTATTATCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTGATTGTGGTGATGAAAATTGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGAt                                                                        >  1:263663/1‑250 (MQ=255)                                                                                                                                                                                     attattTTATCGTTTCAGCACCAATTGCAGCGATGCCTTTTTGCAAGC‑TGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACATGCTCCAACCGCGAGCCAGAACCATGTTGTTATTATCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTGATTGTGGTGATGAAAATCGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGAtt                                                             <  2:427680/251‑1 (MQ=255)                                                                                                                                                                                     aTTACTTTATCGTTTCAGCACCAATTGCAGCGATGCCTTTTTGCAAGC‑TGGGCAAACTAAGTATTTGACCCCGCATAAGGAATAGAACATGCTCCATCCGCGAGCCAGAACCATGTTGTTATTATCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTGATTGTGGTGATGAAAATTGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGAtt                                                             >  2:325017/1‑251 (MQ=255)                                                                                                                                                                                      ttACTTTATCGTTTCAGCACCAATTGCAGCGATGCCTTTTTGCAAGC‑TGGGCAAACTAAGTATTTGACCCCGCATAAGGAATAGAACATGCTCCATCCGCGAGCCAGAACCATGTTGTTATTATCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTGATTGTGGTGATGAAAATTGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAAATCTGGGGATAGCCCAGGATTc                                                            >  1:295574/1‑251 (MQ=255)                                                                                                                                                                                            gaTCGTTTCAGCCCAAATTGCAGCGATGCTTTTTTGCAAGC‑TGGGCAAATTAAGTATTTGACCCACCAAAAGGAATAAAACATGCTCCATCCGCGACCTAGAACCATGTTTTTATTATCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTGATTGTGGTGATGAAAATTGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGAGTCAcccc                                                       <  1:335660/249‑1 (MQ=255)                                                                                                                                                                                              acgtgtcaacACCAATTGCTGCGATTCCTTTTTGCAGGC‑TGGGCAAACTAAGTATTTGACCCCGCATAAGGAATAGAGCTTGCTCCACCCCAGAGCCAGAACCGTTTTGTTATTATCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTGATTGTGGTGATGAAAATTGCCTCGGTATTACAGAATTTGCTCTGGCAACGGCTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCctct                                                    <  1:383688/242‑1 (MQ=255)                                                                                                                                                                                                  ttCAGCACCAATTGCAGCGATGCCTTTTTGCAAGC‑TGGGCAAACTAAGTATTTGACCCCGCATAAGGAATAGAACATGCTCCATCCGCGAGCCAGAACCATGTTGTTATTATCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTGATTGTGGTGATGAAAATTGCCTCGGAATTACAGAATTTGCTCTGG‑AACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGatc                                               >  2:239465/1‑251 (MQ=255)                                                                                                                                                                                                      gcacaaagtgacgcattGCCTTTTTGCAAGC‑GGGCCAAACTAAGTATTTGACCCCGAATAAGTAAAAGAACATGCTCCATCCGCGAGCCAGAACTATGTTGTTATTATCCCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCGTATTCTGATTGTGGTGATGAAAATTGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGatcatc                                            <  2:353493/235‑1 (MQ=255)                                                                                                                                                                                                           aaTTGCAGCGATGCCTTTTTGCAAGC‑TGGGCAAACTAAGAATTTGACCCCGCATAAGGAATAGAACATGCTCCATCCGCGAGCCAGAACCATGTTGTTATTATCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTGATTGTGGTGATGAAAATTGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGgtgt                                       >  1:234885/1‑251 (MQ=255)                                                                                                                                                                                                           aaTTGCAGCGAGGCATTTTTGCAAGC‑TGGGCAAACTAAGTATTTGACCCCGCATAAGGAATAGAACATGCTCCATCCGCGAGCCAGAACCATGATGTTATTATCGATCCACGACGTGGCAATTGGGATTGCGTACAGTCTTATTCTGATTGTGGTGATGAAAATTGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGgtgt                                       >  2:192960/1‑251 (MQ=255)                                                                                                                                                                                                            aTTGCAGCGATGCCTTTTTGCAAGC‑TGGGCAAACTAAGTATTTGACCCCGCATAAGGAATAGAACATGCTCCATCCGCGAGCCAGAACCATGTTGTTATTATCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTGATTGTGGTGATGAAAATTGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGgtgt                                       >  1:368988/1‑250 (MQ=255)                                                                                                                                                                                                                cagcaatcccTTTTTGCAATC‑TGGGCAAACTAAGTATTTGACCCCGCATAAGGAATAGAACATGCCCCATCCGCGAGCCAGAACCATGTTGTTATTATCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTGATTGTGGTGATGAAAATTGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTa                                   <  1:136902/242‑1 (MQ=255)                                                                                                                                                                                                                   cGATGCCTTTTTGCAAGC‑TGGGCAAACTAAGTATTTGACCCCGCATAAGGAATAGAACATGCTCCATCCGCGAGCCAGAACCATGTTGTTATTATCGCTCCCCGCCGTGGCAATTGGGATTGAGTCCAGTCTTATTCTGAATGTGGTGATGAAAATTGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGc                               >  1:50871/1‑251 (MQ=255)                                                                                                                                                                                                                                       gggCAACCTAAGTATTTGCCCCCACTTAAGGAATAGAACATCCTCAATCCGCGAGCCAGAACCGTGTTTTTATTACCCCCCCCCGCCGTGGCAAGTGGGATTGCGTCCAGTCTTATTCTGAATGTGGTGATGAAAATTGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTTGATCGTCGGTGTATTAACGCTAACGGGTATTGCGGTggg            <  2:50871/251‑1 (MQ=255)                                                                                                                                                                                                                                         gCAAACTAAGTATTTGACCCCGCATAAGGAATAGAACATGCTCCATCCGCGAGCCAGAACCATGTTGTTATTATCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTGATTGTGGTGATGAAAATTGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGggt          >  1:96845/1‑251 (MQ=255)                                                                                                                                                                                                                                            aacaaaGTATTTGACCCCGCATAAGAAATAGAACATGCTCCATCCGCGAGCCAGAACCATGTTTTTATTATCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTGATTGTGGTGATGAAAATTGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGggttgg       <  1:53733/247‑1 (MQ=255)                                                                                                                                                                                                                                            aaCTAAGTATTTGACCCCGCATAAGGAATAGAACATGCTCCATCCGCGAGCCAGAACCATGTTGTTATTATCGCTCCCCGCCGTGCCAATTGGGATTGCGTCCAGTCTTATTCTGATTGTGGTGATGAAAATTGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGCACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGggttgg       >  2:374610/1‑251 (MQ=255)                                                                                                                                                                                                                                                aaGTATTTGACCCCGCATAAGGAATAGAACATGCTCCATCCGCGAGCCAGAACCATGTTGTTATTATCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTGATTGTGGTGATGAAAATTGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTAt   >  2:390415/1‑251 (MQ=255)                                                                                                                                                                                                                                                 aGTATTTGACCCCGCATAAGGAATAAAACATGCTCAATCCGCGAGCCAGAACAATTTTGTTATTATCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTGATTGTGGTGATGAAAATTGCCTCGGTATTACAGAATTTGCTCTGGCAATGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATc  <  1:439850/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      ATATCACACAGAGCAAGACGTGCGTTGGTGCCGCCCACATCACCGACTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAGTCAAAATAATTCTTTCTCACACTGTAAATACCTGGGGGCATAACAACAACGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAATTGCAGCGATGCCTTTTTGCAAGC‑TGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACATGCTCCATCCGCGAGCCAGAACCATGTTGTTATTATCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTGATTGTGGTGATGAAAATCGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATC  >  NC_000913/2509362‑2509851

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