Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,701,985	C→A	6.2%	L139M (CTG→ATG) ‡	PP_3267 →	Clp protease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,701,985	0	C	A	6.2%	260.7 / 5.1	81	L139M (CTG→ATG) ‡	PP_3267	Clp protease
Reads supporting (aligned to +/- strand):  ref base C (40/36);  new base A (3/2);  total (43/38)
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.50e-01

CTCTTCCTTGCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATCGACGAGATGCC  >  NC_002947/3701847‑3702127                                                                                                                                           |                                                                                                                                                ctctTCCTTGCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGTATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCa                                                                                                                                       <  8:244399/149‑1 (MQ=255) ctctTCCTTGCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATgcg                                                                                                                                         >  6:134361/1‑147 (MQ=255) ctctTCCTTGCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATgcg                                                                                                                                         >  6:318040/1‑147 (MQ=255) ctctTCCTTGCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCa                                                                                                                                       <  5:237757/149‑1 (MQ=255)  tctTCCTTGCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATgcg                                                                                                                                         >  3:244118/1‑146 (MQ=255)    ttCCTTGCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATgcg                                                                                                                                         >  3:254079/1‑144 (MQ=255)      ccTTGCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATgcgc                                                                                                                                        <  3:318293/143‑1 (MQ=255)      ccTTGCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATgcgc                                                                                                                                        >  4:318293/1‑143 (MQ=255)        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)                                        aGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTAt                                                                                                                                            >  4:279476/1‑105 (MQ=255)                                        aGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGtttt                                                                                                >  2:263364/1‑149 (MQ=255)                                          cAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTg                                                                                               >  3:1804/1‑148 (MQ=255)                                           aaGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGa                                                                                             <  5:318040/149‑1 (MQ=255)                                           aaGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGa                                                                                             <  1:263364/149‑1 (MQ=255)                                                  ttttCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGt                                                                                                                                              <  2:98626/93‑1 (MQ=255)                                                  ttttCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGt                                                                                                                                              >  1:98626/1‑93 (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                                                                                <  3:274891/147‑1 (MQ=255)                                                          ccGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTAt                                                                                >  4:274891/1‑147 (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                                      <  4:177407/149‑1 (MQ=255)                                                                                                 cgGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTGGCGCCGCGGGAGGCGCt                                      <  7:233047/149‑1 (MQ=255)                                                                                                   gATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATCTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTGGCGCCGCGGGAGGCGCTgg                                    <  4:1804/149‑1 (MQ=255)                                                                                                        ccGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTGGCGCCGCGGGAGGCGCTGGAAAAc                               >  7:245398/1‑149 (MQ=255)                                                                                                         cGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACg                              <  7:57256/149‑1 (MQ=255)                                                                                                         cGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACg                              >  2:270413/1‑149 (MQ=255)                                                                                                          gAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGc                                                                        >  3:242293/1‑107 (MQ=255)                                                                                                          gAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGc                                                                        <  4:242293/107‑1 (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                            >  6:299709/1‑144 (MQ=255)                                                                                                                gCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACggg                            <  5:299709/144‑1 (MQ=255)                                                                                                                       gCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTGGCGCCGCGGGAGGCGCTGGAAAAc                               <  1:19058/134‑1 (MQ=255)                                                                                                                       gCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTGGCGCCGCGGGAGGCGCTGGAAAAc                               >  2:19058/1‑134 (MQ=255)                                                                                                                       gCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGa                 >  2:16068/1‑148 (MQ=255)                                                                                                                        ccTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACgg                             >  3:66970/1‑135 (MQ=255)                                                                                                                        ccTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACgg                             <  4:66970/135‑1 (MQ=255)                                                                                                                        ccTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGgatgat                       <  6:270855/141‑1 (MQ=255)                                                                                                                        ccTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGgatgat                       <  4:314963/141‑1 (MQ=255)                                                                                                                        ccTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGgatgat                       >  3:314963/1‑141 (MQ=255)                                                                                                                        ccTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGgatgat                       >  5:270855/1‑141 (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                                      <  6:353215/117‑1 (MQ=255)                                                                                                                                 cGATCTGCGCATACAGGCGCAGGTCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTGGCGCCGCGGGAGGCGCt                                      >  5:353215/1‑117 (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)                                                                                                                                           |                                                                                                                                                CTCTTCCTTGCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCC‑AAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATCGACGAGATGCC  >  NC_002947/3701847‑3702127

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: ‑