Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,863,401:1	+A	5.8%	coding (1045/1068 nt)	PP_1666 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,863,401	1	.	A	5.8%	279.0 / 10.8	69	coding (1045/1068 nt)	PP_1666	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (27/38);  new base A (2/2);  total (29/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.55e-01

ACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTACGCCGTTCCTGGTCGGCATCATGCTGGC  >  NC_002947/1863259‑1863546                                                                                                                                                |                                                                                                                                                  aCAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCtt                                                                                                                                            <  1:770954/150‑1 (MQ=255)   aGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGc                                                                                                                                          >  1:890391/1‑150 (MQ=255)     cTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGt                                                                                                                                        >  1:3489/1‑150 (MQ=255)           gCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGc                                                                                                                                  <  2:561492/150‑1 (MQ=255)            cccAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCt                                                                                                                                 >  2:804630/1‑150 (MQ=255)            cccAGCTTGGCCTCGCCAAGCTGCAGGAAGTGCCCGTCGAGCAGGCACCGGCGGTGCTGGCCACGCCAGTACCGGCAGTCGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCt                                                                                                                                 <  2:265282/150‑1 (MQ=255)             ccAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTg                                                                                                                                <  1:402089/150‑1 (MQ=255)               aGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGt                                                                                                                              <  2:890391/150‑1 (MQ=255)                 cTTGGCCAGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGc                                                                                                                                              <  2:1087847/132‑1 (MQ=255)                 cTTGGCCAGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGc                                                                                                                                              >  1:1087847/1‑132 (MQ=255)                    ggcctggccAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGa                                                                                                                         <  2:969016/150‑1 (MQ=255)                    ggcctggccAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGa                                                                                                                         <  2:434365/150‑1 (MQ=255)                     gcctggccAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAg                                                                                                                        >  2:700356/1‑150 (MQ=255)                        tggccAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAgcgc                                                                                                                     <  2:237048/150‑1 (MQ=255)                             aaGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGcc                                                                                                                <  2:726127/150‑1 (MQ=255)                                 tGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACGCAAGCCG‑TCGAACGGCTTgggt                                                                                                                                        >  2:82422/1‑119 (MQ=255)                                 tGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACGCAAGCCG‑TCGAACGGCTTGGGt                                                                                                                                        <  1:82422/122‑1 (MQ=255)                                 tGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGtga                                                                                                            <  2:994768/150‑1 (MQ=255)                                  gCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCAGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGtgac                                                                                                           >  2:1006767/1‑150 (MQ=255)                                     ggAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGttt                                                                                                                                      >  1:566397/1‑120 (MQ=255)                                     ggAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGttt                                                                                                                                      <  2:566397/120‑1 (MQ=255)                                     ggAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGtgactga                                                                                                        <  1:700356/150‑1 (MQ=255)                                      gAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGtgactgac                                                                                                       >  2:35768/1‑150 (MQ=255)                                       aaGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACa                                                                                                      <  2:246760/150‑1 (MQ=255)                                              ggTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCg                                                                                               >  1:1261982/1‑150 (MQ=255)                                                 cGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTg                                                                                            <  1:700007/150‑1 (MQ=255)                                                  gAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTgg                                                                                           >  1:539285/1‑150 (MQ=255)                                                    gccggCACCGCCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCCCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAGCCCTTGCCGTGACTGACATGCGTCCCTGGAt                                                                                         >  2:661229/4‑150 (MQ=255)                                                    gCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGAt                                                                                         >  1:983333/1‑150 (MQ=255)                                                        gCACCGGCGGTGCCGGCCATGCCAGTACCGTCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTg                                                                                                                                >  1:1534024/1‑107 (MQ=255)                                                        gCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTg                                                                                                                                <  2:1534024/107‑1 (MQ=255)                                                            cggcggTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATctggctgg                                                                                 <  1:499195/150‑1 (MQ=255)                                                                      ggCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCAct                                                                       >  1:907134/1‑150 (MQ=255)                                                                         cATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCActgct                                                                    <  2:539285/150‑1 (MQ=255)                                                                          aTGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCActgctg                                                                   <  2:1261982/150‑1 (MQ=255)                                                                             ccAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATc                                                                >  1:498694/1‑150 (MQ=255)                                                                              cAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTg                                                                                                                  <  1:1094205/99‑1 (MQ=255)                                                                              cAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTg                                                                                                                  >  2:1094205/1‑99 (MQ=255)                                                                                gTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCt                                                             <  1:223682/150‑1 (MQ=255)                                                                                 tACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCtg                                                            <  2:614328/150‑1 (MQ=255)                                                                                   ccGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCtgtg                                                          <  2:3489/150‑1 (MQ=255)                                                                                    cGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGtgc                                                         >  1:420509/1‑150 (MQ=255)                                                                                        aGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACACGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTgctgct                                                     <  1:8041/150‑1 (MQ=255)                                                                                               gccgcAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGc                                              >  2:1120800/1‑150 (MQ=255)                                                                                               gccgcAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGc                                              >  1:602529/1‑150 (MQ=255)                                                                                                 cgcAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTaa                                                                                                                            >  1:437902/1‑70 (MQ=255)                                                                                                 cgcAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTaa                                                                                                                            <  2:437902/70‑1 (MQ=255)                                                                                                 cgcAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCt                                            <  2:498694/150‑1 (MQ=255)                                                                                                  gcAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGAGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTg                                           >  2:279659/1‑150 (MQ=255)                                                                                                        gcgcAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAat                                     <  1:378384/150‑1 (MQ=255)                                                                                                        gcgcAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAat                                     >  1:762221/1‑150 (MQ=255)                                                                                                          gcAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAatat                                   <  1:1120800/150‑1 (MQ=255)                                                                                                               ccGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTa                              >  1:1309968/1‑150 (MQ=255)                                                                                                               ccGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAAGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTa                              >  1:245637/1‑150 (MQ=255)                                                                                                                cgcgCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTAc                             <  2:550206/150‑1 (MQ=255)                                                                                                                   gcCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTACGcc                          <  2:383679/150‑1 (MQ=255)                                                                                                                   gcCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTACGcc                          <  1:1006767/150‑1 (MQ=255)                                                                                                                   gcCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTACGcc                          >  1:975192/1‑150 (MQ=255)                                                                                                                       gcgcAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGAGCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTTCGCCGGTc                      <  2:563817/150‑1 (MQ=255)                                                                                                                        ccatctcagcACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTCCTTCCTGCACAATATTCTTACGCCGTTcc                     <  1:43259/144‑1 (MQ=255)                                                                                                                         gcAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGAACTGCTGTTCGCTGTCCTGCTATATACCCTTGACAATATTCTTACGCCGTTCCt                    <  1:35768/150‑1 (MQ=255)                                                                                                                          cagccagcACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTACGCCGTTCCTg                   <  2:762221/150‑1 (MQ=255)                                                                                                                               agcACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTATTTGCTGTGCTGCTTTATAGCCTGCACAATATTCTTACGCCGTTCCtagatcg              >  2:43259/1‑144 (MQ=255)                                                                                                                                gcACCCAAGCCGTTCG‑ACGGCTTGCGTCTTCTCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTACGCCGTTCCTGGTCGGc             <  1:661229/150‑1 (MQ=255)                                                                                                                                aaacgcAAGCCG‑TCGAACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTACGCCGTTCCTGGTCGGc             >  2:581475/6‑150 (MQ=255)                                                                                                                                  aCCCAAGCCGTTCG‑ACGGCTTGTGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAata                                    >  2:1038698/1‑125 (MQ=255)                                                                                                                                  aCCCAAGCCGTTCG‑ACGGCTTGTGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAata                                    <  1:1038698/125‑1 (MQ=255)                                                                                                                                   cgcAAGCCG‑TCGAACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTACGCCGTTCCTGGTCGGcatc          <  1:581475/148‑1 (MQ=255)                                                                                                                                           gTTCG‑ACGTCGTTCGTTTTCGCCGATCAGCCGCCATTCCCGCGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTACGCCGTTCCTGGTCGGCATCATGCTGGc  <  2:676274/150‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                  ACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCG‑ACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTACGCCGTTCCTGGTCGGCATCATGCTGGC  >  NC_002947/1863259‑1863546

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