Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	262,757:1	+C	6.7%	coding (274/282 nt)	PP_0211 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	262,757	1	.	C	6.7%	235.1 / 11.8	60	coding (274/282 nt)	PP_0211	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (25/31);  new base C (2/2);  total (27/33)
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.79e-01

GGCGCGCCTGTCCGGAGGGATCAGCCTGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTGGGCAACGTGCCACCGATTAC  >  NC_002947/262612‑262904                                                                                                                                                     |                                                                                                                                                    ggcgCGCCTGTCCGGAGGGATCAGCCTGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCa                                                                                                                                                  <  4:43737/149‑1 (MQ=255)  gcgcgcCTGTCCGGAGGGATCAGCCTGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCaa                                                                                                                                                 <  6:98559/149‑1 (MQ=255)  gcgcgcCTGTCCGGAGGGATCAGCCTGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCaa                                                                                                                                                 <  4:144245/149‑1 (MQ=255)   cgcgcCTGTCCGGAGGGATCAGCCTGGTGCGCGACGATGTGAGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAg                                                                                                                                                <  5:68473/149‑1 (MQ=255)       ccTGTCCGGAGGGATCAGCCTGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTg                                                                                                                                              >  5:208281/1‑147 (MQ=255)        cTGTCCGGAGGGATCAGCCTGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAg                                                                                                                                           <  6:159453/149‑1 (MQ=255)        cTGGCCGGAGGGATCAGCCTGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGC‑T‑GACCAGGTCACGCAAGtaagc                                                                                                                                           <  2:366855/149‑5 (MQ=255)          gTCCGGAGGGATCAGCCTGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCt                                                                                                                                         >  6:130514/1‑149 (MQ=255)             cGGAGGGATCAGCCTGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGAt                                                                                                                                      <  6:208281/149‑1 (MQ=255)              ggAGGGATCAGCCTGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATc                                                                                                                                     >  5:369400/1‑149 (MQ=255)               gAGGGATCAGCCTGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGAt                                                                                                                                      <  1:280971/147‑1 (MQ=255)               gAGGGATCAGCCTGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGAt                                                                                                                                      >  2:280971/1‑147 (MQ=255)                aGGGATCAGCCTGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATccc                                                                                                                                   <  7:243829/149‑1 (MQ=255)                          cTGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGcc                                                                                                                         <  5:130514/149‑1 (MQ=255)                           tGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCg                                                                                                                        >  6:291549/1‑149 (MQ=255)                             gTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAAGTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCgg                                                                                                                       <  3:206983/148‑1 (MQ=255)                             gTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCaa                                                                                                                                                 <  6:119217/122‑1 (MQ=255)                             gTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCgg                                                                                                                       >  5:131201/1‑148 (MQ=255)                                cgcgACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAg                                                                                                                    <  6:10576/148‑1 (MQ=255)                                cgcgACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAg                                                                                                                    <  2:72903/148‑1 (MQ=255)                                    acgaTGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGTTCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGt                                                                                                                                               <  8:154158/117‑1 (MQ=255)                                    acgaTGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGt                                                                                                                                               >  7:154158/1‑117 (MQ=255)                                          tgCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGaa                                                                                                          <  2:176355/148‑1 (MQ=255)                                           gcgcGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCGTATGGGCATGAGCGCGGGATCTACCATTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGATGACCCCTTCAGTTGCCGGAAGGGGTGCAg                                                                                                            >  1:282801/1‑146 (MQ=255)                                             gcgTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAg                                                                                                      <  4:316748/149‑1 (MQ=255)                                             gcgTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAg                                                                                                      >  3:303355/1‑149 (MQ=255)                                                 ggtgCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGgac                                                                                                   >  5:47507/1‑148 (MQ=255)                                                      aaaTCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGc                                                                                             >  6:317961/1‑149 (MQ=255)                                                                 cAAGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACg                                                                                   >  4:141982/1‑148 (MQ=255)                                                                 cAAGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACg                                                                                   >  1:148356/1‑148 (MQ=255)                                                                   aGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAAc                                                                                <  6:47507/149‑1 (MQ=255)                                                                       tATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGc                                                                            >  2:44128/1‑149 (MQ=255)                                                                        aTGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGcc                                                                           >  5:207927/1‑149 (MQ=255)                                                                         tGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGcc                                                                           <  6:131201/148‑1 (MQ=255)                                                                           gCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGa                                                                         <  1:192760/148‑1 (MQ=255)                                                                           gCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGa                                                                         <  5:291549/148‑1 (MQ=255)                                                                           gCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAg                                                                        >  2:211170/1‑149 (MQ=255)                                                                            cGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGACCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAgg                                                                       <  2:282801/149‑1 (MQ=255)                                                                             gTACAACTGGTGTGCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGCCCGCCTAGCCAACGAACTGGCCGAggg                                                                      >  4:107636/1‑149 (MQ=255)                                                                              tACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGt                                                                     <  4:303355/149‑1 (MQ=255)                                                                                cAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTaa                                                                   >  4:162045/1‑149 (MQ=255)                                                                                           cAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGTGGACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCttt                                                         >  1:82292/1‑149 (MQ=255)                                                                                             gCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGTTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGc                                                      <  2:82292/149‑1 (MQ=255)                                                                                             gCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGTACTGGCCGAGGGTAACCTGCTCTTTGGc                                                      <  1:211170/149‑1 (MQ=255)                                                                                              cG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCg                                                     <  6:207927/149‑1 (MQ=255)                                                                                                      cATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACAcccc                                              >  6:304312/1‑149 (MQ=255)                                                                                                           gcgcGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAAc                                                                  >  7:219110/1‑124 (MQ=255)                                                                                                           gcgcGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAAc                                                                  <  8:219110/124‑1 (MQ=255)                                                                                                              cgGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTgg                                      <  1:44128/149‑1 (MQ=255)                                                                                                              cgGGAGCTACCCTTGCGCTTACTTGCGT‑GACCTGGGCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTgg                                      <  6:369400/149‑1 (MQ=255)                                                                                                                      aCCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAgg                                                                       >  7:137857/1‑108 (MQ=255)                                                                                                                      aCCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAgg                                                                       <  8:137857/108‑1 (MQ=255)                                                                                                                       cccTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGgc                             >  5:272181/1‑149 (MQ=255)                                                                                                                             gcttcaCTTGC‑T‑GACCAGGTCACGCAAGT‑AAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTg                      >  8:22392/6‑149 (MQ=255)                                                                                                                              gCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTg                      >  8:188966/1‑149 (MQ=255)                                                                                                                              cttcaCTTGC‑T‑GACCAGGTCACGCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTg                      >  5:322592/5‑149 (MQ=255)                                                                                                                                   cTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTGGGCaa                 >  8:39509/1‑149 (MQ=255)                                                                                                                                       cGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTGGGCAACGt              <  7:188966/148‑1 (MQ=255)                                                                                                                                        gT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTGGGCAACGTGc            <  5:304312/149‑1 (MQ=255)                                                                                                                                        ct‑GACCAGGTCACGCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTGGGCAACGt              <  6:322592/147‑1 (MQ=255)                                                                                                                                                    a‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTGGGCAACGTGCCACCGATTAc  <  7:39509/148‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    GGCGCGCCTGTCCGGAGGGATCAGCCTGGTGCGCGACGATGTGCGCGTGGTGCAAATCGAGTCGCAAGGGTATGGCGTACAACTGGTGTTCAGCG‑ATGGGCATGAGCGCGGGATCTACCCTTGGGCTTACTTGCGT‑GACCTGGTCA‑GCAAGTGAAGCTGATCCCTTCAGTTGCCGGAAGGGGTGCAGAAGCAGGACGACCGCCTAGCCAACGAACTGGCCGAGGGTAACCTGCTCTTTGGCGACACCCCGAGTGTGGGTACAGGGCGCAACTGGGCAACGTGCCACCGATTAC  >  NC_002947/262612‑262904

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