Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,238,925	G→T	8.1%	intergenic (+233/‑48)	PP_5701 → / → PP_4616	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,238,925	0	G	T	8.1%	214.3 / 5.0	69	intergenic (+233/‑48)	PP_5701/PP_4616	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (33/30);  new base T (4/2);  total (37/32)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGCCGCCACCCCATTCGCGGGCAAGCCTGTGCCTACAACTGTTTCGGGCAAATGGTTCCTGGCAAAGGCGCTCCCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTA‑AACCGCAACACAACCAGCAAAAAGGCATAGAGGTGATGAATGGCCGTCGAAGTGGTGTACCGCAGCAGCCGCGACCCGGAGCGCTTGTTCATGGATAAGGCCGAAGCAGACCGTCACGACAAGATGCTCGAACT  >  NC_002947/5238779‑5239067                                                                                                                                                   |                                                                                                                                                cAGCCGCCACCCCATTCGCGGGCAAGCCTGTGCCTACAACTGTTTCGGGCAAATGGTTCCTGGCAAAGGCGCTCCCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGa                                                                                                                                                >  3:338100/1‑148 (MQ=255)  aGCCGCCACCCCATTCGCGGGCAAGCCTGTGCCTACAACTGTTTCGGGCAAATGGTTCCTGGCAAAGGCGCTCCCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATa                                                                                                                                              <  2:47893/149‑1 (MQ=255)         aCCCCATTCGCGGGCAAGCCTGTGCCTACAACTGTTTCGGGCAAATGGTTCCTGGCAAAGGCGCTCCCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTa‑aa                                                                                                                                      >  3:137785/1‑149 (MQ=255)          ccccATTCGCGGGCAAGCCTGTGCCTACAACTGTTTCGGGCAAATGGTTCCTGGCAAAGGCGCTCCCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATCTATATACAGTa‑atg                                                                                                                                     >  7:279656/1‑147 (MQ=255)             cATTCGCGGGCAAGCCTGTGCCTACAACTGTTTCGGGCAAATGGTTCCTGGCAAAGGCGCTCCCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTA‑AACCGc                                                                                                                                  <  5:76587/149‑1 (MQ=255)              aTTCGCGGTCAAGGCTGTGCCTACAACTGTTTCGTGCAAATGGTTCCTGGCAAAGGCGCTCCCGCAGCGCTCTGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATCTATATACAGTa‑atggcaa                                                                                                                                 <  5:258638/149‑7 (MQ=255)                    gggCAAGCCTGTGCCTACAACTGTTTCGGGCAAATGGTTCCTGGCAAAGGCGCTCCCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTA‑AACCGCAACACaa                                                                                                                           >  3:19011/1‑149 (MQ=255)                     ggCAAGCCTGTGCCTACAACTGTTTCGGGCAAATGGTTCCTGGCAAAGGCGCTCCCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTA‑AACCGCAACACAAc                                                                                                                          >  2:216842/1‑149 (MQ=255)                     ggCAAGCCGGTGCCTACAACTGTTTCGGGCAAATGGTTCCTGGCAAAGGCGCTCCCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTA‑AACCGCAACACAAc                                                                                                                          <  8:36480/149‑1 (MQ=255)                              gtgCCTACAACTGTTTCGGGCAAATGGTTCCTGGCAAAGGCGCTCCCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTA‑AACCGCAACACAACCa                                                                                                                        >  4:343464/1‑142 (MQ=255)                              gtgCCTACAACTGTTTCGGGCAAATGGTTCCTGGCAAAGGCGCTCCCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTA‑AACCGCAACACAACCa                                                                                                                        <  3:343464/142‑1 (MQ=255)                                  cTACAACTGTTTCGGGCAAATGGTTCCTGGCAAAGGCGCTCCCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATa                                                                                                                                              >  5:223537/1‑117 (MQ=255)                                  cTACAACTGTTTCGGGCAAATGGTTCCTGGCAAAGGCGCTCCCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATa                                                                                                                                              <  6:223537/117‑1 (MQ=255)                                   tACAACTGTTTCGGGCAAATGGTTCCTGGCAAAGGCGCTCCCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATa                                                                                                                                              <  5:173493/116‑1 (MQ=255)                                   tACAACTGTTTCGGGCAAATGGTTCCTGGCAAAGGCGCTCCCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATa                                                                                                                                              >  6:173493/1‑116 (MQ=255)                                    aCAACTGTTTCGGGCAAATGGTTCCTGGCAAAGGCGCTCCCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTA‑AACCGCAACACAACCAGCAAAAAGGCATa                                                                                                           >  7:313377/1‑149 (MQ=255)                                     cAACTGTTTCGGGCAAATGGTTCCTGGCAAAGGCGCTCCCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTA‑AACCGCAACACAACCAGCAAAAAGGCAt                                                                                                            <  7:314465/147‑1 (MQ=255)                                        cTGTTTCGGGCAAATGGTTCCTGGCAAAGGCGCTCCCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTA‑AACCGCAACACAACCAGCAAAAAGGCATAGAgg                                                                                                       >  3:64762/1‑149 (MQ=255)                                           tttCGGGCAAATGGTTCCTGGCAAAGGCGCTCCCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTA‑AACCGCAACACAACCAGCAAAAAGGCATAGAGGtga                                                                                                    >  5:346289/1‑149 (MQ=255)                                           tttCGGGCAAATGGTTCCTGGCAAAGGCGCTCCCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTA‑AACCGCAACACAACCAGCAAAAAGGCATAGAGGtg                                                                                                     <  8:313377/148‑1 (MQ=255)                                                 gCAAATGGTTCCTGGCAAAGGCGCTCCCGCAGCGCTCGGTCTGTGCCCCGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTA‑AACCGCAACACAACCAGCAAAAAGGCATAGAGGTGATGAATg                                                                                              >  7:265826/1‑149 (MQ=255)                                                 gCAAATGGTTCCTGGCAAAGGCGCTCCCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTA‑AACCGCAACACAACCAGCAAAAAGGCATAGAGGTGATGAATg                                                                                              <  4:64762/149‑1 (MQ=255)                                                 gCAAATGGTTCCTGGCAAAGGCGCTCCCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTA‑AACCGCAACACAACCAGCAAAAAGGCATAGAGGTGATGAATg                                                                                              <  6:346289/149‑1 (MQ=255)                                                  cAAATGGTTCCTGGCAAAGGCGCTCCCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTA‑AACCGCAACACAACCAGCAAAAAGGCATAGAGGTGATGAATgg                                                                                             >  3:322460/1‑149 (MQ=255)                                                    aaTGGTTCCTGGCAAAGGCGCTCCCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTA‑AACCGCAACACAACCAGCAAAAAGGCATAGAGGTGATGAATGGcc                                                                                           <  4:338100/149‑1 (MQ=255)                                                          tCCTGGCAAAGGCGCTCCCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATa                                                                                                                                              <  2:13106/93‑1 (MQ=255)                                                          tCCTGGCAAAGGCGCTCCCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATa                                                                                                                                              >  1:13106/1‑93 (MQ=255)                                                            cTGGCAAAGGCGCTCCCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTA‑AACCGCAACACAACCAGCAAAAAGGCATAGAGGTGATGAATGGCCGTCGAAgt                                                                                   >  4:302790/1‑149 (MQ=255)                                                             tGGCAAAGGCGCTCCCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTA‑AACCGCAACACAACCAGCAAAAAGGCATAGAGGTGATGAATGGCCGTCGAAgtg                                                                                  >  4:350163/1‑149 (MQ=255)                                                              ggCAAAGGCGCTCCCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTA‑AACCGCAACACAACCAGCAAAAAGGCATa                                                                                                           >  3:197974/1‑123 (MQ=255)                                                              ggCAAAGGCGCTCCCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTA‑AACCGCAACACAACCAGCAAAAAGGCATa                                                                                                           <  4:197974/123‑1 (MQ=255)                                                               gCAAAGGCGCTCCCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTA‑AACCGCAACACAACCAGCAAAAAGGCATAGAGGTGATGAATGGCCGTCGAAgtggt                                                                                <  5:290710/149‑1 (MQ=255)                                                               gCAAAGGCGCTCCCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTA‑AACCGCAACACAACCAGCAAAAAGGCATAGAGGTGATGAATGGCCGTCGAAgtgg                                                                                 >  2:105191/1‑148 (MQ=255)                                                               gCAAAGGCGCTCCCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTA‑AACCGCAACACAACCAGCAAAAAGGCATAGAGGTGATGAATGGCCGTCGAAgtgg                                                                                 <  1:105191/148‑1 (MQ=255)                                                                cAAAGGCGCTCCCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTGGGGGGGGGGAGTAGGGCCGCAACACAACCAGCAAAAAGGCATAGAGGTGATGAATGGCCGTCGAAgtggt                                                                                >  1:122276/1‑149 (MQ=255)                                                                 aaaGGCGCTCCCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTA‑AACCGCAACACAACCAGCAAAAAGGCATAGAGGTGATGAATGGCCGTCGAAgtg                                                                                  <  2:48701/145‑1 (MQ=255)                                                                 aaaGGCGCTCCCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTA‑AACCGCAACACAACCAGCAAAAAGGCATAGAGGTGATGAATGGCCGTCGAAgtg                                                                                  >  1:48701/1‑145 (MQ=255)                                                                   aGGCGCTCCCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTA‑AACCGCAACACAACCAGCAAAAAGGCATAGAGGTGATGAATGGCCGTCGAAGTGGTGTAc                                                                            >  7:314499/1‑149 (MQ=255)                                                                            cGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTA‑AACCGCAACACAACCAGCAAAAAGGCATAGAGGtgat                                                                                                   >  5:69730/1‑117 (MQ=255)                                                                            cGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTA‑AACCGCAACACAACCAGCAAAAAGGCATAGAGGtgat                                                                                                   <  6:69730/117‑1 (MQ=255)                                                                             gCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTA‑AACCGCAACACAACCAGCAAAAAGGCATAGAGGTGATGAATGGCCGTCGAAGTGGTGTACCGCAGCAGcc                                                                  <  7:116774/149‑1 (MQ=255)                                                                             gCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTA‑AACCGCAACACAACCAGCAAAAAGGCATAGAGGTGATGAATGGCCGTCGAAGTGGTGTACCGCAGCAGcc                                                                  >  8:116774/1‑149 (MQ=255)                                                                                 cgcTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTA‑AACCGc                                                                                                                                  <  6:308620/81‑1 (MQ=255)                                                                                 cgcTCGGTCTGTGCCAAGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTA‑AACCGc                                                                                                                                  >  5:308620/1‑81 (MQ=255)                                                                                            tgCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTA‑AACCGCAACACAACCAGCAAAAAGGCATAGAGGTGATGAATGGCCGTCGAAGTGGTGTACCGCAGCAGCCGCGACCCGGAgcgc                                                    <  3:350163/148‑1 (MQ=255)                                                                                              ccACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTA‑AACCGCAACACAACCAGCAAAAAGGCATAGAGGTGATGAATGGCCGTCGAAGTGGTGTACCGCAGCAGCCGCGACCCGGAGCGCtt                                                  <  3:302790/148‑1 (MQ=255)                                                                                              ccACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTA‑AACCGCAACACAACCAGCAAAAAGGCATAGAGGTGATGAATGGCCGTCGAAGTGGTGTACCGCAGCAGCCGCGACCCGGAGCGCtt                                                  <  8:314499/148‑1 (MQ=255)                                                                                               cACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTA‑AACCGCAACACAACCAGCAAAAAGGCATAGAGGTGATGAATGGCCGTCGAAGTGGTGTACCGCAGCAGCCGCGACCCGGAGCGCTTGt                                                <  4:322460/149‑1 (MQ=255)                                                                                                   cTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTA‑AACCGCAACACAACCAGCAAAAAGGCATAGAGGTGATGAATGGCCGTCGAAGTGGTGTACCGCAGCAGCCGCGACCCGGAGCGCTTGTTCAt                                            <  8:265826/149‑1 (MQ=255)                                                                                                    tACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACa                                                                                                                                            >  2:5471/1‑53 (MQ=255)                                                                                                    tACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACa                                                                                                                                            <  1:5471/53‑1 (MQ=255)                                                                                                    tACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTA‑AACCGCAACACAACCAGCAAAAAGGCATAGAGGtgat                                                                                                   >  8:181532/1‑93 (MQ=255)                                                                                                    tACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTA‑AACCGCAACACAACCAGCAAAAAGGCATAGAGGtgat                                                                                                   <  7:181532/93‑1 (MQ=255)                                                                                                                ttGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTA‑AACCGCAACACAACCAGCAAAAAGGCATAGAGGTGATGAATGGCCGTCGAAGTGGTGTACCGCAGCAGCCGCGACCCGGAGCGCTTGTTCATGGATAAGGCCGaa                               <  4:19011/149‑1 (MQ=255)                                                                                                                  gCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTA‑AACCGCAACACAACCAGCAAAAAGGCATAGAGGTGATGAATGGCCGTCGAAGTGGTGTACCGCAGCAGCCGCGACCCGGAGCGCTTGTTCATGGATAAGGCCGAAg                              <  6:191419/148‑1 (MQ=255)                                                                                                                  gCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTA‑AACCGCAACACAACCAGCAAAAAGGCATAGAGGTGATGAATGGCCGTCGAAGTGGTGTACCGCAGCAGCCGCGACCCGGAGCGCTTGTTCATGGATAAGGCCGAAg                              >  4:10637/1‑148 (MQ=255)                                                                                                                          accaGTATTGCCATTACTGTATATAGATACAGTA‑AACCGCAACACAACCAGCAAAAAGGCATAGAGGTGATGAATGGCCGTCGAAGTGGTGTACCGCAGCAGCCGCGACCCGGAGCGCTTGTTCATGGATAAGGCCGAAGCAGACCGt                      >  3:160123/1‑148 (MQ=255)                                                                                                                              gTATTGCCATTACTGTATATAGATACAGTA‑AACCGCAACACAACCAGCAAAAAGGCATAGAGGTGATGAATGGCCGTCGAAGTGGTGTACCGCAGCAGCCGCGACCCGGAGCGCTTGTTCATGGATAAGGCCGAAGCAGACCGTCACGa                 >  8:261531/1‑149 (MQ=255)                                                                                                                                   cggtttACTGTATCTATATACAGTA‑AACCGCAACACAACCAGCAAAAAGGCATAGAGGTGATGAATGGCCGTCGAAGTGGTGTACCGCAGCAGCCGCGACCCGGAGCGCTTGTTCATGGATAAGGCCGAAGCAGACCGTCACGAc                >  6:378601/5‑145 (MQ=255)                                                                                                                                    ccATTACTGTATATAGATACAGTA‑AACCGCAACACAACCAGCAAAAAGGCATAGAGGTGATGAATGGCCGTCGAAGTGGTGTACCGCAGCAGCCGCGACCCGGAGCGCTTGTTCATGGATAAGGCCGAAGCAGACCGTCACGACAAGAt           >  1:261700/1‑149 (MQ=255)                                                                                                                                     gtttACTGTATCTATATACAGTA‑AACCGCAACACAACCAGCAAAAAGGCATAGAGGTGATGAATGGCCGTCGAAGTGGTGTACCGCAGCAGCCGCGACCCGGAGCGCTTGTTCATGGATAAGGCCGAAGCAGACCGTCACGACAAGATg          >  3:318547/3‑149 (MQ=255)                                                                                                                                      tttACTGTATCTATATACAGTA‑AACCGCAACACAACCAGCAAAAAGGCATAGAGGTGATGAATGGCCGTCGAAGTGGTGTACCGCAGCAGCCGCGACCCGGAGCGCTTGTTCATGGATAAGGCCGAAGCAGACCGTCACGACAAGATGc         >  5:275164/2‑149 (MQ=255)                                                                                                                                      attattGTATATAGATACAGTA‑AACCGCAACACAACCAGCAAAAAGGCATAGAGGTGATGAATGGCCGTCGATGTTGTTTACCGCAGCAGCCGCGACCCTGAGCGCTTTTTCATGGATAAGGCCGAAGCAGACCGTCACGACAAGATGc         <  4:137785/149‑1 (MQ=255)                                                                                                                                      aTTACTGTATATAGATACAGTA‑AACCGCAACACAACCAGCAAAAAGGCATAGAGGTGATGAATGGCCGTCGAAGTGGTGTACCGCAGCAGCCGCGACCCGGAGCGCTTGTTCATGGATAAGGCCGAAGCAGACCGTCACGACAAGATGc         >  3:232521/1‑149 (MQ=255)                                                                                                                                           tGTATATAGATACAGTA‑AACCGCAACACAACCAGCAAAAAGGCATAGAGGTGATGAATGGCCGTCGAAGTGGTGTACCGCAGCAGCCGCGACCCGGAGCGCTTGTTCATGGATAAGGCCGAAGCAGACCGTCACGACAAGATGCTCGaa    >  7:290565/1‑149 (MQ=255)                                                                                                                                             tatctaTATACAGTA‑AACCGCAACACAACCAGCAAAAAGGCATAGAGGTGATGAATGGCCGTCGAAGTGGTGTACCGCAGCAGCCGCGACCCGGAGCGCTTGTTCATGGATAAGGCCGAAGCAGACCGTCACGACAAGATGCTCGaa    <  6:275164/147‑1 (MQ=255)                                                                                                                                             tatataGATACAGTA‑AACCGCAACACAACCAGCAAAAAGGCATAGAGGTGATGAATGGCCGTCGAAGTGGTGTACCGCAGCAGCCGCGACCCGGAGCGCTTGTTCATGGATAAGGCCGAAGCAGACCGTCACGACAAGATGCTCGAACt  >  1:350484/1‑149 (MQ=255)                                                                                                                                                  aGATACAGTA‑AACCGCAACACAACCAGCAAAAAGGCATAGAGGTGATGAATGGCCGTCGAAGTGGTGTACCGCAGCAGCCGCGACCCGGAGCGCTTGTTCATGGATAAGGCCGAAGCAGACCGt                      >  1:364120/1‑124 (MQ=255)                                                                                                                                                  aGATACAGTA‑AACCGCAACACAACCAGCAAAAAGGCATAGAGGTGATGAATGGCCGTCGAAGTGGTGTACCGCAGCAGCCGCGACCCGGAGCGCTTGTTCATGGATAAGGCCGAAGCAGACCGt                      <  2:364120/124‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                CAGCCGCCACCCCATTCGCGGGCAAGCCTGTGCCTACAACTGTTTCGGGCAAATGGTTCCTGGCAAAGGCGCTCCCGCAGCGCTCGGTCTGTGCCACGCTACAATTTCTGCTTGCGAACCAACCAGTATTGCCATTACTGTATATAGATACAGTA‑AACCGCAACACAACCAGCAAAAAGGCATAGAGGTGATGAATGGCCGTCGAAGTGGTGTACCGCAGCAGCCGCGACCCGGAGCGCTTGTTCATGGATAAGGCCGAAGCAGACCGTCACGACAAGATGCTCGAACT  >  NC_002947/5238779‑5239067

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