Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,713,752	C→A	6.6%	intergenic (+21/‑81)	tsaB → / → PP_1508	endopeptidase/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,713,752	0	C	A	6.6%	216.1 / 3.9	70	intergenic (+21/‑81)	tsaB/PP_1508	endopeptidase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (34/30);  new base A (3/2);  total (37/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 = 5.52e-01

CTTGCCCAATGCTCTGGACATTCTCAGCCTGGCCACCTTTGCCTGGGCCCGTGGCGAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACG‑CCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATGCTGCTAAATTGCCATCATTGGTCCTGAGTGCTCATGCCATGCGTATCGACGGTTTCTCATCGCAGTCCTACC  >  NC_002947/1713605‑1713896                                                                                                                                                     |                                                                                                                                                 cTTGCCCAATGCTCTGGACATTCTCAGCCTGGCCACCTTTGCCTGGGCCCGTGGCGAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCCCG‑CCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGcc                                                                                                                                                 >  5:327433/1‑149 (MQ=255) cTTGCCCAATGCTCTGGACATTCTCAGCCTGGCCACCTTTGCCTGGGCCCGTGGCGAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACG‑CCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGcc                                                                                                                                                 >  7:188943/1‑149 (MQ=255)      ccAATGCTCTGGACATTCTCAGCCTGGCCACCTTTGCCTGGGCCCGTGGCGAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACG‑CCCAAGAAGCACTGAGTGCTGTTCGGCGATAATCGACGAaca                                                                                                                                            <  3:293406/149‑1 (MQ=255)       cAATGCTCTGGACATTCTCAGCCTGGCCACCTTTGCCTGGGCCCGTGGCGAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACG‑CCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAacac                                                                                                                                           >  6:98424/1‑149 (MQ=255)        aaTGCTCTGGACATTCTCAGCCTGGCCACCTTTGCCTGGGCCCGTGGCGAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACG‑CCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACg                                                                                                                                          <  6:277919/149‑1 (MQ=255)        aaTGCTCTGGACATTCTCAGCCTGGCCACCTTTGCCTGGGCCCGTGGCGAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACG‑CCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACg                                                                                                                                          >  4:415697/1‑149 (MQ=255)           ccgctGGACCATCTCAGCCGCGCCCCCGTTGCCTGCGCCCGTGGAGAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACG‑CCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTa                                                                                                                                       <  4:336408/146‑1 (MQ=255)            ctctGGACATTCTCAGCCTGGCCACCTTTGCCTGGGCCCGTGGCGAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACG‑CCCAAGAAGCACTGAGTGCTGTTCGGCGATAATCGACGAacagcact                                                                                                                                       <  7:286364/148‑6 (MQ=255)             tctGGACATTCTCAGCCTGGCCACCTTTGCCTGGGCCCGTGGCGAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACG‑CCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAaca                                                                                                                                            >  8:356809/1‑142 (MQ=255)             tctGGACATTCTCAGCCTGGCCACCTTTGCCTGGGCCCGTGGCGAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACG‑CCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAaca                                                                                                                                            <  7:356809/142‑1 (MQ=255)             tctGGACATTCTCAGCCTGGCCACCTTTGCCTGGGCCCGTGGCGAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACG‑CCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTaa                                                                                                                                      >  7:362926/1‑148 (MQ=255)                 gACATTCTCAGCCTGGCCACCTTTGCCTGGGCCCGTGGCGAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACG‑CCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAAcc                                                                                                                                  <  3:397180/148‑1 (MQ=255)                  aCATTCTCAGCCTGGCCACCTTTGCCTGGGCCCGTGGCGAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACG‑CCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCtt                                                                                                                                <  5:380038/149‑1 (MQ=255)                   cATTCTCAGCCTGGCCACCTTTGCCTGGGCCCGTGGCGAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACG‑CCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCttt                                                                                                                               <  1:215339/149‑1 (MQ=255)                       ctcAGCCTGGCCACCTTTGCCTGGGCCCGTGGCGAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACG‑CCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTgcg                                                                                                                           >  1:293553/1‑149 (MQ=255)                       ctcAGCCTGGCCACCTTTGCCTGGGCCCGTGGCGAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACG‑CCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTgcg                                                                                                                           >  3:341110/1‑149 (MQ=255)                        tcAGCCTGGCCACCTTTGCCTGGGCCCGTGGCGAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACG‑CCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTgcgc                                                                                                                          >  1:298359/1‑149 (MQ=255)                        tcAGCCTGGCCACCTTTGCCTGGGCCCGTGGCGAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACG‑CCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTgcgc                                                                                                                          >  1:315808/1‑149 (MQ=255)                         cAGCCTGGCCACCTTTGCCTGGGCCCGTGGCGAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACG‑CCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTgcgc                                                                                                                          >  4:116125/1‑148 (MQ=255)                                gCCACCTTTGCCTGGGCCCGTGGCGAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACG‑CCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCtgtg                                                                                                                  <  4:218181/149‑1 (MQ=255)                                   aCCTTTGCCTGGGCCCGTGGCGAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACG‑CCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTc                                                                                                               <  2:415084/149‑1 (MQ=255)                                        ggCCTGGGGGCGTGGCGAGGGAAGCGTGGCCGAGCAGGCGCATGAGTTTTATCTGCGTGATAATGTAGCCTCT‑CCGATGATGCTCTGAGTGGTGTTCGTTGATTATCGGCATTCTGTCTATAACCTTTTGCGCGATCTGTGGTCCa                                                                                                             <  4:54417/145‑1 (MQ=255)                                         gCCTGGGCCCGTGGCGAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACG‑CCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGtt                                                                                                          <  3:116125/148‑1 (MQ=255)                                          ccTGGGCCCGTGGCGAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACG‑CCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTAt                                                                                                        <  2:298359/149‑1 (MQ=255)                                                 ccGTGGCGAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACG‑CCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCaa                                                                                                 <  3:91455/149‑1 (MQ=255)                                                  cGTGGCGAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACG‑CCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAg                                                                                                <  4:341110/149‑1 (MQ=255)                                                    tGGCGAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACG‑CCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCa                                                                                              >  1:144040/1‑149 (MQ=255)                                                        gAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACG‑CCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAg                                                                                          <  2:290456/149‑1 (MQ=255)                                                         aGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACG‑CCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGc                                                                                         >  5:134587/1‑149 (MQ=255)                                                             aaTCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACG‑CCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACg                                                                                     >  6:89590/1‑149 (MQ=255)                                                              aTCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACG‑CCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACgg                                                                                    <  3:415697/149‑1 (MQ=255)                                                                    gCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACG‑CCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAActct                                                                              >  6:163625/1‑149 (MQ=255)                                                                    gCCGAGCAGGCGCAACCGGTGTATCTGCGCGATAATGTAGCCACG‑CCCAATAAGCACTGAGTGCTGGTCGTCTATTATCGCCGAACACGATAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACAGAActc                                                                               <  1:225932/148‑1 (MQ=255)                                                                     ccGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACGCCCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAActct                                                                              >  1:308084/1‑149 (MQ=255)                                                                     ccGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACG‑CCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAActct                                                                              >  7:412531/1‑148 (MQ=255)                                                                        aGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACG‑CCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGAtg                                                                          >  3:358438/1‑149 (MQ=255)                                                                             gcgcAACCGGTTTATCTGCGCGATAATGTAGCCACG‑CCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGAtgctgc                                                                      >  6:282124/1‑148 (MQ=255)                                                                             gcgcAACCGGTTTATCTGCGCGATAATGTAGCCACG‑CCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGAGATCTTTGGTCCAGTTATCACTCAAGCATTAGCGAAGGAACTCTGAtgctgc                                                                      >  6:323106/1‑148 (MQ=255)                                                                               ggacccGGTTTATCTGCGCGATAATGTAGCGCGG‑CCCAAGACGCACTGAGTGCTGTACCTCGATTATCGCCGCACCCGCTAAAACCTGTTGCGCGCTCTGTGGGCCAGTTATCACTCAGGCATTACCTACGCACCTCGGATGCTGCTa                                                                    >  6:111979/5‑148 (MQ=255)                                                                                   ccGGTTTATCTGCGCGATAATGTAGCCACG‑CCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATGCTGCTAAATTg                                                               <  6:134587/149‑1 (MQ=255)                                                                                     ggTTTATCTGCGCGATAATGTAGCCACG‑CCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATGCTGCTAAATTGcc                                                             <  2:144040/149‑1 (MQ=255)                                                                                      gTTTATCTGCGCGATAATGTAGCCACG‑CCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATGCTGCTAAATTGCca                                                            >  5:310207/1‑149 (MQ=255)                                                                                       tttATCTGCGCGATAATGTAGCCACG‑CCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATGCTGCTAAATTGCca                                                            >  5:152816/1‑148 (MQ=255)                                                                                          aTCTGCGCGATAATGTAGCCACG‑CCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATGCTGCTAAATTGCcatcat                                                        >  5:300/1‑149 (MQ=255)                                                                                           tCTGCGCGATAATGTAGCCACG‑CCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATGCTGCTAAATTGCCATCAtt                                                       >  2:367740/1‑149 (MQ=255)                                                                                               cgcgATAATGTAGCCACG‑CCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATGCTGCTAAATTGCCATCATTg                                                      <  6:327433/146‑1 (MQ=255)                                                                                                         tAGCCACG‑CCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATGCTGCTAAATTGCCATCATTGGTCCTGAGTGCTc                                         <  4:202328/149‑1 (MQ=255)                                                                                                           gCCACG‑CCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATGCTGCTAAATTGCCATCATTGGTCCTGAGTGCTcat                                       <  2:293553/149‑1 (MQ=255)                                                                                                                   ccAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATGCTGCTAAATTGCCATCATTGGTCCTGAGTGCTcatgccatgc                                <  5:163625/149‑1 (MQ=255)                                                                                                                   ccAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATGCTGCTAAATTGCAATCATTGGTCCTGAGTGATcatgccat                                  >  4:384674/1‑147 (MQ=255)                                                                                                                     aagaagCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACAAATTGCGCGATATGTGGTCCAGTTATCACTCAAGCATTAGAGACGGAACTCTGATGCTGCTAAATTGCCATCATTGCTCCTGAGTGCTCATGCCATGCg                               >  2:273110/1‑148 (MQ=255)                                                                                                                     aagaagCACTGAGTGCTGTTCGTCGATTATCGACGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATGCTGCTAAATTGCCATCATTCGTCCTGAGTGCTCATGCCATGCg                               >  2:267007/1‑148 (MQ=255)                                                                                                                          gCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATGCTGCTAAATTGCCATCATTGGTCCTGAGTGCTCATGCCATGCGTATCGa                         >  4:336999/1‑149 (MQ=255)                                                                                                                            aCTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATGCTGCTAAATTGCCATCATTGGTCCTGAGTGCTCATGCCATGCGTATCGACg                       >  6:307651/1‑149 (MQ=255)                                                                                                                               gAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATGCTGCTAAATTGCCATCATTGGTCCTGAGTGCTcatgccat                                  >  5:337140/1‑135 (MQ=255)                                                                                                                               gAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATGCTGCTAAATTGCCATCATTGGTCCTGAGTGCTcatgccat                                  <  6:337140/135‑1 (MQ=255)                                                                                                                                aGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATGCTGCTAAATTGCCATCATTGGTCCTGAGTGCTc                                         <  5:283219/127‑1 (MQ=255)                                                                                                                                aGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATGCTGCTAAATTGCCATCATTGGTCCTGAGTGCTc                                         >  6:283219/1‑127 (MQ=255)                                                                                                                                 gTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATGCTGCTAAATTGCCATCATTGGTCCTGAGTGCTCATGCCATGCGTATCGACGGTTtc                  <  5:89590/149‑1 (MQ=255)                                                                                                                                    gtGTTCGGCGATAATCGACGAACACGCTAAACCCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATCCTGCTAAATTGCCATCATTGGTCCTGAGAGCACCTGCCATGCGTATCGAAGGATTATCa               >  5:434402/2‑149 (MQ=255)                                                                                                                                    gtGTTCGGCGATAATCGACGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATGCTGCTAAATTGCCATCATTGGTCCTGAGTGCTCATGCCATGCGTATCGACGGTTTCTCa               >  1:27635/2‑149 (MQ=255)                                                                                                                                    cTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATGCTGCTAAATTGCCATCATTGGTCCTGAGTGCTCATGCCATGCGTATCGACGGTTTCTCa               <  5:282124/149‑1 (MQ=255)                                                                                                                                     tGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATGCTGCTAAATTGCCATCATTGGTCCTGAGTGCTCATGCCATGCGTATCGACGGTTTCTCAt              <  3:336999/149‑1 (MQ=255)                                                                                                                                      gTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATGCTGCTAAATTGCCATCATTGGTCCTGAGTGCTCATGCCATGCGTATCGACGGTTTCTCATc             >  2:82739/1‑149 (MQ=255)                                                                                                                                          gtcgATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATGCTGCTAAATTGCCATCATTGGTCCTGAGTGCTCATGCCATGCGTATCGACGGTTTCTCATCGCa          <  5:307651/148‑1 (MQ=255)                                                                                                                                           tcgATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATGCTGCTAAATTGCCATCATTGGTc                                                   >  2:325161/1‑106 (MQ=255)                                                                                                                                           tcgATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATGCTGCTAAATTGCCATCATTGGTc                                                   <  1:325161/106‑1 (MQ=255)                                                                                                                                             gATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATGCTGCTAAATTGCCATCATTGGTCCTGAGTGCTCATGCCATGCGTATCGACGGTTTCTCATCGCAGTcc      >  5:101467/1‑149 (MQ=255)                                                                                                                                               ttATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATGCTGCTAAATTGCCATCATTGGTCCTGAGTGCTCATGCCATGCGTATCGACGGTTTCTCATCGCAGTCCTa    >  5:311564/1‑149 (MQ=255)                                                                                                                                                  tCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATGCTGCTAAATTGCCATCATTGGTCCTGAGTGCTCATGCCATGCGTATCGACGGTTTCTCATCGCAGTCCTAcc  <  1:367740/148‑1 (MQ=255)                                                                                                                                                   cGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATGCTGCTAAATTGCCATCATTGGTCCTGAGTGCTCATGCCATGCGTATCGACGGTTTCTCATc             <  6:204759/136‑1 (MQ=255)                                                                                                                                                   cGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATGCTGCTAAATTGCCATCATTGGTCCTGAGTGCTCATGCCATGCGTATCGACGGTTTCTCATc             >  5:204759/1‑136 (MQ=255)                                                                                                                                                     |                                                                                                                                                 CTTGCCCAATGCTCTGGACATTCTCAGCCTGGCCACCTTTGCCTGGGCCCGTGGCGAGGCAATCGTCGCCGAGCAGGCGCAACCGGTTTATCTGCGCGATAATGTAGCCACG‑CCCAAGAAGCACTGAGTGCTGTTCGTCGATTATCGCCGAACACGCTAAAACCTTTTGCGCGATCTGTGGTCCAGTTATCACTCAAGCATTAGCGACGGAACTCTGATGCTGCTAAATTGCCATCATTGGTCCTGAGTGCTCATGCCATGCGTATCGACGGTTTCTCATCGCAGTCCTACC  >  NC_002947/1713605‑1713896

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