Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	924,040	Δ1 bp	coding (777/2277 nt)	clpA →	ATP‑dependent Clp protease ATP‑binding subunit ClpA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	924,040	0	T	.	97.6%	183.8 / ‑3.7	41	coding (777/2277 nt)	clpA	ATP‑dependent Clp protease ATP‑binding subunit ClpA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base . (19/21);  minor base C (1/0);  total (20/21)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GCTTGCGCGCGTGGGCGGAATCGACCCACTGATTGGTCGTGAGAAGGAGCTGGAGCGTGCTATTCAGGTTCTCTGCCGTCGCCGTAAAAACAACCCGCTGCTGGTGGGGGAATCTGGTGTCGGTAAAACCGCGATTGCGGAAGGTCTTGCCTGGCGAATTGTTCAGGGCGATGTGCCGGAAGTGATGGCTGACTGTACGATTTACTCTCTCGATATCGGTTCTCTGTTAGCGGGCACAAAATATCGCGGCGACTTTGAAAAACGTTTTAAAGCGTTGCTCAAGCAGCTGGAGCAGGACACTAACAGCATCCTGTTTATTGATGAGATCCACACCATTATCGGTGCGGGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGAGTTCAGCAACATTTTCGAGAAAGACCGTGCTCTGGCGCGTC  >  NC_000913/923797‑924281                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  gCTTGCGCGCGTGGGCGGAATCGACCCACTGATTGGTCGTGAGAAGGAGCTGGAGCGTGCTATTCAGGTTCTCTGCCGTCGCCGTAAAAACAACCCGCTGCTGGTGGGGGAATCTGGTGTCGGTAAAACCGCGATTGCGGAAGGTCTTGCCTGGCGAATTGTTCAGGGCGATGTGCCGGAAGTGATGGCTGTCTGTACGATTTACTCTCTCGATATCGGTTCTCTGTTAGCGGGCACAAAATA‑Cgcggcg                                                                                                                                                                                                                                            >  1:47416/1‑250 (MQ=255)     gcgcgcGTGGGCGGAATCGACCCACTGATTGGTCGTGAGAAGGAGCTGGAGCGTGCTATTCAGGTTCTCTGCCGTCGCCGTAAAAACAACCCGCTGCTGGTGGGGGAATCTGGTGTCGGTAAAACCGCGATTGCGGAAGGTCTTGCCTGGCGAATTGTTCAGGGCGATGTGCCGGAAGTGATGGCTGACTGTACGATTTACTCTCTCGATATCGGTTCTCTGTTAGCGGGCACAAAATA‑CGCGGCGACttt                                                                                                                                                                                                                                       <  1:72421/251‑1 (MQ=255)        ggcgggggcggactagacacgcTGATTGGTCGTGAGAAGGAGCTGGAGCGGGCTAGTCAGGTTCTCTGGCGTCGCGGTAAAAACAGCCCGCGGCTGGTGGGGGAATCTGGTGTCGGTAAAACCGCGATTGCGGAAGGTCTTGCCTGGCGAATTGTTCAGGGCGATGTGCCGGAAGTGATGGCGGACTGTAAGATTTACTCTCTCGATATCGGTTCTCTGTTAGCGGGCACAAAATA‑CGCGGCGACTTTGaa                                                                                                                                                                                                                                    <  2:115720/230‑1 (MQ=255)                    aTCGACCCACTGATTGGTCGTGAGAAGGAGCTGGAGCGTGCTATTCAGGTTCTCTGCCGTCGCCGTAAAAACAACCCGCTGCTGGTGGTGGAATCTGGTGTCGGTAAAACCGCGATTGCGGAAGGTCTTGCCTGGCGAATTGTTCAGGGCGATGTGCCGGAAGTGATGGCTGACTGTACGATTTACTCTCTCGATATCGGTTCTCTGTTAGCGGGCACAAAATA‑CGCGGCGACTTTGAAAAACGTTTTaaa                                                                                                                                                                                                                        <  2:227664/251‑1 (MQ=255)                      cGACCCACTGATTGGTCGTGAGAAGGAGCTGGAGCGTGCTATTCAGGTTCTCTGCCGTCGCCGTAAAAACAACCCGCTGCTGGTGGGGGAATCTGGTGTCGGTAAAACCGCGATTGCGGAAGGTCTTGCCTGGCGAATTGTTCAGGGCGATGTGCCGGAAGTGATGGCTGACTGTACGATTTACTCTCTCGATATCGGTTCTCTGTTAGCGGGCACAAAATA‑CGCGGCGACTTTGAAAAACGTTTTAAAGc                                                                                                                                                                                                                      <  2:216447/251‑1 (MQ=255)                      cGACCCACTGATTGGTCGTGAGAAGGAGCTGGAGCGTGCTATTCAGGTTCTCTGCCGTCGCCGTAAAAACAACCCGCTGCTGGTGGGGGAATCTGGTGTCGGTAAAACCGCGATTGCGGAAGGTCTTGCCTGGCGAATTGTTCAGGGCGATGTGCCGGAAGTGATGGCTGACTGTACGATTTACTCTCTCGATATCGGTTCTCTGTTAGCGGGCACAAAATA‑CGCGGCGACTTTGAAAAACGTTTTAAAGc                                                                                                                                                                                                                      >  1:392877/1‑251 (MQ=255)                        aCCCACTGCTTGGACGTGAGAAGGAGCTGGAGCGTGATATTCAGGTTCTCTGCCGTCGCCGTAAAAACAACCCGCTGCTGGTGGGGGAATCTGGTGTCGGTAAAACCGCGATTGCGGAAGGTCTTGCCTGGCGAATTGTTCAGGGCGATGTGCCGGAAGTGATGGCTGACTGTACGATTTACTCTCTCGATATCGGTTCTCTGTTAGCGGGCACCAAATA‑CGCGGCGACTTTGACAAACGTTTTAAAGCGt                                                                                                                                                                                                                    >  2:332888/1‑251 (MQ=255)                         cccccTGATTGGTCGTGAGAAGGAGCTGGAGCGTGCTATTCAGGTTCTCTGCCGTCGCCGTAAAAACAACCCGCTGCTGGTGGGGGAATCTGGTGTCGGTAAAACCGCGATTGCGGAAGGTCTTGCCTGGCGAATTGTTCAGGGCGATGTGCCGGAAGTGATGGCTGACTGTACGATTTACTCTCTCGATATCGGTTCTCTGTTAGCGGGCACAAAATA‑CGCGGCGACTTTGAAAAACGTTTTAAAGCGtt                                                                                                                                                                                                                   <  1:413081/247‑1 (MQ=255)                             cTGATTGGTCGTGAGAGGGAGCTGGAGCGTGCTATTCAGGGTCTCTGCCGTCGCCGTAAAAACAACCCGCTGCTGGTGGGGGAATATGGTGTCGGTAAAACCGCGATTGCGGAAGGTCTTGCCTGGCGAATTGTTCAGGGCGATGGGCCGGAAGTGATGGCTCACTGTACGATTTACTCTCTCGCTATCGGTTCTTTGTTAGCTGTCACAAAATA‑CGCGGCGACTTTGACAAACGTTTTAAGGCGTTGCTc                                                                                                                                                                                                               >  2:296799/1‑251 (MQ=255)                               gATTGGGCGTGAGAAGGAGCTGGAGCGTGCTATTCAGGTTCTCTGCCGTCGCCGTAAAAACAACCCGCTGCTGGTGGGGGAATCTGGTGTCGGTAAAACCGCGATTGCGGAAGGTCTTGCCTGGCGAATTGTTCAGGGCGATGTGCCGGAAGTGATGGCTGACTGTACGATTTACTCTCTCGATATCGGTTCTCTGTTAGCGGGCACAAAATA‑CGCGGCGACTTTGAAAAACGTTTTAAAGCGTTGCTCaa                                                                                                                                                                                                             <  1:328724/251‑1 (MQ=255)                                aTTGGTCGTGAGAAGGAGCTGGAGCGTGCTATTCAGGTTCTCTGCCGTCGCCGTAAAAACAACCCGCTGCTGGTGGGGGAATCTGGTGTCGGTAAAACCGCGATTGCGGAAGGTCTTGCCTGGCGAATTGTTCAGGGCGATGTGCCGGAAGTGATGGCTGACTGTACGATTTACTCTCTCGATATCGGTTCTCTGTTAGCGGGCACAAAATA‑CGCGGCGACTTTGAAAAACGTTTTAAAGCGTTGCTCAag                                                                                                                                                                                                            <  2:392522/251‑1 (MQ=255)                                  ttgTCGTGAGAAGGAGCTGGAGCGTGCTATGCAGGTTCTCTGCCGTCGCCGTAAAAACAACCCGCTGCTGGTGGGGGAATCTGGTGTCGGTAAAACCGCGATTGCGGAAGGTCTTGCCTGGCGAATTGTTCAGGGCGATGTGCCGGAAGTGATGGCTGACTGTACGATTTACTCTCTCGATATCGGTTCTCTGTTAGCGGGCACAAAATA‑CGCGGCGACTTTGAAAAACGTTTTAAAGCGTTGCTCAagc                                                                                                                                                                                                           <  2:362588/248‑1 (MQ=255)                                       gTGAGAAGGAGCTGGAGCGTGCTATTCAGGTTCTCTGCCGTCGCCGTAAAAACAACCCGCTGCTGGTGGGGGAATCTGGTGTCGGTAAAACCGCGATTGCGGAAGGTCTTGCCTGGCGAATTGTTCAGGGCGATGTGCCGGAAGTGATGGCTGACTGTACGATTTACTCTCTCGATATCGGTTCTCTGTTAGCGGGCACAAAATA‑CGCGGCGACTTTGAAAAACGTTTTAAAGCGTTGCTCAAGCAGCTgg                                                                                                                                                                                                     >  2:62088/1‑251 (MQ=255)                                               gAGCTGGAGCGTGCTATTCAGGTTCTCTGCCGTCGCCGTAAAAACAACCCGCTGCTGGTGGGGGAATCTGGTGTCGGTAAAACCGCGATTGCGGAAGGTCTTGCCTGGCGAATTGTTCAGGGCGATGTGCCGGAAGTGATGGCTGACTGTACGATTTACTCTCTCGATATCGGTTCTCTGTTAGCGGGCACAAAATA‑CGCGGCGACTTTGAAAAACGTTTTAAAGCGTTGCTCAAGCAGCTGGAGCAGGac                                                                                                                                                                                             <  1:252663/251‑1 (MQ=255)                                                                        tctGCCGTCGCCGTAAAAACAACCCGCTGCTGGTGGGGGAATCTGGTGTCGGTAAAACCGCGATTGCGGAAGGTCTTGCCTGGCGAATTGTTCAGGGCGATGTGCCGGAAGTGATGGCTGACTGTACGATTTACTCTCTCGATATCGGTTCTCTGTTAGCGGGCACAAAATA‑CGCGGCGACTTTGAAAAACGTTTTAAAGCGTTGCTCAAGCAGCTGGAGCAGGACACTAACAGCATCCTGTTTATtgatg                                                                                                                                                                    >  2:22447/1‑251 (MQ=255)                                                                           gCCGTCGCCGTAAAAACAACCCGCTGCTGGTGGGGGAATCTGGTGTCGGTAAAACCGCGATTGCGGAAGGTCTTGCCTGGCGAATTGTTCAGGGCGATGTGCCGGAAGTGATGGCTGACTGTACGATTTACTCTCTCGATATCGGTTCTCTGTTAGCGGGCACAAAATA‑CGCGGCGACTTTGAAAAACGTTTTAAAGCGTTGCTCAAGCAGCTGGAGCAGGACACTAACAGCATCCTGTTTATTGATgag                                                                                                                                                                  <  1:9845/250‑1 (MQ=255)                                                                            ccGTCGCCGTAAAAACAACCCGCTGCTGGTGGGGGAATCTGGTGTCGGTAAAACCGCGATTGCGGAAGGTCTTGCCTGGCGAATTGTTCAGGGCGATGTGCCGGAAGTGATGGCTGACTGTACGATTTACTCTCTCGATATCGGTTCTCTGTTAGCGGGCACAAAATA‑CGCGGCGACTTTGAAAAACGTTTTAAAGCGTTGCTCAAGCAGCTGGAGCAGGACACTAACAGCATCCTGTTTATTGATGAGAt                                                                                                                                                                <  1:293243/251‑1 (MQ=255)                                                                                cGCCGTAAAAACAACCCGCTGCTGGTGGGGGAATCTGGTGTCGGTAAAACCGCGATTGCGGAAGGTCTTGCCTGGCGAATTGTTCAGGGCGATGTGCCGGAAGTGATGGCTGACTGTACGATTTACTCTCTCGATATCGGTTCTCTGTTAGCGGGCACAAAATA‑CGCGGCGACTTTGAAAAACGTTTTAAAGCGTTGCTCAAGCAGCTGGAGCAGGACACTAACAGCATCCTGTTTATTGATGAGATCca                                                                                                                                                             <  1:203476/250‑1 (MQ=255)                                                                                  ccGTAAAAACAACCCGCTGCTGGTGGGGGAATCTGGTGTCGGTAAAACCGCGATTGCGGAAGGTCTTGCCTGGCGAATTGTTCAGGGCGATGTGCCGGAAGTGATGGCTGACTGTACGATTTACTCTCTCGATATCGGTTCTCTGTTAGCGGGCACAAAATA‑CGCGGCGACTTTGAAAAACGTTTTAAAGCGTTGCTCAAGCAGCTGGAGCAGGACACTAACAGCATCCTGTTTATTGATGAGATCCacac                                                                                                                                                          <  1:281078/251‑1 (MQ=255)                                                                                         aacaacCCGCTGCTGGTGGGGGAATCTGGAGTCGGTAAAAAAGCGATTGCGGAAGGTCTTGCCTGGCGAATTGCTCAGGGCGATGTGCCGGCAATGATGGCTGAATGTACCATTTACTCTCTCGATATCGGATCTCTGTTAGCGGGCACACAATA‑CGCGGCGACTTTGAAAAACGTTTTAAAGCGTTGCTCAAGCAGCTGGAGCAGGACACTAACAGCATCCTGTTTATTGATGAGATACACACCATTATc                                                                                                                                                   >  1:393789/1‑251 (MQ=255)                                                                                                              gAATCTGGTGTCGGTAAAACCGCGATTGCGGAAGGTCTTGCCTGGCGAATTGTTCAGGGCGATGTGCCGGAAGTGATGGCTGACTGTACGATTTACTCTCTCGATATCGGTTCTCTGTTAGCGGGCACAAAATA‑CGCGGCGACTTTGAAAAACGTTTTAAAGCGTTGCTCAAGCAGCTGGAGCAGGACACTAACAGCATCCTGTTTATTGATGAGATCCACACCATTATCGGTGCGGGTGCAGCGTCtggt                                                                                                                              <  2:93857/251‑1 (MQ=255)                                                                                                                 tCTGGTGTCGGTAAAACCGCGATTGCGGAAGGTCTTGCCTGGCGAATTGTTCAGGGCGATGTGCCGGAAGTGATGGCTGACTGTACGATTTACTCTCTCGATATCGGTTCTCTGTTAGCGGGCACAAAATA‑CGCGGCGACTTTGAAAAACGTTTTAAAGCGTTGCTCAAGCAGCTGGAGCAGGACACTAACAGCATCCTGTTTATTGATGAGATCCACACCATTATCGGTGCGGGTGCAGCGTCtg                                                                                                                                >  2:127786/1‑246 (MQ=255)                                                                                                                 tCTGGTGTCGGTAAAACCGCGATTGCGGAAGGTCTTGCCTGGCGAATTGTTCAGGGCGATGTGCCGGAAGTGATGGCTGACTGTACGATTTACTCTCTCGATATCGGTTCTCTGTTAGCGGGCACAAAATA‑CGCGGCGACTTTGAAAAACGTTTTAAAGCGTTGCTCAAGCAGCTGGAGCAGGACACTAACAGCATCCTGTTTATTGATGAGATCCACACCATTATCGGTGCGGGTGCAGCGTCtg                                                                                                                                <  1:127786/246‑1 (MQ=255)                                                                                                                                    cgATTGCGGAAGGTCTTGCCTGGCGAATTGTTCAGGGCGATGTGCCGGAAGTGATGGCTGACTGTACGATTTACTCTCTCGATATCGGTTCTCTGTTAGCGGGCACAAAATA‑CGCGGCGACTTTGAAAAACGTTTTAAAGCGTTGCTCAAGCAGCTGGAGCAGGACACTAACAGCATCCTGTTTATTGATGAGATCCACACCATTATCGGTGCGGGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTaaac                                                                                                        <  2:58326/251‑3 (MQ=255)                                                                                                                                                 tCTTGCCTGGCGAATTGTTCAGGGCGATGTGCCGGAAGTGATGGCTGACTGTACGATTTCCTCTCTCGATATCGGTTCTCTGTTAGCGGGCACACAATA‑CGCGGCGACTTTGAAACACGTTTTAAAGCGCTGCTCAAGCAGATGGAGCAGGACACTAACAGCATCCTGTTTATTGATGAGATCCACACCATTATCGGTGCGGGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAATCAAACCGtt                                                                                           >  1:269261/1‑251 (MQ=255)                                                                                                                                                      ccTGGCGAATTGTTCAGGGCGATGTGCCGGAAGTGATGGCTGACTGTACGATTTACTCTCTCGATATCGGTTCTCTGTTAGCGGGCACAAAATA‑CGCGGCGACTTTGAAAAACGTTTTAAAGCGTTGCTCAAGCAGCTGGAGCAGGACACTAACAGCATCCTGTTTATTGATGAGATCCACACCATTATCGGTGCGGGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAATCAAACCGTTGctct                                                                                      >  1:85845/1‑251 (MQ=255)                                                                                                                                                             aaTTGTTCAGGGCGATGTGCCGGAAGTGATGGCTGTCTGTACGATTTACTCTCTCGATATCGGTTCTCTGTTAGCGGGCACAAAATA‑CGCGGCGACTTTGAAAAACGTTTTAAAGCGTTGCTCAAGCAGCTGGAGCAGGACACTAACAGCATCCTGTTTATTGATGAGATCCACACCATTATCGGTGCGGGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAATCAAACCGTTGCTCTCCAGCgg                                                                               <  2:47416/251‑1 (MQ=255)                                                                                                                                                                      gggCGATGTGCCGGAAGAGATGGCTGACTGTACGATTTACTCTCTCGATATCGGTTCTCTGTTAGCGGGCACAAAATA‑CGCGGCGACTTTGAAAAACGTTTTAAAGCGTTGCTCAAGCAGCTGGAGCAGGACACTAACAGCATCCTGTTTATTGATGAGATCCACACCATTATCGGTGCGGGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCg                                                                      >  2:200503/1‑251 (MQ=255)                                                                                                                                                                          gATGTGCCGGAAGTGATGGCTGACTGTACGATTTACTCTCTCGATATCGGTTCTCTGTTAGCGGGCACAAAATA‑CGCGGCGACTTTGAAAAACGTTTTAAAGCGTTGCTCAAGCAGCTGGAGCAGGac                                                                                                                                                                                             >  1:300139/1‑128 (MQ=255)                                                                                                                                                                          gATGTGCCGGAAGTGATGGCTGACTGTACGATTTACTCTCTCGATATCGGTTCTCTGTTAGCGGGCACAAAATA‑CGCGGCGACTTTGAAAAACGTTTTAAAGCGTTGCTCAAGCAGCTGGAGCAGGac                                                                                                                                                                                             <  2:300139/128‑1 (MQ=255)                                                                                                                                                                             gtgCCGGAAGTGATGGCTGACTGTACGATTTACTCTCTCGATATCGGTTCTCTGTTAGCGGGCACAAAATA‑CGCGGCGACTTTGAAAAACGTTTTAAAGCGTTGCTCAAGCAGCTGGAGCAGGACACTAACAGCATCCTGTTTATTGATGAGATCCACACCATTATCGGTGCGGGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAAtt                                                               >  1:87259/1‑251 (MQ=255)                                                                                                                                                                                ccGGAAGTGATGGCTGACTGTACGATTTACTCTCTCGATATCGGTTCTCTGTTAGCGGGCACAAAATA‑CGCGGCGACTTTGAAAAACGTTTTAAAGCGTTGCTCAAGCAGCTGGAGCAGGACACTAACAGCATCCTGTTTATTGATGAGATCCACACCATTATCGGTGCGGGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGt                                                            <  2:151180/251‑1 (MQ=255)                                                                                                                                                                                                aCTGTACGATTTACTCTCTCGATATCGGTTCTCTGTTAGCGGGCACAAAATA‑CGCGGCGACTTTGAAAAACGTTTTAAAGCGTTGCTCAAGCAGCTGGAGCAGGACACTAACAGCATCCTGTTTATTGATGAGATCCACACCATTATCGGTGCGGGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAgg                                            <  2:22017/251‑1 (MQ=255)                                                                                                                                                                                                            aCTCTCTCGATATCGGTTCTCTGTTAGCGGGCACAAAATA‑CGCGGCGACTTTGAAAAACGTTTTAAAGCGTTGCTCAAGCAGCTGGAGCAGGACACTAACAGCATCCTGTTTATTGATGAGATCCACACCATTATCGGTGCGGGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGAGTTCAGCAAc                                 >  2:221098/1‑250 (MQ=255)                                                                                                                                                                                                              tctctcGATATCGGTTCTCTGTTAGCGGGCACAAAATA‑CGCGGCGACTTTGAAAAACGTTTTAAAGCGTTGCTCAAGCAGCTGGAGCAGGACACTAACAg                                                                                                                                                                                     >  2:174306/1‑100 (MQ=255)                                                                                                                                                                                                              tctctcGATATCGGTTCTCTGTTAGCGGGCACAAAATA‑CGCGGCGACTTTGAAAAACGTTTTAAAGCGTTGCTCAAGCAGCTGGAGCAGGACACTAACAg                                                                                                                                                                                     <  1:174306/100‑1 (MQ=255)                                                                                                                                                                                                                 cgcGATATCGGTTATCTGTTAGCGGGCACCAAA‑ACCGCAGCGACTTTGAACAACGTTTTAAAGAGCTGCTCAAGCAGCTGGAGCAGCCCACTACCAGCATCCTGTTTCTTGATGAGATCCACACCATTATCGGTGCGGGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAAACTATCAGGAGTTCAGCAACATTTTc                           >  2:59283/3‑251 (MQ=255)                                                                                                                                                                                                                   cGATATCGGTTCTCTGTTAGCGGGCACAAAATA‑CGCGGCGACTTTGAAAAACGTTTTAAAGCGTTGCTCAAGCAGCTGGAGCAGGACACTAACAGCATCCTGTTTATTGATGAGATCCACACCATTATCGGTGCGGGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGAGTTCAGCAACATTTTCg                          >  2:187058/1‑250 (MQ=255)                                                                                                                                                                                                                     atatCGGTTCTCTGTTAGCGGGCACAAAATA‑CGCGGCGACTTTGAAAAACGTTTTAAAGCGTTGCTCAAGCAGCTGGAGCAGGAAACTAACAGCATCCTGTTTATTGATGAGATCCACACCATTATCGGTGCGGGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGAGTTCAGCAACATTTTCg                          >  2:339031/1‑248 (MQ=255)                                                                                                                                                                                                                         cGGTTCTCTGTTAGCGGGCACAAAATA‑CGCGGCGACTTTGAAAAACGTTTTAAAGCGTTGCTCAAGCAGCTGGAGCAGGACACTAACAGCATCCTGTTTATTGATGAGATCCACACCATTATCGGTGCGGGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAATCAAACCGTTGCTCTCCAGCGGTACAATTCGTGTAATTGGTTCGACAACCTATCAGGAGTTCAGCAACATTTTCGAGaa                      >  1:29166/1‑248 (MQ=255)                                                                                                                                                                                                                           gTTCTCTGTTAGCGGGCACAAAATA‑CGCGGCGACTTTGAAAAACGTTTTAACGCGATGCTCAAGCAGCTGGAGCAGGACACTAACAGCATCCTGTTTATTGATGAGATCCACACCATTATCGGTGAGGGTGCAGCGTCTGGAGGTCAGGTCGATGCGGCTAACCTAATCAAACCGTTGCTCTCCCGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGAGTTCAGCAACATTTTCGAGAAAGAcc                 >  1:403467/1‑251 (MQ=255)                                                                                                                                                                                                                                          gCACAGAATA‑CGCGGCGACTTTGAAAAACGTTTTCAAGCGTTGCTCAAGCAGCTGGAGCAGGACACTAACAGCATCCTGTTTATTGATGAGATCCACACCATTATCGGTGCGGGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGAGTTCAGCAACATTTTCGAGAAAGACCGTGCTCTGGCGCGTc  >  2:64368/1‑251 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  GCTTGCGCGCGTGGGCGGAATCGACCCACTGATTGGTCGTGAGAAGGAGCTGGAGCGTGCTATTCAGGTTCTCTGCCGTCGCCGTAAAAACAACCCGCTGCTGGTGGGGGAATCTGGTGTCGGTAAAACCGCGATTGCGGAAGGTCTTGCCTGGCGAATTGTTCAGGGCGATGTGCCGGAAGTGATGGCTGACTGTACGATTTACTCTCTCGATATCGGTTCTCTGTTAGCGGGCACAAAATATCGCGGCGACTTTGAAAAACGTTTTAAAGCGTTGCTCAAGCAGCTGGAGCAGGACACTAACAGCATCCTGTTTATTGATGAGATCCACACCATTATCGGTGCGGGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGAGTTCAGCAACATTTTCGAGAAAGACCGTGCTCTGGCGCGTC  >  NC_000913/923797‑924281

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