Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,173,343	T→C	6.7%	T146A (ACG→GCG)	PP_3672 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,173,343	0	T	C	6.7%	174.4 / 2.0	61	T146A (ACG→GCG)	PP_3672	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (28/28);  new base C (2/2);  total (31/30)
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.98e-01

CGCCGCGCAAGTCATTGTTCTTGGGTGATGACCTGCTTCAACCGCTCGCTCGCCGGTGCCGGCCGCCCGTAAAAGTAGCCCTGGAAGTGCGAACAGCCTTCCGCCGCGAGCCGGCGCATCTGCTCCTCCGACT‑CCACGCCTTCAGCCGTGATCTTGATCATCAGGCTGTTGGACAGCTCGGTAATGGCGCGAATGATCGACATGGCCTCGCGGCTGTCGCACATGTCATGCACGAACGACTTGTCGATCTTGATCCGGTCGAACGGGAACGAACGCAGGTAGCCCAACGAT  >  NC_002947/4173195‑4173485                                                                                                                                                      |                                                                                                                                               cgccgcGCAAGTCATTGTTCTTGGGTGATGACCTGCTTCAACCGCTCGCTCGCCGGTGCCGGCCGCCCGTAAAAGTAGCCCTGGAAGTGCGAACAGCCTTCCGCCGCGAGCCGGCGCATCTGCTCCTCCGACT‑CCACGCCTTCAGCCGt                                                                                                                                                >  3:167355/1‑149 (MQ=255)  gccgcGCAAGTCATTGTTCTTGGGTGATGACCTGCTTCAACCGCTCGCTCGCCGGTGCCGGCCGCCCGTAAAAGTAGCCCTGGAAGTGCGAACAGCCTTCCGCCGCGAGCCGGCGCATCTGCTCCTCCGACT‑CCACGCCTTCAGCCGt                                                                                                                                                >  4:60751/1‑148 (MQ=255)    cgcgCAAGTCATTGTTCTTGGGTGATGACCTGCTTCAACCGCTCGCTCGCCGGTGCCGGCCGCCCGTAAAAGTAGCCCTGGAAGTGCGAACAGCCTTCCGCCGCGAGCCGGCGCATCTGCTCCTCCGACT‑CCACGCCTTCAGCCGTGAt                                                                                                                                             <  6:236588/149‑1 (MQ=255)       gcAAGTCATTGTTCTTGGGTGATGACCTGCTTCAACCGCTCGCTCGCCGGTGCCGGCCGCCCGTAAAAGTAGCCCTGGAAGTGCGAACAGCCTTCCGCCGCGAGCCGGCGCATCTGCTCCTCCGACT‑CCACGCCTTCAGCCGTGATCtt                                                                                                                                          >  3:102184/1‑149 (MQ=255)           gTCATTGTTCTTGGGTGATGACCTGCTTCAACCGCTCGCTCGCCGGTGCCGGCCGCCCGTAAAAGTAGCCCTGGAAGTGCGAACAGCCTTCCGCCGCGAGCCGGCGCATCTGCTCCTCCGACT‑CCACGCCTTCAGCCGTGATCTTGat                                                                                                                                       >  6:92150/1‑148 (MQ=255)             cATTGTTCTTGGGTGATGACCTGCTTCAACCGCTCGCTCGCCGGTGCCGGCCGCCCGTAAAAGTAGCCCTGGAAGTGCGAACAGCCTTCCGCCGCGAGCCGGCGCATCTGCTCCTCCGACT‑CCACGCCTTCAGCCGTGATCTTGatcct                                                                                                                                    >  4:126620/1‑147 (MQ=255)                  ttCTTGGGTGATGACCTGCTTCAACCGCTCGCTCGCCGGTGCCGGCCGCCCGTAAAAGTAGCCCTGGAAGTGCGAACAGCCTTCCGCCGCGAGCCGGCGCATCTGCTCCTCCGACT‑CCACGCCTTCAGCCGTGATCTTGATCATCAGGc                                                                                                                               <  4:167355/149‑1 (MQ=255)                    cTTGGGTGATGACCTGCTTCAACCGCTCGCTCGCCGGTGCCGGCCGCCCGTAAAAGTAGCCCTGGAAGTGCGAACAGCCTTCCGCCGCGAGCCGGCGCATCTGCTCCTCCGACT‑CCACGCCTTCAGCCGTGATCTTGATCATCAGGCTg                                                                                                                             >  5:170772/1‑149 (MQ=255)                       gggTGATGACCTGCTTCAACCGCTCGCTCGCCTGTGCCGGCCGCCCGTAAAAGTAGCCCTGGAAGTGCGAACAGCCTTCCGCCGCGAGCCGGCGCATCTGCTCCTCCGACT‑CCACGCCTTCAGCCGTGATCTTGATCATCAGGCTGtt                                                                                                                           <  8:262162/148‑1 (MQ=255)                        ggTGATGACCTGCTTCAACCGCTCGCTCGCCGGTGCCGGCCGCCCGTAAAAGTAGCCCTGGAAGTGCGAACAGCCTTCCGCCGCGAGCCGGCGCATCTGCTCCTCCGACT‑CCACGCCTTCAGCCGTGATCTTGATCATCAGGCTGTTg                                                                                                                          <  3:60751/148‑1 (MQ=255)                          tgatgaCCTGCTTCAACCGCTCGCTCGCCGGTGCCGGCCGCCCGTAAAAGTAGCCCTGGAAGTGCGAACAGCCTTCCGCCGCGAGCCGGCGCATCTGCTCCTCCGACT‑CCACGCCTTCAGCCGTGATCTTGATCATCAGGCTGTTGGa                                                                                                                        >  2:247809/1‑148 (MQ=255)                            aagaCCTGCTTCAACAGCTCGCTCGCCGGTGCCGGCCGCCCGTAAAAGAACCCCTGGAAGTCCGAAAAGCCTTCCGCCACGAGCCGGCGCATCTGCTCCTCCGACT‑CCACTCCTTCAGCCGAGATCTTGATCATCAGGCTGTTGTACAg                                                                                                                     >  6:309337/3‑149 (MQ=255)                                ccTGCTTCAACCGCTCGCTCGCCGGTGCCGGCCGCCCGTAAAAGTAGCCCTGGAAGTGCGAACAGCCTTCCGCCGCGAGCCGGCGCATCTGCTCCTCCGACT‑CCACGCCTTCAGCCGTGATCTTGATCATCAGGCTGTTGGACAGCTCg                                                                                                                 <  4:258378/149‑1 (MQ=255)                                   gCTTCAACCGCTCGCTCGCCGGTGCCGGCCGCCCGTAAAAGTAGCCCTGGAAGTGCGAACAGCCTTCCGCCGCGAGCCGGCGCATCTGCTCCTCCGACT‑CCACGCCTTCAGCCGTGATCTTGAtcatca                                                                                                                                  <  6:136684/129‑1 (MQ=255)                                   gCTTCAACCGCTCGCTCGCCGGTGCCGGCCGCCCGTAAAAGTAGCCCTGGAAGTGCGAACAGCCTTCCGCCGCGAGCCGGCGCATCTGCTCCTCCGACT‑CCACGCCTTCAGCCGTGATCTTGAtcatca                                                                                                                                  >  5:136684/1‑129 (MQ=255)                                     ttCAACCGCTCGCTCGCCGGTGCCGGCCGCCCGTAAAAGTAGCCCTGGAAGTGCGAACAGCCTTCCGCCGCGAGCCGGCGCATCTGCTCCTCCGACT‑CCACGCCTTCGGCCGTGATCTTGATCATCAGGCTGTTGGACAGCTCGGTAAt                                                                                                            >  5:276050/1‑149 (MQ=255)                                       cAACCGCTCGCTCGCCGGTGCCGGCCGCCCGTAAAAGTAGCCCTGGAAGTGCGAACAGCCTTCCGCCGCGAGCCGGCGCATCTGCTCCTCCGACT‑CCACGCCTTCAGCCGTGATCTTGat                                                                                                                                       <  8:173224/120‑1 (MQ=255)                                         aCCGCTCGCTCGCCGGTGCCGGCCGCCCGTAAAAGTAGCCCTGGAAGTGCGAACAGCCTTCCGCCGCGAGCCGGCGCATCTGCTCCTCCGACTCCCACGCCTTCAGCCGTGATCTTGATCATCAGGCTGTTGGACAGCTCGGTAATGgc                                                                                                         <  1:247809/149‑1 (MQ=255)                                          ccGCTCGCTCGCCGGTGCCGGCCGCCCGTAAAAGTAGCCCTGGAAGTGCGAACAGCCTTCCGCCGCGAGCCGGCGCATCTGCTCCTCCGACT‑CCACGCCTTCAGCCGTGATCTTGATCATCAGGCTGTTGGACAGCTCGGTAATGgcgc                                                                                                       >  7:63698/1‑149 (MQ=255)                                           cgctcgctCGCCGGTGCCGGCCGCCCGTAAAAGTAGCCCTGGAAGTGCGAACAGCCTTCCGCCGCGAGCCGGCGCATCTGCTCCTCCGACT‑CCACGCCTTCAGCCGTGATCTTGATCATCAGGCTGTTGGACAGCTCGGTAATgg                                                                                                          <  8:58121/145‑1 (MQ=255)                                           cgctcgctCGCCGGTGCCGGCCGCCCGTAAAAGTAGCCCTGGAAGTGCGAACAGCCTTCCGCCGCGAGCCGGCGCATCTGCTCCTCCGACT‑CCACGCCTTCAGCCGTGATCTTGATCATCAGGCTGTTGGACAGCTCGGTAATgg                                                                                                          >  7:58121/1‑145 (MQ=255)                                             ctcgctcgCCGGTGCCGGCCGCCCGTAAAAGTAGCCCTGGAAGTGCGAACAGCCTTCCGCCGCGAGCCGGCGCATCTGCTCCTCCGACT‑CCACGCCTTCAGCCGTGATCTTGATCATCAGGCTGTTGGACAGCTCGGTAATGGCGCGaa                                                                                                    >  7:28907/1‑149 (MQ=255)                                                    gcCGGTGCCGGCCGCCCGTAAAAGTAGCCCTGGAAGTGCGAACAGCCTTCCGCCGCGAGCCGGCGCATCTGCTCCTCCGACT‑CCACGCCTTCAGCCGTGATCTTGATCATCAGGCTGTTGGACAGCTCGGTAATGGCGCGAATGATCGa                                                                                             >  6:831/1‑149 (MQ=255)                                                    gcCGGTGCCGGCCGCCCGTAAAAGTAGCCCTGGAAGTGCGAACAGCCTTCCGCCGCGAGCCGGCGCATCTGCTCCTCCGACT‑CCACGCCTTCAGCCGTGATCTTGATCATCAGGCTGTTGGACAGCTCGGTAATGGCGCGAATGATCGa                                                                                             >  4:169997/1‑149 (MQ=255)                                                     ccGGTGCCGGCCGCCCGTAAAAGTAGCCCTGGAAGTGCGAACAGCCTTCCGCCGCGAGACGGCGCATCTGCTCCTCCGACT‑CCACGCCTTCAGCCGTGATCTTGATCATCAGGCTGTTGGACAGCTCGGTAATGGCGCGAATGATCGAc                                                                                            <  3:302261/149‑1 (MQ=255)                                                      cGGTGCCGGCCGCCCGTAAAAGTAGCCCTGGAAGTGCGAACAGCCTTCCGCCGCGAGCCGGCGCATCTGCTCCTCCGACT‑CCACGCCTTCAGCCGTGATCTTGATCATCAGGCTGTTGGACAGCTCGGTAATGGCGCGAATGATCGAc                                                                                            >  4:130707/1‑148 (MQ=255)                                                      cGGTGCCGGCCGCCCGTAAAAGTAGCCCTGGAAGTGCGAACAGCCTTCCGCCGCGAGCCGGCGCATCTGCTCCTCCGACT‑CCACGCCTTCAGCCGTGATCTTGATCATCAGGCTGTTGGACAGCTCGGTAATGGCGCGAATGATCGAc                                                                                            <  3:130707/148‑1 (MQ=255)                                                        gTGCCGGCCGCCCGTAAAAGTAGCCCTGGAAGTGCGAACAGCCTTCCGCCGCGAGCCGGCGCATCTGCTCCTCCGACT‑CCACGCCTTCAGCCGTGATCTTGATCATCAGGCTGTTGGACAGCTCGGTAATGGCGCGAATGATCGACATg                                                                                         <  6:170772/149‑1 (MQ=255)                                                                   ccGTAAAAGTAGCCCTGGAAGTGCGAACAGCCTTCCGCCGCGAGCCGGCGCATCTGCTCCTCCGACT‑CCACGCCTTCAGCCGTGATCTTGATCATCAGGCTGTTGGACAGCTCGGTAATGGCGCGAATGATCGACATGGCCTCGCGGCt                                                                              >  5:95909/1‑149 (MQ=255)                                                                    cGTAAAAGTAGCCCTGGAAGTGCGAACAGCCTTCCGCCGCGAGCCGGCGCATCTGCTCCTCCGACT‑CCACGCCTTCAGCCGTGATCTTGATCATCAGGCTGTTGGACAGCTCGGTAATGGCGCGAATGATCGACATGGCCTCGCGGCTg                                                                             >  5:327337/1‑149 (MQ=255)                                                                    cGTAAAAGTAGCCCTGGAAGTGCGAACAGCCTTCCGCCGCGAGCCGGCGCATCTGCTCCTCCGACT‑CCACGCCTTCAGCCGCGATCTTGATCATCAGGCTGTTGGACAGCTCGGTAATGGCGCGAATGATCGACATGGCCTCGCGGCTg                                                                             >  6:316613/1‑149 (MQ=255)                                                                    cGTAAAAGTACCCCTGGAAGTGCGAACAGCCTTCCGCCGCGAGCCGGCGCATCTGCTCCTCCGACT‑CCACGCCTTCAGCCGCGATCTTCATCATCAGGCTGTTGGACAGCTCGGTAATGGCGCGAATGATCGACATGGCCTCGCGGCt                                                                              >  4:318439/1‑148 (MQ=255)                                                                      tAAAAGTAGCCCTGGAAGTGCGAACAGCCTTCCGCCGCGAGCCGGCGCATCTGCTCCTCCGACT‑CCACGCCTTCAGCCGTGATCTTGATCATCAGGCTGTTGGACAGCTCGGTAATGGCGCGAATGATCGACATGGCCTCGCGGCTGTc                                                                           >  2:200183/1‑149 (MQ=255)                                                                               cccTGGAAGTGCGAACAGCCTTCCGCCGCGAGCCGGCGCATCTGCTCCTCCGACT‑CCACGCCTTCAGCCGTGATCTTGATCATCAGGCTGTTGGACAGCTCGGTAATGGCGCGAATGATCGACATGGCCTCGCGGCTGTCGCACATGt                                                                   >  2:47076/1‑148 (MQ=255)                                                                                 cTGGAAGTGCGAACAGCCTTCCGCCGCGAGCCGGCGCATCTGCTCCTCCGACT‑CCACGCCTTCAGCCGTGATCTTGATCATCAGGCTGTTGGACAGCTCGGTAATGGCGCGAATGATCGACATGGCCTCGCGGCTGTCGca                                                                        >  7:43285/1‑141 (MQ=255)                                                                                 cTGGAAGTGCGAACAGCCTTCCGCCGCGAGCCGGCGCATCTGCTCCTCCGACT‑CCACGCCTTCAGCCGTGATCTTGATCATCAGGCTGTTGGACAGCTCGGTAATGGCGCGAATGATCGACATGGCCTCGCGGCTGTCGca                                                                        <  8:43285/141‑1 (MQ=255)                                                                                 cTGGAAGTGCGAACAGCCTTCCGCCGCGAGCCGGCGCATCTGCTCCTCCGACT‑CCACGCCTTCAGCCGTGATCTTGATCATCAGGCTGTTGGACAGCTCGGTAATGGCGCGAATGATCGACATGGCCTCGCGGCTGTCGCACATGTCAt                                                                >  3:60057/1‑149 (MQ=255)                                                                                       gTGCGAACAGCCTTCCGCCGCGAGCCGGCGCATCTGCTCCTCCGACT‑CCACGCCTTCAGCCGTGATCTTGATCATCAGGCTGTTGGACAGCTCGGTAATGGCGCGAATGATCGACATGGCCTCGCGGCTGTCGCACATGTCATGCacg                                                           >  8:73424/1‑148 (MQ=255)                                                                                        tGCGAACAGCCTTCCGCCGCGAGCCGGCGCATCTGCTCCTCCGACT‑CCACGCCTTCAGCCGTGATCTTGATCATCAGGCTGTTGGACAGCTCGGTAATGGCGCGAATGATCGACATGGCCTCGCGGCTGTCGCACATGTCATGCacgaa                                                         <  5:831/149‑1 (MQ=255)                                                                                                  cTTCCGCCGCGAGCCGGCGCATCTGCTCCTCCGACT‑CCACGCCTTCAGCCGTGATCTTGATCATCAGGCTCTTGGACAGCTCGCTAATGGCTCGAATGATCGACATGGCCTAGCGGCTGTCGCACATGACATGCACGAACGAcatctc                                                >  8:301927/1‑143 (MQ=255)                                                                                                   ttCCGCCGCGAGCCGGCGCATCTGCTCCTCCGACT‑CCACGCCTTCAGCCGTGATCTTGATCATCAGGCTGTTGGACAGCTCGGTAATGGCGCGAATGATCGACATGGCCTCGCGGCTGTCGCACATGTCATGCACGAACGACTTg                                                  >  1:222656/1‑145 (MQ=255)                                                                                                   ttCCGCCGCGAGCCGGCGCATCTGCTCCTCCGACT‑CCACGCCTTCAGCCGTGATCTTGATCATCAGGCTGTTGGACAGCTCGGTAATGGCGCGAATGATCGACATGGCCTCGCGGCTGTCGCACATGTCATGCACGAACGACTTg                                                  <  2:222656/145‑1 (MQ=255)                                                                                                    tCCGCCGCGAGCCGGCGCATCTGCTCCTCCGACT‑CCACGCCTTCAGCCGTGATCTTGATCATCAGGCTGTTGGACAGCTCGGTAATGGCGCGAATGATCGACATGGCCTCGCGGCTGTCGCACATGTCATGCACGAACGACTTGTCGAt                                             >  5:127125/1‑149 (MQ=255)                                                                                                        ccgcGAGCCGGCGCATCTGCTCCTCCGACT‑CCACGCCTTCAGCCGTGATCTTGATCATCAGGCTGTTGGACAGCTCGGTAATGGCGCGAATGATCGACATGGCCTCGCGGCTGTCGCACATGTCATGCACGAACGACTTGTCGATCTTg                                         <  1:47076/149‑1 (MQ=255)                                                                                                        ccgcGAGCCGGCGCATCTGCTCCTCCGACT‑CCACGCCTTCAGCCGTGATCTTGATCATCAGGCTGTTGGACAGCTCGGTAATGGCGCGAATGATCGACATGGCCTCGCGGCTGTCGCACATGTCATGCACGAACGACTTGTCGATCTTg                                         <  4:102184/149‑1 (MQ=255)                                                                                                             aGCCGGCGCATCTGCTCCTCCGACT‑CCACGCCTTCAGCCGTGATCTTGATCATCAGGCTGTTGGACAGCTCGGTAATGGCGCGAATGATCGACATGGCCTCGCGGCTGTCGCACATGTCATGCACGAACGACTTGTCGATCTTGATCCg                                    <  6:127125/149‑1 (MQ=255)                                                                                                                  gcgcATCTGCTCCTCCGACT‑CCACGCCTTCAGCCGTGATCTTGATCATCAGGCTGTTGGACAGCTCGGTAATGGCGCGAATGATCGACATGGCCTCGCGGCTGTCGCACATGTCATGCACGAACGACTTGTCGATCTTGATCCGGTCg                                >  4:67032/1‑148 (MQ=255)                                                                                                                   cgcATCTGCTCCTCCGACT‑CCACGCCTTCAGCCGTGATCTTGATCATCAGGCTGTTGGACAGCTCGGTAATGGCGCGAATGATCGACATGGCCTCGCGGCTGTCGCACATGTCATGCACGAACGACTTGTCGATCTTGATCCGGGCGaa                              <  8:28907/149‑1 (MQ=255)                                                                                                                     cATCTGCTCCTCCGACT‑CCACGCCTTCAGCCGTGATCTTGATCATCAGGCTGTTGGACAGCTCGGTAATGGCGCGAATGATCGACATGGCCTCGCGGCTGTCGCACATGTCATGCACGAACGACTTGTCGATCTTGATCCGGTCGAACg                            <  1:200183/149‑1 (MQ=255)                                                                                                                      aTCTGCTCCTCCGACT‑CCACGCCTTCAGCCGTGATCTTGATCATCAGGCTGTTGGACAGCTCGGTAATGGCGCGAATGATCGACATGGCCTCGCGGCTGTCGCACTTGTCATGCACGAACGACTTGTCGATCTTGATCCGGTCGAACgg                           <  7:301927/149‑1 (MQ=255)                                                                                                                      aTCTGCTCCTCCGACT‑CCACGCCTTCAGCCGTGATCTTGATCATCAGGCTGTTGGACAGCTCGGTAATGGCGCGAATGATCGACATGGCCTCGCGGCTGTCGCACATGTCATGCACAAACGACTTGTCGATCTTGATCCGGTCGAACgg                           <  4:60057/149‑1 (MQ=255)                                                                                                                           ctcctcCGACT‑CCACGCCTTCAGCCGTGATCTTGATCATCAGGCTGTTGGACAGCTCGGTAATGGCGCGAATGATCGACATGGCCTCGCGGCTGTCGCACATGTCATGCACGAACGACTTGTCGATCTTGATCCGGTCGAACggg                          <  6:327337/145‑1 (MQ=255)                                                                                                                           ctcctcCGACT‑CCACGCCTTCAGCCGTGATCTTGATCATCAGGCTGTTGGACAGCTCGGTAATGGCGCGAATGATCGACATGGCCTCGCGGCTGTCGCACATGTCATGCACGAACGACTTGTCGATCTTGATCCGGTCGAACGGgaac                       >  5:149374/1‑148 (MQ=255)                                                                                                                           ctcctcCGACT‑CCACGCCTTCAGCCGTGATCTTGATCATCAGGCTGTTGGACAGCTCGGTAATGGCGCGAATGATCGACATGGCCTCGCGGCTGTCGCACATGTCATGCACGAACGACTTGTCGATCTTGATCCGGTCGAACGGgaac                       <  7:73424/148‑1 (MQ=255)                                                                                                                           ctcctcCGACT‑CCACGCCTTCAGCCGTGATCTTGATCATCAGGCTGTTGGACAGCTCGGTAATGGCGCGAATGATCGACATGGCCTCGCGGCTGTCGCACATGTCATGCACGAACGACTTGTCGATCTTGATCCGGTCGAACGGgaa                        <  6:95909/147‑1 (MQ=255)                                                                                                                              ctccGACT‑CCACGCCTTCAGCCGTGATCTTGATCATCAGGCTGTTGGACAGCTCGGTAATGGCGCGAATGATCGACATGGCCTCGCGGCTGTCGCACATGTCATGCACGAACGACTTGTCGATCTTGATCCGGTCGAACGGgaacga                     <  5:92150/147‑1 (MQ=255)                                                                                                                               tccGACT‑CCACGCCTTCAGCCGTGATCTTGATCATCAGGCTGTTGGACAGCTCGGTAATGGCGCGAATGATCGACATGGCCTCGCGGCTGTCGCACATGTCATGCACGAACGACTTGTCGATCTTGATCCGGTCGAACGGGaacgaacg                  >  5:192374/1‑149 (MQ=255)                                                                                                                                       ccACGCCTTCAGCCGTGATCTTGATCATCAGGCTGTTGGACAGCTCGGTAATGGCGCGAATGATCGACATGGCCTCGCGGCTGTCGCACATGTCATGCACGAACGACTTGTCGATCTTGATCCGGTCGAACGGGAACGAACGCAGGTAg           >  1:250655/1‑149 (MQ=255)                                                                                                                                        cACGCCTTCAGCCGTGATCTTGATCATCAGGCTGTTGGACAGCTCGGTAATGGCGCGAATGATCGACATGGCCTCGCGGCTGTCGCACATGTCATGCACGAACGACTTGTCGATCTTGATCCGGTCGAACGGGAACGAACGCAGGTAg           <  6:192374/148‑1 (MQ=255)                                                                                                                                          cGCCTTCAGCCGTGATCTTGATCATCAGGCTGTTGGACAGCTCGGTAATGGCGCGAATGATCGACATGGCCTCGCGGCTGTCGCACATGTCATGCACGAACGACTTGTCGATCTTGATCCGGTCGAACGGGAACGAACGCAGGTAGcc         >  2:136723/1‑148 (MQ=255)                                                                                                                                                 aGCCGTGATCTTGATCATCAGGCTGTTGGACAGCTCGGTAATGGCGCTAATGATCGACATGGCCTCGCGGCTGTCGCACATGTCATGCACGAACGACTTGTCGATCTTGATCCGGTCGAACGGGAACGAACGCAGGTAGCCCAACGAt  <  5:309337/148‑1 (MQ=255)                                                                                                                                                 aGCCGCGATCTTGATCATCAGGCTGTTGGACAGCTCGGTAATGGCGCGAATGATCGACATGGCCTCGCGGCTGTCGCACATGTCATGCACGAACGACTTGTCGATCTTGATCCGGTCGAACGGGAACGAACGCAGGTAGCCCAACGAt  <  3:318439/148‑1 (MQ=255)                                                                                                                                                 aGCCGCGATCTTGATCATCAGGCTGTTGGACAGCTCGGTAATGGCGCGAATGATCGACATGGCCTCGCGGCTGTCGCACATGTCATGCACGAACGACTTGTCGATCTTGATCCGGTCGAACGGGAACGAACGCAGGTAGCCCAACGAt  <  5:316613/148‑1 (MQ=255)                                                                                                                                                      tGATCTTGATCATCAGGCTGTTGGACAGCTCGGTTATGGCGCGAATGATCGACATGGCCTCGCGGCTGTCGCACATGTCATGCACGAACGACTTGTCGATCTTGATCCGGTCGAACGGga                         <  8:39238/120‑1 (MQ=255)                                                                                                                                                      tGATCTTGATCATCAGGCTGTTGGACAGCTCGGTAATGGCGCGAATGATCGACATGGCCTCGCGGCTGTCGCACATGTCATGCACGAACGACTTGTCGATCTTGATCCGGTCGAACGGga                         >  7:39238/1‑120 (MQ=255)                                                                                                                                                      |                                                                                                                                               CGCCGCGCAAGTCATTGTTCTTGGGTGATGACCTGCTTCAACCGCTCGCTCGCCGGTGCCGGCCGCCCGTAAAAGTAGCCCTGGAAGTGCGAACAGCCTTCCGCCGCGAGCCGGCGCATCTGCTCCTCCGACT‑CCACGCCTTCAGCCGTGATCTTGATCATCAGGCTGTTGGACAGCTCGGTAATGGCGCGAATGATCGACATGGCCTCGCGGCTGTCGCACATGTCATGCACGAACGACTTGTCGATCTTGATCCGGTCGAACGGGAACGAACGCAGGTAGCCCAACGAT  >  NC_002947/4173195‑4173485

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