Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	473,702	T→C	8.5%	Q428R (CAG→CGG)	dnaG ←	DNA primase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	473,702	0	T	C	8.5%	163.9 / 2.9	58	Q428R (CAG→CGG)	dnaG	DNA primase
Reads supporting (aligned to +/- strand):  ref base T (25/28);  new base C (2/3);  total (27/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTCCTGATGAACGGGGGCATAGTCGTGCTGCGGCATGTCGCCGTAGTCGGGGGTGTAGCTGGCCATGGCGTCGTAATCGTAGCCAGGGTCGTAGTCCGGCATGCTGCTGGTCGCCGGCGCTTGTTGCGCCAGCTGCTCGACCTGCTGCGGGTCGAGGCC‑GGTGATTTCCTTCAGGCGGTTGCGCATCAGCTGACGCAGGTTGGCGCCGGGGATCTTTTCGATCAACGGTGCGGCCAAGGTCGCCATATGCGCCTTGCCTTCCAGCGAGCGCGGGTCGGCTTCTACGCCAA  >  NC_002947/473559‑473849                                                                                                                                                |                                                                                                                                                     gCTCCTGATGAACGGGGGCATAGTCGTGCTGCGGCATGTCGCCGTAGGCGGGGGGGTAGCTGGCCATGGCGTCGTAATCGTAGCCAGGGTCGTAGTCCGGCATGCTGCTGGTCGCCGGCGCTTGTTGCGCCAGCTGCTCGACCtgctgc                                                                                                                                                 <  5:221214/149‑1 (MQ=255)  cTCCTGATGAACGGGGGCATAGTCGTGCTGCGGCATGTCGCCGTAGTCGGGGGTGTAGCTGGCCATGGCGTCGTAATCGTAGCCAGGGTCGTAGTCCGGCATGCTGCTGGTCGCCGGCGCTTGTTGCGCCAGCTGCTCGACCtgctgc                                                                                                                                                 <  4:296095/148‑1 (MQ=255)   tCCTGATGAACGGGGGCATAGTCGTGCTGCGGCATGTCGCCGTAGTCGGGGGTGTAGCTGGCCATGGCGTCGTAATCGTAGCCAGGGTCGTAGTCCGGCATGCTGCTGGTCGCCGGCGCTTGTTGCGCCAGCTGCTCGACCTGCTGCgg                                                                                                                                               >  2:22495/1‑149 (MQ=255)      tgatgaACGGGGGCATAGTCGTGCTGCGGCATGTCGCCGTAGTCGGGGGTGTAGCTGGCCATGGCGTCGTAATCGTAGCCAGGGTCGTAGTCCGGCATGCTGCTGGTCGCCGGCGCTTGTTGCGCCAGCTGCTCGACCTGCTGCGGGTc                                                                                                                                            >  3:192624/1‑149 (MQ=255)       gatgaACGGGGGCATAGTCGTGCTGCGGCATGTCGCCGTAGTCGGGGGTGTAGCTGGCCATGGCGTCGTAATCGTAGCCAGGGTCGTAGTCCGGCATGCTGCTGGTCGCCGGCGCTTGTTGCGCCAGCTGCTCGACCtgctgc                                                                                                                                                 <  1:187911/143‑1 (MQ=255)       gatgaACGGGGGCATAGTCGTGCTGCGGCATGTCGCCGTAGTCGGGGGTGTAGCTGGCCATGGCGTCGTAATCGTAGCCAGGGTCGTAGTCCGGCATGCTGCTGGTCGCCGGCGCTTGTTGCGCCAGCTGCTCGACCtgctgc                                                                                                                                                 >  2:187911/1‑143 (MQ=255)        atgaACGGGGGCATAGTCGTGCTGCGGCATGTCGCCGTAGTCGGGGGTGTAGCTGGCCATGGCGTCGTAATCGTAGCCAGGGTCGTAGTCCGGCATGCTGCTGGTCGCCGGCGCTTGTTGCGCCAGCTGCTCGACCCGCAGCAGGTCGa                                                                                                                                          <  7:176575/149‑1 (MQ=255)          gaACGGGGGCATAGTCGTGCTGCGGCATGTCGCCGTAGTCGGGGGTGTAGCTGGCCATGGCGTCGTAATCGTAGCCAGGGTCGTAGTCCGGCATGCTGCTGGTCGCCGGCGCTTGTTGCGCCAGCTGCTCGACCTGCTGCGGGTCGAgg                                                                                                                                        >  4:42957/1‑149 (MQ=255)            aCGGGGGCATAGTCGTGCTGCGGCATGTCGCCGTAGTCGGGGGTGTAGCTGGCCATGGCGTCGTAATCGTAGCCAGGGTCGTAGTCCGGCATGCTGCTGGTCGCCGGCGCTTGTTGCGCCAGCTGCTCGACCTGCTGCGGGTCGAGGc                                                                                                                                       <  4:172925/148‑1 (MQ=255)             cGGGGGCATAGTCGTGCTGCGGCATGTCGCCGTAGTCGGGGGTGTAGCTGGCCATGGCGTCGTAATCGTAGCCAGGGTCGTAGTCCGGCATGCTGCTGGTCGCCGGCGCTTGTTGCGCCAGCTGCTCGACCTGCTGCGGGTCGAGGCC‑g                                                                                                                                    <  5:11439/149‑1 (MQ=255)              gggggCATAGTCGTGCTGCGGCATGTCGCCGTAGTCGGGGGTGTAGCTGGCCATGGCGTCGTAATCGTAGCCAGGGTCGTAGTCCGGCATGCTGCTGGTCGCCGGCGCTTGTTGCGCCAGCTGCTCGACCTGCTGCGGGTCGAGGCC‑gg                                                                                                                                   >  2:131486/1‑149 (MQ=255)                 ggCATAGTCGTGCTGCGGCATGTCGCCGTAGTCGGGGGTGTAGCTGGCCATGGCGTCGTAATCGTAGCCAGGGTCGTAGTCCGGCATGCTGCTGGTCGCCGGCGCTTGTTGCGCCAGCTGCTCGACCTGCTGCGGGTCGAGGCC‑GGTGa                                                                                                                                <  2:72608/149‑1 (MQ=255)                   cATAGTCGTGCTGCGGCATGTCGCCGTAGTCGGGGGTGTAGCGGGCCATGGCGTCGTAATCGTAGCCAGGGTCGTAGTCCGGCATGCTGCTGGTCGCCGGCGCTTGTTGCGCCAGCTGCTCGACCTGCTGCGGGTCGAGGCC‑GGTGAtt                                                                                                                              <  8:100747/149‑1 (MQ=255)                             ctgcGGCATGGCGCCGTTGTCGGGGGTGTAGCTGGCCATGGCGTCGTAATCGTAGCCAGGGTCGTAGTCCGGCATGCTGCTGGTCGCCGGCGCTTGTTGCGCCAGCTGCTCGACCTGCTGCGGGTCGAGGCC‑GGTGATTTCCTTCAGGc                                                                                                                    <  7:228813/149‑1 (MQ=255)                                  gCATGTCGCCGTAGTCGGGGGTGTAGCTGGCCATGGCGTCGTAATCGTAGCCAGGGTCGTAGTCCGGCATGCTGCTGGTCGCCGGCGCTTGTTGCGCCAGCTGCTCGACCTGCTGCGGGTCGAGGCC‑GGTGATTTCCTTCAGGCGGTTg                                                                                                               >  5:124373/1‑149 (MQ=255)                                   cATGTCGCCGTAGTCGGGGGTGTAGCTGGCCATGGCGTCGTAATCGTAGCCAGGGTCGTAGTCCGGCATGCTGCTGGTCGCCGGCGCTTGTTGCGCCAGCTGCTCGACCTGCTGCGGGTCGAGGCC‑GGTGATTTCCTTCAGGCGGTTgc                                                                                                              <  5:121834/149‑1 (MQ=255)                                   cATGTCGCCGTAGTCGGGGGTGTAGCTGGCCATGGCGTCGTAATCGTAGCCAGGGTCGTAGTCCGGCATGCTGCTGGTCGCCGGCGCTTGTTGCGCCAGCTGCTCGACCTGCTGCGGGTCGAGGCC‑GGTGATTTCCTTCAGGCGGTTgc                                                                                                              <  6:210323/149‑1 (MQ=255)                                   cATGTCGCCGTAGTCGGGGGTGTAGCTGGCCATGGCGTCGTAATCGTAGCCAGGGTCGTAGTCCGGCATGCTGCTGGTCGCCGGCGCTTGTTGCGCCAGCTGCTCGACCTGCTGCGGGTCGAGGCC‑GGTGATTTCCTTCAGGCGGTTgc                                                                                                              >  6:121834/1‑149 (MQ=255)                                    aTGTCGCCGTAGTCGGGGGTGTAGCTGGCCATGGCGTCGTAATCGTAGCCAGGGTCGTAGTCCGGCATGCTGCTGGTCGCCGGCGCTTGTTGCGCCAGCTGCTCGACCTGCTGCGGGTCGAGGCC‑GGTGATTTCCTTCAGGCGGTTgcg                                                                                                             <  4:192624/149‑1 (MQ=255)                                         gCCGTAGTCGGGGGTGTAGCTGGCCATGGCGTCGTAATCGTAGCCAGGGTCGTAGTCCGGCATGCTGCTGGTCGCCGGCGCTTGTTGCGCCAGCTGCTCGACCTGCTGCGGGTCGAgg                                                                                                                                        >  3:1528/1‑118 (MQ=255)                                         gCCGTAGTCGGGGGTGTAGCTGGCCATGGCGTCGTAATCGTAGCCAGGGTCGTAGTCCGGCATGCTGCTGGTCGCCGGCGCTTGTTGCGCCAGCTGCTCGACCTGCTGCGGGTCGAgg                                                                                                                                        <  4:1528/118‑1 (MQ=255)                                         gCCGTAGTCGGGGGTGTAGCTGGCCATGGCGTCGTAATCGTAGCCAGGGTCGTAGTCCGGCATGCTGCTGGTCGCCGGCGCTTGTTGCGCCAGCTGCTCGACCTGCTGCGGGTCGAGGCC‑GGTGATTTCCTTCAGGCGGTTGCGCAt                                                                                                          >  1:40257/1‑147 (MQ=255)                                         gCCGTAGTCGGGGGTGTAGCTGGCCATGGCGTCGTAATCGTAGCCAGGGTCGTAGTCCGGCATGCTGCTGGTCGCCGGCGCTTGTTGCGCCAGCTGCTCGACCTGCTGCGGGTCGAGGCC‑GGTGATTTCCTTCAGGCGGTTGCGCATc                                                                                                         <  6:124373/148‑1 (MQ=255)                                           cGTAGTCGGGGGTGTAGCTGGCCATGGCGTCGTAATCGTAGCCAGGGTCGTAGTCCGGCATGCTGCTGGTCGCCGGCGCTTGTTGCGCCAGCTGCTCGACCTGCTGCGGGTCGAGGCC‑GGTGATTTCCTTCAGGCGGTTGCGCATCAGc                                                                                                      <  3:241173/149‑1 (MQ=255)                                           cGTAGTCGGGGGTGTAGCTGGCCATGGCGTCGTAATCGTAGCCAGGGTCGTAGTCCGGCATGCTGCTGGTCGCCGGCGCTTGTTGCGCCAGCTGCTCGACCTGCTGCGGGTCGAGGCC‑GGTGATTTCCTTCAGGCGGTTGCGCATCAGc                                                                                                      >  5:79554/1‑149 (MQ=255)                                                  gggggTGTAGCTGGCCATGGCGTCGTAATCGTAGCCAGGGTCGTAGTCCGGCATGCTGCTGGTCGCCGGCGCTTGTTGCGCCAGCTGCTCGACCTGCTGCGGGTCGAGGCC‑GGTGATTTCCTTCAGGCGGTTGCGCATCAGCTGACg                                                                                                 >  5:336163/1‑147 (MQ=255)                                                  gggggTGTAGCTGGCCATGGCGTCGTAATCGTAGCCAGGGTCGTAGTCCGGCATGCTGCTGGTCGCCGGCGCTTGTTGCGCCAGCTGCTCGACCTGCTGCGGGTCGAGGCC‑GGTGATTTCCTTCAGGCGGTTGCGCATCAGCTGACGc                                                                                                <  1:131486/148‑1 (MQ=255)                                                    gggTGTAGCTGGCCATGGCGTCGTAATCGTAGCCAGGGTCGTAGTCCGGCATGCTGCTGGTCGCCGGCGCTTGTTGCGCCAGCTGCTCGACCTGCTGCGGGTCGAGGCC‑GGTGATTTCCTTCAGGCGGTTGCGCATCAGCTGACGCAgg                                                                                             <  6:79554/149‑1 (MQ=255)                                                         tAGCTGGCCATGGCGTCGTAATCGTAGCCAGGGTCGTAGTCCGGCATGCTGCTGGTCGCCGGCGCTTGTTGCGCCAGCTGCTCGACCTGCTGCGGGTCGAGGCC‑GGTGATTTCCTTCAGGCGGTTGCGCATCAGCTGACGCAGGTTGgc                                                                                        <  2:182299/149‑1 (MQ=255)                                                             tGGCCATGGCGTCGTAATCGTAGCCAGGGTCGTAGTCCGGCATGCTGCTGGTCGCCGGCGCTTGTTGCGCCAGCTGCTCGACCTGCTGCGGGTCGAGGCC‑GGTGATTTCCTTCAGGCGGTTGCGCATCAGCTGACGCAGGTTGgcg                                                                                       >  8:135719/1‑146 (MQ=255)                                                             tGGCCATGGCGTCGTAATCGTAGCCAGGGTCGTAGTCCGGCATGCTGCTGGTCGCCGGCGCTTGTTGCGCCAGCTGCTCGACCTGCTGCGGGTCGAGGCC‑GGTGATTTCCTTCAGGCGGTTGCGCATCAGCTGACGCAGGTTGgcg                                                                                       <  7:135719/146‑1 (MQ=255)                                                                 cATGGCGTCGTAATCGTAGCCAGGGTCGTAGTCCGGCATGCTGCTGGTCGCCGGCGCTTGTTGCGCCAGCTGCTCGACCTGCTGCGGGTCGAGGCC‑GGTGATTTCCTTCAGGCGGTTGCGCATCAGCTGACGCAGGTTGGCGCCGGGGa                                                                                >  5:355950/1‑149 (MQ=255)                                                                 cATGGCGTCGTAATCGTAGCCAGGGTCGTAGTCCGGCATGCTGCTGGTCGCCGGCGCTTGTTGCGCCAGCTGCTCGACCTGCTGCGGGTCGAGGCC‑GGTGATTTCCTTCAGGCGGTTGCGCATCAGCTGACGCAGGTTGGCGCCGGGGa                                                                                <  3:259216/149‑1 (MQ=255)                                                                   tGGCGTCGTAATCGTAGCCAGGGTCGTAGTCCGGCATGCTGCTGGTCGCCGGCGCTTGTTGCGCCAGCTGCTCGACCTGCTGCGGGTCGAGGCC‑GGTGATTTCCTTCAGGCGGTTGCGCATCAGCTGACGCAGGTTGGCGCCGGGGATc                                                                              <  6:355950/149‑1 (MQ=255)                                                                       gtcgtAATCGTAGCCAGGGTCGTAGTCCGGCATGCTGCTGGTCGCCGACGCTTGTTGCGCCAGCTGCTCGACCTGCTGCGGGTCGAGGCC‑GGTGATTTCCTTCAGGCGGTTGCGCATCAGCTGACGCAGGTTGGCGCCGGGGATCtttt                                                                          >  1:79307/1‑149 (MQ=255)                                                                        tcgtAATCGTAGCCAGGGTCGTAGTCCGGCATGCTGCTGGTCGCCGGCGCTTGTTGCGCCAGCTGCTCGACCTGCTGCGGGTCGAGGCC‑GGTGATTTCCTTCAGGCGGTTGCGCATCAGCTGACGCAGGTTGGCGCCGGGGATCTTTTc                                                                         >  8:125955/1‑149 (MQ=255)                                                                          gtAATCGTAGCCAGGGTCGTAGTCCGGCATGCTGCTGGTCGCCGGCGCTTGTTGCGCCAGCTGCTCGACCTGCTGCGGGTCGAGGCC‑GGTGATTTCCTTCAGGCGGTTGCGCATCAGCTGACGCAGGTTGGCGCCGGGGATCTTTTCGa                                                                       <  3:42957/149‑1 (MQ=255)                                                                            aaTCGTAGCCAGGGTCGTAGTCCGGCATGCTGCTGGTCGCCGGCGCTTGTTGCGCCAGCTGCTCGACCTGCTGCGGGTCGAGGCC‑GGTGATTTCCTTCAGGCGGTTgcg                                                                                                             <  1:90608/109‑1 (MQ=255)                                                                            aaTCGTAGCCAGGGTCGTAGTCCGGCATGCTGCTGGTCGCCGGCGCTTGTTGCGCCAGCTGCTCGACCTGCTGCGGGTCGAGGCC‑GGTGATTTCCTTCAGGCGGTTgcg                                                                                                             >  2:90608/1‑109 (MQ=255)                                                                                    ccAGGGTCGTAGTCCGGCATGCTGCTGGTCGCCGGCGCTTGTTGCGCCAGCTGCTCGACCTGCTGCGGGTCGAGGCC‑GGTGATTTCCTTCAGGCGGTTGCGCATCAGCTGACGCAGGTTGGCGCCGGGGATCTTTTCGATCAACGGTGc                                                             >  4:335324/1‑149 (MQ=255)                                                                                          tCGTAGTCCGGCATGCTGCTGGTCGCCGGCGCTTGTTGCGCCAGCTGCTCGACCTGCTGCGGGTCGAGGCC‑GGTGATTTCCTTCAGGCGGTTGCGCATCAGCTGACGCAGGTTGGCGCCGGGGATCTTTTCGATCAACGGTGCGGcc                                                         >  6:187744/1‑147 (MQ=255)                                                                                          tCGTAGTCCGGCATGCTGCTGGTCGCCGGCGCTTGTTGCGCCAGCTGCTCGACCTGCTGCGGGTCGAGGCC‑GGTGATTTCCTTCAGGCGGTTGCGCATCAGCTGACGCAGGTTGGCGCCGGGGATCTTTTCGATCAACGGTGCGGcc                                                         >  7:280568/1‑147 (MQ=255)                                                                                            gTAGTCCGGCATGCTGCTGGTCGCCGGCGCTTGTTGCGCCAGCTGCTCGACCTGCTGCGGGTCGAGGCC‑GGTGATTTCCTTCAGGCGGTTGCGCATCAGCTGACGCAGGTTGGCGCCGGGTATCTTTTCGATCAACGGTGCGGCCAAgg                                                     <  1:22495/149‑1 (MQ=255)                                                                                               gTCCGGCATGCTGCTGGTCGCCGGCGCTTGTTGCGCCAGCTGCTCGACCTGCTGCGGGTCGAGGCC‑GGTGATTTCCTTCAGGCGGTTGCGCATCAGCTGACGCAGGTTGGCGCCGGGGATCTTTTCGATCAACGGTGCGGCCAAGGTc                                                   <  2:40257/148‑1 (MQ=255)                                                                                                          tgctgGTCGCCGGCGCTTGTTGCGCCAGCTGCTCGACCTGCTGCGGGTCGAGGCC‑GGTGATTTCCTTCAGGCGGTTGCGCATCAGCTGACGCAGGTTGGCGCCGGGGATCTTTTCGATCAACGGTGCGGCCAAGGTCGCCATATGCgcc                                       >  5:363954/1‑149 (MQ=255)                                                                                                            ctgGTCGCCGGCGCTTGTTGCGCCAGCTGCTCGACCTGCTGCGGGTCGAGGCC‑GGTGATTTCCTTCAGGCGGTTGCGCATCAGCTGACGCAGGTTGGCGCCGGGGATCTTTTCGATCAACGGTGCGGCCAAGGTCGCCATATGCgcctt                                     >  5:216024/1‑149 (MQ=255)                                                                                                            ctgGTCGCCGGCGCTTGTTGCGCCAGCTGCTCGACCTGCTGCGGGTCGAGGCC‑GGTGATTTCCTTCAGGCGGTTGCGCATCAGCTGACGCAGGTTGGCGCCGGGGATCTTTTCGATCAACGGTGCGGCCAAGGTCGCCATATGCgcct                                      <  7:125955/148‑1 (MQ=255)                                                                                                             tgGTCGCCGGCGCTTGTTGCGCCAGCTGCTCGACCTGCTGCGGGTCGAGGCC‑GGTGATTTCCTTCAGGCGGTTGCGCATCAGCTGACGCAGGTTGGCGCCGGGGATCTTTTCGATCAACGGTGCGGCCAAGGTCGCCATATGCgccttg                                    <  2:79307/149‑1 (MQ=255)                                                                                                                tCGCCGGCGCTTGTTGCGCCAGCTGCTCGACCTGCTGCGGGTCGAGGCCGGGTGATTTCCTTCAGGCGGTTGCGCATCAGCTGACGCAGGTTGGCGCCGGGGATCTTTTCGATCAACGGTGCGGCCAAGGTCGCCATATGCgccttgcc                                  >  1:272004/1‑149 (MQ=255)                                                                                                                  gCCGGCGCTTGTTGCGCCAGCTGCTCGACCTGCTGCGGGTCGAGGCC‑GGTGATTTCCTTCAGGCGGTTGCGCATCAGCTGACGCAGGTTGGCGCCGGGGATCTTTTCGATCAACGGTGCGGCCAAGGTCGCCATATGCGCCTTGCCTTc                               <  6:216024/149‑1 (MQ=255)                                                                                                                  gCCGGCGCTTGTTGCGCCAGCTGCTCGACCTGCTGCGGGTCGAGGCC‑GGTGATTTCCTTCAGGCGGTTGCGCATCAGCTGACGCAGGTTGGCGCCGGGGATCTTTTCGATCAACGGTGCGGCCAAGGTCGCCATATGCGCCTTGCCTTc                               >  4:250133/1‑149 (MQ=255)                                                                                                                        gcTTGTTGCGCCAGCTGCTCGACCTGCTGCGGGTCGAGGCC‑GGTGATTTCCTTCAGGCGGTTGCGCATCAGCTGACGCAGGTTGGCGCCGGGGATCTTTTCGATCAACGGTGCGGCCAAGGTCGCCATATGCGCCTTGCCTTc                               <  4:327046/143‑1 (MQ=255)                                                                                                                        gcTTGTTGCGCCAGCTGCTCGACCTGCTGCGGGTCGAGGCC‑GGTGATTTCCTTCAGGCGGTTGCGCATCAGCTGACGCAGGTTGGCGCCGGGGATCTTTTCGATCAACGGTGCGGCCAAGGTCGCCATATGCGCCTTGCCTTc                               >  3:327046/1‑143 (MQ=255)                                                                                                                             ttgCGCCAGCTGCTCGACCTGCTGCGGGTCGAGGCC‑GGTGATTTCCTTCAGGCGGTTGCGCATCAGCTGACGCAGGTTGGCGCCGGGGATCTTTTCGATCAACGGTGCGGCCAAGGTCGCCATATGCGCCTTGCCTTCCAGCGAgcgc                     <  6:336163/148‑1 (MQ=255)                                                                                                                                 gcCAGCTGCTCGACCTGCTGCGGGTCGAGGCC‑GGTGATTTCCTTCAGGCGGTTGCGCATCAGCTGACGCAGGTTGGCGCCGGGGATCTTTTCGATCAACGGTGCGGCCAAGGTCGCCATATGCGCCTTGCCTTCCAGCGAGCGCGGGTc                <  6:363954/149‑1 (MQ=255)                                                                                                                                    ccggcctCGACCCGCAGCAGGTCGA‑GCC‑GGTGATTTCCTTCAGGCGGTTGCGCATCAGCTGACGCAGGTTGGCGCCGGGGATCTTTTCGATCAACGGTGCGGCCAAGGTCGCCATATGCGCCTTGCCTTCCAGCGAgcgc                     >  6:227904/6‑140 (MQ=255)                                                                                                                                    ccggcctCGACCCGCAGCAGGTCGA‑GCC‑GGTGATTTCCTTCAGGCGGTTGCGCATCAGCTGACGCAGGTTGGCGCCGGGGATCTTTTCGATCAACGGTGCGGCCAAGGTCGCCATATGCGCCTTGCCTTCCAGCGAgcgc                     <  5:227904/135‑1 (MQ=255)                                                                                                                                        cctCGACCCGCAGCAGGTCGAGGCC‑GGTGATTTCCTTCAGGCGGTTGCGCATCAGCTGACGCAGGTTGGCGCCGGGGATCTTTTCGATCAACGGTGCGGCCAAGGTCGCCATATGCGCCTTGCCTTCCAGCGAGCGCGGGTCg               <  4:277299/142‑1 (MQ=255)                                                                                                                                        cctCGACCCGCAGCAGGTCGAGGCC‑GGTGATTTCCTTCAGGCGGTTGCGCATCAGCTGACGCAGGTTGGCGCCGGGGATCTTTTCGATCAACGGTGCGGCCAAGGTCGCCATATGCGCCTTGCCTTCCAGCGAGCGCGGGTCg               >  3:277299/2‑143 (MQ=255)                                                                                                                                          tCGACCTGCTGCGGGTCGAGGCC‑GGTGATTTCCTTCAGGCGGTTGCGCATCAGCTGACGCAGGTTGGCGCCGGGGATCTTTTCGATCAACGGTGCGGCCAAGGTCGCCATATGCGCCTTGCCTTCCAGCGAGCGCGGGTCGGCTTCTAc       >  6:16720/1‑149 (MQ=255)                                                                                                                                           cGACCTGCTGCGGGTCGAGGCC‑GGTGATTTCCTTCAGGCGGTTGCGCATCAGCTGACGCAGGTTGGCGCCGGGGATCTTTTCGATCAACGGTGCGGCCAAGGTCGCCATATGCGCCTTGCCTTCCAGCGAGCGCGGGTCGGCTTCTACg      >  8:33545/1‑149 (MQ=255)                                                                                                                                            gACCTGCTGCGGGTCGAGGGC‑GGTGATTTCCTTCAGGCGGTTGCGCATCAGATGACGCAGGTTGGCGCCGGTGATCTTTTCGATCAACGGTGCTGCCAAGGTCGCCATATGCGCCTTGCCTTCCAGCGAGCGCGGGTCGGCTTCTAc       <  6:95905/147‑1 (MQ=255)                                                                                                                                               ctgctgCGGGTCGAGGCC‑GGTGATTTCCTTCAGGCGGTTGCGCATCAGCTGACGCAGGTTGGCGCCGGGGATCTTTTCGATCAACGGTGCGGCCAAGGTCGCCATATGCGCCTTGCCTTCCAGCGAGCGCGGGTCGGCTTCTACGCCaa  <  5:127935/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                     GCTCCTGATGAACGGGGGCATAGTCGTGCTGCGGCATGTCGCCGTAGTCGGGGGTGTAGCTGGCCATGGCGTCGTAATCGTAGCCAGGGTCGTAGTCCGGCATGCTGCTGGTCGCCGGCGCTTGTTGCGCCAGCTGCTCGACCTGCTGCGGGTCGAGGCC‑GGTGATTTCCTTCAGGCGGTTGCGCATCAGCTGACGCAGGTTGGCGCCGGGGATCTTTTCGATCAACGGTGCGGCCAAGGTCGCCATATGCGCCTTGCCTTCCAGCGAGCGCGGGTCGGCTTCTACGCCAA  >  NC_002947/473559‑473849

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