Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,096,754	G→A	G25D (GGT→GAT)	rdgB →	dITP/XTP pyrophosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,096,754	0	G	A	92.3%	152.0 / NA	53	G25D (GGT→GAT)	rdgB	dITP/XTP pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  major base A (25/23);  minor base C (4/0);  total (30/23)
Fisher's exact test for biased strand distribution p-value = 1.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

CGGCCAGGCCAACAGTCATCTGGTGAAGTTTCTCGGTAAGCAATTCCGGGTTGCCAAAAGCCAGGTGGTGATTGAAAAAGGCGAACTTGGCCGCCACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTG  >  NC_000913/3096514‑3096989                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                            aggccaggccaACAGTCATCTGGTGAAGTTTCTCGGTAAGCACTTCCGGGTTGCCAAAAGCCAGGTGGTGATTGAAAAAGGCGAACTTGGCCGCCACAAACAAATTACAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTCTGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTTGCGTCGCTGCTTAGCACCTTCGATCTTGatctc                                                                                                                                                                                                                                   >  2:105664/2‑248 (MQ=255)       ggccaACAGTAATCTGGTGAAGTTTCTCGGTAAGCAATTCCGGGTTGCCAACAGCCAGGAGGTGATTGAAAAAGGCGAACTGGGCCGCCACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCGTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGAGCGTGAGCTGGCGTCGTTGCTTTGCGGCTTCGATCTTGTTATAGTGGcc                                                                                                                                                                                                                             >  2:483584/1‑251 (MQ=255)                      ggTGAAGTTTCACGGTAAGCAATTCCGGGTTGCCAAAAGCCAGGTGGTGATTGAAAAAGGCGAACTTGGCCGCCACAAACAAATTAAAATCATTAATCAGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAACAAGTTGTCCTCGCAACCGGTAATGTCGGTAAAGTGCGTGAGCTGGCGTCGATGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCgcc                                                                                                                                                                                                              >  2:265006/1‑249 (MQ=255)                                          aaTTCCGGGTTGCCAAAAGCCAGGTGGTGATTGAAAAAGGCGAACTTGGCCGCCACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAAc                                                                                                                                                                                          >  1:348013/1‑251 (MQ=255)                                          aaTTCCGGGTTGCCAAAAGCCAGGTGGTGATTGAAAAAGGCGAACTGGGCCGCCACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAAc                                                                                                                                                                                          <  1:282528/251‑1 (MQ=255)                                                          aaGCCAGGTGGTGATTGAAAAAGGCGAACTTGGCCGCCACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGATGATTCCGCTGAAGAAACCGGCCTGACCTTTATc                                                                                                                                                                          >  2:193589/1‑251 (MQ=255)                                                           aGCCAGGTGGTGATTGAAAAAGGCGAACTGGGCCGCCACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCg                                                                                                                                                                         <  1:116528/251‑1 (MQ=255)                                                                     cgagtaaaaaaGGCGAACTCGGCCGCCACAAAAAAATTAAAATCATTAAACCGCAACAAAGCCCGCCAGAAATCGCGGCGTTAAGTAATTATGTATCCTATGCAAAAAGTTGTCCTCGCAAACGGCAATGTCGGTAAAGTGCTCGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGTCCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAAcgcgt                                                                                                                                                                <  2:41653/244‑2 (MQ=255)                                                                           aaaaaGGCGAACTTGGCCGCCACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGa                                                                                                                                                          >  2:234123/1‑250 (MQ=255)                                                                            aaaaGGCGAACTTGGCCGCCACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACCGACCTCAGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCAAAACCGCGATTCTGaa                                                                                                                                                         >  2:70305/1‑250 (MQ=255)                                                                              aaGGCGAACGTCGCCGCCAAAAACAAATTAAAATCATTATGCCGCACCAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAgc                                                                                                                                                      <  2:1891/251‑1 (MQ=255)                                                                               aGGCGAACTTGGCCGCCACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTCAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAgcg                                                                                                                                                     >  2:282291/1‑251 (MQ=255)                                                                               aGGCGAACTTGGCCGCCAAAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATAGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAgcg                                                                                                                                                     <  2:159633/251‑1 (MQ=255)                                                                                  cGAACTTGACCGCCACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGCTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTAACGAATACGCGATTCTGAAAacgcgc                                                                                                                                                  >  2:469260/1‑251 (MQ=255)                                                                                          gccgccACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGc                                                                                                                                          <  1:234123/251‑1 (MQ=255)                                                                                           ccgccACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCACAACGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGCTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGACGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGACAGCGCGCCATGCGGc                                                                                                                                          >  2:61016/1‑250 (MQ=255)                                                                                               cacaAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCaaa                                                                                                                                       >  2:412794/1‑249 (MQ=255)                                                                                                caaaacaaaTTAAAATCATTAATCCGCCACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTAGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGACGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCATAAGt                                                                                                                                    <  2:9816/249‑1 (MQ=255)                                                                                                acaaacaaATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTgaagaa                                                                                                                                                                                            <  2:252044/195‑1 (MQ=255)                                                                                                acaaacaaATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTgaagaa                                                                                                                                                                                            >  1:252044/1‑195 (MQ=255)                                                                                                   aacaaaTTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGAc                                                                                                                                 >  1:32058/1‑251 (MQ=255)                                                                                                     caaaTTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGATGTCATCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGAGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGAGTTGATTCCGCTGAAGAAACCGGCCTGATCTTTATAGAAAACGCGATTCTGACAACGCGCCATGCGGCAAAAGTGACCg                                                                                                                               >  1:230107/1‑251 (MQ=255)                                                                                                      aaaTTAAAATCATTAATCCGCAACAAATCCCGCCAGAAAGCGCGGCGTTAATTAATTATGTAACCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGGCGGTAAAGTGCGTTAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCTCCATGCGGCAAAAGTGTCCGc                                                                                                                              <  2:273024/251‑1 (MQ=255)                                                                                                           aaaaTCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTa                                                                                                                          <  1:427359/250‑1 (MQ=255)                                                                                                            aaaTCATTAATCCGCAACAAATCCCGCCAAAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTAcc                                                                                                                        >  1:327555/1‑251 (MQ=255)                                                                                                                 aTTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAACACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCa                                                                                                                    >  2:242070/1‑250 (MQ=255)                                                                                                                     aTCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGc                                                                                                               >  1:394535/1‑251 (MQ=255)                                                                                                                       ccGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTCCCGGCAATTGCcg                                                                                                             >  2:272555/1‑251 (MQ=255)                                                                                                                           aacaaATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCcgacg                                                                                                          >  1:219186/1‑250 (MQ=255)                                                                                                                             caaATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGCTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAAtt                                                                                                                 >  2:366008/1‑241 (MQ=255)                                                                                                                             caaATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCAGGCGGCAAAAGTGACCGCTTTACCGGCAAtt                                                                                                                 <  1:366008/241‑1 (MQ=255)                                                                                                                                aTCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACtctg                                                                                                    >  1:75390/1‑251 (MQ=255)                                                                                                                                  cccGCCAGAAAACGCGGCGTTAATTAATTAGGTATCCTAGGCAAAAAGTTGTCCTCGAACCCGGGAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCGGCTGAGCGACTTCGATCTTGATAGCGCGGCGCAAACAGAACTCGGCGTTGATGCCGCTGAAGAAACCGGCCTGACCTTGATCGAAAACGCGATTCTGAAAGCGCGCTATGCGTCATAAGTGACCGCGTTTCCGGCATTTGCCTAAGACtctggt                                                                                                  <  2:56136/251‑1 (MQ=255)                                                                                                                                    cGCCATAAATCGCGTCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGTCAATGGCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACtctggtct                                                                                                <  2:455255/251‑1 (MQ=255)                                                                                                                                        aGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGCTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCAGCCATTGCCGACGACtctggtctgg                                                                                              >  1:408985/1‑249 (MQ=255)                                                                                                                                        aGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCg                                                                                            <  1:1693/251‑1 (MQ=255)                                                                                                                                            aTCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGAGATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCGTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGTCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAg                                                                                        <  1:56303/251‑1 (MQ=255)                                                                                                                                                 ggcgTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCGTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGAtgtg                                                                                   <  1:30839/251‑1 (MQ=255)                                                                                                                                                     ttaattaaTTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTg                                                                               >  2:232624/1‑251 (MQ=255)                                                                                                                                                     ttaattaaTTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTg                                                                               >  1:238537/1‑251 (MQ=255)                                                                                                                                                               aGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGcc                                                                     <  2:408985/251‑1 (MQ=255)                                                                                                                                                                 gTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTGGGTGGGCTTGGCGGCGCGCCgg                                                                   >  1:81098/1‑251 (MQ=255)                                                                                                                                                                   aTCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGc                                                                 <  2:394535/251‑1 (MQ=255)                                                                                                                                                                    tCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGcc                                                                >  1:464463/1‑250 (MQ=255)                                                                                                                                                                        aTGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTAGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCAtt                                                              <  1:353760/249‑1 (MQ=255)                                                                                                                                                                           cAAAAAGTTGGCCTCGCAACCGGCAATGTCGGTAAAGGGCGTGAGCTGTCGTCGCTGCTTAGCGACTTCGATCTTGATAGCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTc                                                         <  2:214343/251‑1 (MQ=255)                                                                                                                                                                              aaaGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCcgc                                                      <  2:219186/251‑1 (MQ=255)                                                                                                                                                                                             aCCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACCGACCTCGGCGTTGATTCCGCCGAAAACACCGGCCTGACCTCTATCGCAAAAGCGATTCTAAAAGCGCGCCATGCGGCAAACGTCACCGCTTTACCGGCAATTGCCGCCGACTCTGGTCTGGCGGCAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTgg                                       >  1:440493/1‑250 (MQ=255)                                                                                                                                                                                                        cggtaaatgacGTGGGCTGGGGTCGCTGCTTAGCTACTTCGATATTGATAGCGTGGACCAAACAGACCTCGGCGTTGAGGCCGCGAAAGAAACCGGCCTGAACTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGTCTCTGGTCTGGCGGTATATGTGCTTTGCGGCGCGCGGTGCAGTTACTCCGCGCGTTATGCCTGTGAAgtctcgagt                             <  2:292436/240‑9 (MQ=255)                                                                                                                                                                                                                   gTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAAt                 >  1:242076/1‑251 (MQ=255)                                                                                                                                                                                                                          ggCGTCGCTGCTTAGCGTCTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCa            <  1:70305/249‑1 (MQ=255)                                                                                                                                                                                                                               cGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAActg     >  1:171207/1‑251 (MQ=255)                                                                                                                                                                                                                                  tgctTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAACGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCACAAActgctg  >  2:45160/1‑251 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                            CGGCCAGGCCAACAGTCATCTGGTGAAGTTTCTCGGTAAGCAATTCCGGGTTGCCAAAAGCCAGGTGGTGATTGAAAAAGGCGAACTTGGCCGCCACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTG  >  NC_000913/3096514‑3096989

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