Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,193,707	A→G	G166G (GGT→GGC)	thiE ←	thiamine phosphate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,193,707	0	A	G	96.2%	181.1 / NA	52	G166G (GGT→GGC)	thiE	thiamine phosphate synthase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (34/16);  minor base T (1/1);  total (35/17)
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.86e-01

CCGCCCAGACCGATAATCAGCACCTGGCTGTCGAGCAGTTTTTGCTGCCCGTCCAGAGCGATATCGTCGAGCAGGATTTGGCGGCTATAACGCATAAAGTCACGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATACCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGT‑CGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCC  >  NC_000913/4193464‑4193950                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     ccgcccAGTCCGATAAGCAGCGCCTGGCTGTCGAGCAGTTTTTGCTGCCCGTCTAGAGCGATATCGTAGAGCAGGCTTTGGCGGCTATAACTCATAAAGTCACGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGATTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCaa                                                                                                                                                                                                                                               <  2:270783/251‑1 (MQ=255)        gACCGATAATCAGCGCCTGGCTGTCGAGCAGTTTTTGCGGCCCGTCCAGTGCGATATCGGCGAGCAGGATGTGGCGGCTATAACGCATAAAGTCACGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACGGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATTCCGCCAAGCGCCAc                                                                                                                                                                                                                                        <  2:293518/251‑1 (MQ=255)          ccGATAATCAGCACCTGGCTGTCGAGCAGTTTTTGCTGCCCGTCCAGAGCGATATCGTCGAGCAGGATTTGGCGGCTATAACGCATAAAGTCACGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACgg                                                                                                                                                                                                                                      >  1:269493/1‑251 (MQ=255)                    gCACCTGGCTGTCGAGCAGTTTTTGCTGCCCGTCCAGAGCGATATCGTCGAGCAGGATTTGGCGGCTATAACGCATAAAGTCACGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTCATGGCGCTGACGACGGCGATACTGCCGACACCCGTTTCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCCCCACGGGGGGATaagc                                                                                                                                                                                                                            >  2:260555/1‑249 (MQ=255)                     cACCTGGCTGTCGAGCAGTTTTTGCTGCCCGTCCAGAGCGATATCGTCGAGCAGGATTTGGCGGCTATAACGCATAAAGTCACGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATcc                                                                                                                                                                                                                           >  2:159964/1‑251 (MQ=255)                      aCCTGGCTGTCGAGCAGTTTTTGCTGCCCGTCCAGAGCGATATCGTCGAGCAGGATTTGGCGGCTATAACGCATAAAGTCACGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCg                                                                                                                                                                                                                          <  1:274556/251‑1 (MQ=255)                                 gAGCAGTTTTTACTGCCCGTCCAGAGCGATATCGTCGTGCAGGATTTGGCGGCTATAACGCATAAAGTCACGGTCATTCATCGCCAGCTCCTGCAATTGCCATCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCTACACCCGTTGCTATCACCGCAGGCGCTCGTGCCAGACTGAGGCCGCCAATCGCCTCCGTGGGATCAGCCGTTAGTCGCTCa                                                                                                                                                                                                               <  2:291529/251‑1 (MQ=255)                                            gctgcCCGTCCAGAGCGATATCGTCGAGCAGGATTTGGCGGCTATAACGCATAAAGTCACGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACCCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGc                                                                                                                                                                                                    >  1:127433/1‑251 (MQ=255)                                               gcCCGTCCAGAGCGATATCGTCGAGCAGGATTTGGCGGCTATAACGCATAAAGTCACGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCa                                                                                                                                                                                                  <  1:215463/250‑1 (MQ=255)                                                cccGTCCAGAGCGATATCGTCGAGCAGGATTTGGCGGCTATAACGCATAAAGTCACGGTCATTCATCGCGAACTCCGGCAATTTCCAGCAACTGTGCCTTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGGTAATCCGCCAGTCGCTCAACATGCCGTGCCAGc                                                                                                                                                                                                <  2:14086/251‑1 (MQ=255)                                                                  gtcgAGCAGGATTTGGCGGCTATAACGCATAAAGTCACGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGt                                                                                                                                                                              >  1:154242/1‑251 (MQ=255)                                                                       gCAGGATTTGGCGGATATAACGCATAAAGTCACGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGc                                                                                                                                                                            >  1:251053/1‑248 (MQ=255)                                                                            aTTTGGCGGCTATAACGCATAAAGTCACGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGc                                                                                                                                                                     >  1:306926/1‑250 (MQ=255)                                                                                     cTATAACGCATAAAGTCACGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTgtttgg                                                                                                                                                           >  1:249877/1‑251 (MQ=255)                                                                                      tataACGCATAAAGTCACGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTgtttggt                                                                                                                                                          >  2:202062/1‑251 (MQ=255)                                                                                       ataACGCATAAAGTCACGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTgtttggtt                                                                                                                                                         >  2:18725/1‑251 (MQ=255)                                                                                       ataACGCATAAAGTCACGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTgtttggtt                                                                                                                                                         <  2:249877/251‑1 (MQ=255)                                                                                        taACGCATAAAGTCACGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTgtttggttt                                                                                                                                                        >  1:302848/1‑251 (MQ=255)                                                                                           cGCATAAAGTCACGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCg                                                                                                                                                     >  1:25607/1‑251 (MQ=255)                                                                                                aaaGTCACGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCggg                                                                                                                                                >  1:67671/1‑251 (MQ=255)                                                                                                  ggtcccGGTCATTCATCCCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCAAAACGCCAGTCTGCGCCTTGAGTAGTGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCGAGACTGATGCCGCCAATCGCCCCGGTGGGATAATCCGCCAGTCGCTCAACATGCCCACCCAGCGGTTCCCGGCCCTGCGGCGCAGAAGGCATCAGTTTGGATTGCGTCgggc                                                                                                                                               >  2:198503/6‑249 (MQ=255)                                                                                                        ggTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATg                                                                                                                                        >  1:138491/1‑251 (MQ=255)                                                                                                           cattcatCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGCTGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCCGCTGTTCCCGCCCCTGCGGAGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACCCATGTct                                                                                                                                     >  2:163950/1‑250 (MQ=255)                                                                                                                atCGCCAACTCCTGCAATTTCCAGCAACTGTGCCTTGGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGCCGGCGATACTGCCGACACCCGTTGCTATCCCCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGGGGGATAATCCGCCAGGCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAgcg                                                                                                                                <  2:96702/251‑1 (MQ=255)                                                                                                                atCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGACGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACCTGTCTCAgcg                                                                                                                                >  1:184678/1‑251 (MQ=255)                                                                                                                     cAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATCCTGACGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGat                                                                                                                            >  1:305161/1‑250 (MQ=255)                                                                                                                        cTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGAtata                                                                                                                         >  2:39603/1‑250 (MQ=255)                                                                                                                          ccTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGAtata                                                                                                                         >  2:76824/1‑248 (MQ=255)                                                                                                                          ccTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGCGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGCAGGCATCTGTTTGGTTTGCGTCGGGCACACATGTCTCCGCGCGagctaaa                                                                                                                       >  2:186962/1‑244 (MQ=255)                                                                                                                          ccTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAAAGGCGCTGAAGACGGCGATCCTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCCGACTGATGACGCCAATCGCCACGGTGGGGTAATCCGCCAGTCGCTAAGCATGGCGTGCCAGCTGTTACAGCCCCTGAGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGAtata                                                                                                                         >  1:26007/1‑248 (MQ=255)                                                                                                                                    ccAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGCCTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCCGCCCCCGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACCCATGTCTCAGCGCGATATAAGAGAGGCGtgct                                                                                                            >  2:86795/1‑251 (MQ=255)                                                                                                                                          aCTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTGGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAgc                                                                                                      <  1:39603/251‑1 (MQ=255)                                                                                                                                             gtgCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGCCGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGcgcg                                                                                                   >  1:108678/1‑251 (MQ=255)                                                                                                                                              tgCCGTTGCCAAACGCCAGTCCGCGGCTTAAGTAATGGAGCTGACGACGACGATAATGCCGACACCCGTTGCTATCACCGCAGGCGCACGTGCCAGACTGATTCCGCCAATCGCCACGGTGGGATAATCCGCCAGCCCCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGCGCAGAAGGCATCTGTTCGGTTTGCGTCGGGAACACATGTCTCAGCCCGATAAAAACGGGGCGcgctgacaacccg                                                                                                   >  2:280752/1‑242 (MQ=255)                                                                                                                                                  gTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTc                                                                                              <  1:189432/251‑1 (MQ=255)                                                                                                                                                   ttGCCAACCGCCAGTCTGCGGCTTGAGTAATGGCGCTGCCGACGGCGATACTTCCAACACCCGTTCCTATCCCCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGTGAACACATGTCGCATCGCGATATAAGAGGGGCGTTCTGCCAGCGCGACGTc                                                                                              <  2:306926/250‑1 (MQ=255)                                                                                                                                                    tGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCg                                                                                             <  1:202062/250‑1 (MQ=255)                                                                                                                                                       cAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGAttt                                                                                         >  1:69343/1‑251 (MQ=255)                                                                                                                                                          aCGCCAGTCTGCGGCTTGCGTAATGGCGCTGACGACGAAGATACTGCCGCCACACGTTGCTATCACCGCAGGCGCGCGTACCACACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGCGCCAGCTGTTCCAGCCCCTGCGGCGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGCGCTGCCAGCGCGACGTCGATTTccc                                                                                      >  1:201318/1‑250 (MQ=255)                                                                                                                                                                gTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGt                                                                                >  1:41364/1‑251 (MQ=255)                                                                                                                                                                             gTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGT‑CGAAAc                                                                  <  1:293306/251‑1 (MQ=255)                                                                                                                                                                                     gcTGACGCCGGCGCTACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATCTAAGAGTGGCGTGCTGCCAACGCGACGTCGATTTCCATATCGTCATG‑GTGCGAAACGCCCAGcc                                                          >  2:273903/1‑251 (MQ=255)                                                                                                                                                                                       tGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGT‑CGAAACGCCCAGCCGc                                                        >  1:243159/1‑251 (MQ=255)                                                                                                                                                                                        gacgacGGCGATACTGCCGACCCCCGTTGCTATCACCCCAGACGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGT‑CGAAACGCCCAGcctcc                                                       >  2:132385/1‑248 (MQ=255)                                                                                                                                                                                                    aCTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGT‑CGAAACGCCCAGCCGCAGGCCTGCcgcgc                                           >  1:218008/1‑251 (MQ=255)                                                                                                                                                                                                         cGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGT‑CGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATg                                      <  2:277719/251‑1 (MQ=255)                                                                                                                                                                                                          gACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGT‑CGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATgg                                     <  1:263616/251‑1 (MQ=255)                                                                                                                                                                                                                   ggCTATCACCGCAGGCGCGCGTGACAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGATTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGT‑CGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTgaga                             <  2:211321/249‑1 (MQ=255)                                                                                                                                                                                                                    gCTATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGT‑CGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTg                                >  2:281326/1‑246 (MQ=255)                                                                                                                                                                                                                      tATCACCGCAGGCGCGCGTGCCAGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGT‑CGAAACGCCCAGCCGCAGGCCTGCCGCGCGGCTGGCCTTg                                >  2:61407/1‑244 (MQ=255)                                                                                                                                                                                                                                            aGACTGATGCCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGCTTTCCATATCGTCATGTGT‑CGAAACGCCCAGCCGCAGGCCTGCCGCGCGGCTGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGccc   >  2:74942/1‑251 (MQ=255)                                                                                                                                                                                                                                             gACTGATGCCGCCAACCGCCCCGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGT‑CGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGcccc  <  1:163950/251‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     CCGCCCAGACCGATAATCAGCACCTGGCTGTCGAGCAGTTTTTGCTGCCCGTCCAGAGCGATATCGTCGAGCAGGATTTGGCGGCTATAACGCATAAAGTCACGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATACCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGT‑CGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCC  >  NC_000913/4193464‑4193950

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