Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,639,005	G→C	P276P (CCG→CCC)	creD →	putative inner membrane protein CreD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,639,005	0	G	C	98.0%	163.1 / NA	51	P276P (CCG→CCC)	creD	putative inner membrane protein CreD
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (27/23);  minor base A (1/0);  total (28/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTCTACCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATACCGGTTTTACCGTGGCATGGATAATCG  >  NC_000913/4638758‑4639243                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                               aCCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTT‑CAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCTGCTGGGAAAACTTCCCCGCGt                                                                                                                                                                                                                                            >  1:350631/1‑251 (MQ=255)      aaGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAg                                                                                                                                                                                                                                        <  1:142449/251‑1 (MQ=255)                        aaTTTAAGCGATACCGGCGATCTTACTGGGGGGCCTGCCGGGAGTACTAGCGAAATGACCTTAACCAGCAATTGGACGCATCCCAGTTTTTTAGGCGATTTTTTGCCAGCCCAACGGGAAGTTAGCGAGTCCGGTTTTCACGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGAGTGCGTTTTGCTTCAGGCAATGCTACCGGCTGGGCAAACTTCCCAGCGTTTCGCGTCCCAGtggcggcgcc                                                                                                                                                                                                                      >  2:553631/1‑242 (MQ=255)                         aTTTAAGCGGTACCGGTGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTGTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTGCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCa                                                                                                                                                                                                                     <  2:431291/251‑1 (MQ=255)                             aaGCGGTACCGGCGATCTTTCTGGGGGGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACGGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGcc                                                                                                                                                                                                                  <  2:274633/250‑1 (MQ=255)                                         cGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCCGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAACCTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAAt                                                                                                                                                                                                     >  2:593921/1‑251 (MQ=255)                                             atgtCTGTGGTGCCTGCCGTGCGTAATAGCGAAATGACCGTAACCAGCAACTGGACGCATCCCATTTTTTTAGGTGATGTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGTCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGACACCGGCTGGGAAAACTTCCCCGCGTTTAGGTTCGCAGTGACGACGCCAGCCGATCAAGACCAATTaa                                                                                                                                                                                                  <  2:450811/247‑1 (MQ=255)                                               ttCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTaa                                                                                                                                                                                                  <  1:593921/248‑1 (MQ=255)                                                   gtggtgCCTGCCGGGCGTAATAGCGAAAAGGCCTTACCCAGCAACTGGCCGCATCCCAGTTTTTTAGTGGATTTTTTGCCGGCCAAGCGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCg                                                                                                                                                                                           <  2:135426/251‑1 (MQ=255)                                                     ggtgCCTGCCGGGCGTAATAGCGAAAGGACCTTAACCAGCAACGGGCCGCATCCCAGTTTTGTAGGTGATTTTTTGCCAGCCAAACGGTAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAAAACCAATTAACGGACCggg                                                                                                                                                                                         <  1:435282/251‑1 (MQ=255)                                                            gccgggcgTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTa                                                                                                                                                                                   >  1:460339/3‑250 (MQ=255)                                                               gggcgTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTCGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGAACGGGCGACTAAGt                                                                                                                                                                                >  1:311968/1‑250 (MQ=255)                                                                             aaggACCGTAACCAGCAACTGGCGGCAACCCAGTTTTTGAGGTGAGTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAAGGGCAAAGCAGCTGGTGGGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAAGGATACCGGCTGGGAAAACTCCCCCGCGTTTAGCGCCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTTCGCCATTCTGCTGa                                                                                                                                                                 <  2:229075/248‑1 (MQ=255)                                                                              aTGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGAt                                                                                                                                                                <  2:160885/251‑1 (MQ=255)                                                                                     tAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACt                                                                                                                                                          <  1:155259/250‑1 (MQ=255)                                                                                               cTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCCCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGCCGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTCAGTACGCCATTCTGTTGATTGCACTGACTTTTATgg                                                                                                                                               >  2:192898/1‑251 (MQ=255)                                                                                                                   tAGGTGATTTTTTGCCAGCCAAACGGGACGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCCGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCACTACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACg                                                                                                                           >  2:309403/1‑251 (MQ=255)                                                                                                                   tAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCCGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGCCGACTAAGTACGCCCTTCTGCTGATTTCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACg                                                                                                                           >  2:126720/1‑251 (MQ=255)                                                                                                                          ttttttgccAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACcg                                                                                                                    >  1:189617/8‑251 (MQ=255)                                                                                                                               tgccAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCa                                                                                                                >  1:482898/3‑250 (MQ=255)                                                                                                                                        aaCGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTacac                                                                                                      >  2:220106/1‑251 (MQ=255)                                                                                                                                         aCGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACAcc                                                                                                     <  1:338137/251‑1 (MQ=255)                                                                                                                                         aCGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACAcc                                                                                                     >  1:262924/1‑251 (MQ=255)                                                                                                                                            ggAAGTTAGCGAGTCAGCATTTCAGGCGCATTCGCAAAGCCGCTGGCTTGCTAATAATCTCGCTGCGCGTTTTGCTTCAGGCAATGATACCGGCTGGGTACACTTCCCCGCGTTTAGCGTCGCAGCGACGCCGCCTGCCGATCACTACCACTTAACGGACCGGGCGACTACGTTCGCCATTCTGCTGATTGCACTTACTTTTCTGGCGTTCTTTGTTTTTTAAACGCTACCCGCGCCACGTTTatccccc                                                                                                   >  2:293662/1‑244 (MQ=255)                                                                                                                                              aaGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTAGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATg                                                                                                >  1:280468/1‑251 (MQ=255)                                                                                                                                                       aGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGTCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTg                                                                                       <  1:110076/251‑1 (MQ=255)                                                                                                                                                            ggTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCACCGAGTTTAGCCTCCCACTGACGCCGCCAGCCGATCACTACCATTTATCGGACCCGGCGACTCCGTAAGCCATTCTGCTGATTGCACTGCCTTTTATGGAGTTCTTTGTTTTTGCAACGCTCACCGCGCAACGTTTACACCACATGCAATCTTTGCtggt                                                                                  >  1:299634/1‑251 (MQ=255)                                                                                                                                                                    ggCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCttt                                                                          >  2:267907/1‑251 (MQ=255)                                                                                                                                                                     gcgcAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTc                                                                         <  1:126720/251‑1 (MQ=255)                                                                                                                                                                              cAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCACGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGAt                                                                >  2:44797/1‑251 (MQ=255)                                                                                                                                                                                aaGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGTCTGGGAAAACTTCCCCGCTTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCTGGCGACTAAGTACGCCATTCTGCTGATTGCACGGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGTTGGTGCTTTCATTGGGGATGt                                                              <  2:280468/251‑1 (MQ=255)                                                                                                                                                                                aaGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGt                                                              <  1:98825/251‑1 (MQ=255)                                                                                                                                                                                    agcTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCAcc                                                                                                                     >  2:500760/1‑192 (MQ=255)                                                                                                                                                                                    agcTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCAcc                                                                                                                     <  1:500760/192‑1 (MQ=255)                                                                                                                                                                                                   aatctagacgAGCGTTTTGCTTCAGGCACTGATACCGCCTGGGAAAACTTCCCCGCCTTTAGCATCGAAGTGACGACGCCAGCCGCTCAATCACAATTAACGGACCGGGCGACTAAAAACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTggg                                                                               >  2:368556/10‑215 (MQ=255)                                                                                                                                                                                                   aatCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGGCGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCGATGCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTggg                                                                               <  1:368556/215‑1 (MQ=255)                                                                                                                                                                                                   aatCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGgc                                           >  2:441824/1‑251 (MQ=255)                                                                                                                                                                                                   aatCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGgc                                           >  1:459231/1‑251 (MQ=255)                                                                                                                                                                                                        cGGTGAGCGTGTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGgcg                                          <  1:320696/247‑1 (MQ=255)                                                                                                                                                                                                         ggagAGCGTTTTGCTTCAGGCAATGATCCCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGAAGACGCCAGCCGCTCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCCTTGGTGATGTTTTATTTGCTCTTGCTGGCGCtttt                                     >  2:55628/4‑250 (MQ=255)                                                                                                                                                                                                                         cAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGttt                      <  2:236351/250‑1 (MQ=255)                                                                                                                                                                                                                             cAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCg                 <  1:55628/251‑1 (MQ=255)                                                                                                                                                                                                                              aaTGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACcg                                                                                                                    >  1:535222/1‑151 (MQ=255)                                                                                                                                                                                                                              aaTGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACcg                                                                                                                    <  2:535222/151‑1 (MQ=255)                                                                                                                                                                                                                                 gATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGCCGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTAACGTGGc             >  2:336067/1‑251 (MQ=255)                                                                                                                                                                                                                                  aTACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCCGCCGATCAATACCAATTAACGGACCGGGCGCCTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGc             >  1:598626/1‑250 (MQ=255)                                                                                                                                                                                                                                    aCCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATg          >  1:410600/1‑251 (MQ=255)                                                                                                                                                                                                                                    aCCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATg          >  1:543647/1‑251 (MQ=255)                                                                                                                                                                                                                                    aCCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCGGGTGGTGCTTTCATTGGTGATGTTTTATGGGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATg          <  1:95029/251‑1 (MQ=255)                                                                                                                                                                                                                                      cGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCCGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATgg         >  2:379892/1‑250 (MQ=255)                                                                                                                                                                                                                                            ggAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCGTTTTTTTTGAAACGCTCACCGCTCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGTGCTTGCTGGCGCTTTCGGAACATATCGGTTTTACCGTGGCATGGATAATCg  <  1:390201/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                               ACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTCTACCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATACCGGTTTTACCGTGGCATGGATAATCG  >  NC_000913/4638758‑4639243

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