Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,593,764	T→C	8.4%	V267A (GTG→GCG)	PP_3171 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,593,764	0	T	C	8.4%	198.1 / 7.5	71	V267A (GTG→GCG)	PP_3171	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (31/34);  new base C (2/4);  total (33/38)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.81e-01

GTTCGCCGCTGCCAGCGATGCTGAGGCGGACTACCTGGCCAGCTCCCACCGCAAGTGGATGCTGGACCTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCC  >  NC_002947/3593618‑3593912                                                                                                                                                   |                                                                                                                                                     gTTCGCCGCTGCCAGCGATGCTGAGGCGGACTACCTGGCCAGCTCCCACCGCAAGTGGATGCTGGACCTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCACGgcgcgc                                                                                                                                                    >  5:382582/1‑149 (MQ=255)  ttCGCCGCTGCCAGCGATGCTGAGGCGGACTACCTGGCCAGCTCCCACCGCAAGTGGATGCTGGACCTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGc                                                                                                                                                    <  5:386098/148‑1 (MQ=255)  ttCGCCGCTGCCAGCGATGCTGAGGCGGACTACCTGGCCAGCTCCCACCGCAAGTGGATGCTGGACCTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCt                                                                                                                                                   <  7:302821/149‑1 (MQ=255)  ttCGCCGCTGCCAGCGATGCTGAGGCGGACTACCTGGCCAGCTCCCAACGCAAGTTGATGCAGGACCGGGAAGGGGGGCGGGTTGGGTTTCAGCCCAAGCCGCAAGAGGGGTTCCTGGAACGGGTGCCGGACCACGCGCGCGCCGTGCt                                                                                                                                                   >  1:83041/1‑149 (MQ=255)    cgccgcTGCCAGCGATGCTGAGGCGGACTACCTGGCCAGCTCCCACCGCAAGTGGATGCTGGACCTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTc                                                                                                                                                  >  4:364710/1‑148 (MQ=255)    cgccgcTGCCAGCGATGCTGAGGCGGACTACCTGGCCAGCTCCCACCGCAAGTGGATGCTGGACCTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTc                                                                                                                                                  <  6:259526/148‑1 (MQ=255)     gccgcTGCCAGCGATGCTGAGGCGGACTACCTGGCCAGCTCCCACCGTAAGTGGATGCTGGACCTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCACGGCGCGCTCgt                                                                                                                                                <  8:220520/149‑2 (MQ=255)     gccgcTGCCAGCGATGCTGAGGCGGACTACCTGGCCAGCTCCCACCGCAAGTGGATGCTGGACCTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCACGGCGCGCTCgt                                                                                                                                                <  6:382582/149‑2 (MQ=255)     gccgcTGCCAGCGATGCTGAGGCGGACTACCTGGCCAGCTCCCACCGCAAGTGGATGCTGGACCTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCACGGCGCGCTCg                                                                                                                                                 <  5:243158/148‑1 (MQ=255)           gCCAGCGATGCTGAGGCGGACTACCTGGCCAGCTCCCACCGCAAGTGGATGCTGGACCTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCACGGCGCGCTCgtggtccg                                                                                                                                          <  5:255844/149‑8 (MQ=255)                cGATGCTGAGGCGGACTACCTGGCCAGCTCCCACCGCAAGTGGATGCTGGACCTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATgg                                                                                                                                     >  1:14437/1‑149 (MQ=255)                      tGAGGCGGACTACCTGGCCAGCTCCCACCGCAAGTGGATGCTGGACCTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCt                                                                                                                               >  4:226901/1‑149 (MQ=255)                                  ccTGGCCAGCTCCCACCGCAAGTGGATGCTGGACCTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCgg                                                                                                                     <  2:14437/147‑1 (MQ=255)                                   cTGGCCAGCTCCCACCGCAAGTGGATGCTGGACCTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAgg                                                                                                                  >  6:24264/1‑149 (MQ=255)                                        cAGCTCCCACCGCAAGTGGATGCTGGACCTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGccc                                                                                                              >  7:297344/1‑148 (MQ=255)                                        cAGCTCCCACCGCAAGTGGATGCTGGACCTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGccc                                                                                                              >  3:72523/1‑148 (MQ=255)                                            tCCCACCGCAAGTGGATGCTGGACCTGCATGTCGGGCGCTTTGGCTTGCTGCCGAAGCCGCAAGAGGGTTACCTCGAACGCTTGACGGACCAAGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTCGCCGGAGGGCCCGaagg                                                                                                         >  3:235385/1‑147 (MQ=255)                                            tCCCACCGCAAGTGGATGCTGGACCTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACg                                                                                                         >  3:184422/1‑149 (MQ=255)                                                 ccGCAAGTGGATGCTGGACCTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCGCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGa                                                                                                    <  5:225799/149‑1 (MQ=255)                                                  cGCAAGTGGATGCTGGACCTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAg                                                                                                   >  5:189896/1‑149 (MQ=255)                                                      aGTGGATGCTGGACCTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGaa                                                                                               <  4:291842/149‑1 (MQ=255)                                                       gTGGATGCTGGACCTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAg                                                                                              <  5:208897/149‑1 (MQ=255)                                                             gCTGGACCTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGt                                                                                                      <  3:30107/135‑1 (MQ=255)                                                             gCTGGACCTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGt                                                                                                      >  4:30107/1‑135 (MQ=255)                                                                ggACCTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGcgct                                                                                     >  8:305040/1‑149 (MQ=255)                                                                   ccTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGcgctcgc                                                                                  <  2:42301/149‑1 (MQ=255)                                                                   ccTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGcgctcgc                                                                                  <  4:184422/149‑1 (MQ=255)                                                                    cTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGcgctcgct                                                                                 >  5:25721/1‑149 (MQ=255)                                                                      gCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCg                                                                                                             <  4:296827/119‑1 (MQ=255)                                                                      gCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCg                                                                                                             >  3:296827/1‑119 (MQ=255)                                                                      gCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAggg                                                                                                 <  4:256605/131‑1 (MQ=255)                                                                      gCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAggg                                                                                                 >  3:256605/1‑131 (MQ=255)                                                                      gCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGa                                                                               >  7:310057/1‑149 (MQ=255)                                                                      gCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGa                                                                               >  6:129778/1‑149 (MQ=255)                                                                       cATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGAt                                                                              <  3:226901/149‑1 (MQ=255)                                                                        aTGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATc                                                                             <  5:166420/149‑1 (MQ=255)                                                                         tGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCg                                                                            <  5:24264/149‑1 (MQ=255)                                                                           tCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGa                                                                           <  4:72523/148‑1 (MQ=255)                                                                           tCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGa                                                                           <  4:197941/148‑1 (MQ=255)                                                                               gcgcTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGa                                                                      >  2:323716/1‑149 (MQ=255)                                                                                 gcTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACGTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGAcc                                                                    >  1:38125/1‑149 (MQ=255)                                                                                            ctgcCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGGCTGCTCCGTAGCCGGTGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGa                                                         <  8:310057/149‑1 (MQ=255)                                                                                            ctgcCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGa                                                         <  6:189896/149‑1 (MQ=255)                                                                                             tgcCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAg                                                        >  5:281633/1‑149 (MQ=255)                                                                                              gcCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGc                                                       >  8:99585/1‑149 (MQ=255)                                                                                                ccAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCtg                                                     <  7:305040/149‑1 (MQ=255)                                                                                                 cAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGCGCAGACCGGGGCGGACGAGCtga                                                    >  2:284048/1‑149 (MQ=255)                                                                                                     ccGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTgatgat                                                <  5:129778/149‑1 (MQ=255)                                                                                                     caggaaGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTgatgat                                                <  1:325429/145‑1 (MQ=255)                                                                                                      cGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATc                                               >  3:280477/1‑149 (MQ=255)                                                                                                                    cTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCa                                 >  4:200829/1‑149 (MQ=255)                                                                                                                    cTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCa                                 >  3:341285/1‑149 (MQ=255)                                                                                                                       gAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGa                              <  2:38125/149‑1 (MQ=255)                                                                                                                       gAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGa                              >  6:259355/1‑149 (MQ=255)                                                                                                                       gAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGAGCGACTGCCGCATCCAt                                <  3:364710/147‑1 (MQ=255)                                                                                                                             ttGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGa                              <  6:386316/143‑1 (MQ=255)                                                                                                                             ttGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGa                              >  5:386316/1‑143 (MQ=255)                                                                                                                               gCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTg                       <  2:83041/148‑1 (MQ=255)                                                                                                                                   ctgatcGAGCACGGCGCGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGcgcg                  >  3:4062/6‑149 (MQ=255)                                                                                                                                    tgatcGAGCACGGCGCGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCgg                 >  5:319500/5‑149 (MQ=255)                                                                                                                                    aCCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCgg                 <  4:341285/149‑1 (MQ=255)                                                                                                                                    aCCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCgg                 <  4:280477/149‑1 (MQ=255)                                                                                                                                    aCCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCgg                 <  7:99585/149‑1 (MQ=255)                                                                                                                                       aCGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCtga                                                    <  7:165088/111‑1 (MQ=255)                                                                                                                                       aCGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCtga                                                    >  8:165088/1‑111 (MQ=255)                                                                                                                                        cGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGcc                                        <  8:280119/122‑1 (MQ=255)                                                                                                                                        cGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGcc                                        >  7:280119/1‑122 (MQ=255)                                                                                                                                         gAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCg            <  5:259355/149‑1 (MQ=255)                                                                                                                                          aGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAg                                                                                                                   >  5:337729/1‑45 (MQ=255)                                                                                                                                          aGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAg                                                                                                                   <  6:337729/45‑1 (MQ=255)                                                                                                                                                 cGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTc         >  1:272338/1‑144 (MQ=255)                                                                                                                                                 cGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTc         <  2:272338/144‑1 (MQ=255)                                                                                                                                                   tGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGcc  >  5:132503/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                     GTTCGCCGCTGCCAGCGATGCTGAGGCGGACTACCTGGCCAGCTCCCACCGCAAGTGGATGCTGGACCTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCC  >  NC_002947/3593618‑3593912

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