Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	336,125	+T	coding (210/237 nt)	PP_0278 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	336,124	1	.	T	98.4%	215.5 / ‑3.4	63	coding (209/237 nt)	PP_0278	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base T (32/30);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.05e-01

TTAACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATA‑TGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑A‑TACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTC  >  NC_002947/335977‑336260                                                                                                                                                       |                                                                                                                                         ttAACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTGGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGTGCGATA‑TGGAGGACGACCATAAATTCCTTGGGGCAGCCG‑A‑t                                                                                                                                         <  3:297928/149‑1 (MQ=255)    aCCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATA‑TGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑A‑ttac                                                                                                                                      <  5:86016/149‑4 (MQ=255)        aTCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATA‑TGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCtgt                                                                                                                                   >  8:41959/1‑149 (MQ=255)          cccTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGGCGACGCGGCCTACGGGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATA‑TGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGc                                                                                                                                 <  8:98535/149‑1 (MQ=255)           ccTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGATGGCGATATTGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGc                                                                                                                                 <  1:76231/149‑1 (MQ=255)               aaTCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATA‑TGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGAc                                                                                                                              <  6:196756/147‑1 (MQ=255)               aaTCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATA‑TGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGAc                                                                                                                              <  7:196588/147‑1 (MQ=255)                     cTACTCGTAGAGCCCGTCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTTCTCGAGGGCGATA‑TGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATgg                                                                                                                          <  6:330520/145‑1 (MQ=255)                     cTACGCGTAGAGCCCGGCGTCCCCATCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATA‑TGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCaa                                                                                                                      <  3:179403/149‑1 (MQ=255)                        cgcgTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATA‑TGGAGGACGTCCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAg                                                                                                                     <  1:81436/147‑1 (MQ=255)                        cgcgTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATA‑TGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAgg                                                                                                                    <  8:105375/148‑1 (MQ=255)                          cgTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATA‑TGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGgcgc                                                                                                                 <  4:83085/149‑1 (MQ=255)                          cgTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATA‑TGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGgcgc                                                                                                                 >  4:105629/1‑149 (MQ=255)                           gTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGCTCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATA‑TGGAGGACGACCATAAATAGTAGGGGGCCGCCG‑ATTAACGGTGCGACATGGCCAAGGCGCt                                                                                                                >  1:308790/1‑149 (MQ=255)                                        tCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATA‑TGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATc                                                                                                   <  4:12760/149‑1 (MQ=255)                                           ccGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATA‑TGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCtgat                                                                                                             >  7:342766/1‑136 (MQ=255)                                           ccGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATA‑TGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCtgat                                                                                                             <  8:342766/136‑1 (MQ=255)                                           ccGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATA‑TGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATc                                                                                                   >  6:2005/1‑146 (MQ=255)                                             gTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATA‑TGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATc                                                                                                   >  3:34658/1‑144 (MQ=255)                                                 aCGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATA‑TGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGc                                                                                           >  3:41095/1‑148 (MQ=255)                                                        cTACGAGCATGTCTCGATCCAGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATA‑TGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGt                                                                                    >  2:125684/1‑148 (MQ=255)                                                          aCGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATA‑TGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAg                                                                                 >  2:255188/1‑149 (MQ=255)                                                           cGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATA‑TGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAg                                                                                 <  3:105629/148‑1 (MQ=255)                                                            gAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATA‑TGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGgc                                                                               >  6:367684/1‑149 (MQ=255)                                                              gCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATA‑TGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGgcg                                                                              >  4:225133/1‑148 (MQ=255)                                                                    ctcGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATA‑TGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAg                                                                        <  4:147058/148‑1 (MQ=255)                                                                             gATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATA‑TGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGAcgc                                                              >  6:99819/1‑149 (MQ=255)                                                                                gggCTACATCAAGCACCTGGTGCGCGAGGGCGATA‑TGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCAttagatcgg                                                            >  3:297561/1‑141 (MQ=255)                                                                                gggCTACATCAAGCACCTGGTGCGCGAGGGCGATA‑TGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCAt                                                                    <  4:297561/140‑1 (MQ=255)                                                                                 ggCTACATCAAGCACCTGGTGCGCGAGGGCGATA‑TGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGt                                                           >  3:159222/1‑148 (MQ=255)                                                                                 ggCTACATCAAGCACCTGGTGCGCGAGGGCGATA‑TGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTc                                                          <  1:125684/149‑1 (MQ=255)                                                                                   cTACATCAAGCACCTGGTGCGCGAGGGCGATA‑TGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCAt                                                        >  7:161807/1‑149 (MQ=255)                                                                                    tACATCAAGCACCTGGTGCGCGAGGGCGATA‑TGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCAtt                                                       >  5:355512/1‑149 (MQ=255)                                                                                      cATCAAGCACCTGGTGCGCGAGGGCGATA‑TGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTc                                                          <  2:249492/144‑1 (MQ=255)                                                                                      cATCAAGCACCTGGTGCGCGAGGGCGATA‑TGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTc                                                          >  1:249492/1‑144 (MQ=255)                                                                                      cATCAAGCACCTGGTGCGAGAGGGCGATA‑TGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCTACATGGACAATGCGCTGATGAAAGATATAGATATCGCCAAGGGTAAGTCTCATTAGTCATTGCCGAGTTCAtttt                                                     >  4:230344/1‑148 (MQ=255)                                                                                       aTCAAGCACCTGGTGCGCGAGGGCGATA‑TGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAg                                                                        <  7:170029/129‑1 (MQ=255)                                                                                       aTCAAGCACCTGGTGCGCGAGGGCGATA‑TGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAg                                                                        >  8:170029/1‑129 (MQ=255)                                                                                       aTCAAGCACCTGGTGCGCGAGGGCGATA‑TGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTtgt                                                    >  7:249962/1‑149 (MQ=255)                                                                                       aTCAAGCACCTGGTGCGCGAGGGCGATA‑TGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTtgt                                                    >  5:32314/1‑149 (MQ=255)                                                                                          aaGCACCTGGTGCGCGAGGGCGATA‑TGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTgg                                                  >  8:352580/1‑148 (MQ=255)                                                                                                cTGGTGCGCGAGGGCGATA‑TGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCgg                                                            <  6:95235/132‑1 (MQ=255)                                                                                                cTGGTGCGCGAGGGCGATA‑TGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCgg                                                            >  5:95235/1‑132 (MQ=255)                                                                                                cTGGTGCGCGAGGGCGATA‑TGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGc                                           >  2:94672/1‑149 (MQ=255)                                                                                                   gTGCGCGAGGGCGATA‑TGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCtt                                        <  4:41095/149‑1 (MQ=255)                                                                                                     gcgcGAGGGCGATA‑TGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGc                                      >  6:71925/1‑149 (MQ=255)                                                                                                      cgcgAGGGCGATA‑TGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGc                                      <  4:159222/148‑1 (MQ=255)                                                                                                       gcgAGGGCGATA‑TGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGt                                    >  4:295346/1‑149 (MQ=255)                                                                                                        cgAGGGCGATA‑TGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTc                                   <  7:41959/149‑1 (MQ=255)                                                                                                           gggCGATA‑TGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTcgc                                 >  4:31775/1‑148 (MQ=255)                                                                                                            ggCGATA‑TGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTcgcg                                >  8:98997/1‑148 (MQ=255)                                                                                                                ata‑tGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATgg                           <  2:308790/149‑1 (MQ=255)                                                                                                                    tGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGc                         <  5:99819/148‑1 (MQ=255)                                                                                                                    tGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAAAGGGcc                        <  7:352580/149‑2 (MQ=255)                                                                                                                     ggaggaCGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACAAGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGACCGGCCTTGCGTCTCGAATGGGCTg                       >  3:99729/1‑149 (MQ=255)                                                                                                                      gaggaCGACCATAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTg                       <  5:2005/148‑1 (MQ=255)                                                                                                                                aTAAATTCCTGGGGGCAGCCGTT‑TACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTtg                                                     <  1:361575/108‑1 (MQ=255)                                                                                                                                aTAAATTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTtg                                                     >  2:361575/1‑108 (MQ=255)                                                                                                                                    aTTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCTTTGCAGcccattc        <  8:161807/149‑5 (MQ=39)                                                                                                                                    aTTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCTTTGCAGcccattc        <  1:94672/149‑5 (MQ=39)                                                                                                                                    aTTCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGc        <  6:32314/149‑1 (MQ=255)                                                                                                                                     ttCCTGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATTGGCTGGCAAGCAGCCCCGGCg       >  1:265407/1‑149 (MQ=255)                                                                                                                                        cTGGGGGCAGCCT‑ATTACCTGTGCGACTTGGCCAAGGCTGTTATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGAt     <  3:225133/148‑1 (MQ=38)                                                                                                                                         tGGGGGCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGAttt   >  5:123138/1‑149 (MQ=255)                                                                                                                                           ggggCAGCCG‑ATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTc  <  5:367684/148‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                         TTAACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATA‑TGGAGGACGACCATAAATTCCTGGGGGCAGCCG‑A‑TACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTC  >  NC_002947/335977‑336260

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