Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,615,343	G→A	7.1%	V572M (GTG→ATG) ‡	PP_4085 →	rhs family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,615,343	0	G	A	7.1%	205.3 / 4.2	70	V572M (GTG→ATG) ‡	PP_4085	rhs family protein
Reads supporting (aligned to +/- strand):  ref base G (36/29);  new base A (2/3);  total (38/32)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

CATCAACTGGCTAGCGGTGCTGTCCGCGAATACGAATGGGCGCGCTTCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCACTACTACTGGGGGCCACTGTACGAGGTCTATAAGCACATCGACCC  >  NC_002947/4615199‑4615488                                                                                                                                                 |                                                                                                                                                  cATCAACTGGCTAGCGGTGCTGTCCGCGAATACGAATGGGCGCGCTTCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGAt                                                                                                                                               >  4:25012/1‑149 (MQ=255)   taaaCTGGCTAGCGTTGCTGTCCGCTAATACGAATGGGCGCGCTTCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTc                                                                                                                                             <  3:189386/147‑1 (MQ=255)          gCTAGCGGTGCTGTCCGCGAATACGAATGGGCGCGCTTCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAATCACTCGCCat                                                                                                                                        <  6:2555/147‑2 (MQ=255)           cTAGCGGTGCTGTCCGCGAATACGAATGGGCGCGCTTCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTgg                                                                                                                                     <  8:154843/149‑1 (MQ=255)           cTAGCGGTGCTGTCCGCGAATACGAATGGGCGCGCTTCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTgg                                                                                                                                     >  3:299360/1‑149 (MQ=255)              gCGGTGCTGTCCGCGAATACGAATGGGCGCGCTTCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTgggg                                                                                                                                   >  6:148068/1‑148 (MQ=255)                   gCTGTCCGCGAATACGAATGGGCGCGCTTCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGa                                                                                                                             <  7:57851/149‑1 (MQ=255)                                  gAATGGGCGCGCTTCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGAt                                                                                                                            <  1:270177/135‑1 (MQ=255)                                  gAATGGGCGCGCTTCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGAt                                                                                                                            >  2:270177/1‑135 (MQ=255)                                  gAATGGGCGCGCTTCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGt                                                                                                                         <  1:198180/138‑1 (MQ=255)                                  gAATGGGCGCGCTTCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGt                                                                                                                         >  2:198180/1‑138 (MQ=255)                                       ggCGCGCTTCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTc                                                                                                          <  1:108450/148‑1 (MQ=255)                                              ttCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACaa                                                                                                  >  1:297335/1‑149 (MQ=255)                                              ttCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAATCACTCGCCattc                                                                                                                                      <  4:158966/113‑4 (MQ=255)                                              ttCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAATCACTCGCCattc                                                                                                                                      >  3:158966/1‑110 (MQ=255)                                                 gCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCt                                                                                               >  3:51966/1‑149 (MQ=255)                                                  ccATTCCTGAAACGCGCCCTAGTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTg                                                                                              >  5:328201/1‑149 (MQ=255)                                                   cATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTg                                                                                              >  3:14857/1‑148 (MQ=255)                                                      tCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATg                                                                                                                           >  4:254873/1‑116 (MQ=255)                                                      tCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATg                                                                                                                           <  3:254873/116‑1 (MQ=255)                                                         tGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTg                                                                                              >  6:24506/1‑142 (MQ=255)                                                         tGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTg                                                                                              <  5:24506/142‑1 (MQ=255)                                                         tGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCg                                                                                       >  4:176949/1‑149 (MQ=255)                                                          gAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTACACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCgg                                                                                      >  7:6050/1‑149 (MQ=255)                                                           aagcgcgCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGt                                                                                     <  8:110381/146‑1 (MQ=255)                                                           aaaCGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGt                                                                                     <  4:299360/149‑1 (MQ=255)                                                            aaCGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTg                                                                                    <  4:14857/149‑1 (MQ=255)                                                                 gcCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAAcc                                                                               <  4:228528/149‑1 (MQ=255)                                                                  cccTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTgggg                                                                                                                                   >  2:185345/1‑96 (MQ=255)                                                                  cccTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAAccc                                                                              <  3:176949/149‑1 (MQ=255)                                                                    cTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTc                                                                            <  3:25012/149‑1 (MQ=255)                                                                    cTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCCGGAAGCCGGTGAAACCCTc                                                                            >  5:108331/1‑149 (MQ=255)                                                                      aCTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGt                                                                          >  5:178376/1‑149 (MQ=255)                                                                          ccAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCACCCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGGGAATAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGATACCCTCGTccag                                                                      >  4:66113/1‑146 (MQ=255)                                                                          ccAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAg                                                                                                                                                   >  3:133161/1‑72 (MQ=255)                                                                          ccAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAg                                                                                                                                                   <  4:133161/72‑1 (MQ=255)                                                                          ccAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCAc                                                                       <  5:148068/148‑1 (MQ=255)                                                                          ccAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACg                                                                      >  2:72447/1‑149 (MQ=255)                                                                          ccAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACg                                                                      >  2:69036/1‑149 (MQ=255)                                                                                   tGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTc                                                             >  3:269658/1‑149 (MQ=255)                                                                                     gACGGCACGCCCTACCGCTTGCCTCCACTTCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCg                                                            <  5:205275/148‑1 (MQ=255)                                                                                     gACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCg                                                            <  2:185701/148‑1 (MQ=255)                                                                                       cGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGgcc                                                         <  6:328201/149‑1 (MQ=255)                                                                                         ccACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAAAGGCGAGAGAATCGCCACTGGGGAAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTCAAACCCTCGTAACGGACAACCTCGgccgc                                                       >  3:285901/2‑149 (MQ=255)                                                                                                cTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGAt                                                 >  7:236936/1‑148 (MQ=255)                                                                                                     gCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCactac                                           >  2:114810/1‑149 (MQ=255)                                                                                                     gCTTGCCTCCACTACTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCactac                                           <  1:114810/149‑1 (MQ=255)                                                                                                         gCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTGCGAGTACAACAGGGACGCCGGTGAACCCAGCGGCACGGACAACCACGGCCGACTTGAGAactactact                                       >  6:57843/1‑149 (MQ=255)                                                                                                           cTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGTCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGAt                                                 <  4:76534/137‑1 (MQ=255)                                                                                                           cTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGTCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGAt                                                 >  3:76534/1‑137 (MQ=255)                                                                                                              cACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCACTACTACTgggg                                   <  6:178376/148‑1 (MQ=255)                                                                                                                   cagcaGCCCCAAGCCGATCACGACTGGCTAGTGATTGGCCACTGGGCCAGCGATGGTCAAGAGTACCGGTTCGAGTACAACCTTGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCACTACTACTGGGGGCGAc                              <  8:270403/144‑1 (MQ=255)                                                                                                                    tgGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCACTACTACTGGGGGCCACTg                            >  7:170477/1‑149 (MQ=255)                                                                                                                     ggAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGTTGAAACCCTCGTCACGGACAACCTCGCCCGCCGTGATCACTACTACTGGGGGCCACTg                            >  4:172894/1‑148 (MQ=255)                                                                                                                     ggAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCACTACTACTGGGGGCCACTg                            >  1:8471/1‑148 (MQ=255)                                                                                                                      gAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCACTACTACTGGGGGCCACTGt                           <  8:236936/148‑1 (MQ=255)                                                                                                                        gCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCACTACTACTGGGGGCCACTGTACg                        >  8:230311/1‑149 (MQ=255)                                                                                                                          cccAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCACTACTACTGGGGGCCACTGTACGAg                      <  2:68824/149‑1 (MQ=255)                                                                                                                           caaaCCCGATAACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCACTACTACTGGGGGCCACTGTACGAgg                     >  1:83932/3‑149 (MQ=255)                                                                                                                            cAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCACTACTACTGGGGGCCACTGTACGAGGt                    >  3:60884/1‑149 (MQ=255)                                                                                                                                  gATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCACTACTACTGGGGGCCACTGTACGAgg                     >  8:290388/1‑142 (MQ=255)                                                                                                                                  gATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCACTACTACTGGGGGCCACTGTACGAgg                     <  7:290388/142‑1 (MQ=255)                                                                                                                                   aTCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCACTACTACTGGGGGCCACTGTAc                         <  4:176558/137‑1 (MQ=255)                                                                                                                                   aTCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCACTACTAATGGGGGCCACTGTAc                         >  3:176558/1‑137 (MQ=255)                                                                                                                                     cACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCACTACTACTGGGGGCCACTGTACGAGGTCTATAAGca           >  7:300268/1‑149 (MQ=255)                                                                                                                                       ccagtGGCGAATCACTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCACTACTACTGGGGGCCACTGTACGAg                      >  2:279075/5‑136 (MQ=255)                                                                                                                                       ccagtGGCGAATCACTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCACTACTACTGGGGGCCACTGTACGAg                      <  1:279075/132‑1 (MQ=255)                                                                                                                                       cGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCACTACTACTGGGGGCCACTGTACGAGGTCTATAAGcac          <  2:218987/148‑1 (MQ=255)                                                                                                                                       cGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCACTACTACTGGGGGCCACTGTACGAGGTCTATAAGcac          >  1:218987/1‑148 (MQ=255)                                                                                                                                        gAATGTCGAGTGATTCGGCACTCGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAATCCGGTGAATCCCTCGTCACGGAGATCCTCGGCCGCCGTGATCACTACTACTGGCGGCCACTGTACGAGGTGTATAAGCACAt        <  4:285901/149‑1 (MQ=255)                                                                                                                                        gAATGGCGAGTTATTCGGCAGTGTGTCAGCGTTGGTTAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCACTACTACTGGGGGCCACTGTACGAGGTCTATAAGCACAt        <  2:83932/149‑1 (MQ=255)                                                                                                                                              cGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCACTAATACTGGGGGCCACTGTACGATGTCTTTAAGCACATCGAgcc  >  6:249276/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                  CATCAACTGGCTAGCGGTGCTGTCCGCGAATACGAATGGGCGCGCTTCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCACTACTACTGGGGGCCACTGTACGAGGTCTATAAGCACATCGACCC  >  NC_002947/4615199‑4615488

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