Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,593,758	G→A	9.7%	R265H (CGC→CAC)	PP_3171 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,593,758	0	G	A	9.7%	199.7 / 9.3	72	R265H (CGC→CAC)	PP_3171	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (31/34);  new base A (3/4);  total (34/38)
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.65e-01

ATGAACCTGTTCGCCGCTGCCAGCGATGCTGAGGCGGACTACCTGGCCAGCTCCCACCGCAAGTGGATGCTGGACCTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCG  >  NC_002947/3593610‑3593902                                                                                                                                                     |                                                                                                                                                 aTGAACCTGTTCGCCGCTGCCAGCGATGCTGAGGCGGACTACCTGGCCAGCTCCCACCGCAAGTGGATGCTGGACCTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGCCCACGAgcg                                                                                                                                                  >  6:225799/1‑149 (MQ=255)     aCCTGTTCGCCGCTGCCAGCGATGCTGAGGCGGACTACCTGGCCAGCTCCCACCGCAAGTGGATGCTGGACCTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAgcacggc                                                                                                                                              >  6:243158/1‑147 (MQ=255)         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                                                                                                                                        <  6:259526/148‑1 (MQ=255)            cgccgcTGCCAGCGATGCTGAGGCGGACTACCTGGCCAGCTCCCACCGCAAGTGGATGCTGGACCTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTc                                                                                                                                        >  4:364710/1‑148 (MQ=255)             gccgcTGCCAGCGATGCTGAGGCGGACTACCTGGCCAGCTCCCACCGTAAGTGGATGCTGGACCTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCACGGCGCGCTCgt                                                                                                                                      <  8:220520/149‑2 (MQ=255)             gccgcTGCCAGCGATGCTGAGGCGGACTACCTGGCCAGCTCCCACCGCAAGTGGATGCTGGACCTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCACGGCGCGCTCgt                                                                                                                                      <  6:382582/149‑2 (MQ=255)             gccgcTGCCAGCGATGCTGAGGCGGACTACCTGGCCAGCTCCCACCGCAAGTGGATGCTGGACCTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCACGGCGCGCTCg                                                                                                                                       <  5:243158/148‑1 (MQ=255)               cgcTGCCAGCGATGCTGAGGCGGACTACCTGGCCAGCTCCCACCGCAAGTGGATGCTGGACCTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGcc                                                                                                                                              <  4:288090/139‑1 (MQ=255)               cgcTGCCAGCGATGCTGAGGCGGACTACCTGGCCAGCTCCCACCGCAAGTGGATGCTGGACCTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGcc                                                                                                                                              >  3:288090/1‑139 (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                                                                                            >  4:30107/1‑135 (MQ=255)                                                                     gCTGGACCTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGt                                                                                            <  3:30107/135‑1 (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:197941/148‑1 (MQ=255)                                                                                   tCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGa                                                                 <  4:72523/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                    >  5:386316/1‑143 (MQ=255)                                                                                                                                     ttGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGa                    <  6:386316/143‑1 (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:280477/149‑1 (MQ=255)                                                                                                                                            aCCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCgg       <  4:341285/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                              >  7:280119/1‑122 (MQ=255)                                                                                                                                                cGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGcc                              <  8:280119/122‑1 (MQ=255)                                                                                                                                                 gAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCg  <  5:259355/149‑1 (MQ=255)                                                                                                                                                  aGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAg                                                                                                         <  6:337729/45‑1 (MQ=255)                                                                                                                                                  aGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAg                                                                                                         >  5:337729/1‑45 (MQ=255)                                                                                                                                                     |                                                                                                                                                 ATGAACCTGTTCGCCGCTGCCAGCGATGCTGAGGCGGACTACCTGGCCAGCTCCCACCGCAAGTGGATGCTGGACCTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCG  >  NC_002947/3593610‑3593902

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