Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	251,629	C→T	G244G (GGC→GGT)	dinB →	DNA polymerase IV

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	251,629	0	C	T	95.8%	144.7 / NA	48	G244G (GGC→GGT)	dinB	DNA polymerase IV
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (18/28);  minor base G (0/2);  total (18/30)
Fisher's exact test for biased strand distribution p-value = 5.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGAAGTTCCGGCATTTTTACAAACCTTACCGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAAGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGATTGCTCGCCAGGGGGTGAAATTAAAGTTCGACGATTTTCAGCAAACCACCCAGGAGCACGTCTGGCCGCGGCTGAATAAAGCTGATCTAATCGCCACCGCGCGTAAAACCTGGGATGAACGCCGCGGCG  >  NC_000913/251398‑251873                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     aGAAGGTCCGGCATTTTTACAAACCTTCCCGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAAGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCgg                                                                                                                                                                                                                                   <  1:6657/251‑1 (MQ=255)              tttttACAAACCTTACCGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAAGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGCTGGCGGAAGATATtcatc                                                                                                                                                                                                                       >  2:314440/1‑250 (MQ=255)               ttttACAAACCTTACCGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAAGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCAc                                                                                                                                                                                                                     >  1:292309/1‑251 (MQ=255)                     aaaCCTGACCGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAAGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCt                                                                                                                                                                                                               <  2:110566/251‑1 (MQ=255)                                  ggCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAAGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGa                                                                                                                                                                                                  <  1:185974/251‑1 (MQ=255)                                    cAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAAGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGAtt                                                                                                                                                                                                <  1:64261/251‑1 (MQ=255)                                     aaaaaTCCCCGTCGGCGGCAAAGTCTCAGCGGCAAAACGGGAAGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATGTTGTGGGAGCGTAGTCAGGGGATTGACTAACGCGATGTTAACAGCGAACGTTTGCGAAAATCCGTCGGTGTGGAACGCACTATGGCGTAAGATATTCATCACTGGTCTGAATGTGAAGCGATTa                                                                                                                                                                                               <  2:208040/251‑1 (MQ=255)                                     aaaaaTCCCCGGCGGCGGCAAAGGCTCAGCGGGAAAACTGGAAGCGATGGTGCTGCGGACCTGCGGTTATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTTAGGGGCTTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTa                                                                                                                                                                                               <  1:219804/251‑1 (MQ=255)                                       aaaTCCCCGGGGTCGGCAAAGGCTCAGCGGCAAAACTGGAAGCGAGGGGGCTGCGGACCTGCGTTGATGTACAAAAGTGTGATCTGGTGAGGCTGCTTAAACGCGTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAAGCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGACTATc                                                                                                                                                                                             <  2:156697/251‑1 (MQ=255)                                       aaaTCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAAGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATc                                                                                                                                                                                             >  1:57397/1‑251 (MQ=255)                                       aaaTCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAAGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATc                                                                                                                                                                                             >  1:281221/1‑251 (MQ=255)                                        aaTCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAAGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCg                                                                                                                                                                                            <  1:20868/251‑1 (MQ=255)                                            cccGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAAGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGTCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTTTAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCg                                                                                                                                                                                        <  1:239212/251‑1 (MQ=255)                                            cccGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAAGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCg                                                                                                                                                                                        >  1:350247/1‑251 (MQ=255)                                              cGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAAGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGc                                                                                                                                                                                      >  1:188401/1‑251 (MQ=255)                                                gCGTCGGCAAAGTCTCAGCGGCAAAACTGGAAGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTg                                                                                                                                                                                    <  1:355285/251‑1 (MQ=255)                                                        aaaGTCTCAGCGGCAAAACTGGAAGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCgg                                                                                                                                                                             >  1:248630/1‑250 (MQ=255)                                                                aGCGGCAAAACTGGAAGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAAc                                                                                                                                                                    <  2:350247/251‑1 (MQ=255)                                                                        aaCTGGAAGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTg                                                                                                                                                            <  1:374192/251‑1 (MQ=255)                                                                         aCTGGAAGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGGGTGATCTGGTGATGCTGCTGAAACGCTTTGGCAAATTTGGCCGCATCTTGTGGGTGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGTAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTgg                                                                                                                                                           <  1:189491/251‑1 (MQ=255)                                                                                             ggtccTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTa                                                                                                                                       <  2:281221/248‑1 (MQ=255)                                                                                                    cGGTGATGTACAAAAGTTTGATCTGGTGATGCGGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGAGTGACGAACGCGATGTTAACAGCGAACGTTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGGGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGATTg                                                                                                                                <  1:224186/251‑1 (MQ=255)                                                                                                        gATGTACAAAAGGGTGATCTGGTGATGCTGCTTAAACGCTTTTGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCTGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGATTGCTc                                                                                                                             <  1:213096/250‑1 (MQ=255)                                                                                                          tGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAg                                                                                                                                                                                                     <  1:93159/176‑1 (MQ=255)                                                                                                          tGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAg                                                                                                                                                                                                     >  2:93159/1‑176 (MQ=255)                                                                                                             aCAAAAGGGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATGTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAACCCTGATTTACTGATTGCTCGCCAgg                                                                                                                       <  2:57397/251‑1 (MQ=255)                                                                                                              cAAAAGGGTGAACTGGTGCTCATGCTTAAACGGTTTGGCAAATTTGGGCGAATTTTGTGGGAGCGTTGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGGGTGGACCGCACTATGGCGGAAGAAATTCATCACGGTTCTGAATGTGAAGCGATTCTCGAGCGGCTGTATCCTGAACTTGAACGCCGTCTGGCAAAGGTAAAGCCTGATTTACTGATTGCCCGCCAggg                                                                                                                      <  2:226150/251‑1 (MQ=255)                                                                                                                       gATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGATTGCTCGCCAGGGGGTGAAAt                                                                                                              <  1:90705/250‑1 (MQ=255)                                                                                                                              tGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGATTGCTCGCCAGGGGGTGAAATTAAAGTTc                                                                                                      >  2:332316/1‑251 (MQ=255)                                                                                                                               gATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGATTGCTCGCCAGGGGGTGAAATTAAAGTTcg                                                                                                     <  1:238645/251‑1 (MQ=255)                                                                                                                                aTGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGATTGCTCGCCAGGGGGTGAAATTAAAGTTcga                                                                                                    <  1:248508/251‑1 (MQ=255)                                                                                                                                    tgctTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGATTGCTCGCCAGGGGGTGAAATTAAAGTTCGACGAt                                                                                                >  1:351362/1‑251 (MQ=255)                                                                                                                                                      aaTTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGGTTGCTCGCCAGGGGGTGAAATTAAAATTCGACGATTTTCAGCACACCAccc                                                                                >  2:255394/1‑249 (MQ=255)                                                                                                                                                         ttGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGATTGCTCGCCAGGGGGTGAAATTAAAGTTCGACGATTTTCAGCAAACCACCCAgg                                                                             >  1:289600/1‑249 (MQ=255)                                                                                                                                                                cATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGATTGCTCGCCAGGGGGTGAAATTAAAGTTCGACGATTTTCAGCAAACCACCCAGGAGCACGTCt                                                                    >  1:363816/1‑251 (MQ=255)                                                                                                                                                                          gAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCTATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGATTGCTCGCCAGGGGGTGAAATTAAAGTGCGACGATGTTCAGCAAACCACCCAGTAGCACGTCTGGCCGCGGCt                                                          <  2:218750/251‑1 (MQ=255)                                                                                                                                                                              gTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGATTGCTCGCCAGGGGGTGAAATTAAAGTTCGACGATTTTCAGCAAACCACCCAGGAGCACGTCTGGCCGCGGCTGAAt                                                      >  1:155738/1‑251 (MQ=255)                                                                                                                                                                               tAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTATGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAAGGAAAACCTGTTTTACTGATTGCTCGCCAGGAGGTGAAATTAAAGTTCGACGATTTTCAGCAAACCACCCACGCGCACGTCGGGCCGCGGCTgcat                                                      >  2:27544/1‑247 (MQ=255)                                                                                                                                                                                   catgggaattgcGAACGCGATGTGAACAGCGAACGGTTGCGAAAATCCGTCGGTGGTGAACGCACGATGGCGGAATATATTCAGCCCTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGATTGCTCTCCAGGGGGTGAAATTAAAGTTCGACGATTTTCAGCAAACCACCCAGGAGCACGTCTGGCCGCGGCTGAATAAAGc                                                 <  2:229893/240‑1 (MQ=255)                                                                                                                                                                                   cAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGATTGCTCGCCAGGGGGTGAAATTAAAGTTCGACGATTTGCAGCAAACGACCCAGGAGCACGTCTGGCCGCGGCTGAATAAAg                                                  <  1:98623/250‑1 (MQ=255)                                                                                                                                                                                     ggggATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGATTGCTCGCCAGGGGGTGAAATTAAAGTTCGACGATTTTCAGCAAACCACCCAGGAGCACGTCTGGCCGCGGCTGAATAAAGCTg                                               >  1:310393/1‑251 (MQ=255)                                                                                                                                                                                            gACGAACGCGATGATAACAGCTAACGGTTGCGAAAAACCGTCGGTGTGGAACGCACGATGGCGGAAGATAGTCATCACTGGTCTGAATGTGAAGCGATTACCGAGCGTCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGAGTTACAGATTGCTCGCCAGGGGGTGAAATTAAAGTTCGACGAGTTTCAGCAAACCACCCAGGAGCACGTCTGGCCGCGGCTGAATAAAGCTGATCTAAt                                        <  1:353514/251‑1 (MQ=255)                                                                                                                                                                                                           aaCAGCGAACGGTTGTGAAAATCCTGCGGGGTGGAACGCACGATGGCGGAAGATATTCAACACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAAATTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGATGGCTCGCCAGGGGGTGAAATTAAAGTTAGACGATTTTCAGCAAACCACCCAGGAGCACGTCTGGCCGCGGCTGAATAAAGCTGATCTAATCGCCACCGCGCGTaa                         <  2:203418/251‑1 (MQ=255)                                                                                                                                                                                                            aCAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGATTGCTCGCCAGGGGGTGAAATTAAAGTTCGACGATTTTCAGCAAACCACCCAGGAGCACGTCTGGCCGCGGCTGAATAAAGCTGATCTAATCGACACCGCGCGTaa                         >  1:280868/1‑250 (MQ=255)                                                                                                                                                                                                                  aaCGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGATTGCTCGCCAGGGGGTGAAATTAAAGTTCGACGATTTTCAGCAAACCACCCAGGAGCACGTCTGGCCGCGGCTGAATAAAGCTGATCTAATCGCCACCGCGCGTAAAACCTgg                  <  1:262620/251‑1 (MQ=255)                                                                                                                                                                                                                            aaaaTCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGATTGCTCGCCAGGGGGTGAAATTAAAGTTCGACGATTTTCAGCAAACCACCCAGGAGCACGTCTGGCCGCGGCTGAATAAAGCTGATCTAATCGCCACCGCGCGTAAAACCTGGGATGAAcgcc        <  2:155738/251‑1 (MQ=255)                                                                                                                                                                                                                             aaaTCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGATTGCTCGCCAGGGGGTGAAATTAAAGTTCGACGATTTTCAGCAAACCACCCAGGAGCACGTCTGGCCGCGGCTGAATAAAGCTGATCTAATCGCCACCGCGCGTAAAACCTGGGATGAAcgccg       <  1:159227/251‑1 (MQ=255)                                                                                                                                                                                                                                  cGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGATTGCTCGCCAGGGGGTGAAATTAAAGTTCGACGATTTTCAGCAAACCACCCAGGAGCACGTCTGGCCGCGGCTGAATAAAGCTGATCTAATCGCCACCGCGCGTAAAACCTGGGATGAACGCCgcggcg  >  1:179772/1‑251 (MQ=255)                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     AGAAGTTCCGGCATTTTTACAAACCTTACCGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAAGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGATTGCTCGCCAGGGGGTGAAATTAAAGTTCGACGATTTTCAGCAAACCACCCAGGAGCACGTCTGGCCGCGGCTGAATAAAGCTGATCTAATCGCCACCGCGCGTAAAACCTGGGATGAACGCCGCGGCG  >  NC_000913/251398‑251873

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