Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,701,997	A→G	6.8%	I143V (ATC→GTC)	PP_3267 →	Clp protease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,701,997	0	A	G	6.8%	205.7 / 4.5	73	I143V (ATC→GTC)	PP_3267	Clp protease
Reads supporting (aligned to +/- strand):  ref base A (35/33);  new base G (3/2);  total (38/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.18e-01

TTGCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATCGACGAGATGCCGATCGCGTCAGACGCAA  >  NC_002947/3701854‑3702144                                                                                                                                                |                                                                                                                                                     ttGCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTg                                                                                                                                                 >  7:12122/1‑149 (MQ=255)   gCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTgg                                                                                                                                                <  4:254079/148‑1 (MQ=255)   gCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGc                                                                                                                                               >  1:317884/1‑149 (MQ=255)    ccGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCg                                                                                                                                              <  8:12122/149‑1 (MQ=255)    ccGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCg                                                                                                                                              <  5:134361/149‑1 (MQ=255)           gCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCa                                                                                                                                       <  3:122626/149‑1 (MQ=255)                       ttGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCgg                                                                                                                           >  8:233047/1‑149 (MQ=255)                             tttCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTgg                                                                                                                                                >  1:26232/1‑122 (MQ=255)                             tttCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTgg                                                                                                                                                <  2:26232/122‑1 (MQ=255)                                 aGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGtttt                                                                                                                 >  2:263364/1‑149 (MQ=255)                                   cAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTg                                                                                                                >  3:1804/1‑148 (MQ=255)                                    aaGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGa                                                                                                              <  1:263364/149‑1 (MQ=255)                                    aaGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGa                                                                                                              <  5:318040/149‑1 (MQ=255)                                             ttCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACcgc                                                                                                      <  4:244118/148‑1 (MQ=255)                                                cACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGtt                                                                                                   <  2:317884/148‑1 (MQ=255)                                                cACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTa                                                                                                  >  8:92007/1‑149 (MQ=255)                                                   ccGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTg                                                                                                                >  5:13881/1‑132 (MQ=255)                                                   ccGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTg                                                                                                                <  6:13881/132‑1 (MQ=255)                                                   ccGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTAt                                                                                                 >  4:274891/1‑147 (MQ=255)                                                   ccGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTAt                                                                                                 <  3:274891/147‑1 (MQ=255)                                                           ttGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGcc                                                                                        >  7:95080/1‑148 (MQ=255)                                                                  ttGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATcc                                                                                   >  3:246297/1‑146 (MQ=255)                                                                    gCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑aa                                                                                <  4:246297/146‑1 (MQ=255)                                                                    gCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAgg                                                                              <  8:95080/148‑1 (MQ=255)                                                                        ggTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGAGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAGGATACt                                                                         >  3:54971/1‑149 (MQ=255)                                                                         gTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTg                                                                        >  8:80732/1‑149 (MQ=255)                                                                          ttCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTg                                                                                                                <  6:77094/109‑1 (MQ=255)                                                                            cGGTGGTCCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTGgc                                                                      <  4:54971/148‑1 (MQ=255)                                                                              gTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTGGcgcc                                                                   >  4:88931/1‑149 (MQ=255)                                                                              gTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTGGcgcc                                                                   >  3:177407/1‑149 (MQ=255)                                                                                 gACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTGGCGCcgcg                                                                <  3:88931/149‑1 (MQ=255)                                                                                     aTTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTGGCGCCGCGGGAg                                                            >  8:57256/1‑149 (MQ=255)                                                                                         gcgGATGTCCGAGTATGCGATGAGACTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTGGCGCCGCGGGAGgcgc                                                        >  8:75118/1‑149 (MQ=255)                                                                                          cgGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTGGCGCCGCGGGAGGCGCt                                                       <  7:233047/149‑1 (MQ=255)                                                                                          cgGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTGGCGCCGCGGGAGGCGCt                                                       <  4:177407/149‑1 (MQ=255)                                                                                            gATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATCTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTGGCGCCGCGGGAGGCGCTgg                                                     <  4:1804/149‑1 (MQ=255)                                                                                                 ccGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTGGCGCCGCGGGAGGCGCTGGAAAAc                                                >  7:245398/1‑149 (MQ=255)                                                                                                  cGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACg                                               >  2:270413/1‑149 (MQ=255)                                                                                                  cGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACg                                               <  7:57256/149‑1 (MQ=255)                                                                                                   gAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGc                                                                                         <  4:242293/107‑1 (MQ=255)                                                                                                   gAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGc                                                                                         >  3:242293/1‑107 (MQ=255)                                                                                                   gAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACgg                                              <  7:75118/149‑1 (MQ=255)                                                                                                         gCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATg                                                                                                            >  3:192790/1‑82 (MQ=255)                                                                                                         gCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATg                                                                                                            <  4:192790/82‑1 (MQ=255)                                                                                                         gCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACggg                                             <  5:299709/144‑1 (MQ=255)                                                                                                         gCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACggg                                             >  6:299709/1‑144 (MQ=255)                                                                                                                gCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTGGCGCCGCGGGAGGCGCTGGAAAAc                                                >  2:19058/1‑134 (MQ=255)                                                                                                                gCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTGGCGCCGCGGGAGGCGCTGGAAAAc                                                <  1:19058/134‑1 (MQ=255)                                                                                                                gCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGa                                  >  2:16068/1‑148 (MQ=255)                                                                                                                 ccTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACgg                                              <  4:66970/135‑1 (MQ=255)                                                                                                                 ccTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACgg                                              >  3:66970/1‑135 (MQ=255)                                                                                                                 ccTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGgatgat                                        >  3:314963/1‑141 (MQ=255)                                                                                                                 ccTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGgatgat                                        <  6:270855/141‑1 (MQ=255)                                                                                                                 ccTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGgatgat                                        >  5:270855/1‑141 (MQ=255)                                                                                                                 ccTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGgatgat                                        <  4:314963/141‑1 (MQ=255)                                                                                                                 ccTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGaaa                                >  6:21869/1‑149 (MQ=255)                                                                                                                 ccTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGaaa                                >  4:282453/1‑149 (MQ=255)                                                                                                                   tGGATTACGACCTGCGCCTGTATGCGCAGATAGTGGCGGAGCGCACGAAAGGCGCGGCCGAGCTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGCATGATCCATGAAATc                              >  2:142269/1‑149 (MQ=255)                                                                                                                     gATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATcga                            <  5:21869/149‑1 (MQ=255)                                                                                                                      aTTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATcgac                           <  8:245398/149‑1 (MQ=255)                                                                                                                       ccaCGATCTGCGCATACAGGCGCAGGTCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATcgac                           >  7:324329/3‑148 (MQ=255)                                                                                                                          cGATCTGCGCATACAGGCGCAGGTCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTGGCGCCGCGGGAGGCGCt                                                       >  5:353215/1‑117 (MQ=255)                                                                                                                          cGATCTGCGCATACAGGCGCAGGTCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTGGCGCCGCGGGAGGCGCt                                                       <  6:353215/117‑1 (MQ=255)                                                                                                                          cGATCTGCGCATACAGGCGCAGGTCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACg                                               <  8:276019/125‑1 (MQ=255)                                                                                                                          cGATCTGCGCATACAGGCGCAGGTCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACg                                               >  7:276019/1‑125 (MQ=255)                                                                                                                          cGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCTCAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATcgacga                         <  1:270413/148‑1 (MQ=255)                                                                                                                          cGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATcgac                           >  2:16461/1‑145 (MQ=255)                                                                                                                              cTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGGAATCGACGAGATGcc                   >  3:59403/1‑149 (MQ=255)                                                                                                                               tGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATCGACGAGATGcc                   >  1:96135/1‑148 (MQ=255)                                                                                                                                  gcCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTAt                                                                                                 <  6:53193/68‑1 (MQ=255)                                                                                                                                        aTGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATCGACGAGATGCCGATCGCGt           >  2:87697/1‑147 (MQ=255)                                                                                                                                            gcAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATCGACGAGATGCCGATCGCGTCAGacg     >  3:342694/1‑149 (MQ=255)                                                                                                                                               gATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGTTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATCGACGAGATGCCGATCGCGTCAGacgcaa  <  1:142269/149‑1 (MQ=255)                                                                                                                                                aTCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATCGACGAGATGc                    <  1:306189/130‑1 (MQ=255)                                                                                                                                                aTCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATCGACGAGATGc                    >  2:306189/1‑130 (MQ=255)                                                                                                                                                |                                                                                                                                                     TTGCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATCGACGAGATGCCGATCGCGTCAGACGCAA  >  NC_002947/3701854‑3702144

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