Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,329,005	G→C	5.9%	G70A (GGC→GCC)	PP_3799 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,329,005	0	G	C	5.9%	209.3 / 2.9	68	G70A (GGC→GCC)	PP_3799	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (32/31);  new base C (2/2);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

CACAACCCGATGTGCGAATGCGAGGAACTGGCCAACGGCCAGGTCAAGTGCACCGGCGGTTTCTCGGACGGCAGCGGTGCCCCCGGGGTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTATGTGCTGTTCGATGCCGGCCCCGGGCACGTGGTTGAAGTGGACTACGCCGATATCGCCCAACCATGAGCCGCGCCCAAGTCATTCGCCCGGCCGGTGCCGG  >  NC_002947/4328863‑4329151                                                                                                                                               |                                                                                                                                                   cacaACCCGATGTGCGAATGCGAGGAACTGGCCAACGGCCAGGTCAAGTGCACCGGCGGTTTCTCGGACGGCAGCGGTGCCCCCGGGGTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCgacgac                                                                                                                                             <  2:668458/150‑1 (MQ=255) cacaACCCGATGTGCGAATGCGAGGAACTGGCCAACGGCCAGGTCAAGTGCACCGGCGGTTTCTCGGACGGCAGCGGTGCCCCCGGGGTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCgacgac                                                                                                                                             <  1:164365/150‑1 (MQ=255)  acaACCCGATGTGCGAATGCGAGGAACTGGCCAACGGCCAGGTCAAGTGCACCGGCGGTTTCTCGGACGGCAGCGGTGCCCCCGGGGTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTAGTCGTCGCCAAGCTTGCCgccgacc                                                                                                                                            <  2:1341008/150‑2 (MQ=255)  acaACCCGATGTGCGAATGCGAGGAACTGGCCAACGGCCAGGTCAAGTGCACCGGCGGTTTCTCGGACGGCAGCGGTGCCCCCGGGGTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTAGTCGTCGCCAAGCTTGCCgccgacc                                                                                                                                            >  1:1341008/1‑149 (MQ=255)   caACCCGATGTGCGAATGCAAGGACCTGGCCAACGGCCAGGTCAAGTGCACCGGCGGTTTCTCGGACGGCAGCGGTGCCCCCGGGGGGACCTTCGACGTGATCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTc                                                                                                                                           <  1:712140/150‑1 (MQ=255)    aaCCCGATGTGCGAATGCGAGGAACTGGCCAACGGCCAGGTCAAGTGCACCGGCGGTTTCTCGGACGGCAGCGGTGCCCCCGGGGTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCg                                                                                                                                          >  2:916147/1‑150 (MQ=255)    aaCCCGATGTGCGAATGCGAGGAACTGGCCAACGGCCAGGTCAAGTGCACCGGCGGTTTCTCGGACGGCAGCGGTGCCCCCGGGGTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCg                                                                                                                                          >  1:1413072/1‑150 (MQ=255)      cccGATGTGCGAATGCGAGGAACTGGCCAACGGCCAGGTCAAGTGCACCGGCGGTTTCTCGGACGGCAGCGGTGCCCCCGGGGTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGTCGCCAAGCTTGCCgccgaccag                                                                                                                                          >  2:58445/1‑145 (MQ=255)      cccGATGTGCGAATGCGAGGAACTGGCCAACGGCCAGGTCAAGTGCACCGGCGGTTTCTCGGACGGCAGCGGTGCCCCCGGGGTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGTCGCCAAGCTTGCCgccgaccag                                                                                                                                          <  1:58445/148‑4 (MQ=255)        cGATGTGCGAATGCGAGGAACTGGCCAACGGCCAGGTCAAGTGCACCGGCGGTTTCTCGGACGGCAGCGGTGCCCCCGGGGTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGc                                                                                                                                      >  1:585346/1‑150 (MQ=255)            gtgCGAATGCGAGGAACTGTCCAACGGCCAGTTGAAGTGCATCGGTGGTTTCGCGGACGTCAGCGGTGCCCCCGGTGTTACCCTCGACGTGATCGGTTATGACGAACAGGGGCTGGTCCGCGGCAAGCTTGGCGACGACTCGACGCTGAc                                                                                                                                  <  2:253064/150‑1 (MQ=255)            gtgCGAATGCGAGGAACTGGCCAACGGCCAGGTCAAGTGCACCGGCGGTTTCTCGGACGGCAGCGGTGCCCCCGGGGTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGAc                                                                                                                                  <  2:1337325/150‑1 (MQ=255)               cGAATGCGAGGAACTGGCCAACGGCCAGGTCAAGTGCACCGGCGGTTTCTCGGACGGCAGCGGTGCCCCCGGGGTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGtt                                                                                                                               >  1:1888205/1‑150 (MQ=255)                      gAGGAACTGGCCAACGGCCAGGTCAAGTGCACCGGCGGTTTCTCGGACGGCAGCGGTGCCCCCGGGGTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGa                                                                                                                                         >  2:1851477/1‑133 (MQ=255)                      gAGGAACTGGCCAACGGCCAGGTCAAGTGCACCGGCGGTTTCTCGGACGGCAGCGGTGCCCCCGGGGTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGa                                                                                                                                         <  1:1851477/133‑1 (MQ=255)                          aaCTGGCCAACGGCCAGGTCAAGTGCACCGGCGGTTTCTCGGACGGCAGCGGTGCCCCCGGGGTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCgg                                                                                                                    >  1:1017523/1‑150 (MQ=255)                               gCCAACGGCCAGGTCAAGTGCACCGGCGGTTTCTCGGACGGCAGCGGTGCCCCCGGGGTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGc                                                                                                               >  1:1415383/1‑150 (MQ=255)                                 cAACGGCCAGGTCAAGTGCACCGGCGGTTTCTCGGACGGCAGCGGTGCCCCCGGGGTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGc                                                                                                                                                   >  2:1153285/1‑112 (MQ=255)                                 cAACGGCCAGGTCAAGTGCACCGGCGGTTTCTCGGACGGCAGCGGTGCCCCCGGGGTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGc                                                                                                                                                   <  1:1153285/112‑1 (MQ=255)                                 cAACGGCCAGGTCAAGTGCACCGGCGGTTTCTCGGACGGCAGCGGTGCCCCCGGGGTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGa                                                                                                             >  1:1846442/1‑150 (MQ=255)                                  ccgctctaGGTCAAGGGCAACGGCGGTTTCTGGGACGGCAGCGGTGCCCCCGGGGTGACCGTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCTGGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGaa                                                                                                            <  2:726595/143‑1 (MQ=255)                                   aCGGCCAGGTCAAGTGCACCGGCGGTTTCTCGGAAGGCCGCGTTGCCCCCGGGGTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAAt                                                                                                           <  1:1557303/150‑1 (MQ=255)                                    cGGCCAGGTCAAGTGCACCGGCGGTTTCTCGGACGGCAGCGGTGCCCCCGGGGTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAAtt                                                                                                          <  2:1415383/150‑1 (MQ=255)                                       ccAGGTCAAGTGCACCGGCGGTTTCTCGGACGGCAGCGGTGCCCCCGGGGTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTa                                                                                                       >  2:312435/1‑150 (MQ=255)                                        cAGGTCAAGTGCACCGGCGGTTTCTCGGACGGCAGCGGTGCCCCCGGGGTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTAt                                                                                                      <  2:1888205/150‑1 (MQ=255)                                        cAGGTCAAGTGCACCCGCGGTTGCTCGGACGGCAGCGGTGCCCCCGGGGTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGGCCGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTAt                                                                                                      <  2:1255877/150‑1 (MQ=255)                                          ggTCAAGTGCACCGGCGGTTTCTCGGACGGCAGCGGTGCCCCCAGGGTGACCCTAAACGTGATCGGTTATTACTTTCAGGGGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAAtagatcgg                                                                                                    >  1:726595/1‑143 (MQ=255)                                              aaGTGCACCGGCGGTTTCTCGGACGGCAGCGGTGCCCCCGGGGTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTATGTGCTg                                                                                                >  2:377321/1‑150 (MQ=255)                                               aGTGCACCGGCGGTTTCTCGGACGGCAGCGGTGCCCCCGGGGTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGGCAGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTATGTGCTGt                                                                                               <  2:1696847/150‑1 (MQ=255)                                                      cggcggTTTCTCGGACGGCAGCGGTGCCCCCGGGGTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTATGTGCTGTTCGATGc                                                                                        <  1:916147/150‑1 (MQ=255)                                                        gcggTTTCTCGGACGGCAGCGGTGCCCCCGGGGTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTATGTGCTGTTCGATGCCg                                                                                      >  1:178618/1‑150 (MQ=255)                                                           gTTTCTCGGACGGCAGCGGTGCCCCCGGGGTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTATGTGCTGTTCGATGCCGGcc                                                                                   >  2:1340032/1‑150 (MQ=255)                                                                tcGGACGGCAGCGGTGCCCCCGGGGTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTATGTGCTGTTCGATGCCGGCCCCggg                                                                              >  1:880146/1‑150 (MQ=255)                                                                  ggACGGCAGCGGTGCCCCCGGGGTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGc                                                                                                               >  2:856512/1‑115 (MQ=255)                                                                  ggACGGCAGCGGTGCCCCCGGGGTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGc                                                                                                               <  1:856512/115‑1 (MQ=255)                                                                  ggACGGCAGCGGTGCCCCCGGGGTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTATGTGCTGTTCGATGCCGGCCCCGGGCa                                                                            >  1:1568684/1‑150 (MQ=255)                                                                   gACGGCAGCGGTGCCCCCGGGGTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTATGTGCTGTTCGATGCCGGCCCCGGGCAc                                                                           <  2:1846442/150‑1 (MQ=255)                                                                      ggCAGCGGTGCCCCCGGGGTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATGCTATGTGCTGTTCGATGCCGGCCCCGGGCACGTg                                                                        <  2:585346/150‑1 (MQ=255)                                                                        cAGCGGTGCCCCCGGGGTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTATGTGCTGTTCGATGCCGGCCCCGGGCACGTGGt                                                                      <  2:1017523/150‑1 (MQ=255)                                                                           cGGTGCCCCCGGGGTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTATGTGCTGTTCGATGCCGGCCCCGGACACGTGGTTGa                                                                   >  1:1161419/1‑150 (MQ=255)                                                                               gCCCCCGGGGTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTATGTGCTGTTCGATGCCGGCCCCGGGCACGTGGTTGAAGTg                                                               <  1:77207/150‑1 (MQ=255)                                                                               cccccgGGGGTGACCCTGGACGTGATCGGTTATGACGAACCGGTGCTGGTCGGCGGCAAGCTTGTCAACGACTCGACGCTGACGTTCAAGCGCCCGGACGCCGAATTCTATGTGCTGTTCGATGCCGGCCCCGGGCACGTGGTTGCAGTg                                                               <  2:708518/149‑1 (MQ=255)                                                                                       ggTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTATGTGCTGTTCGATGCCGGCCCCGGGCACGTGGTTGAAGTGGACTACGc                                                       <  2:1576487/150‑1 (MQ=255)                                                                                        gTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTATGTGCTGTTCGATGCCGGCCCCGGGCACGTGGTTGAAGTGGACTACGcc                                                      <  2:1568684/150‑1 (MQ=255)                                                                                         tGACCCTGGACGTGGTCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTATGTGCTGTTCGAt                                                                                          <  1:524418/113‑1 (MQ=255)                                                                                         tGACCCTGGACGTGGTCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTATGTGCTGTTCGAt                                                                                          >  2:524418/1‑113 (MQ=255)                                                                                            cccTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTATGTGCTGTTCGATGCCGGCCCCGGGCACGTGGTTGAAGTGGACTACGCCGata                                                  >  2:456691/1‑150 (MQ=255)                                                                                            cccTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTATGTGCTGTTCGATGCCGGCCCCGGGCACGTGGTTGAAGTGGACTACGCCGata                                                  >  2:1557912/1‑150 (MQ=255)                                                                                              cTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTATGTGCTGTTCGATGCCGGCCCCGGGCACGTGGTTGAAGTGGACTACGCCGATATc                                                >  1:252532/1‑150 (MQ=255)                                                                                              cTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTATGTGCGGTTCGATGCCGGCCCCGGGCACGTGGTTGAAGTGGACTACGCCGATATc                                                >  1:1078693/1‑150 (MQ=255)                                                                                                 gACGTGATCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTATGTGCTGTTCGATGCCGGCCCCGGGCACGTGGTTGAAGTGGACTACGCCGATATCGcc                                             >  2:1138435/1‑150 (MQ=255)                                                                                                    gTGATCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTATGTGCTGTTCGATGCCGGCCCCGGGCACGTGGTTGAAGTGGACTACGCCGATATCGCCCaa                                          >  1:1023150/1‑150 (MQ=255)                                                                                                            ttATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTATGTGCTGTTCGATGCCGGCCCCGGGCACGTGGTTGAAGTGGACTACGCCGATATCGCCCAACCATGAGc                                  >  2:259551/1‑150 (MQ=255)                                                                                                            ttATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTATGTGCTGTTCGATGCCGGCCCCGGGCACGTGGTTGAAGTGGACTACGCCGATATCGCCCAACCATGAGc                                  <  1:377321/150‑1 (MQ=255)                                                                                                             tATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTATGTGCTGTTCGATGCCGGCCCCGGGCACGTGGTTGAAGTGGACTACGCCGATATCGCCCAACCATGAGcc                                 <  1:602750/150‑1 (MQ=255)                                                                                                                 aCGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTATGTGCTGTTCGATGCCGGCCCCGGGCACGTGGTTg                                                                    >  1:989441/1‑111 (MQ=255)                                                                                                                 aCGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTATGTGCTGTTCGATGCCGGCCCCGGGCACGTGGTTg                                                                    <  2:989441/111‑1 (MQ=255)                                                                                                                 aCGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTATGTGCTGTTCGATGCCGGCCCCGGGCACGTGGTTGAAGTGGACTACGcc                                                      >  2:1729795/1‑125 (MQ=255)                                                                                                                 aCGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTATGTGCTGTTCGATGCCGGCCCCGGGCACGTGGTTGAAGTGGACTACGcc                                                      <  1:1729795/125‑1 (MQ=255)                                                                                                                     aCAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTATGTGCTGTTCGATGCCGGCCCCGGGCACGTGGTTGAAGTGGACTACGCCGATATCGCCCAACCATGAGCCGCGCCCaa                         <  2:252532/150‑1 (MQ=255)                                                                                                                          tGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTATGTGCTGTTCGATGCCGGCCCCGGGCACGTGGTTGAAGTGGACTACGCCGATATCGCCCAACCATGAGCCGCGCCCAAGTCAt                    >  1:1609193/1‑150 (MQ=255)                                                                                                                          tGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTATGTGCTGTTCGATGCCGGCCCCGGGCACGTGGTTGAAGTGGACTACGCCGATATCGCCCAACCATGAGCCGCGCCCAAGTCAt                    >  1:1348877/1‑150 (MQ=255)                                                                                                                             tGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTATGTGCTGTTCGATGCCGGCCCCGGGCACGTGGTTGAAGTGGACTACGCCGATATCGCCCAACCATGAGCCGCGCCCAAGTCATTCg                 >  2:580761/1‑150 (MQ=255)                                                                                                                                tCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTATGTGCTGTTCGATGCCGGCCCCGGGCACGTGGTTGAAGTGGACTACGCCGATATCGCCCAACCATGAGCCGCGCCCAAGTCATTCGccc              <  2:1348877/150‑1 (MQ=255)                                                                                                                                   gcggcAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTATGTGCTGTTCGATGCCGGCCCCGGGCACGTGGTTGAAGTGGACTACGCCGATATCGCCCAACCATGAGCCGCGCCCAAGTCATTCGCccggc           <  1:312435/150‑1 (MQ=255)                                                                                                                                      gcAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTATGTGCTGTTCGATGCCGGCCCCGGGCACGTGGTTGAAGTGGACTACGCCGATATCGCCCAACCATGAGCCGCGCCCAAGTCATTCGCccggccgg        >  2:721993/1‑150 (MQ=255)                                                                                                                                           cTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTATGTGCTGTTCGATGCCGGCCCCGGGCACGTGGTTGAAGTGGACTACGCCGATATCGCCCAACCATGAGCCGCGCCCAAGTCATTCGCCCGGCCGGTGCCg   <  2:1609193/150‑1 (MQ=255)                                                                                                                                            ttGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTATGTGCTGTTCGATGCCGGCCCCGGGCACGTGGTTGAAGTGGACTACGCCGATATCGACCAAACCTGAGCCGCGCCCAAGTCATTCGCCCGGCCGGTGCCgg  <  1:1340032/150‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                   CACAACCCGATGTGCGAATGCGAGGAACTGGCCAACGGCCAGGTCAAGTGCACCGGCGGTTTCTCGGACGGCAGCGGTGCCCCCGGGGTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTATGTGCTGTTCGATGCCGGCCCCGGGCACGTGGTTGAAGTGGACTACGCCGATATCGCCCAACCATGAGCCGCGCCCAAGTCATTCGCCCGGCCGGTGCCGG  >  NC_002947/4328863‑4329151

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