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	90.2%	98.9 / 1.1	41	G25D (GGT→GAT)	rdgB	dITP/XTP pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base A (17/20);  total (19/22)
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.61e-01

CGGCCAGGCCAACAGTCATCTGGTGAAGTTTCTCGGTAAGCAATTCCGGGTTGCCAAAAGCCAGGTGGTGATTGAAAAAGGCGAACTTGGCCGCCACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGA  >  NC_000913/3096514‑3096999                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                      cGGCCAGGCCAACAGTCATCTGGTGAATTTTCTCGGTAAGCAATTCCGGGTTGCCAAAAGCCAGGTGGTGATTGAAAAAGGCGAACTTGGCCGCCACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATc                                                                                                                                                                                                                                             <  1:357922/251‑1 (MQ=255)        gccaACAGTCATCTGGTGAAGTTTCTCGGTAAGCAATTCCGGGTTGCCAAAAGCCAGGTGGTGATTGAAAAAGGCGAACTTGGCCGCCACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGccc                                                                                                                                                                                                                                      <  1:229049/251‑1 (MQ=255)                               cctcGGTAAGAAATTCCGGGTGGCAAAAAGCCAGTGGGTATTTAAAAAAGGCAAACTGGGCCGCCACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTcc                                                                                                                                                                                                               <  2:130099/250‑1 (MQ=255)                                       aGCAATCCCGGGTTGCAAAAAGCCAGGTGGTGATTGAAAAAGGCGAACTTGGCCGCCACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTgaaga                                                                                                                                                                                                       <  1:433628/251‑1 (MQ=255)                                           aTTCCGGGTTGCAAAAAGCCAGGTGGTGATTGAAAAAGGCGAACTTGGCCGCCACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAAcc                                                                                                                                                                                                   <  1:114491/251‑1 (MQ=255)                                            ttCCGGGTTGCCAAAAGCCAGGTGGTGATTGAAAAAGGCGAACTTGGCCGCCACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAACCAGACCTCGGCGTTGATTCCGCTGAAGAAACCg                                                                                                                                                                                                  >  1:352174/1‑251 (MQ=255)                                                                      gATTGAAAAAGGCGAACTTGGCCGCCACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGCAGAAACCGGCCTGACCTTTATCGAAAAAGCGAtt                                                                                                                                                                        >  2:244114/1‑251 (MQ=255)                                                                       aTTGAAAAAGGCGAACTTGGCCGCCACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTc                                                                                                                                                                       >  1:497150/1‑251 (MQ=255)                                                                        ttGAAAAAGGCGAATTGGGCCGCCACAAACAAATTAAAATCATTAATCCGCAACAAAACCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTGGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTc                                                                                                                                                                       <  2:82663/250‑1 (MQ=255)                                                                            aaaagactatcttccccgccACAAAGAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAACCGAGGCGTTAATTAATGAGGTATCCTATGCAAAAATTGGTCCTCGCAACCGGCAATGTCGGTAAAGTGCTTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATACCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAATAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGaaa                                                                                                                                                                  <  2:163064/236‑1 (MQ=255)                                                                                        tGGCCGCCACAAACAAATTAAAAACATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCTTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGGCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATg                                                                                                                                                        <  1:890/249‑1 (MQ=255)                                                                                        gggCCGCCACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCg                                                                                                                                                      <  1:13160/250‑1 (MQ=255)                                                                                                    acaaaTTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGAcc                                                                                                                                          >  1:448219/1‑251 (MQ=255)                                                                                                       aaTTAAATTCATTAATCCGCAATAAACCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCTTGGGCCAAACAGACCGCTGCGTTGATTCCGCTGAAGAAAGCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCt                                                                                                                                       <  1:254554/251‑1 (MQ=255)                                                                                                          aaaaaTCATTATCCCGAAACAAATCCCTCCATAAATCTCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAAGGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTa                                                                                                                                    <  2:266843/250‑1 (MQ=255)                                                                                                                              aaaTCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCAGCAATGTCGGTAAAGTGCGTGAGCTGGAGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAACCAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGCCCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGActc                                                                                                                >  2:311535/1‑251 (MQ=255)                                                                                                                                aTCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACtctg                                                                                                              >  1:135727/1‑251 (MQ=255)                                                                                                                                                   cgTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCt                                                                                           <  1:334331/251‑1 (MQ=255)                                                                                                                                                    gTTAATTAATTAGGTATCCTATGCAAAAAGTTGTACTCGCAAACGGCAATGTCGGAAAAGTGCGTGAGCTGGCGTGGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGACAACGCGATTCTGAAAGCGCGCCATACGGCAAACGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGtgcgtgtgttt                                                                                          >  2:72901/1‑241 (MQ=255)                                                                                                                                                          taattaGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGAGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTggcggc                                                                                    >  1:327716/1‑251 (MQ=255)                                                                                                                                                             ttaGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGAATTCGATCTTGATATCGTGGACCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGgcgcg                                                                                 >  1:29857/1‑251 (MQ=255)                                                                                                                                                              taGGAATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGACAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGGAGATGTGCTTGGCGgcgccc                                                                                >  2:135228/1‑249 (MQ=255)                                                                                                                                                                     ccTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGGTAGCTGGCGGCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCa                                                                          <  2:436327/250‑1 (MQ=255)                                                                                                                                                                         tGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAACGTGCGTGAGCTGGCGTAGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCACGGCGTTGATTCCGCTGAAGACACCGGGCTGCCCTTTATCGAAAACGCGATTCTGCAACAGCGCCATGAGGCACAAGTGACCGCTTTCACGGCAATGGCCGACGACTCTGGTCTGGCGGGAGAtgagacacggggcgcgccgtgcgattc                                                                      >  2:423896/1‑225 (MQ=255)                                                                                                                                                                           cAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCg                                                                                                                                         <  1:131236/181‑1 (MQ=255)                                                                                                                                                                           cAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCg                                                                                                                                         >  2:131236/1‑181 (MQ=255)                                                                                                                                                                            gaaacgTTGTCCTCGCAACCGGCAATGTCGGTCAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAACCGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGAACGCTTTACCGGCAATTGCCTCCGACTCTTGTCTGGCGGGGGATGTTCTTGGGTGCGCCCCGGGCACTTACTcc                                                                  >  2:369841/6‑251 (MQ=255)                                                                                                                                                                                aGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGAAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCcgcgc                                                              >  1:927/1‑251 (MQ=255)                                                                                                                                                                                    gTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGAc                                                                                                                                           >  1:15998/1‑170 (MQ=255)                                                                                                                                                                                    gTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTGCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGAc                                                                                                                                           <  2:15998/170‑1 (MQ=255)                                                                                                                                                                                         cgaacccGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGTCGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTcc                                                      <  2:387827/245‑1 (MQ=255)                                                                                                                                                                                          gCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACTGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCgg                                                    >  1:203715/1‑251 (MQ=255)                                                                                                                                                                                           cAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTTCTCCTCGCGTTATTCCggg                                                   >  2:184573/1‑250 (MQ=255)                                                                                                                                                                                                  cAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGAcg                                            >  1:471331/1‑251 (MQ=255)                                                                                                                                                                                                               gTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCTATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAAGTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCaaaa                               <  1:311535/251‑1 (MQ=255)                                                                                                                                                                                                                      aGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTTCTCCGCGCGTTATTCCGGTGAAGGCGCGACAGGTCAAAAAAATCTg                        >  2:492918/1‑251 (MQ=255)                                                                                                                                                                                                                      aGCTGGCGTCGCGGCTTAGCGACTGCGATCTTGATATCGTGGACCAAACAGACCTCGGCGTTGAGTCCGCTGAAGAAACCGGCCTGTCCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCTGGCATGTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTg                        <  2:478904/251‑1 (MQ=255)                                                                                                                                                                                                                       ggtagcgtctgtgctTAGCGACTGCGGTCTTGATATCGTGCCACAATCAGACCCCAGCGTTGAGTCGGCTGGAGAAACAGGCATGACCGTTATCTAAAACGTGATTCTGAAAGCGCGCCATGCAGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTTGTCTGGGGGTAGATGTGCTTGGCGGCTCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGc                       <  2:429282/240‑1 (MQ=255)                                                                                                                                                                                                                               cGCTGCTTAGCGACTTCGATCTTGATAACGTGGCACAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAAct                <  2:144648/250‑1 (MQ=255)                                                                                                                                                                                                                                        gCGACTCCGATCTTGATATCGTGGCCCAATCAGACCTCGGCGTTGATTCCGCTGCAGAAACCGGCCTGACGTTTTTCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCTGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACa      <  2:189254/251‑1 (MQ=255)                                                                                                                                                                                                                                              tCGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGa  <  1:244114/249‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                      CGGCCAGGCCAACAGTCATCTGGTGAAGTTTCTCGGTAAGCAATTCCGGGTTGCCAAAAGCCAGGTGGTGATTGAAAAAGGCGAACTTGGCCGCCACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGA  >  NC_000913/3096514‑3096999

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