Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	495,622	C→T	T168I (ACC→ATC)	htpG →	chaperone protein HtpG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	495,622	0	C	T	80.0%	102.6 / 25.0	51	T168I (ACC→ATC)	htpG	chaperone protein HtpG
Reads supporting (aligned to +/- strand):  ref base C (4/6);  new base T (19/21);  total (23/28)
Fisher's exact test for biased strand distribution p-value = 7.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.52e-01

CGCGACGAAGTGATTGACCATCTGGGGACTATCGCTAAATCCGGTACCAAATCATTCCTCGAATCCCTGGGTTCTGACCAGGCGAAAGACAGCCAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCACCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTCCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACTCGTAACAAGTCGGAAATCACCGATGAAGAGTACAAAGAGTTCT  >  NC_000913/495378‑495867                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      cgcgACGAAGTGATTGACCATCTGGGGACTATCGCTAAATCCGGTACCAAATCATTCCTCGAATCCCTGGGTTCTGACCAGGCGAAAGACAGCCAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCATCAaaga                                                                                                                                                                                                                                                 <  2:399759/251‑1 (MQ=255)        aaGTTATTGACCATCTGGGGACTATCGCTAAATCCGGTACCAAATCATTCCTCGAATCCCTGGGTTCTGACCAGGCGAAAGACAGCCAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCATCAAAGAAGATCGt                                                                                                                                                                                                                                          <  1:230218/251‑1 (MQ=255)        aaGTGATTGACCATCTGGGGACTATCGCTAAATCCGGTACCAAATCATTCCTCGAATCCCTGGGTTCTGACCAGGCGAAAGACAGCCAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCATCAaaga                                                                                                                                                                                                                                                 >  2:65082/1‑244 (MQ=255)        aaGTGATTGACCATCTGGGGACTATCGCTAAATCCGGTACCAAATCATTCCTCGAATCCCTGGGTTCTGACCAGGCGAAAGACAGCCAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCATCAAAGAAGATCGt                                                                                                                                                                                                                                          <  2:337218/251‑1 (MQ=255)            gATTGACCATCTGGGGACTATCGCTAAATCCGGTTCCAAATCATTCCTCGAACCCGTGGGTTCTGACCAGGCTAAAGACAGCCAGCTGATCGGTCAGTGTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGCCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCATCAAAGAAGATCGTGGt                                                                                                                                                                                                                                       <  1:207043/250‑1 (MQ=255)                   cATCTGGGGACTATCGCTAAATCCGGTACCAAATCATTCCTCTAATCCCTGGGTTCTGACCAGTCGAAAGACAGCCAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAAt                                                                                                                                                                                                                               <  1:650/251‑1 (MQ=255)                         gggACTATCGCTAAATCCGGTACCAAATCATTCCTCGAATCCCTGGGTTCTGACCAGGCGAAAGACAGCCAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGGGAATACAC‑CGGTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCt                                                                                                                                                                                                                         >  1:452075/1‑251 (MQ=255)                          ggACTATCGCTAAATCCGGTACCAAATCATTCCTCGAATCCCTTGGTTCTGACCAGGCGAAAGACAGCCAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTGCTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCACCAAAGAAGATCGTGGTACTGAAATCACCCTg                                                                                                                                                                                                                        <  2:126123/251‑1 (MQ=255)                          ggACTATCGCTAAATCCGGTACCAAATCATTCCTCGAATCCCTGGGTTCTGACCAGGCGAAAGACAGCCAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCACCAAAGAAGATCGTGGTACTGAAATCACCCTg                                                                                                                                                                                                                        <  1:233955/251‑1 (MQ=255)                                      aaTCCGGTACCAAATCATTCCTCGAATCCCTGGGTTCTGACCAGGCGAAAGACAGCCAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGaa                                                                                                                                                                                                            <  2:396899/251‑1 (MQ=255)                                         ccGGTACCAAATCATTCCTCGAATCCCTGGGTTCTGACCAGGCGAAAGACAGCCAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGc                                                                                                                                                                                                         >  1:54138/1‑251 (MQ=255)                                                       ttCCTCGAATCCCTGGGTTCTGACCAGGCGAAAGACAGCCAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGtttct                                                                                                                                                                                           >  2:120558/1‑248 (MQ=255)                                                              aaTCCCTGGGTTCTGACCAGGCGAAAGACAGCCAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCACCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTCCTCGATGAc                                                                                                                                                                                    <  1:145781/251‑1 (MQ=255)                                                               aTCCCTGGGTTCTGACCAGGCGAAAGACAGCCAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACt                                                                                                                                                                                   <  1:173384/251‑1 (MQ=255)                                                                 cccTGGGTTCTGACCAGGCGAAAGACAGCCAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTgg                                                                                                                                                                                 <  1:472362/251‑1 (MQ=255)                                                                  ccTGGGTTCTGACCAGGCGAAAGACAGCCAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCACCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTCCTCGATGACTGgc                                                                                                                                                                                <  1:283799/251‑1 (MQ=255)                                                                    tGGGTTCTGACCAGGCGAAAGACAGCCAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTGgcgc                                                                                                                                                                              >  1:376148/1‑251 (MQ=255)                                                                         tCTGACCAGGCGAAAGACAGCCAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCACCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTCCTCGATGACTGGCgcgtgc                                                                                                                                                                          <  1:249656/250‑1 (MQ=255)                                                                               cAGGCGAAAGACAGCCAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGCCAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAACGCGAAGACGAGTTTCTCGATGACTGGCGCGTGCGTTCcat                                                                                                                                                                   >  2:54351/1‑251 (MQ=255)                                                                                    gAAAGACAGCCAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTGGCGCGTGCGTTCCatcatc                                                                                                                                                               >  2:391683/1‑250 (MQ=255)                                                                                    gAAAGACAGCCAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTGGCGCGTGCGTTCCAtcatca                                                                                                                                                              <  2:331723/251‑1 (MQ=255)                                                                                          cagccagcTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAAt                                                                                                                                                        >  2:242512/1‑251 (MQ=255)                                                                                                       ggTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCatat                                                                                                                                           <  1:235149/251‑1 (MQ=255)                                                                                                          cAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATcgc                                                                                                                                        <  1:166940/251‑1 (MQ=255)                                                                                                               tGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGc                                                                                                                                   >  2:463503/1‑251 (MQ=255)                                                                                                                ggTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCACCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTCCTCGATGACTGGCGCGTGCGGTCCATCATCAGCAAATACTCCGACCATATCGCGCTGcc                                                                                                                                  >  2:266311/1‑251 (MQ=255)                                                                                                                   gtTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTGGGGCGTGCGTGCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGt                                                                                                                               <  1:24618/251‑1 (MQ=255)                                                                                                                      ggTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCATCAAAGAAGATCTTGGTACTGAAATCACCCTGCATCTTCGTGAAGGCGAAGACGAGTTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTag                                                                                                                             <  1:65655/250‑1 (MQ=255)                                                                                                                             aCTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATc                                                                                                                        <  1:54351/248‑1 (MQ=255)                                                                                                                                  gCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAAc                                                                                                                 <  1:463503/250‑1 (MQ=255)                                                                                                                                      ttATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTGGCgcgtgcgt                                                                                                                                                                        <  2:421157/191‑1 (MQ=255)                                                                                                                                      ttATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTGGCgcgtgcgt                                                                                                                                                                        >  1:421157/1‑191 (MQ=255)                                                                                                                                            tGGCCGACAAAGTGACCGTGCGTACTCGCGCGCCAGGCGAAAAACCAGAAAATCGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACgcgaagcgaa                                                                                                       >  1:28733/1‑250 (MQ=255)                                                                                                                                               ccGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCACCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTCCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGAc                                                                                                   <  1:131573/251‑1 (MQ=255)                                                                                                                                                 gacaacgtgaccccgctaaCTCGCGCGGCAGGCGAAAAACCCGAGAATGGCGTCTTATGGGCATCGGCTGGCGAAGGAGAATACACACATT‑CCGACATCACCCAAGAAGATCGTGGTACTGAAATCACCCTGCCTTTGCGTATAGGCGAAGACGAGTTACTCGATGAATTGCCCTTGCGTTCCCTCAACCGCAAACACTCCGACCATATAGCGCTGCCGGTAGAGGGCGAAAAACGCGCAGAGAAAGACg                                                                                                  >  2:266003/19‑250 (MQ=255)                                                                                                                                                             gTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCATCAAAGAAGATCGTGTTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTAt                                                                                     <  2:201177/251‑1 (MQ=255)                                                                                                                                                                              aGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGACGGTGAATACAC‑CGTTGCCGACATCACCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGGTCCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGAGAAAACGTTATCTCCTGGGAGaaactc                                                                     >  2:359107/1‑247 (MQ=255)                                                                                                                                                                              aGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATc                                                                     >  1:244332/1‑250 (MQ=255)                                                                                                                                                                                      aaCCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAgcg                                                            >  1:490870/1‑251 (MQ=255)                                                                                                                                                                                                    gTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCACCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTCCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGAc                                              >  2:327393/1‑251 (MQ=255)                                                                                                                                                                                                           gggAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACTCGTaa                                        <  1:120558/250‑1 (MQ=255)                                                                                                                                                                                                                      ggCGAAGGTGAATACAC‑CGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTGATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACTCGTAACAAGTCGGAAAt                            <  2:490870/251‑1 (MQ=255)                                                                                                                                                                                                                          aaGGTGAATACAC‑CGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACTCGTAACAAGTCGGAAATCAcc                        >  2:108079/1‑251 (MQ=255)                                                                                                                                                                                                                           aGGTGAATACAC‑CGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGAAGGCGCTGTGGACTCGTAACAAGTCGGAAATCACCg                       >  2:347552/1‑251 (MQ=255)                                                                                                                                                                                                                               gAATACAC‑CGTTGCCGCCATCATCAAAGAAGATCGTGGTACTGACATCAACCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGCGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACt                                             >  2:269562/1‑225 (MQ=255)                                                                                                                                                                                                                               gAATACAC‑CGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACAGCGAAACCGTTATCTCCTGGGAGAAAATCAGCAAAGCGCAGGCGCTGTGGACt                                             <  1:269562/225‑1 (MQ=255)                                                                                                                                                                                                                                 attcac‑CGTTGCCGACATCAGCAAAGAAGATCGGGGTACTGAAATCACCCTGCATCTGCGTGAAGGCAAAGACGAGTTTCTCGATGACTGGCGCCTGCGTTCCATCATCAGCAAATACTCCGACCATAGCGCGCTGCCGGTATATATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTAGCTCCTGGGAGAAAAGCAACAAAGCGCAGGCGCTGTTGTCTCTTAACAAGTCGGAAATCACCGATGAag                 <  2:462394/248‑1 (MQ=255)                                                                                                                                                                                                                                  tACAC‑CGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACTCGTAACAAGTCGGAAATCACCGATGa                   >  2:114839/1‑248 (MQ=255)                                                                                                                                                                                                                                   acac‑CGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACTCGTAACAAGTCGGAAATCACCGATGAag                 >  1:3104/1‑249 (MQ=255)                                                                                                                                                                                                                                          ttGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCg                                                                                                                                 <  1:158360/131‑1 (MQ=255)                                                                                                                                                                                                                                          ttGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCg                                                                                                                                 >  2:158360/1‑131 (MQ=255)                                                                                                                                                                                                                                                 catcatCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACTCGTAACAAGTCGGAAATCACCGATGAAGAGTACAAAGAGTTCt  <  1:65082/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      CGCGACGAAGTGATTGACCATCTGGGGACTATCGCTAAATCCGGTACCAAATCATTCCTCGAATCCCTGGGTTCTGACCAGGCGAAAGACAGCCAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCACCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTCCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACTCGTAACAAGTCGGAAATCACCGATGAAGAGTACAAAGAGTTCT  >  NC_000913/495378‑495867

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