Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,151,298	G→A	G210D (GGC→GAC)	fabG →	3‑oxoacyl‑[acyl‑carrier‑protein] reductase FabG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,151,298	0	G	A	95.5%	142.9 / NA	45	G210D (GGC→GAC)	fabG	3‑oxoacyl‑[acyl‑carrier‑protein] reductase FabG
Reads supporting (aligned to +/- strand):  ref base G (1/0);  major base A (15/27);  minor base C (2/0);  total (18/27)
Fisher's exact test for biased strand distribution p-value = 1.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCGTCATGGTCGTATTATCACTATCGGTTCTGTGGTTGGTACCATGGGAAATGGCGGTCAGGCCAACTACGCTGCGGCGAAAGCGGGCTTGATCGGCTTCAGTAAATCACTGGCGCGCGAAGTTGCGTCACGCGGTATTACTGTAAACGTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACA‑TGGTCTGACCGCGATTTGCACAAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCGGGATTT  >  NC_000913/1151052‑1151539                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   aGCGTCATGGTCGTATTATCACTATCGGTTCTTTGGTTGGTACCATGGGAAATGGCGGTCAGGCCAACTACGCTGCGGCGAAAGCGGGCTTGATCGGCTTCAGTAAATCACTGGCGCGCGAAGTTGCGTCACGCGGTATTACTGTAAACGTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCgacggc                                                                                                                                                                                                                                                <  1:231859/251‑1 (MQ=255)     tCATGGTCGTATTATCACTATCGGTTCTGTGGTTGGTACCATGGGAAATGGCGGTCAGGCCAACTACGCTGCGGCGAAAGCGGGCTTGATTGGCTTCAGTAAATCACTGGCGCGCGAAGTTGCGTCACGCGGTATTACTGTAAACGTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGGGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGACGGCGcac                                                                                                                                                                                                                                            >  2:398710/1‑251 (MQ=255)                       tATCGGTTCTGTGGTTGGTACCATGGGAAATGGCGGTCAGGCCAACTACGCTGCGGCGAAAGCGGGCTTGATCGGCTTCAGTAAATCACTGGCGCGCGAAGTTGCGTCACGCGGTATTACTGTAAACGTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGACGGCGCACAGGAAATCGCCAACGCgg                                                                                                                                                                                                                          >  1:453968/1‑251 (MQ=255)                                 gtgGTTGGTACCATGGGAAATGGCGGTCAGGCCAACTACGCTGCGGCGAAAGCGGGCTTGATCGGCTTCAGTAAATCACTGGCGCGCGAAGTTGCGTCACGCGGTATTACTGTAAACGTTGTTGCTCCGGGCTTTATTGAAACGGACATGGCACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGACGGCGCACAGGAAATCGCCaa                                                                                                                                                                                                                               <  1:71/236‑1 (MQ=255)                                 gtgGTTGGTACCATGGGAAATGGCGGTCAGGCCAACTACGCTGCGGCGAAAGCGGGCTTGATCGGCTTCAGTAAATCACTGGCGCGCGAAGGTGCGTCACGCGGTATTACTGTAAACGTTGTTGCTCCGGGCTTTATTGAAACGGACATGGCACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGGTCCTGCGGGTCGCCTCGGCGGCGCACCGGAAATCGCCaa                                                                                                                                                                                                                               >  2:71/1‑236 (MQ=255)                                      tggtACCATGGGAAATGGCGGTCAGGCCAACTACGCTGCGGCGAAAGCGGGCTTGATCGGCTTCAGTAAATCACTGGCGCGCGAAGTTGCGTCACGCGGTATTACTGTAAACGTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGACGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCa                                                                                                                                                                                                            <  1:248686/250‑1 (MQ=255)                                        gtACCATGGGAAATGGCGGTCAGGCCAACTACGCTGCGGCGAAAGCGGGCTTGATCGGCTTCAGTAAATCACTGGCGCGCGAAGTTGCGTCACGCGGTATTACTGTAAACGTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGACGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATcc                                                                                                                                                                                                         <  1:202594/251‑1 (MQ=255)                                                 gAAATGGGGGTCAGGACAACTACGCTGCGGCGAAAGCGGGATTGATCGGCTTCAGTAAATCACTGGCGCGCGAAGTTGCGTCACGCGGTATTACTGTAAACGCTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCATCGTGCGGGTATGCTGGCGCAGGTTCCTGCGGGTCGCCTCGACGGCGCACAGGAAATCGCCAACGCGGTTGCATTTCTGGCATCCGACGAagca                                                                                                                                                                                                <  2:19157/251‑1 (MQ=255)                                                  aaaTGGCGGTCAGGCCAACTACGCTGCGGCGAAAGCGGGCTTGATCGGCTTCAGTAAATCACTGGCGCGCGAAGTTGCGTCACGCGGTATTACTGTAAACGTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGACGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAagcag                                                                                                                                                                                               <  1:400841/251‑1 (MQ=255)                                                         ggTCAGGCCAACTACGCTGCGGCGAAAGCGGGCTTGATCGGCTTCAGTAAATCACTGGCGCGCGAAGTTGCGTCACGCGGTATTACTGTAAACGTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGACGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACAt                                                                                                                                                                                        >  1:4410/1‑251 (MQ=255)                                                              ggacAACTACGCTGCGGCGAAAGCGGGCTGGATCGGTTTCAGTAAATCACTGGCGCGTGAAGTTGCGTCACGCGGTATTGCTGTAAACGTGGTTGCTCCGGGATTTAGTGAAACGGACAGGACACTTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCGGGTTCCTGCGGGTCGCCTCGACGGCGCACACGAAATCGCCAACGCGGTTGCATTCCTGGCTTCCGACTAATCCGCTTACagcccgt                                                                                                                                                                                   <  2:290400/248‑7 (MQ=255)                                                                 cAACTACGGTGAGGCGAAAGCGGGCGTGATCGGCTTCAGTAAATCACTGGCGCGCGAAGTTGCGTCACGCGGTATTACTGTAAACGTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGCGGTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGACGGCGCACAGGAAATCGCGAACGCGGTTGCAGTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTg                                                                                                                                                                                <  2:219094/251‑1 (MQ=255)                                                                       cGCTGCGGCGAAAGCGGGCTTGATCGGCTTCAGTAAATCACTGGCGCGCGAAGTTGCGTCACGCGGTATTACTGTAAACGTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGACGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACtt                                                                                                                                                                          >  1:137596/1‑251 (MQ=255)                                                                        gCTGCGGCGAAAGCGGGCTTGATCGGCTTCAGTAAATCACTGGCGCGCGAAGTTGCGTCACGCGGTATTACTGTAAACGTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGACGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACttt                                                                                                                                                                         <  1:351416/251‑1 (MQ=255)                                                                             ggcgAAAGCGGGCTTGATCGGCTTCAGTAAATCACTGGCGCGCGAAGTTGCGTCACGCGGTATTACTGTAAACGTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGACGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCAtg                                                                                                                                                                    >  2:206420/1‑251 (MQ=255)                                                                              gcgAAAGCGGGCTTGATCGGCTTCAGTAAATCACTGGCGCGCGAAGTTGCGTCACGCGGTATTACTGTAAACGTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGACGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCAtgt                                                                                                                                                                   <  1:306831/251‑1 (MQ=255)                                                                                      gggCTTGATCGGCTTCAGTAAATCACTGGCGCGCGAAGTTGCGTCACGCGGTATTACTGTAAACGTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGACGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAAcggcg                                                                                                                                                           <  1:121606/251‑1 (MQ=255)                                                                                                    tCAGTAAATCACTGGCGCGCGAAGTTGCGTCACGCGGTATTACTGTAAACGTTGTTGCTCCGGGCTTTATTGAAAAGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGAAGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAACCTTTGCATGTGAAAGGCGGGATGTAAA‑TGGTc                                                                                                                                            >  2:264834/1‑251 (MQ=255)                                                                                                     cAGTAAATCACTGGCGCGCGAAGTTGCGTCACGCGGTATTACTGTAAACGTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGACGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACA‑TGGTc                                                                                                                                            <  1:206420/250‑1 (MQ=255)                                                                                                         aaaTCACTGGCGCGCGAAGTTGCGTCACGCGGTATTACTGTAAACGTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGACGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACA‑TGGTCTGAcc                                                                                                                                       <  2:182054/251‑1 (MQ=255)                                                                                                                 ggCGCGCGAAGTTGCGTCACGCGGTATTACTGTAAACGTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGACGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACA‑TGGTCTGACCGCGATTTg                                                                                                                               <  1:361278/251‑1 (MQ=255)                                                                                                                  gcgagcgAAGTTGCGTCACGCGGTATTACTGTAAACGTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTGCCTGCGGGTCGCCTCGACGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGTTGAAACTTTGCATGTGAACGGCGGGATGTACA‑TGGTCTGACCGCGATTTGc                                                                                                                              <  2:279485/247‑1 (MQ=255)                                                                                                                    gcgcgAAGTTGCGTCACGCGGTATTACTGTAAACGTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGACGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACA‑TGGTCTGACCGCGATTTGcac                                                                                                                            <  2:4410/251‑1 (MQ=255)                                                                                                                        gAAGTTGCGTCACGCGGTATTACTGTAAACGTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGACGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACA‑TGGTCTGACCGCGATTTTCACaaaa                                                                                                                        <  2:215540/251‑1 (MQ=255)                                                                                                                          aGTTGCGTCACGCGTTATTACTGTAAACGTTGTTGCTCCGGGCGTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGACGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATGCGACGAAGCAGCTTACATCACGGGTGAACCTTTGCATGGGAACGGCTGGATGTACA‑GGGTCTGACCGCGATGTGCACAAAATg                                                                                                                      <  2:213609/251‑1 (MQ=255)                                                                                                                               catcaagcggaataaCTGTAAACGTAGCTCCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTCGCGCAGGTTCCTGCGGGTCGCCTCGCCGGCGCACCGGAACTCGCCAACGCGGTCGCATTCCTGGCATCCGCCGAAGCAGCCTACATCACGGGTGAAACTTTGCATGTCAACGGCGGGATGCACAATGTTC‑GACCGCGGTTTGCACAAAATGCTcct                                                                                                                 >  2:65109/15‑249 (MQ=255)                                                                                                                                        gTATTACTGTAAACGTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGACGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACA‑TGGTCTGACCGCGATTTGCACAAAATGCTCATGTTgcgcgc                                                                                                        >  1:254251/1‑251 (MQ=255)                                                                                                                                                          ttgCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGACGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGGGAAACTTTGCATGTGAACGGCGGGATGTACA‑TGGTCTGACCGCGATTTGCACAAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATg                                                                                        >  2:127937/1‑249 (MQ=255)                                                                                                                                                           tgCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGACGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACA‑TGGTCTGACCGCGATTTGCACAAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAg                                                                                      <  1:127937/250‑1 (MQ=255)                                                                                                                                                               ccGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGACGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACA‑TGGTCTGACCGCGATTTGCACAAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGtaata                                                                                 <  2:24686/251‑1 (MQ=255)                                                                                                                                                                cGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGACGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGGGAAACTTTGCATGTGAACGGCGGGATGTACA‑TGGTCTGACCGCGATTTGCACAAAATGCTCATGTTTCGCGCAGTCTGCGGGGTTCTGAGtaata                                                                                 >  2:198309/1‑250 (MQ=255)                                                                                                                                                                 ggtcTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGAGTCGCCTCGACGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAATCTTTGCATGTGAACGTCGGGATGTACA‑TGGTCTGACCGCGATTTGCACAAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGAaataat                                                                               <  2:436368/248‑1 (MQ=255)                                                                                                                                                                 gggCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGACGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACA‑TGGTCTGACCGCGATTTGCACAAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTaataat                                                                               <  2:453968/251‑1 (MQ=255)                                                                                                                                                                 gggCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGACGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACA‑TGGTCTGACCGCGATTTGCACAAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTaataat                                                                               >  1:262299/1‑251 (MQ=255)                                                                                                                                                                    cTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGGACCTGCGGGTCGCCTCGACGGCGCACAGGAAATCGCCAACGCGGATGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACA‑TGGTCTGACCGCGATTTGCACAAAATGCTCATGTTGCGCGCAGTCTGCGAGGTTATGAGTAATAAtt                                                                              >  1:277691/1‑249 (MQ=255)                                                                                                                                                                                 gACATGCCACGTGCGCAGAGCGATCACCAGCGTGCGGGTATCGCGGCGCAGGATTCTGCGGGTCGCCTCGACGGCGCACAGGAAATCGCCGACGCGGTTGCATTCCTGGCATCCGACGCCGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGGTCTACA‑TGGTCTGACCGCGATTTGCACAAAATGCTCATGTTGCGCGCAGGCTGCGTCGGTATGAGTCATAACTAGTGCAAAATGAttt                                                               >  2:265856/1‑251 (MQ=255)                                                                                                                                                                                           gTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGACGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACA‑TGGTCTGACCGCGATTTGCACAAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTgg                                                     <  1:213188/251‑1 (MQ=255)                                                                                                                                                                                                    gCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGACGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACA‑TGGTCTGACCGCGATTTGCACAAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGTAAgg                                            >  2:181294/1‑251 (MQ=255)                                                                                                                                                                                                    gCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGACGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACA‑TGGTCTGACCGCGATTTGCACAAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGTAAgg                                            <  1:198309/251‑1 (MQ=255)                                                                                                                                                                                                     cGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGACGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACA‑TGGTCTGACCGCGATTTGCACAAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGTAAGGc                                           >  2:92801/1‑251 (MQ=255)                                                                                                                                                                                                                   cGGGTATCCTGGCGCAGGTTCCTGCGGGACGCCTCGACGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACA‑TGGTCTGACCGCGATTTGCACAAAATTCTCATGTTGCGCGCAGTCTTCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGTAAGGCCTCAAAATAACGt                              >  1:51077/1‑250 (MQ=255)                                                                                                                                                                                                                              gcgcAGGTTCCTGCGGGTCGCCTCGACGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACA‑TGGTCTGACCGCGATTTGCACAAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTa                  <  2:210531/251‑1 (MQ=255)                                                                                                                                                                                                                                      tCCTGCGGGACGCGTCGCCGGCGCACAGGAAATCGCCACCGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACA‑TGGTCTGACCGCGATTTGCACAAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGAAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGCAAGGCCTCAAAATAACGTAAAATCGTGGTCAGACCTGc          >  2:292901/1‑251 (MQ=255)                                                                                                                                                                                                                                          gCGGGTCGCCTCGACGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACA‑TGGTCTGACCGCGATTTGCACAAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCggg      <  2:172091/251‑1 (MQ=255)                                                                                                                                                                                                                                          gCGGGTCGCCTCGACGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACA‑TGGTCTGACCGCGATTTGCACAAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCggg      <  1:274760/251‑1 (MQ=255)                                                                                                                                                                                                                                              gTCGGCTCGACGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACA‑TGGTCTGACCGCGATTTGCACAAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGTGTGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACTTGCCGGGAttt  <  1:40016/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   AGCGTCATGGTCGTATTATCACTATCGGTTCTGTGGTTGGTACCATGGGAAATGGCGGTCAGGCCAACTACGCTGCGGCGAAAGCGGGCTTGATCGGCTTCAGTAAATCACTGGCGCGCGAAGTTGCGTCACGCGGTATTACTGTAAACGTTGTTGCTCCGGGCTTTATTGAAACGGACATGACACGTGCGCTGAGCGATGACCAGCGTGCGGGTATCCTGGCGCAGGTTCCTGCGGGTCGCCTCGGCGGCGCACAGGAAATCGCCAACGCGGTTGCATTCCTGGCATCCGACGAAGCAGCTTACATCACGGGTGAAACTTTGCATGTGAACGGCGGGATGTACA‑TGGTCTGACCGCGATTTGCACAAAATGCTCATGTTGCGCGCAGTCTGCGTGGTTATGAGTAATAATTAGTGCAAAATGATTTGCGTTATTGGGGGGTAAGGCCTCAAAATAACGTAAAATCGTGGTAAGACCTGCCGGGATTT  >  NC_000913/1151052‑1151539

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