Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,615,347	T→C	7.1%	I573T (ATT→ACT)	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,347	0	T	C	7.1%	197.2 / 4.4	71	I573T (ATT→ACT)	PP_4085	rhs family protein
Reads supporting (aligned to +/- strand):  ref base T (35/30);  new base C (2/3);  total (38/33)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

CATCAACTGGCTAGCGGTGCTGTCCGCGAATACGAATGGGCGCGCTTCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCACTACTACTGGGGGCCACTGTACGAGGTCTATAAGCACATCGACCCTCA  >  NC_002947/4615199‑4615491                                                                                                                                                     |                                                                                                                                                 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                                                                                                                                         >  3:158966/1‑110 (MQ=255)                                              ttCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAATCACTCGCCattc                                                                                                                                         <  4:158966/113‑4 (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                                                                                                 <  5:24506/142‑1 (MQ=255)                                                         tGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTg                                                                                                 >  6:24506/1‑142 (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)                                                                          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                        <  7:290388/142‑1 (MQ=255)                                                                                                                                  gATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCACTACTACTGGGGGCCACTGTACGAgg                        >  8:290388/1‑142 (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)                                                                                                                                                  tGATTCGCCACTGGGGCATCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCACTACTACTGGGGGCCACTGTACGAGGTCTATAAGCACATCGACCCt    >  5:40686/1‑146 (MQ=255)                                                                                                                                                  tGATTCGCCACTGGGGCATCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCACTACTACTGGGGGCCACTGTACGAGGTCTATAAGCACATCGACCCTCa  <  6:40686/148‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                 CATCAACTGGCTAGCGGTGCTGTCCGCGAATACGAATGGGCGCGCTTCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAGCCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCACTACTACTGGGGGCCACTGTACGAGGTCTATAAGCACATCGACCCTCA  >  NC_002947/4615199‑4615491

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