Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,389,596	T→G	8.5%	Q28P (CAA→CCA) ‡	PP_2991 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,389,596	0	T	G	8.5%	177.1 / 7.8	59	Q28P (CAA→CCA) ‡	PP_2991	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (28/26);  new base G (2/3);  total (30/29)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01

GCAGCTCTGCTGCAATCACCGACAAGGGCGGGTTCGAGCGATAACGGTCGTA‑CAGGCCGCCTTTGACGACCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGGGTGGGAACAACTTGAGGGTTCAGCGAGCTACGTTGTCGCCAGAGAGCACCGCCATCACCGCTTCAGCTACGGC  >  NC_002947/3389450‑3389741                                                                                                                                                    |                                                                                                                                                  gCAGCTCTGCTGCAATCACCGACAAGGGCGGGTTCGAGCGATAACGGTCGTAGCAGGCCGCCTTTGACGACCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCttg                                                                                                                                                  >  1:84638/1‑149 (MQ=255)      tctGCTGCAATCACCGACAAGGGCGGGTTCGAGCGATAACGGTCGTA‑CAGGCCGCCTTTGACGACCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCAAGGGTTGGc                                                                                                                                            <  4:191982/149‑1 (MQ=255)          ctgcAATCACCGACAAGGGCGGGTTCGAGCGATAACGGTCGTA‑CAGGCCGCCTTTGACGACCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAAtt                                                                                                                                        <  2:84638/149‑1 (MQ=255)          ctgcAATCACCGACAAGGGCGGGTTCGAGCGATAACGGTCGTA‑CAGGCCGCCTTTGACGACCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCAAGGGTTGGcg                                                                                                                                           <  1:180126/146‑2 (MQ=255)              aaTCACCGACAAGGGCGGGTTCGAGCGATAACGGTCGTA‑CAGGCCGCCTTTGACGACCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGc                                                                                                                                    <  5:11018/149‑1 (MQ=255)              aaTCACCGACAAGGGCGGGTTCGAGCGATAACGGTCGTA‑CAGGCCGCCTTTGACGACCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGc                                                                                                                                    <  8:216043/149‑1 (MQ=255)                tCACCGACAAGGGCGGGTTCGAGCGATAACGGTCGTA‑CAGGCCGCCTTTGACGACCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCt                                                                                                                                   >  3:316641/1‑148 (MQ=255)                tCACCGACAAGGGCGGGTTCGAGCGATAACGGTCGTA‑CAGGCCGCCTTTGACGACCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCAAGGGTTGGcg                                                                                                                                           >  6:336535/1‑139 (MQ=255)                tCACCGACAAGGGCGGGTTCGAGCGATAACGGTCGTA‑CAGGCCGCCTTTGACGACCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCAAGGGTTGGcg                                                                                                                                           <  5:336535/140‑2 (MQ=255)                   ccGACAAGGGCGGGTTCGAGCGATAACGGTCGTA‑CAGGCCGCCTTTGACGACCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGaa                                                                                                                                >  5:127441/1‑148 (MQ=255)                    cGACAAGGGCGGGTTCGTGCGATAACGGTCGTA‑CAGGCCTCCTTTGACGACCGGCAGCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCGGCCAACCCTTGGTTGGCAATTGAGCTGAACa                                                                                                                              <  5:152711/149‑1 (MQ=255)                            gCGGGTTCGAGCGATAACGGTCGTA‑CAGGCCGCCTTTGACGACCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCt                                                                                                                      <  3:310710/149‑1 (MQ=255)                                   cGAGCGATAACGGTCGTA‑CAGGCCGCCTTTGACGACCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAAt                                                                                                                <  1:277447/148‑1 (MQ=255)                                    gAGCGATAACGGTCGTA‑CAGGCCGCCTTTGACGACCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGc                                                                                                              >  7:70006/1‑149 (MQ=255)                                     aGCGATAACGGTCGTA‑CAGGCCGCCTTTGACGACCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCg                                                                                                             >  5:71622/1‑149 (MQ=255)                                              ggTCGTA‑CAGGCCGCCTTTGACGACCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGt                                                                                                                    <  1:247833/133‑1 (MQ=255)                                              ggTCGTA‑CAGGCCGCCTTTGACGACCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGt                                                                                                                    >  2:247833/1‑133 (MQ=255)                                              ggTCGTA‑CAGGCCGCCTTTGACGACCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGc                                                                                                     >  5:84791/1‑148 (MQ=255)                                               gTCGTA‑CAGGCCGCCTTTGACGACCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCaa                                                                                                   <  5:8833/149‑1 (MQ=255)                                                   tA‑CAGGCCGCCTTTGACGACCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGttt                                                                                               >  1:39860/1‑149 (MQ=255)                                                   tA‑CAGGCCGCCTTTGACGACCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGtt                                                                                                >  8:195310/1‑148 (MQ=255)                                                    a‑CAGGCCGCCTTTGACGACCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATg                                                                                                               >  2:83568/1‑132 (MQ=255)                                                    a‑CAGGCCGCCTTTGACGACCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATg                                                                                                               <  1:83568/132‑1 (MQ=255)                                                      cAGGCCGCCTTTGACGACCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTc                                                                                              <  4:234741/148‑1 (MQ=255)                                                       aGGCCGCCTTTGACGACCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCaa                                                                                            <  8:70006/149‑1 (MQ=255)                                                         gccgccTTTGACGACCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATc                                                                                          <  7:195310/149‑1 (MQ=255)                                                           cgccTTTGACGACCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGtt                                                                                                                   <  8:239720/122‑1 (MQ=255)                                                           cgccTTTGACGACCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGtt                                                                                                                   >  7:239720/1‑122 (MQ=255)                                                             ccTTTGACGACCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCttggtt                                                                                                                                               >  6:318348/1‑92 (MQ=255)                                                                   acgacCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCATGCAAGTTTCAATCGGTCTTtctc                                                                                >  7:269159/1‑149 (MQ=255)                                                                        cGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGAGGTGCTTGCAAGTTTCAATCGGTCTTTCTCcgatg                                                                           >  7:188605/1‑146 (MQ=255)                                                                         ggCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCggg                                                                          <  2:39860/149‑1 (MQ=255)                                                                             ccGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGGGTggg                                                                      <  6:127441/149‑1 (MQ=255)                                                                                      cTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGGGTGGGAACAACtt                                                              >  2:21095/1‑148 (MQ=255)                                                                                      cTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGGGTGGGAACAACtt                                                              >  3:305577/1‑148 (MQ=255)                                                                                        gCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTACAAGTTTCAATCGGTCTTTCTCCGCGGGTGGGAACAACTTGAg                                                           >  6:125905/1‑149 (MQ=255)                                                                                         cGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGGGTGGGAACAACTTGAgg                                                          <  4:316641/149‑1 (MQ=255)                                                                                         cGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGGGTGGGAACAACTTGAg                                                           <  6:84791/148‑1 (MQ=255)                                                                                          gATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGGGTGGGAACAACTTGAggg                                                         >  6:79127/1‑149 (MQ=255)                                                                                           aTTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGGGTGGGAACAACTTGAGGGt                                                        >  3:165272/1‑149 (MQ=255)                                                                                                gATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGGGTGGGAACAACTTGAGGGTTCAGc                                                   >  6:58627/1‑149 (MQ=255)                                                                                                gATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGGGTGGGAACAACTTGAGGGTTCAGc                                                   >  2:103923/1‑149 (MQ=255)                                                                                                       ttGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGGGTGGGAACAACTTGAGGGTTCAGCTAGCTAc                                            >  7:120207/1‑149 (MQ=255)                                                                                                         ggTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGGGTGGGAACAACTTGAGGGTTCAGCGAGCTACGt                                          >  1:334488/1‑149 (MQ=255)                                                                                                              cAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTACAAGTTTCAATCGGTCTTTCTCCGCGGGTGGGAACAACTTGAGGGTTCAGCGAGCTACGTTGTc                                      <  5:125905/148‑1 (MQ=255)                                                                                                                 tCGACGTGAGCATTTCAGGCGTCCCCGCCATCCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGGGTGGGAACAACTTGAGGGTTCAGCGAGCTACGTTGTCGc                                    >  1:28926/1‑147 (MQ=255)                                                                                                                     cGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGGGTGGGAACAACATGAGGGTTCAGCGAGCTACGTTGTCGCCaga                                >  3:176746/1‑147 (MQ=255)                                                                                                                       tGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGGGTGGGAACAACTTGAGGGTTCAGCGAGCTACGTTGTCGCCAGAGAGc                            <  6:71622/149‑1 (MQ=255)                                                                                                                          gCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGGGTGGGAACAACTTGAGGGTTCAGc                                                   <  8:136205/123‑1 (MQ=255)                                                                                                                            aTTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGGGTGGGAACAACCTGAGGGTTCAGCGAGCTACGTTGTCGCCAGAGAGCACCGc                       <  4:176746/149‑1 (MQ=255)                                                                                                                             tttCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGGGTGGGAACAACTTGAGGGTTCAGCGAGCTACGTTGTCGCCAGAGAGCACCGc                       >  4:177601/1‑148 (MQ=255)                                                                                                                                cAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCTGGTGGGAACAACTTGAGGGTTCAGCGAGCTACGTTGTCGCCAGAGAGCACCGCCAt                    <  5:154952/148‑1 (MQ=255)                                                                                                                                 aGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGGGTGGGAACAACTTGAGGGTTCAGCGAGCTACGTTGTCGCCAGAGAGCACCGCCATCa                  >  5:298592/1‑149 (MQ=255)                                                                                                                                      tCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGGGTGGGAACAACTTg                                                             <  2:66680/101‑1 (MQ=255)                                                                                                                                      tCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGGGTGGGAACAACTTg                                                             >  1:66680/1‑101 (MQ=255)                                                                                                                                        cccGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGGGTGGGAACAACTTGAGGGTTCAGCGAGCTACGTTGTCGCCagaga                              >  3:266643/1‑130 (MQ=255)                                                                                                                                        cccGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGGGTGGGAACAACTTGAGGGTTCAGCGAGCTACGTTGTCGCCagaga                              <  4:266643/130‑1 (MQ=255)                                                                                                                                        cccGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGGGTGGGAACAACTTGAGGGTTCAGCGAGCTACGTTGTCGCCAGAGAGCACCGCCATCACCGCtt            <  2:28926/148‑1 (MQ=255)                                                                                                                                         ttgCCAACCAAGGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGGGTGGGAACAACTTGAGGGTTCAGCGAGCTACGTTGTCGCCAGAGAGCACCGCCATCACCGCtt            >  5:81922/3‑147 (MQ=255)                                                                                                                                             cAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGGGTGGGAACAACTTGAGGGTTCAGCGACCTACGTTGTCGCCAGAGAGCACCGCCATCACCGCTTCAGCTa      >  1:114988/1‑149 (MQ=255)                                                                                                                                                 ccTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCATGCAAGTTTCAATCGGTCTTTCTCCGCGGGTGGGAACAACTTGAGGGTTCAGCGAGCTACGTTGTCGCCAGAGAGCACCGCCATCACCGCTTCAGCTACGGc  <  8:269159/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                  GCAGCTCTGCTGCAATCACCGACAAGGGCGGGTTCGAGCGATAACGGTCGTA‑CAGGCCGCCTTTGACGACCGGCACCGGGTGGCCTGCGATTTCGATTGTGTTGGTTTCAGTCGACGTGAGCATTTCAGGCGTCCCCGCCAACCCTTGGTTGGCAATTGAGCTGAACAGCAATCCTGTTAATGCGGTGCTTGCAAGTTTCAATCGGTCTTTCTCCGCGGGTGGGAACAACTTGAGGGTTCAGCGAGCTACGTTGTCGCCAGAGAGCACCGCCATCACCGCTTCAGCTACGGC  >  NC_002947/3389450‑3389741

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