Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,883,834	G→A	6.5%	intergenic (‑167/‑784)	PP_1691 ← / → PP_1692	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,883,834	0	G	A	6.5%	177.5 / 2.8	61	intergenic (‑167/‑784)	PP_1691/PP_1692	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (28/29);  new base A (2/2);  total (30/31)
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.21e-01

TTTAAACACGGGGGTTATGATCCGTGCATGATAGAAACGTTCCAGTGACGCCTTTATTAAGGGACGAACGAACCTGACTGGAACAGGGCAGCCAGAGGGGGCGCAATCAAATGTAGGCGCTTGGGCTGACTGGTAGGTACTTCGGAAGGCAGGCAGCACGCCGGGGCCAAGGTTCTAGGCCGTCGAGGTGGGATC‑CTGATCAATCTGTCATGTGATTCGTGCTTGCGCAGCGTAAGGCTGGCGAGGGAATCATTGAATGGGCCGGTCCTTGTTTCAGCAGGGACGCTCTTGGGGG  >  NC_002947/1883686‑1883980                                                                                                                                                     |                                                                                                                                                    tttAAACACGGGGGTTATGATCCGTGCATGATAGAAACGTTCCAGTGACGCCTTTATTAAGGGACGAACGAACCTGACTGGAACAGGGCAGCCAGAGGGGGCGCAATCAAATGTAGGCGCTTGGGCTGACTGGTAGGTACTTCGGAagg                                                                                                                                                     >  8:148642/1‑149 (MQ=255)   tAAACACGGGGGTTATGATCCGTGCATGATAGAAACGTTCCAGTGACGCCTTTATTAAGGGACGAACGAACCTGACTGGAACAGGGCAGCCAGAGGGGGCGCAATCAAATGTAGGCGCTTGGGCTGACTGGTAGGTACTTCGGAaggc                                                                                                                                                    <  7:195867/148‑1 (MQ=255)    aaaCACGGGGGTTATGATCCGTGCATGATAGAAACGTTCCAGTGACGCCTTTATTAAGGGACGAACGAACCTGACTGGAACAGGGCAGCCAGAGGGGGCGCAATCAAATGTAGGCGCTTGGGCTGACTGGTAGGTACTTCGGAaggcag                                                                                                                                                  >  8:228852/1‑149 (MQ=255)      acacGGGGGTTATGATCCGTGCATGATAGAAACGTTCCAGTGACGCCTTTATTAAGGGACGAACGAACCTGACTGGAACAGGGCAGCCAGAGGGGGCGCAATCAAATGTAGGCGCTTGGGCTGACTGGTAGGTACTTCGGAaggcagg                                                                                                                                                 >  4:258711/1‑148 (MQ=255)               ttATGATCCGTGCATGATAGAAACGTTCCAGTGACGCCTTTATTAAGGGACGAACGAACCTGACTGGAACAGGGCAGCCAGAGGGGGCGCAATCAAATGTAGGCGCTTGGGCTGACTGGTAGGTACTTCGGAAGGCAGGCAGCACGCCg                                                                                                                                       <  2:85839/149‑1 (MQ=255)                    aTCCGTGCATGATAGAAACGTTCCAGTGACGCCTTTATTAAGGGACGAACGAACCTGACTGGAACAGGGCAGCCAGAGGGGGCGCAATCAAATGTAGGCGCTTGGGCTGACTGGTAGGTACTTCGGAAGGCAGGCAGCACGCCGGGGc                                                                                                                                   <  7:281377/148‑1 (MQ=255)                           cATGATAGAACCGTTCCAGTGACGCCTTTATTAAGGGACGAACGAACCTGACTGGAACAGGGCAGCCAGAGGGGGCGCAATCAAATGTAGGCGCTTGGGCTGACTGGTAGGTACTTCGGAAGGCAGGCAGCACGCCGGGGCCAAGGtt                                                                                                                            <  6:305450/148‑1 (MQ=255)                             tGATAGAAACGTTCCAGTGACGCATTTATTAAGGGACGAACGAACCTGACTGGAACAGGGCAGCCAGAGGGGGCGCAATCAAATGTAGGCGCTTGGGCTGACTGGTAGGTACTTCGGAAGGCAGGCAGCACGCCGGGGCCAAGGTTCTa                                                                                                                         <  8:219061/149‑1 (MQ=255)                                 aGAAACGTTCCAGTGACGCCTTTATTAAGGGACGAACGAACCTGACTGGAACAGGGCAGCCAGAGGGGGCGCAATCAAATGTAGGCGCTTGGGCTGACTGGTAGGTACTTCGGAAGGCAGGCAGCACGCCGGGGCCAAGGTTCTAGGcc                                                                                                                     >  1:48997/1‑149 (MQ=255)                                 aGAAACGTTCCAGTGACGCCTTTATTAAGGGACGAACGAACCTGACTGGAACAGGGCAGCCAGAGGGGGCGCAATCAAATGTAGGCGCTTGGGCTGACTGGTAGGTACTTCGGAAGGCAGGCAGCACGCCGGGGCCAAGGTTCTAGGc                                                                                                                      >  3:163795/1‑148 (MQ=255)                                 aGAAACGTTCCAGTGACGCCTTTATTAAGGGACGAACGAACCTGACTGGAACAGGGCAGCCAGAGGGGGCGCAATCAAATGTAGGCGCTTGGGCTGACTGGTAGGTACTTCGGAAGGCAGGCAGCACGCCGGGGCCAAGGTTCTAGGc                                                                                                                      <  4:163795/148‑1 (MQ=255)                                   aaaCGTTCCAGTGACGCCTTTATTAAGGGACGAACGAACCTGACTGGAACAGGGCAGCCAGAGGGGGCGCAATCAAATGTAGGCGCTTTGGCTGACTGGTAGGTACTTCGGAAGGCAGGCAGCACGCCGGGGCCAAGGTTCTAGGCCGt                                                                                                                   <  5:253676/149‑1 (MQ=255)                                   aaaCGTTCCAGTGACGCCTTTATTAAGGGACGAACGAACCTGACTGGAACAGGGCAGCCACAGGGGGCGCAATCAAATGTAGGCGCTTGTGCTGACTGGTAGGTACTTCGGAAGGCAGGCAGCACGCCGGCGCCAAGGATCTAGGCCga                                                                                                                   >  2:246384/1‑148 (MQ=255)                                     aCGTTCCAGTGACGCCTTTATTAAGGGACGAACGAACCTGACTGGAACAGGGCAGCCAGAGGGGGCGCAATCAAATGTAGGCGCTTGGGCTGACTGGTAGGTACTTCGGAAGGCAGGCAGCACGCCCGGGCCAAGGTTCTAGGCCGTCg                                                                                                                 <  7:228852/149‑1 (MQ=255)                                                aCGCCTTTATTAAGGGACGAACGAACCTGACTGGAACAGGGCAGCCAGAGGGGGCGCAATCAAATGTAGGCGCTTGGGCTGACTGGTAGGTACTTCGGAAGGCAGGCAGCACGCCGGGGCCAAGGTTCTAGGCCGTCGAGGTGGGATc‑c                                                                                                     >  2:313239/1‑149 (MQ=255)                                                    cTTTATTAAGGGACGAACGAACCTGACTGGAACAGGGCAGCCAGAGGGGGCGCAATCAAATGTAGGCGCTTGGGCTGACTGGTAGGTACTTCGGAAGGCAGGCAGCACGCCGGGGCCAAGGTTCTAGGCCGTCGAGGTGGGATC‑CTGAt                                                                                                 <  3:258711/149‑1 (MQ=255)                                                    cTTTATTAAGGGACGAACGAACCTGACTGGAACAGGGCAGCCAGAGGGGGCGCAATCAAATGTAGGCGCTTGGGCTGACTGGTAGGTACTTCGGAAGGCAGGCAGCACGCCGGGGCCAAGGTTCTAGGCCGTCGAGGTGGGATC‑CTGAt                                                                                                 >  4:102386/1‑149 (MQ=255)                                                    cTTTATTAAGGGACGAACGAACCTGACTGGAACAGGGCAGCCAGAGGGGGCGCAATCAAATGTAGGCGCTTGGGCTGACTGGTAGGTACTTCGGAAGGCAGGCAGCACGCCGGGGCCAAGGTTCTAGGCCGTCGAGGTGGGATC‑CTGAt                                                                                                 <  3:271109/149‑1 (MQ=255)                                                    cTTTATTAAGGGACGAACGAACCTGACTGGAACAGGGCAGCCAGAGGGGGCGCAATCAAATGTAGGCGCTTGGGCTGACTGGTAGGTACTTCGGAAGGCAGGCAGCACGCCGGGGCCAAGGTTCTAGGCCGTCGAGGTGGGATC‑CTGAt                                                                                                 <  1:92640/149‑1 (MQ=255)                                                               gACGAACGAACCTGACTGGAACAGGGCAGCCAGAGGGGGCGCAATCAAATGTAGGCGCTTGGGCTGACTGGTAGGTACTTCGGAAGGCAGGCAGCACGCCGGGGCCAAGGTTCTAGGCCGTCGAGGTGGGATC‑CTGATCAATCTGTCAt                                                                                      >  6:313725/1‑149 (MQ=255)                                                                   aacgaacCTGACTGGAACAGGGCAGCCAGAGGGGGCGCAATCAAATGTAGGCGCTTGGGCTGACTGGTAGGTACTTCGGAAGGCAGGCAGCACGCCGGGGCCAAGGTTCTAGGCCGTCGAGGTGGGATC‑CTGATCAATCTGTCATGTGa                                                                                  >  7:61570/1‑149 (MQ=255)                                                                                 gAACAGGGCAGCCAGAGGGGGCGCAATCAAATGTAGGCGCTTGGGCTGACTGGTAGGTACTTCGGAAGGCAGGCAGCACGCCGGGGCCAAGGTTCTAGGCCGTCGAGGTGGGATC‑CTGATCAATCTGTCATGTGATTCGTGCTTgcg                                                                      >  4:176714/1‑147 (MQ=255)                                                                                  aaCAGGGCAGCCAGAGGGGGCGCAATCAAATGTAGGCGCTTGGGCTGACTGGTAGGTACTTCGGGAGACAGGCAGCACGCCGGGTCCAAGGTTCTAGGCCGTCGAGGTGGGATC‑CTGATCAATCTGTCATGTGATTCGTGCTTGCGAAg                                                                   >  8:256557/1‑149 (MQ=255)                                                                                  aaCAGGGCAGCCAGAGGGGGCGCAATCAAATGTAGGCGCTTGGGCTGACTGGTAGGTACTTCGGGAGACAGGCAGCACGCCGGGGCCAAGGTTCTAGGCCGTCGAGGTGGGATC‑CTGATCAATCTGTCATGTGATTCGTGCTTGCGCAg                                                                   >  8:261641/1‑149 (MQ=255)                                                                                       ggCAGCCAGAGGGGGCGCAATCAAATGTAGGCGCTTGGGCTGACTGGTAGGTACTTCGGAAGGCAGGCAGCACGCCGGGGCCAAGGTTCTAGGCCGTCGAGGTGGGATC‑CTGATCAATCTGTCATGTGATTCGTGCTTGCGCAGCGt                                                                >  8:239868/1‑147 (MQ=255)                                                                                       ggCAGCCAGAGGGGGCGCAATCAAATGTAGGCGCTTGGGCTGACTGGTAGGTACTTCGGAAGGCAGGCAGCACGCCGGGGCCAAGGTTCTAGGCCGTCGAGGTGGGATC‑CTGATCAATCTGTCATGTGATTCGTGCTTGCGCAGCGt                                                                >  1:129815/1‑147 (MQ=255)                                                                                       ggCAGCCAGAGGGGGCGCAATCAAATGTAGGCGCTTGGGCTGACTGGTAGGTACTTCGGAAGGCAGGCAGCACGCCGGGGCCAAGGTTCTAGGCCGTCGAGGTGGGATC‑CTGATCAATCTGTCATGTGATTCGTGCTTGCGCAGCGTaa                                                              <  8:61570/149‑1 (MQ=255)                                                                                       ggCAGCCAGAGGGGGCGCAATCAAATGTAGGCGCTTGGGCTGACTGGTAGGTACTTCGGAAGGCAGGCAGCACGCCGGGGCCAAGGTTCTAGGCCGTCGAGGTGGGATC‑CTGATCAATCTGTCATGTGATTCGTGCTTGCGCAGCGTaa                                                              <  7:239868/149‑1 (MQ=255)                                                                                              agagGGGGCGCAATCAAATGTAGGCGCTTGGGCTGACTGGTAGGTACTTCGGAAGGCAGGCAGCACGCCGGGGCCAAGGTTCTAGGCCGTCGAGGTg                                                                                                            >  2:57800/1‑97 (MQ=255)                                                                                              agagGGGGCGCAATCAAATGTAGGCGCTTGGGCTGACTGGTAGGTACTTCGGAAGGCAGGCAGCACGCCGGGGCCAAGGTTCTAGGCCGTCGAGGTGGGATC‑CTGATCAATCTGTCATGTGATTCGTGCTTGCGCAGCGTAAGGCTGGc                                                       <  2:48997/149‑1 (MQ=255)                                                                                                agGGGGCGCAATCAAATGTAGGCGCTTGGGCTGACTGGTAGGTACTTCGGAAGGCAGGCAGCACGCCGGGGCCAAGGTTCTAGGCCGTCGAGGTGGGATC‑CTGATCAATCTGTCATGTGATTCGTGCTTGCGCAGCGTAAGGCTGGCGa                                                     <  7:148642/149‑1 (MQ=255)                                                                                                   gggCGCAATCAAATGTAGGCGCTTGGGCTGACTGGTAGGTACTTCGGAAGGCAGGCAGCACGCCGGTGCCAAGGTTCTAGGCCGTCGAGGTGGGATC‑CTGATCAATCTGTCATGTGATTCGTGCTTGCGCAGCGTAAGGCTGGCGAgg                                                   >  8:37176/1‑148 (MQ=255)                                                                                                        cAATCAAATGTAGGCGCTTGGGCTGACTGGTAGGTACTTCGGAAGGCAGGCAGCACGCCGGGGCCAAGGTTCTAGGCCGTCGAGGTGGGATC‑CTGATCAATCTGTCATGTGATTCGTGCTTGCGCAGCGTAAGGCTGGCGAGGGa                                                 <  3:28286/145‑1 (MQ=255)                                                                                                        cAATCAAATGTAGGCGCTTGGGCTGACTGGTAGGTACTTCGGAAGGCAGGCAGCACGCCGGGGCCAAGGTTCTAGGCCGTCGAGGTGGGATC‑CTGATCAATCTGTCATGTGATTCGTGCTTGCGCAGCGTAAGGCTGGCGAGGGa                                                 >  4:28286/1‑145 (MQ=255)                                                                                                        cAATCAAATGTAGGCGCTTGGGCTGACTGGTAGGTACTTCGGAAGGCAGGCAGCACGCCGGGGCCAAGGTTCTAGGCCGTCGAGGTGGGATC‑CTGATCAATCTGTCATGTGATTCGTGCTTGCGCAGCGTAAGGCTGGCGAGGGAATCa                                             <  1:313239/149‑1 (MQ=255)                                                                                                        cAATCAAATGTAGGCGCTTGGGCGGACTGGTAGGTACTTCGGAAGGCAGGCAGCACGCCGGGGCCAAGGTTCTAGGCCGTCGAGGTGGGATCTCTGATCAATCTGTCATGTGATTCGTGCTTGCGCAGCGTAAGGCTGGCGAGGGAATc                                              <  1:280716/149‑1 (MQ=255)                                                                                                         aaTCATATGTAGGCGCTTGGGCTGACTGGTAGCTACTTCGGAAGGCAGGCAGCACGCCTGGGCCAAGGTTCTAGGCCGTCGAGGTGGGATC‑GTGATCATTCTGTCATGTGAGTCGTGCTTGCGCAGCGTAAGGCTGGCGAGGGAATCAt                                            <  1:272384/149‑1 (MQ=255)                                                                                                         aaTCAAATGTAGGCGCTTGGGCTGACTGGTAGGTACTTCGGAAGGCAGGCAGCACGGCGGGGCCAAGGTTCTAGGCCGTCGAGGTGGGATC‑CTGATCAATCTGTCATGTGATTCGTGCTTGCGCAGCGTAAGGCTGGCGAGGGAATCAt                                            <  1:246384/149‑1 (MQ=255)                                                                                                           tCAAATGTAGGCGCTTGGGCTGACTGGTAGGTACTTCGGAAGGCAGGCAGCACGCCGGGGCCAAGGTTCTAGGCCGTCGAGGTGGGATC‑CTGATCAATCTGTCATGTGATTCGTGCTTGCGCAGCGTAAgg                                                            <  3:20482/131‑1 (MQ=255)                                                                                                           tCAAATGTAGGCGCTTGGGCTGACTGGTAGGTACTTCGGAAGGCAGGCAGCACGCCGGGGCCAAGGTTCTAGGCCGTCGAGGTGGGATC‑CTGATCAATCTGTCATGTGATTCGTGCTTGCGCAGCGTAAgg                                                            >  4:20482/1‑131 (MQ=255)                                                                                                           tCAAATGTAGGCGCTTGGGCTGACTGGTAGGTACTTCGGAAGGCAGGCAGCACGCCGGGGCCAAGGTTCTAGGCCGTCGAGGTGGGATC‑CTGATCAATCTGTCATGTGATTCGTGCTTGCGCAGCGTAAGGCTGGCGAGGGAATCAtt                                           >  2:337891/1‑148 (MQ=255)                                                                                                           tCAAATGTAGGCGCTTGGGCTGACTGGTAGGTACTTCGGAAGGCAGGCAGCACGCCGGGGCCAAGGTTCTAGGCCGTCGAGGTGGGATC‑CTGATCAATCTGTCATGTGATTCGTGCTTGCGCAGCGTAAGGCTGGCGAGGGAATCAtt                                           >  6:214838/1‑148 (MQ=255)                                                                                                               aTGTAGGCGCTTGGGCTGACTGGTAGGTACTTCGGAAGGCAGGCAGCACGCCGGGGCCAAGGTTCTAGGCCGTCGAGGTGGGATC‑CTGATCAATCTGTCATGTGATTCGTGCTTGCGCAGCGTAAGGCTGGCGAGGGAATCATTGAATg                                      >  6:366057/1‑149 (MQ=255)                                                                                                                tGTAGGCGCTTGGGCTGACTGGTAGGTACTTCGGAAGGCAGGCAGCACGCCGGGGCCAAGGTTCTAGGCCGTCGAGGTGGGATC‑CTGATCAATCTGTCATGTGATTCGTGCTTGCGCAGCGTAAGTCTGGCGAGGGAATCATTGAATgg                                     <  5:313725/149‑1 (MQ=255)                                                                                                                  tAGGCGCTTGGGCTGACTGGTAGGTACTTCGGAAGGCAGGCAGCACGCCGGGGCCAAGGTTCTAGTCCGTCGAGGTGGGATC‑CTGATCAATCTGTCATGTGATTCGTGCTTGCGCAGCGTAAGGCTGGCGAGGTAATCATTGAATGGGc                                   <  7:37176/149‑1 (MQ=255)                                                                                                                  tAGGCGCTTGGGCTGACTGGTAGGTACTTCGGAAGGCAGGCAGCACGCCGGGGCCAAGGTTCTAGGCCGTCGAGGTGGGATC‑CTGATCAATCTGTCATGTGATTCGTGCTTGCGCAGCGTAAGGCTGGCGAGGGAATCATTGAATGGGc                                   >  4:342457/1‑149 (MQ=255)                                                                                                                  tAGGCGCTTGGGCTGACTGGTAGGTACTTCGGAAGGCAGGCAGCACGCCGGGGCCAAGGTTCTAGGCCGTCGAGGTGGGATC‑CTGATCAATCTGTCATGTGATTCGTGCTTGCGCAGCGTAAGGCTGGCGAGGGAATCATTGAATGGGc                                   >  1:121959/1‑149 (MQ=255)                                                                                                                  tAGGCGCTTGGGCTGACTGGTAGGTACTTCGGAAGGCAGGCAGCACGCCGGGGCCAAGGTTCTAGGCCGTCGAGGTGGGATC‑CTGATCAATCTGTCATGTGATTCGTGCTTGCGCAGCGTAAGGCTGGCGAGGGAATCATTGAATGGGc                                   >  1:130829/1‑149 (MQ=255)                                                                                                                    ggCGCTTGGGCTGACTGGTAGGTACTTCGGAAGGCAGGCAGCACGCCGGGGCCAAGGTTCTAGGCCGTCGAGGTGGGATC‑CTGATCAATCTGTCATGTGATTCGTGCTTGCGCAGCGTAAGGCTGGCGAGGGAATCATTGAATGGGCCg                                 <  1:337891/149‑1 (MQ=255)                                                                                                                                 aCTGGTAGGTACTTCGGAAGGCAGGCAGCACGCCGGGGCCAAGGTTCTAGGCCGTCGAGGTGGGATC‑CTGATCAATCTGTCATGTGATTCGTGCTTGCGCAGCGTAAGGCTGGCGAGGGAATCATTGAATGGGCCGGTCCTTGTTTca                     <  5:366057/148‑1 (MQ=255)                                                                                                                                     gtaggtaCTTCGGAAGGCAGGCAGCACGCCGGGGCCAAGGTTCTAGGCCGTCGAGGTGGGATC‑CTGATCAATCTGTCATGTGATTCGTGCTTGCGCAGCGTAAGGCTGGCGAGGGAATCATTGAATGGGCCGGTCCTTGTTTcagcag                 <  3:102386/148‑1 (MQ=255)                                                                                                                                      taggtaCTTCGGAAGGCAGGCAGCACGCCGGGGCCAAGGTTCTAGGCCGTCGAGGTGGGATC‑CTGATCAATCTGTCATGTGATTCGTGCTTGCGCAGCGTAAGGCTGGCGAGGGAATCATTGAATGGGCCGGTCCTTGTTTCAGCAggg               >  6:345086/1‑149 (MQ=255)                                                                                                                                        ggtaCTTCGGGAGACAGGCAGCACGCCGGGGCCAAGGTTCTAGGCCGTCGAGGTGGGATC‑CTGATCAATCTGTCATGTGATTCGTGCTTGCGCAGCGTAAGGCTGGCGAGGGAATCATTGAATGGGCCGGTCCTTGTTTCAGCAGGGAc             <  7:256557/149‑1 (MQ=255)                                                                                                                                        ggtaCTTCGGGAGACAGGCAGCACGCCGGGGCCAAGGTTCTAGGCCGTCGAGGTGGGATC‑CTGATCAATCTGTCATGTGATTCGTGCTTGCGCAGCGTAAGGCTGGCGAGGGAATCATTGAATGGGCCGGTCCTTGTTTCAGCAGGGAc             <  7:261641/149‑1 (MQ=255)                                                                                                                                         gtaCTTCGGAAGGCAGGCAGCACGCCGGGGCCAAGGTTCTAGGCCGTCGAGGTGGGATC‑CTGATCAATCTGTCATGTGATTCGTGCTTGCGCAGc                                                                  <  7:295435/95‑1 (MQ=255)                                                                                                                                         gtaCTTCGGAAGGCAGGCAGCACGCCGGGGCCAAGGTTCTAGGCCGTCGAGGTGGGATC‑CTGATCAATCTGTCATGTGATTCGTGCTTGCGCAGc                                                                  >  8:295435/1‑95 (MQ=255)                                                                                                                                         gtaCTTCGGAAGGCAGGCAGCACGCCGGGGCCAAGGTTCTAGGCCGTCGAGGTGGGATC‑CTGATCAATCTGTCATGTGATTCGTGCTTGCGCAGCGTAAGGCTGGCGAGGGAATCATTGAATGGGCCGGTCCTTGTTTCAGCAGGGACg            >  1:195283/1‑149 (MQ=255)                                                                                                                                           aCTTCGGAAGGCAGGCAGCACGCCGGGGCCAAGGTTCTAGGCCGTCGAGGTGGGATC‑CTGATCAATCTGTCATGTGATTCGTGCTTGCGCAGCGTAAGGCTGGCGAGGGAATCATTGAATGGGCCGGTCCTTGTTTCAGCAGGGACGc           >  2:277530/1‑148 (MQ=255)                                                                                                                                            cTTCGGAAGGCAGGCAGCACGCCGGGGCCAAGGTTCTAGGCCGTCGAGGTGGGATC‑CTGATCAATCTGTCATGTGATTCGTGCTTGCGCAGCGTAAGGCTGGCGAGGGAATCATTGAATGGGCCGGTCCTTGTTTCAGCAGGGACGctc         >  3:157291/1‑149 (MQ=255)                                                                                                                                             ttCGGAAGGCAGGCAGCACGCCGGGGCCAAGGTTCTAGGCCGTCGAGGTGGGATC‑CTGATCAATCTGTCATGTGATTCGTGCTTGCGCAGCGTAAGGCTGGCGAGGGAATCATTGAATGGGCCGGTCCTTGTTTCAGCAGGGACGctct        >  6:104421/1‑149 (MQ=255)                                                                                                                                                 gAAGGCAGGCAGCACGCCGGGGCCAAGGTTCTAGGCCGTCGAGGTGGGATC‑CTGATCAATCTGTCATGTGATTCGTGCTTGCGCAGCGTAAGGCTGGCGAGGGAATCATTGAATGGGCCGGTCCTTGTTTCAGCAGGGACGCTCTTggg    <  3:342457/149‑1 (MQ=255)                                                                                                                                                 gAAGGCAGGCAGCACGCCGGGGCCAAGGTTCTAGGCCGTCGAGGTGGGATC‑CTGATCAATCTGTCATGTGATTCGTGCTTGCGCAGCGTAAGGCTGGCGAGGGAATCATTGAATGGGCCGGTCCTTGTTTCAGCAGGGACGCTCTTggg    >  6:137019/1‑149 (MQ=255)                                                                                                                                                 gAAGGCAGGCAGCACGCCGGGGCCAAGGTTCTAGGCCGTCGAGGTGGGATC‑CTGATCAATCTGTCATGTGATTCGTGCTTGCGCAGCGTAAGGCTGGCGAGGGAATCATTGAATGGGCCGGTCCTTGTTTCAGCAGGGACGCTCTTg      <  2:262690/147‑1 (MQ=255)                                                                                                                                                 gAAGGCAGGCAGCACGCCGGGGCCAAGGTTCTAGGCCGTCGAGGTGAGATC‑CTGATCAATCTGTCATGTGATTCGTGCTTGCGCAGCGTAAGGCTGGCGAGGGAATCATTGAATGGGCCGGTCCTTGTTTCAGCAGGGACGCTCTTg      >  1:262690/1‑147 (MQ=255)                                                                                                                                                  aaGGCAGGCAGCACGCCGGGGCCAAGGTTCTAGGCCGTCGAGGTGGGATC‑CTGATCAATCTGTCATGTGATTCGTGCTTGCGCAGCGTAAGGCTGGCGAGGGAATCATTGAATGGGCCGGTCCTTGTTTCAGCAGGGACGCTCTTgggg   <  4:312206/149‑1 (MQ=255)                                                                                                                                                   aggcaggcAGCACGCCGGGGCCAAGGTTCTAGGCCGTCGAGGTGGGATC‑CTGATCAATCTGTCATGTGATTCGTGCTTGCGCAGCGTAAGGCTGGCGAGGGAATCATTGAATGGGCCGGTCCTTGTTTCAGCAGGGACGCTCTTggggg  <  2:130829/149‑1 (MQ=255)                                                                                                                                                   aggcaggcAGCACGCCGGGGCCAAGGTTCTAGGCCGTCGAGGTGGGATC‑CTGATCAATCTGTCATGTGATTCGTGCTTGCGAAGCGTAAGGCTGGCGAGGGAATCATTGAATGGGCCGGTCCTTGTTTCAGCAGGGACGCTCTTggggg  >  7:214872/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                    TTTAAACACGGGGGTTATGATCCGTGCATGATAGAAACGTTCCAGTGACGCCTTTATTAAGGGACGAACGAACCTGACTGGAACAGGGCAGCCAGAGGGGGCGCAATCAAATGTAGGCGCTTGGGCTGACTGGTAGGTACTTCGGAAGGCAGGCAGCACGCCGGGGCCAAGGTTCTAGGCCGTCGAGGTGGGATC‑CTGATCAATCTGTCATGTGATTCGTGCTTGCGCAGCGTAAGGCTGGCGAGGGAATCATTGAATGGGCCGGTCCTTGTTTCAGCAGGGACGCTCTTGGGGG  >  NC_002947/1883686‑1883980

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