Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	473,708	G→A	8.2%	P426L (CCG→CTG)	dnaG ←	DNA primase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	473,708	0	G	A	8.2%	172.0 / 5.0	61	P426L (CCG→CTG)	dnaG	DNA primase
Reads supporting (aligned to +/- strand):  ref base G (26/30);  new base A (2/3);  total (28/33)
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.54e-01

TCCTGATGAACGGGGGCATAGTCGTGCTGCGGCATGTCGCCGTAGTCGGGGGTGTAGCTGGCCATGGCGTCGTAATCGTAGCCAGGGTCGTAGTCCGGCATGCTGCTGGTCGCCGGCGCTTGTTGCGCCAGCTGCTCGACCTGCTGCGGGTCGAGGCC‑GGTGATTTCCTTCAGGCGGTTGCGCATCAGCTGACGCAGGTTGGCGCCGGGGATCTTTTCGATCAACGGTGCGGCCAAGGTCGCCATATGCGCCTTGCCTTCCAGCGAGCGCGGGTCGGCTTCTACGCCAAGTTGCT  >  NC_002947/473561‑473855                                                                                                                                                    |                                                                                                                                                     tCCTGATGAACGGGGGCATAGTCGTGCTGCGGCATGTCGCCGTAGTCGGGGGTGTAGCTGGCCATGGCGTCGTAATCGTAGCCAGGGTCGTAGTCCGGCATGCTGCTGGTCGCCGGCGCTTGTTGCGCCAGCTGCTCGACCTGCTGCgg                                                                                                                                                     >  2:22495/1‑149 (MQ=255)    tgatgaACGGGGGCATAGTCGTGCTGCGGCATGTCGCCGTAGTCGGGGGTGTAGCTGGCCATGGCGTCGTAATCGTAGCCAGGGTCGTAGTCCGGCATGCTGCTGGTCGCCGGCGCTTGTTGCGCCAGCTGCTCGACCTGCTGCGGGTc                                                                                                                                                  >  3:192624/1‑149 (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                                                                                                                                              <  4:1528/118‑1 (MQ=255)                                       gCCGTAGTCGGGGGTGTAGCTGGCCATGGCGTCGTAATCGTAGCCAGGGTCGTAGTCCGGCATGCTGCTGGTCGCCGGCGCTTGTTGCGCCAGCTGCTCGACCTGCTGCGGGTCGAgg                                                                                                                                              >  3:1528/1‑118 (MQ=255)                                       gCCGTAGTCGGGGGTGTAGCTGGCCATGGCGTCGTAATCGTAGCCAGGGTCGTAGTCCGGCATGCTGCTGGTCGCCGGCGCTTGTTGCGCCAGCTGCTCGACCTGCTGCGGGTCGAGGCC‑GGTGATTTCCTTCAGGCGGTTGCGCAt                                                                                                                >  1:40257/1‑147 (MQ=255)                                       gCCGTAGTCGGGGGTGTAGCTGGCCATGGCGTCGTAATCGTAGCCAGGGTCGTAGTCCGGCATGCTGCTGGTCGCCGGCGCTTGTTGCGCCAGCTGCTCGACCTGCTGCGGGTCGAGGCC‑GGTGATTTCCTTCAGGCGGTTGCGCATc                                                                                                               <  6:124373/148‑1 (MQ=255)                                         cGTAGTCGGGGGTGTAGCTGGCCATGGCGTCGTAATCGTAGCCAGGGTCGTAGTCCGGCATGCTGCTGGTCGCCGGCGCTTGTTGCGCCAGCTGCTCGACCTGCTGCGGGTCGAGGCC‑GGTGATTTCCTTCAGGCGGTTGCGCATCAGc                                                                                                            >  5:79554/1‑149 (MQ=255)                                         cGTAGTCGGGGGTGTAGCTGGCCATGGCGTCGTAATCGTAGCCAGGGTCGTAGTCCGGCATGCTGCTGGTCGCCGGCGCTTGTTGCGCCAGCTGCTCGACCTGCTGCGGGTCGAGGCC‑GGTGATTTCCTTCAGGCGGTTGCGCATCAGc                                                                                                            <  3:241173/149‑1 (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                                                                                             <  7:135719/146‑1 (MQ=255)                                                           tGGCCATGGCGTCGTAATCGTAGCCAGGGTCGTAGTCCGGCATGCTGCTGGTCGCCGGCGCTTGTTGCGCCAGCTGCTCGACCTGCTGCGGGTCGAGGCC‑GGTGATTTCCTTCAGGCGGTTGCGCATCAGCTGACGCAGGTTGgcg                                                                                             >  8:135719/1‑146 (MQ=255)                                                               cATGGCGTCGTAATCGTAGCCAGGGTCGTAGTCCGGCATGCTGCTGGTCGCCGGCGCTTGTTGCGCCAGCTGCTCGACCTGCTGCGGGTCGAGGCC‑GGTGATTTCCTTCAGGCGGTTGCGCATCAGCTGACGCAGGTTGGCGCCGGGGa                                                                                      <  3:259216/149‑1 (MQ=255)                                                               cATGGCGTCGTAATCGTAGCCAGGGTCGTAGTCCGGCATGCTGCTGGTCGCCGGCGCTTGTTGCGCCAGCTGCTCGACCTGCTGCGGGTCGAGGCC‑GGTGATTTCCTTCAGGCGGTTGCGCATCAGCTGACGCAGGTTGGCGCCGGGGa                                                                                      >  5:355950/1‑149 (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                           <  5:227904/135‑1 (MQ=255)                                                                                                                                  ccggcctCGACCCGCAGCAGGTCGA‑GCC‑GGTGATTTCCTTCAGGCGGTTGCGCATCAGCTGACGCAGGTTGGCGCCGGGGATCTTTTCGATCAACGGTGCGGCCAAGGTCGCCATATGCGCCTTGCCTTCCAGCGAgcgc                           >  6:227904/6‑140 (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)                                                                                                                                                ctgcGGGTCGAGGCC‑GGTGATTTCCTTCAGGCGGTTGCGCATCAGCTGACGCAGGTTGGCGCCGGGGATCTTTTCGATCAACGGTGCGGCCAAGGTCGCCATATGCGCCTTGCCTTCCAGCGAGCGCGGGTCGGCTTCTACGCCAAGtt     >  2:159303/1‑149 (MQ=255)                                                                                                                                                ctgcGGGTCGAGGCC‑GGTGATTTCCTTCAGGCGGTTGCGCATCAGCTGACGCAGGTTGGCGCCGGGGATCTTTTCGATCAACGGTGCGGCCAAGGTCGCCATATGCGCCTTGCCTTCCAGCGAGCGCGGGTCGGCTTCTACGCCAAGt      <  7:33545/148‑1 (MQ=255)                                                                                                                                                   cGGGTCGAGGCC‑GGTGATTTCCTTCAGGCGGTTGCGCATCAGCTAACGCAGGTTGGCGCCGGGGATCTTTTCGATCAACGGTGCGGCCAAGGTCGCCATATGCGCCTTGCCTTCCAGCGAGCGCGGGTCGGCTTCTACGCCAAGTTGCt  >  1:256663/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                     TCCTGATGAACGGGGGCATAGTCGTGCTGCGGCATGTCGCCGTAGTCGGGGGTGTAGCTGGCCATGGCGTCGTAATCGTAGCCAGGGTCGTAGTCCGGCATGCTGCTGGTCGCCGGCGCTTGTTGCGCCAGCTGCTCGACCTGCTGCGGGTCGAGGCC‑GGTGATTTCCTTCAGGCGGTTGCGCATCAGCTGACGCAGGTTGGCGCCGGGGATCTTTTCGATCAACGGTGCGGCCAAGGTCGCCATATGCGCCTTGCCTTCCAGCGAGCGCGGGTCGGCTTCTACGCCAAGTTGCT  >  NC_002947/473561‑473855

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