Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,604,205	C→G	9.5%	P12R (CCG→CGG) ‡	PP_4076 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,604,205	0	C	G	9.5%	183.9 / 10.4	63	P12R (CCG→CGG) ‡	PP_4076	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (28/29);  new base G (3/3);  total (31/32)
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.03e-01

GACGTGGCGGGTGATCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTG‑CTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACTTCAACGGGCAT  >  NC_002947/4604061‑4604350                                                                                                                                                 |                                                                                                                                                   gACGTGGCGGGTGATCCTGGCTGGTTCCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCGTGGTAGGCCGGCTCGAGTCGCAATGAGTGGTTTCTTCGATGGCTTGATGGCGGACGGGTcc                                                                                                                                                <  6:321300/149‑1 (MQ=18)       gCGGGTGATACGGGATGGTACCAGGTCTCGCTAGCCGTTCGCCCGCACTTAAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACGGGTCCGCCATc                                                                                                                                          <  8:168582/149‑1 (MQ=18)        cGGGTGAGCCCGGCTCGGACCAGGCCACGCTCGCCGTGCGCCCGCACGTCAAGTACCTGGGCGCGAACTTCGAAATTACCCTGCTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATcc                                                                                                                                         <  8:290140/149‑1 (MQ=255)         gggTGATCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCg                                                                                                                                        >  8:9866/1‑149 (MQ=255)          ggTGATCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACGGGTCCGc                                                                                                                                              <  3:235530/142‑1 (MQ=255)          ggTGATCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACGGGTCCGc                                                                                                                                              >  4:235530/1‑142 (MQ=255)               tCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGc                                                                                                                                  <  4:311059/149‑1 (MQ=255)                  gggCTGGTACTTGGTCTGGCTCGCCGTGCGGGGGGAGTTGTAGTTCATGGGGGCGAACTTGGTATTTTGGGTGGTTGGGCGGCTGGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGAGCCGTCCGCGATCCGGGAGGCGtc                                                                                                                               <  4:46221/149‑1 (MQ=255)                   ggCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGtct                                                                                                                              <  8:52898/149‑1 (MQ=255)                    gCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGtctc                                                                                                                             <  2:220804/149‑1 (MQ=255)                          aCCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAg                                                                                                                        >  8:127789/1‑148 (MQ=255)                           ccAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAgg                                                                                                                       >  3:127187/1‑148 (MQ=255)                                tctcGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAgcg                                                                                                                 >  2:145869/1‑149 (MQ=255)                                tctcGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAgcg                                                                                                                 <  4:32267/149‑1 (MQ=255)                                  tcgctcgcCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGCGAGGCGTCTCGCCAGGAAAGCGCt                                                                                                               >  2:152579/1‑149 (MQ=255)                                    gctcgcCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTc                                                                                                             >  1:134820/1‑149 (MQ=255)                                      tcgcCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGa                                                                                                                                     >  1:307535/1‑123 (MQ=255)                                      tcgcCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGa                                                                                                                                     <  2:307535/123‑1 (MQ=255)                                      tcgcCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCg                                                                                                            >  6:136349/1‑148 (MQ=255)                                      tcgcCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCg                                                                                                            >  2:207873/1‑148 (MQ=255)                                       cgcCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGca                                                                                                          >  4:128355/1‑149 (MQ=255)                                        gcCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGCGAGGCGTCTCGCCAGGAAAGCGCTTCGc                                                                                                           <  1:152579/147‑1 (MQ=255)                                         ccGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGa                                                                                                                      <  3:156069/135‑1 (MQ=255)                                         ccGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGa                                                                                                                      >  4:156069/1‑135 (MQ=255)                                          cGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACaa                                                                                                       <  2:159064/149‑1 (MQ=255)                                               gcCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTcg                                                                                                  >  2:223300/1‑149 (MQ=255)                                                 ccGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAg                                                                                                                   <  7:221845/130‑1 (MQ=255)                                                 ccGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAg                                                                                                                   >  8:221845/1‑130 (MQ=255)                                                 ccGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTcgcg                                                                                                >  2:160840/1‑149 (MQ=255)                                                 ccGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTcgcg                                                                                                >  2:190189/1‑149 (MQ=255)                                                  cGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATCGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCgg                                                                                               >  4:35881/1‑149 (MQ=255)                                                        tCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATc                                                                                          >  7:4079/1‑148 (MQ=255)                                                         cAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATc                                                                                          >  1:319628/1‑147 (MQ=255)                                                           aGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGt                                                                                      <  1:160840/149‑1 (MQ=255)                                                            gTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGt                                                                                      >  6:25886/1‑148 (MQ=255)                                                            gTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGt                                                                                      <  5:25886/148‑1 (MQ=255)                                                               cATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCa                                                                                  <  4:127187/149‑1 (MQ=255)                                                                aTGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCAc                                                                                 <  8:129221/149‑1 (MQ=255)                                                                        gAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAAc                                                                         >  8:172107/1‑149 (MQ=255)                                                                            ttCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGaa                                                                           <  4:203885/143‑1 (MQ=255)                                                                            ttCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGaa                                                                           >  3:203885/1‑143 (MQ=255)                                                                              cGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTg                                                                     <  8:4079/147‑1 (MQ=255)                                                                                      cccTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTc                                                           >  1:137527/1‑149 (MQ=255)                                                                                       ccTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACTTTGCTGAACTCTCg                                                          <  2:134820/149‑1 (MQ=255)                                                                                            tAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAg                                                     <  1:207873/149‑1 (MQ=255)                                                                                              ggCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAgg                                                    <  7:9866/148‑1 (MQ=255)                                                                                                 cGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTG‑c                                               <  7:127789/149‑1 (MQ=255)                                                                                                       gACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTG‑CTCCCag                                         <  2:319628/149‑1 (MQ=255)                                                                                                       gACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGAGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTG‑CTcc                                            >  8:17028/1‑146 (MQ=255)                                                                                                                tGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTG‑CTCCCAGAGCAGcccc                                >  7:21840/1‑149 (MQ=255)                                                                                                                    tGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTG‑CTCCCAGAGCAGCCCCGAGc                            >  6:245712/1‑149 (MQ=255)                                                                                                                         tcttcGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTG‑CTCCCAGAGCAGCCCCGAGCTGGGc                       >  3:250511/1‑149 (MQ=255)                                                                                                                          cttcGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTG‑CTCCCAGAGCAGCCCCGAGCTGGGc                       <  5:136349/148‑1 (MQ=255)                                                                                                                              gATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCCTCCCAGAGCAGCCCCGAGCTGGGCCTgc                   <  1:190189/149‑1 (MQ=255)                                                                                                                               aTCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTG‑CTCCCAGAGCAGCCCCGAGCTGGGCCTgcgc                 >  2:197818/1‑149 (MQ=255)                                                                                                                                  gCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTG‑CTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGAc              <  1:145869/149‑1 (MQ=255)                                                                                                                                      gATGGCGGACGGGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTG‑CTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACTTc           >  2:60782/1‑148 (MQ=255)                                                                                                                                      gATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTG‑CTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACTTc           >  6:53924/1‑148 (MQ=255)                                                                                                                                       aTGGCGGACGGGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTG‑CTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACTTCaa         >  8:147937/1‑149 (MQ=255)                                                                                                                                        tGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTG‑CTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACTTCAAc        <  1:197818/149‑1 (MQ=255)                                                                                                                                        tGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTA‑CTCCCAGAGCAGCCCGGAGCTGGGCCTGCGCGACTTCAAc        <  5:194130/149‑1 (MQ=255)                                                                                                                                           cGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTG‑CTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACTTCAACggg     <  5:245712/149‑1 (MQ=255)                                                                                                                                              aCCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTG‑CTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACTTCAACGGGCAt  <  2:137527/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                   GACGTGGCGGGTGATCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTG‑CTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACTTCAACGGGCAT  >  NC_002947/4604061‑4604350

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