Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,181,660	C→T	L189L (CTG→TTG)	yegV →	putative sugar kinase YegV

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,181,660	0	C	T	97.8%	147.6 / NA	45	L189L (CTG→TTG)	yegV	putative sugar kinase YegV
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (20/24);  minor base G (0/1);  total (20/25)
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.99e-01

GGTGTCTGGCGCTGGTTGAGCCGGATGGCGAACGCACTTTTATGTCATTCAGCGGTGTTGAAAATCAGTGGAATCGCCAGTGGCTGGCGCGATTAACCGTTGCGCCTGGCAGCCTGCTCTATTTTTCCGGTTATCAACTGGCCTCGCCCTGCGGCGAATTGTTAGTGGAATGGCTGGAAGAGCTGCAAGACGTG‑ACGCCGTTTATCGATTTTGGCCCACGTATTGGCGATATCCCGGATGCATTACTGGCGCGGATCATGGCCTGTCGACCTTTAGTGTCGCTCAATCGTCAAGAGGCTGAGATTGCCGCCGAACGTTTTGCTTTATCCGCAGAGATAACAACACTTGGCAAGCAATGGCAGGAGAAATTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCATGGTATTTCAGCAACGACGCTTCTGGCTGCATTCCGGCATTTCCAACGCAAGTTGTCGACACCATTGGG  >  NC_000913/2181415‑2181893                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                          ggTGTCTGGCGCTGGTTGAGCCGGATGGCGAACGCACTTTTATGTCATTCAGCGGTGTTGAAAATCAGTGTAATCGCCAGTGGCTGGCGCGATTAACCGTTGCGCCTGGCAGCCTGCTCTATTTTTCCGGTTATCAACTGGCCTCGCCCTGCGGCGAATTGTTAGTGGAATTGCGGGAAGAGCTGCAAGACGTG‑ACGCCGTTTATCGATTTTTGCCCACGTATGGGCTATATCCCGGAGGCATTATTGgc                                                                                                                                                                                                                                       <  2:183483/250‑1 (MQ=255)    gtCTGGCGCTGGTTGAGCCGGATGGCGAACGCACTTTTATGTCATGCAGCGGTGTTGAAAATCAGTGGAATCGCCAGTGGCTGGCGCGATTAACCGTTGCGCCTGGCAGCCTGCTCTATTTTTCCGGTTATCAACTGGCCTCGCCCTGCGGCGAATTGTTAGTGGAATGGCTGGAAGAGCTGCAAGACGTG‑ACGCCGTTTATCGATTTTGGCCCACGTATTGGCGATATCCCGGATGCATTATTGGcgcg                                                                                                                                                                                                                                    <  2:123506/250‑1 (MQ=255)      cTGGCGCTGGTTGAGCCGGATGGCGAACGCACTTTTATGTCATTCAGCGGTGTTGAAAATCAGTGGAATCGCCAGTGGCTGGCGCGATTAACCGTTGCGCCTGGCAGCCTGCTCTATTTTTCCGGTTATCAACTGGCCTCGCCCTGCGGCGAATTGTTAGTGGAATGGCTGGAAGAGCTGCAAGACGTG‑ACGCCGTTTATCGATTTTGGCCCACGTATTGGCGATATCCCGGATGCATTATTGGCGCGGAt                                                                                                                                                                                                                                 >  1:138388/1‑251 (MQ=255)          cgcTGGTTGAGCCGGATGGCGAACGCACTTTTATGTCATTCAGCGGTGTTGAAAATCAGTGGAATCGCCAGTGGCTGGCGCGATTAACCGTTGCGCCTGGCAGCCTGCTCTATTTTTCCGGTTATCAACTGGCCTCGCCCTGCGGCGAATTGTTAGTGGAATGGCTGGAAGAGCTGCAAGACGTG‑ACGCCGTTTATCGATTTTGGCCTACGTATTGGCGATATCCCGGATGCATTATTGGCGCTGATCATg                                                                                                                                                                                                                             <  2:194835/251‑1 (MQ=255)          cgcTGGTTGAGCCGGATGGCGAACGCACTTTTATGTCATTCAGCGGTGTTGAAAATCAGTGGAATCGCCAGTGGCTGGCGCGATTAACCGTTGCGCCTGGCAGCCTGCTCTATTTTTCCGGTTATCAACTGGCCTCGCCCTGCGGCGAATTGTTAGTGGAATGGCTGGAAGAGCTGCAAGACGTG‑ACGCCGTTTATCGATTTTGGCCCACGTATTGGCGATATCCCGGATGCATTATTGGCGCGGATCATg                                                                                                                                                                                                                             <  2:98698/251‑1 (MQ=255)                   aGCCGGATGGCGAACGCACTTTTATGTCATTCAGCGGTGTTGAAAATCAGTGGAATCGCCAGTGGCTGGCGCGATTAACCGTTGCGCCTGGCAGCCTGCTCTATTTTTCCGGTTATCAACTGGCCTCGCCCTGCGGCGAATTGTTAGTGGAATGGCTGGAAGAGCTGCAAGACGTG‑ACGCCGTTTATCGATTTTGGCCCACGTATTGGCGATATCCCGGATGCATTATTGGCGCGGATCATGGCCTGTCg                                                                                                                                                                                                                     >  1:98189/1‑250 (MQ=255)                    gCCGGATGGCGAACGCACTTTTATGTCAATCAGCGGTGTTGAAAAGCAGTGGAATCGCCAGTGGCTGGCTCGATTAACCGTTGCGCCTGGCAGCCTGCTCTATTTTTCCGGTTATCAACTGGCCTCGCCCTGCGGCGAATTGTTAGTGGAATGGCTGGAAGAGCTGCAAGACGTG‑ACGCCGTTTATCGATTTTGGCCCACGTATTGGCGATATCCCGGATGCATTATTGGCGCGGATCATGGCCTGTCGAc                                                                                                                                                                                                                   <  2:236711/251‑1 (MQ=255)                       ggATGGCGACCGCACTTTTATGTCATTCAGCGGTGTTGAAAATCAGTGGAATCGCCAGTGGCTGGCGCGATTAACCGTTGCGCCTGGCCGCCTGCTCTATTTTTCCGGTTATCAACTGGCCTCGCCCTGCGGCGAATTGTTAGTGGAATGGCTGGAAGAGCTGCACGACGTG‑ACGCCGTTTATCGATTTTGGCACACGTATTGGCGATATCCCGGATGCATTATTGGCGCGGATCATGGCCTGTCGACCtt                                                                                                                                                                                                                >  1:314891/1‑251 (MQ=255)                                        ttATGTCATTCAGCGGTGTTGAAAATCAGTGGAATCGCCAGTGGCTGGCGCGATTAACCGTTGCGCCTGGCAGCCTGCTCTATTTTTCCGGTTATCAACTGGCCTCGCCCTGCGGCGAATTGTTAGTGGAATGGCTGGAAGAGCTGCAAGACGTG‑ACGCCGTTTATCGATTTTGGCCCACGTATTGGCGATATCCCGGATGCATTATTGGCGCGGATCATGGCCTGTCGACCTTTAGTGTCGCTCAATCGt                                                                                                                                                                                               >  1:84114/1‑251 (MQ=255)                                                  cAGCGGTGTTGAAAATCAGTGGAATCGCCAGTGGCTGGCGCGATTAACCGTTGCGCCTGGCAGCCTGCTCTATTTTTCCGGTTATCAACTGGCCTCGCCCTGCGGCGAATTGTTAGTGGAATGGCTGGAAGAGCTGCAAGACGTG‑ACGCCGTTTATCGATTTTGGCCCACGTATTGGCGATATCCCGGATGCATTATTGGCGCGGATCATGGCCTGTCGACCTTTAGTGTCGCTCAATCGTCAAGAGGCTg                                                                                                                                                                                     >  2:224852/1‑251 (MQ=255)                                                             aaaatcatgcaaaacGCCCGTGGCTGGGGCGCTGAACCCTTGCGCCTGGCGGCGTGCTCTATTTTTCCGGTTATCAACTGGCCTCGCTCTGCGGCGAATTTTTAGTGGAATGGCTGGAAGAGCTGCAAGTCGGG‑ACGCCGGTTGACGATTTTGGCCCACGTATTGGCGATATCCCGGATGCGTTATTGGCGCGGATCATGGCCTGTCGACCTTTTGTGTCGCTCAATCGTCAAGAGGCTGGGATTgccgcc                                                                                                                                                                          <  2:241880/237‑1 (MQ=255)                                                                  cAGTGGAAGCGCCAGTGGCTGGCGCGATTAACCGTTGCGCCTGGCAGCCTGCTCTATTTTTCCGGTTATCAACTGGCCTCGCCCTGCGGCGAATTGTTAGTGGAATGGCTGGAAGAGCTGCAAGACGTG‑ACGCCGTTTATCGATTTTGGCCCACGTATTGGCGATATCCCGGATGCATTATTGGCGCGGATCATGGCCTGTCGACCTTTAGTGTCGCTCAATCGTCAAGAGGCTGAGATTGCCGCCGAAc                                                                                                                                                                      <  1:422986/250‑1 (MQ=255)                                                                        aagcGCCAGTGGCTGGCGCGATTAACCGTTGCGCCTGGCAGCCTGCTCTATTTTTCCGGTTATCAACTGGCCTCGCCCTGCGGCGAATTGTTAGTGGAATGGCTGGAAGAGCTGCAAGACGTG‑ACGCCGTTTATCGATTTTGGCCCACGTATTGGCGATATCCCGGATGCATTATTGGCGCGGATCATGGCCTGTCGACCTTTAGTGTCGCTCAATCGTCAAGAGGCTGAGATTGCCGCCGAACGTTTTGc                                                                                                                                                               <  1:430426/248‑1 (MQ=255)                                                                                          ggagtAACCGGTGCGCCTGGCAGCCTGCTCTATTTTTCCGGTTATCAACTGGCCTCGCCCTGCGGCGAATTGTTAGTGGAATGGCTGGAAGAGCTGCAAGACGTG‑ACGCCGTTTATCGATTTTGGCCCACGTATTGGCGATATCCCGGATGCATTATTGGCGCGGATCATGGCCTGTCGACCTTTAGTGTCGCTCAATCGTCAAGAGGCTGAGATTGCCGCCGAACGTTTTGCTTTATCCGCAGAGATaac                                                                                                                                             <  2:15043/247‑1 (MQ=255)                                                                                           gATTAACCGTTGCGCCTGGCAGCCTGCTCTATTTTTCCGGTTATCAACTGGCCTCGCCCTGCGGCGAATTGTTAGTGGAATGGCTGGAAGAGCTGCAAGACGTG‑ACGCCGTTTATCGATTTTGGCCCACGTATTGGCGATATCCCGGATGCATTATTGGCGCGGATCATGGCCTGTCGACCTTTAGTGTCGCTCAATCGTCAAGAGGCTGAGATTGCCGCCGAACGTTTTGCTTTATCCGCAGAGATaac                                                                                                                                             >  1:403821/1‑250 (MQ=255)                                                                                                   gTTGCGCCTGGCAGCCTGCTCTATTTTTCCGGTTATCAACTGGCCTCGCCCTGCGGCGAATTGTTAGTGGAATGGCTGGAAGAGCTGCAAGACGTG‑ACGCCGTTTATCGATTTTGGCCCACGTATTGGCGATATCCCGGATGCATTATTGGCGCGGATCATGGCCTGTCGACCTTTAGTGTCGCTCAATCGTCAAGAGGCTGAGATTGCCGCCGAACGTTTTGCTTTATCCGCAGAGATAACAACACTTgg                                                                                                                                    <  1:109967/251‑1 (MQ=255)                                                                                                       cgcCTGGCAGCCTGCTCTATTTGTCCGGTTATCAACTGGCCTCGCCCTGCGGCGAATTGTTAGTGGAATGGCTGGAAGAGCTGCAAGACGTG‑ACGCCGTTTATCGATTTTGGCCCACGTATTGGCGATATCCCGGATGCATTATTGGCGCGGATCATGGCCTGTCGACCTTTAGTGTCGCTCAATCGTCAAGAGGCTGAGATTGCCGCCGAACGTTTTGCTTTATCCGCAGAGATAACAACACTTGgcaag                                                                                                                                <  2:249615/251‑1 (MQ=255)                                                                                                            ggCAGCCTGCTCTATTTTTCCGGTTATCAACTGGCCTCGCCCTGCGGCGAATTGTTAGTGGAATGGCTGGAAGAGCTGCAAGACGTG‑ACGCCGTTTATCGATTTTGGCCCACGTATTGGCGATATCCCGGATGCATTATTGGCGCGGATCATGGCCTGTCGACCTTTAGTGTCGCTCAATCGTCAAGAGGCTGAGATTGCCGCCGAACGTTTTGCTTTATCCGCAGAGATAACAACACTTGgcaagcaa                                                                                                                             <  1:373944/249‑1 (MQ=255)                                                                                                            ggCAGCCTGCTCTATTTTTCCGGTTATCAACTGGCCTCGCCCTGCGGCGAATTGTTAGTGGAATGGCTGGAAGAGCTGCAAGACGTG‑ACGCCGTTTATCGATTTTGGCCCACGTATTGGCGATATCCCGGATGCATTATTGGCGCGGATCATGGCCTGTCGACCTTTAGTGTCGCTCAATCGTCAAGAGGCTGAGATTGCCGCCGAACGTTTTGCTTTATCCGCAGAGATAACAACACTTGGCAAGCAATg                                                                                                                           <  1:77567/251‑1 (MQ=255)                                                                                                                 ccTGCTCTATTTTTCCGGTTATCAACTGGCCTCGCCCTGCGGCGAATTGTTAGTGGAATGGCTGGAAGAGCTGCAAGACGTG‑ACGCCGTTTATCGATTTTGGCCCACGTATTGGCGATATCCCGGATGCATTATTGGCGCGGATCATGGCCTGTCGACCTTTAGTGTCGCTCAATCGTCAAGAGGCTGAGATTGCCGCCGAACGTTTTGCTTTATCCGCAGAGATAACAACACTTGGCAAGCAATGGCAgg                                                                                                                      <  2:316047/251‑1 (MQ=255)                                                                                                                  cTGCTCTATTTTTCCGGTTATCAACTGGCCTCGCCCTGCGGCGAATTGTTCGTGGAATGGCTGGAAGAGCTGCAAGACGTG‑ACGCCGTTTATCGATTTTGGCCCACGTATTGGCGATATCCCGGATGCATTATTGGCGCGGATCATGGCCTGTCGACCTTTAGTGTCGCTCCATCGTCAAGAGGCTGAGATTGCCGCCGAACGTTTTGCTTTATCCGCAGCGATAACAACACTTGGCAAGCAATGGCAggg                                                                                                                     >  1:181611/1‑250 (MQ=255)                                                                                                                     ctctATTTTTCCGGTTATCAACTGGCCTCGCCCTGCGGCGACTTGATAGTGGAATGGCTGGAAGAGCTGCAAGACGTG‑ACGCCGTTTATCGATTTTGGCCCACGTATTGGCGATATCCAGGATGCATTATTGGCGCGGATCATGGCCTGTCGACCTCTAGTGTCGCTCAATCGTCAAGAGGCTGAGATTGCCGCCGAACGTTTTGCTTTATCCGCAGAGATAACAACACTTGGCAAGCAATGGCCGgaga                                                                                                                   >  1:443026/1‑250 (MQ=255)                                                                                                                         aTTTTTCCGGTTATCAACTGGCCTCGCCCTGCGGCGAATTGTTAGTGGAATGGCTGGAAGAGCTGCAAGACGTG‑ACGCCGTTTATCGATTTTGGCCCACGTATTGGCGATATCCCGGATGCATTATTGGCGCGGATCATg                                                                                                                                                                                                                             <  1:84009/140‑1 (MQ=255)                                                                                                                             ttCCGGTTATCAACTGGCCTCGCCCTGCGGCCAATTGTTAGTGGAATGGCTGGAAGAGCTGCAAGACGTG‑ACGCCGTTTATCGATTTTGGCCCACGTATTGGCGATATCCCGGATGCATTATTGGCGCGGATCATGGCCTGTCGACCTTTAGTGTCGCTCAATCGTCAAGAGGCTGAGATTGCCGCCGAACGTTTTGCTTTATCCGCAGAGATAACAACACTTGGCAAGCAATGGCAGGAGAAATTTGCcg                                                                                                          >  1:144828/1‑251 (MQ=255)                                                                                                                                cGGTTATCAAATGGCCTCGCCCTTCGGCGAATTGTTAGTGGAATGGCTGGAAGAGCTGCAAGACGTG‑ACGCCGTTTATCGATTTTGGCCCACGTATTGGCGATCCCCCGGATGCATTATTGGCGCGGATCATGGCCTGTCGACCTTTAGTGTCGCTCAATCGTCAAGAGGCTGAGCTTGCCGCCGAACGTTTTGCTTTATCCGCCGAGATAACAACACTTGGCAAGCACTGGCAGGAGAAATTTGCCgcgc                                                                                                       >  2:106874/1‑251 (MQ=255)                                                                                                                                 ggTTATCAACTGGCCTCGCCCTGCGGCTAATTGTTAGTGGAATGGCTGGAAGAGCTGCAAGACGTG‑ACGCCGTTTATCGATTTTGGCCCACGTATTGGCGATATCCCGGATGCATTATTGGCGCGGATCATGGCCTGTCGACCTTTAGTGTCGCTCAATCGTCAAGAGGCTGAGATTGCCGCCGAACGtt                                                                                                                                                                   >  1:458362/1‑190 (MQ=255)                                                                                                                                 ggTTATCAACTGGCCTCGCCCTGCGGCTAATTGTTAGTGGAATGGCTGGAAGAGCTGCAAGACGTG‑ACGCCGTTTATCGATTTTGGCCCACGTATTGGCGATATCCCGGATGCATTATTGGCGCGGATCATGGCCTGTCGACCTTTAGTGTCGCTCAATCGTCAAGAGGCTGAGATTGCCGCCGAACGtt                                                                                                                                                                   <  2:458362/190‑1 (MQ=255)                                                                                                                                        aaCTGGCCTCGCCCTGCGGCGAATTGTTAGTGGAATGGCTGGAAGAGCTGCAAGACGTG‑ACGCCGTTTATCGATTTTGGCCCACGTATTGGCGATATCCCGGATGCATTATTGGCGCGGATCATGGCCTGTCGACCTTTAGTGTCGCTCAATCGTCAAGAGGCTGAGATTGCCGCCGAACGTTTTGCTTTATCCGCAGAGATAACAACACTTGGCAAGCAATGGCAGGAGAAATTTGCCGCGCCGTTGAtc                                                                                               <  2:403821/251‑1 (MQ=255)                                                                                                                                        aaCTGCCCTCGCCCTGCCGCGAATTGTTAGTGGAGTGGTTGGAAGAGCTGCAAGCCGTGCATGC‑GTTTATCGATTTTGGCCCACGTATTTGCGATATCCCGGATGGGTTGTTGGCGCGTATCATGGCCTGTCGACCTTTAGTGTCGCTCAATCGTCAAGAGGCTGAGATTGCCGCTCAACGTTTTGCTTTATCCGCAGAGAGAACAACACTGGACAAGCAATGGCAGGAGAGATTTGCCGCTCCGTTGAtc                                                                                               <  2:456220/251‑1 (MQ=255)                                                                                                                                         aatGGCCCCGCCCTGCGGCGAATTGTTAGTGGAATGGCTGGAAGAGCTGCAAGACGTG‑ACGCCGTTTATCGATTTTGGCCCACGTATTGGCGATATCCCGGATGCATTATTGGCGCGGATCATGGCCTGACGACCTTTAGTGTCGCTCAATCGTCAAGAGGCTGAGATTGCCGCCGAACGTTTTAGTTTGTCCGCAGAGATAACAACAATTGGCAAGCAAGGGCAGGAGCAATTTGCCTGGCCGTTGAtcg                                                                                              >  2:246960/3‑251 (MQ=255)                                                                                                                                         aCTGGCCTCGCCCTGCGGCGAATTGTTAGTGGAATGGCTGGAAGAGCTGCAAGACGTG‑ACGCCGTTTATCGATTTTGGCCCACGTATTGGCGATATCCCGGATGCATTATTGGCGCGGATCATGGCCTGTCGACCTTTAGTGTCGCTCAATCGTCAAGAGGCTGAGATTGCCGCCGAACGTTTTGCTTTATCCGCAGAGATAACAACACTTGGCAAGCAATGGCAGGAGAAATTTGCCGCGCCGTTGAtcg                                                                                              <  1:13134/251‑1 (MQ=255)                                                                                                                                            ggCCTCGCCCTGCGGCGAATTGTTAGTGGAATGGCTGGAAGAGCTGCAAGACGTG‑TCGCCGTTTATCGATTTTGGCCCACGTATTGGCGATATCCCGGATGCATTAGTGGCTCGGATCATGGCCTGTCGACCTTTAGTGTCGCTCAATCGTCAAGAGGCTGAGATTGCCGCCGAACGTTTTGCTGTATCCGCAGAGATAACAACACTTGGCAAGCAATGGCAGGATAAATTTGCCGCGCCGTTGAtcgttc                                                                                           <  2:89760/251‑1 (MQ=255)                                                                                                                                              ccTCGCCCTGCGGCGAATTGTTAGTGGAATGGCTGGAAGTGCTGCAAGACGTG‑ACGCCGTTTATCGATTTTGGCCCACGTATTGGCGATATCCCGGATGCATTATTGGCGCGGATCATGGCCTGTCGACCTTTAGTGTCGCTCAATCGTCAAGAGGCTGAGATTGCCGCCGAACGTTTTGCTTTATCCGCAGAGATAACAACACTTGGCAAGCAATGGCAGGAGAAATTTGCCGCGCCTTTGAtcgttcgt                                                                                         <  1:34306/251‑1 (MQ=255)                                                                                                                                                  gCCCTGCGGCGAATTGTTAGTGGAATGGCTGGAAGAGCTGCAAGACGTG‑ACGCCGTTTATCGATTTTGGCCCACGTATTGGCGATATCCCGGATGCATTATTGGCGCGGATCATGGCCTGTCGACCTTTAGTGTCGCTCAATCGTCAAGAGGCTGAGATTGCCGCCGCACGTTTTGCTTTATCCGCAGAGATAACAACACTTGGCAAGCAATGGCAGGAGAAATTTGCCGCGCCGTTGATCGTTCGTCTCg                                                                                     >  2:121770/1‑251 (MQ=255)                                                                                                                                                                       ggAATGGCTGGAAGAGCTGCAAGACGTG‑ACGCCGTTTATCGATTTTGGCCCACGTATTGGCGATATCCCGGATGCATTATTGGCGCGGATCATGGCCTGTCGACCTTTAGTGTCGCTCAATCGTCAAGAGGCTGAGATTGCCGCCGAACGTTTTGCTTTATCCGCAGAGATAACAACACTTGGCAAGCAATGGCAGGAGAAATTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCATGGTAtt                                                                <  1:45514/251‑1 (MQ=255)                                                                                                                                                                             gCTGGAAGAGCTGCAAGACGTG‑ACGCCGTTTATCGATTTTGGCCCGCGTATTGGCGATATCCCGGATGCATTATTGGCGCGGATCATGGCCTGTCGACCTTTAGTGTCGCTCAATCGTCAAGAGGCTGAGATTGCCGCCGAACGTTTTGCTTTATCCGCAGAGATAACAACACTTGGCAAGCAATGGCAGGAGAAATTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCATGGTATTTCAGCa                                                          <  2:345266/251‑1 (MQ=255)                                                                                                                                                                                ggtagagCTGCACGCCGTG‑TCGCCGTTTATCGATTTTCGCCCACGCATTGGCGACACCCCGGATGCATTATTGTCCCGCATCATGGTCTGTCGCACTTTAATGTCGCTCAATCGTCAAACGGCTGAGATTGCCGCCGAACGTTTTGCTTTATCCGCAGGGATTAAAACACTTGCCAAGCAATGGCAGGCGAAGTTTGACGCGCCGCTGATCGATCGTCTTGTTACAGAAGCCGCATGGCATTTCAGCTacg                                                       >  2:194688/4‑251 (MQ=255)                                                                                                                                                                                  aaGAGCTGCAAGCCGTG‑ACGCCGTTTATCGATTTTGGCCCACGTATTGGCGATATCCCGGATGCATTATTGGCGCGGATCATGGCCTGTCGACCTTTAGTGTCGCTCAATCGTCAAGAGGCTGAGATTGCCGCCGAACGTTTTGCTTTATCCGCAGAGATAACAACACTTGGCAAGCAATGGCAGGAGAAATTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCATGGTATTTCAGCAacgac                                                     <  2:181611/251‑1 (MQ=255)                                                                                                                                                                                                     cGCCGTTTATCGATTTTGGCCCACGTATTGGCGATATCCCGGATGCATTATTGGCGCGGATCATGGCCTGTCGACCTTTAGTGTCGCTCAATCGTCAAGAGGCTGAGATTGCCGCCGAACGTTTTGCTTTATCCGCAGAGATAACAACACTTGGCAAGCAATGGCAGGAGAAATTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCATGGTATTTCAGCAACGACGCTTCTGGCTGCATTCCg                                   >  1:166851/1‑251 (MQ=255)                                                                                                                                                                                                            tATCGGTTTTGGCCCATGTATCGGCGATCTCCCGGATGCATTATTGGCTCGGATCATGGCCTGTCGACCTTTAGTGTCGCTCAATCGTCAAGAGGCTGAGATTGCCGCCGAACGGTTTGCTTTATCCGCAGAGATAACACCACTTGGCAAGCAATGGCAGGAGAAATTTGCCGCGCCGTTGATCGTTCGTCTCGTTAAAGAAGGCGCCTGGTATTTCTGCAAACACGCTTCTGGCTGCATTCCGGCCTTTc                            >  2:321076/1‑251 (MQ=255)                                                                                                                                                                                                               cGATTTTGGCCCACGCATTGGCGACATCCCGGATGCATTATTGGCGCGCATCATGGCCTGTCGACCTTTATCGTCGCTGAATCGACAACAGGCTGAGTGTGCCGCCGAACGCTTTGCTTTATCCGCAGAGATAACAACACTTGGCACGCAATGGCAGGAGATAATTGCCGCGCCGTTGAACGTTCGTCTCTTTAATGCAGGCGCATGGTATCTTAGCTACGACGCTCCTGGCTGCATTCCGGCATTTcta                          >  2:330605/1‑248 (MQ=255)                                                                                                                                                                                                                          cACGTATTGGCGATATCCCGGATGCATTATTGGCGCGGATCATGGCCTGTCGACCTTTAGTGTCGCTCAATCGTCAAGAGGCTGAGATTGCCGCCGAACGTTTTGCTTTATCCGCAGAGATAACAACACTTGGCAAGCAATGGCAGGAGAAATTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCATGGGATTTCAGCAACGACGCTTCTGGCTGCATTCCGGCATTTCCAACGCAAGTTGTc              >  2:77067/1‑251 (MQ=255)                                                                                                                                                                                                                                 tGGCGATATCCCGGATGCATTATTGGCGCGGATCATGGCCTGTCGACCTTTAGTGTCGCTCAATCGTCAAGAGGCTGAGATTGCCGCCGAACGTTTTGCTTTATCCGCAGAGATAACAACACTTGGCAAGCAATGGCAGGAGAAATTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCATGGTATTTCAGCAACGACGCTTCTGGCTGCATTCCGGCATTTCCAACGCAAGTTGTCGACAcc        >  1:259634/1‑250 (MQ=255)                                                                                                                                                                                                                                      atatCCCGGATGCATTATTGGCGCGGATCATGGCCTGTCGACCTTTAGTGTCGCTCAATCGTCAAGAGGCTGAGATTGCCGCCGAACGTTTTGCTTTATCCGCAGAGATAACAACACTTGGCAAGCAATGGCAGGAGAAATTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCATGGTATTTCAGCAACGACGCTTCTGGCTGCATTCCGGCATTTCCAACGCAAGTTGTCGACACCATTggg  <  1:121770/251‑1 (MQ=255)                                                                                                                                                                                                                                      atatCACGGATGCATTATTGGCGCGGATCATGGCCTGTCGACTTTTAGTGTCGCTCAAGCGTCAAGAGGCTGAGATTGCCGCCGAACGTTTTGCTTTATCCGCAGAGATAACAACACTTGGCAAGCAATGGCAGGAGCACTTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCCTGGTATTTCAGCAACGACGCTTCTGGCTGCATTCCGGCATTTCCAACGCAAGTTGTCGACACCCTTggg  >  2:227778/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                          GGTGTCTGGCGCTGGTTGAGCCGGATGGCGAACGCACTTTTATGTCATTCAGCGGTGTTGAAAATCAGTGGAATCGCCAGTGGCTGGCGCGATTAACCGTTGCGCCTGGCAGCCTGCTCTATTTTTCCGGTTATCAACTGGCCTCGCCCTGCGGCGAATTGTTAGTGGAATGGCTGGAAGAGCTGCAAGACGTG‑ACGCCGTTTATCGATTTTGGCCCACGTATTGGCGATATCCCGGATGCATTACTGGCGCGGATCATGGCCTGTCGACCTTTAGTGTCGCTCAATCGTCAAGAGGCTGAGATTGCCGCCGAACGTTTTGCTTTATCCGCAGAGATAACAACACTTGGCAAGCAATGGCAGGAGAAATTTGCCGCGCCGTTGATCGTTCGTCTCGATAAAGAAGGCGCATGGTATTTCAGCAACGACGCTTCTGGCTGCATTCCGGCATTTCCAACGCAAGTTGTCGACACCATTGGG  >  NC_000913/2181415‑2181893

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