Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,556,984	C→A	8.9%	intergenic (+36/‑57)	purU‑II → / → PP_5468	formyltetrahydrofolate hydrolase/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,556,984	0	C	A	8.9%	171.4 / 5.2	56	intergenic (+36/‑57)	purU‑II/PP_5468	formyltetrahydrofolate hydrolase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (25/26);  new base A (3/2);  total (28/28)
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.69e-01

ATGGTCCGCTTCGGCCGTGATGTGGAGAAGATGGTGCTGGCCCGTGGCCTGCGTTATCACCTGGAGGACCGCGTGCTGGTGCATGGTAACAAGACTGTGGTGTTCGACTGACACCG‑TCCGGCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAACCCCGTCGCCTCACCTCCACCGATCCTTGGCGAAGGCTGCCTGGCCCGCTACGACCCTGAAGCCCTGGATGA  >  NC_002947/1556838‑1557132                                                                                                                                                    |                                                                                                                                                     aTGGTCCGCTTCGGCCGTGATGTGGAGAAGATGGTGCTGGCCCGTGGCCTGCGTTATCACCTGGAGGACCGCGTGCTGGTGCATGGTAACAAGACTGTGGTGTTCGACTGACACCG‑TCCGGCTTGCTGTGAGATTCCGGGCTGCTgcg                                                                                                                                                     >  3:113813/1‑148 (MQ=255)    gTCCGCTTCGGCCGTGATGTGGAGAAGATGGTGCTGGCCCGTGGCCTGCGTTATCACCTGGAGGACCGCGTGCTGGTGCATGGTAACAAGACTGTGGTGTTCGACTGACACCG‑TCCGGCTTGCTGTGAGATTCCGGGCTGCTGCGCAGc                                                                                                                                                 >  1:83629/1‑149 (MQ=255)       cGCTTCGGCCGTGATGTGGAGAAGATGGTGCTGGCCCGTGGCCTGCGTTATCACCTGGATGACCGCGTGCTGGTGCATGGTAACAAGACTGTGGTGTTCGACTGACACCGTTCCGGCTTGCTGTGAGATTCCGGGCTGCGCAGCAGccc                                                                                                                                               <  1:65206/149‑1 (MQ=255)         cTTCGGCCGTGATGTGGAGAAGATGGTGCTGGCCCGTGGCCTGCGTTATCACCTGGAGGACCGCGTGCTGGTGCATGGTAACAAGACTGTGGTGTTCGACTGACACCG‑TCCGGCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCCCTGt                                                                                                                                            <  4:140811/149‑1 (MQ=255)         cTTCGGCCGTGATGTGGAGAAGATGGTGCTGGCCCGTGGCCTGCGTTATCACCTGGAGGACCGCGTGCTGGTGCATGGTAACAAGACTGTGGTGTTCGACTGACACCG‑TCCGGCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCCCTGt                                                                                                                                            >  4:335211/1‑149 (MQ=255)            cGGCCGTGATGTGGAGAAGATGGTGCTGGCCCGTGGCCTGCGTTATCACCTGGAGGACCGCGTGCTGGTGCATGGTAACAAGACTGTGGTGTTCGACTGACACCG‑TCCGGCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCCCTGTCAc                                                                                                                                         >  4:161751/1‑149 (MQ=255)              gCCGTGATGTGGAGAAGATGGTGCTGGCCCGTGGCCTGCGTTATCACCTTGAGGACCGCGTGCTGGTGCATGGTAACAAGACTGTGTTGTTCGACTGACACGG‑TCCGGCTTGCTGTGAGATTCCGGGCtgctgc                                                                                                                                                      <  3:163503/134‑1 (MQ=255)              gCCGTGATGTGGAGAAGATGGAGCTGGACCGTGGACTGCGTTATCACATGGAGGACCGCGTGCTGGTGCATGGTAACAAGACTGTGGTGTTCGACTGACACCG‑TCCGGCTTCCTGTGAGATTCCGGGCtgatgc                                                                                                                                                      >  4:163503/1‑134 (MQ=255)               ccGTGATGTGGAGAAGATGGTGCTGGCCCGTGGCCTGCGTTATCACCTGGAGGACCGCGTGCTGGTGCATGGTAACAAGACTGTGGTGTTCGACTGACACCG‑TCCGGCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCCCTGTCACGa                                                                                                                                       >  5:248795/1‑148 (MQ=255)               ccGTGATGTGGAGAAGATGGTGCTGGCCCGTGGCCTGCGTTATCACCTGGAGGACCGCGTGCTGGTGCATGGTAACAAGACTGTGGTGTTCGACTGACACCG‑TCCGGCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCCCTGTCACGAc                                                                                                                                      <  1:207028/149‑1 (MQ=255)                    aTGTGGAGAAGATGGTGCTGGCCCGTGGCCTGCGTTATCACCTGGAGGACCGCGTGCTGGTGCATGGTAACAAGACTGTGGTGTTCGACTGACACCG‑TCCGGCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAggg                                                                                                                                  <  8:218781/148‑1 (MQ=255)                       tgGAGAAGATGGTGCTGGCCCGTGGCCTGCGTTATCACCTGGAGGACCGCGTGCTGGTGCATGGTAACAAGACTGTGGTGTTCGACTGACACCG‑TGCGGCTTGCTGTGAGATTGCGGGCTGCTGCGCAGCCCTGTCACGACAGGGgcgc                                                                                                                              <  1:316947/149‑1 (MQ=255)                         gagaAGATGGTGCTGGCCCGTGGCCTGCGTTATCACCTGGAGGACCGCGTGCTGGTGCATGGTAACAAGACTGTGGTGTTCGACTGACACCG‑TCCGGCTTGCTGTGGGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGgcg                                                                                                                               >  7:334278/1‑146 (MQ=255)                         gagaAGATGGTGCTGGCCCGTGGCCTGCGTTATCACCTGGAGGACCGCGTGCTGGTGCATGGTAACAAGACTGTGGTGTTCGACTGACACCG‑TCCGGCTTGCTGTGAGATTCCGGGCTGCGCAGCAGcccgg                                                                                                                                             >  4:219376/1‑130 (MQ=255)                         gagaAGATGGTGCTGGCCCGTGGCCTGCGTTATCACCTGGAGGACCGCGTGCTGGTGCATGGTAACAAGACTGTGGTGTTCGACTGACACCG‑TCCGGCTTGCTGTGAGATTCCGGGCTGCGCAGCAGcccgg                                                                                                                                             <  3:219376/132‑3 (MQ=255)                            aagaTGGTGCTGGCCCGTGGCCTGCGTTATCACCTGGAGGACCGCGTGCTGGTGCATGGTAACAAGACTGTGGTGTTCGACTGACACCG‑TCCGGCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAgc                                                                                                                         >  4:124318/1‑149 (MQ=255)                            aagaTGGTGCTGGCCCGTGGCCTGCGTTATCACCTCGAGGACCGCGTGCTGGTGCATGGTAACAACACTGTGGTGTTCGACTGACACCG‑TCCGGCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAgc                                                                                                                         >  1:357695/1‑149 (MQ=255)                                  gTGCTGGCCCGTGGCCTGCGTTATCACCTGGAGGACCGCGTGCTGGTGCATGGTAACAAGACTGTGGTGTTCGACTGACACCG‑TCCGGCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGcccc                                                                                                                    <  3:335211/148‑1 (MQ=255)                                   tGCTGGCCCGTGGCCTGCGTTATCACCTGGAGGACCGCGTGCTGGTGCATGGTAACAAGACTGTGGTGTTCGACTGACACCG‑TCCGGCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGccccag                                                                                                                  <  3:161751/149‑3 (MQ=255)                                     cTGGCCCGTGGCCTGCGTTATCACCTGGAGGACCGCGTGCTGGTGCATGGTAACAAGACTGTGGTGTTCGACTGACACCG‑TCCGGCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCTCTGCGcccctgtc                                                                                                                <  6:248795/149‑5 (MQ=255)                                      tGGCCCGTGGCCTGCGTTATCAGCTGGAGGACCGCGTGCTGGTGCAGGGTAGCAAGCCTGTGGTGTGCGACTGACACCG‑GCCCGCTGGCTGTGAGCTTCCGGGCTGCTGCGCAGCCCGGTGGCGACAGGGGCGCGGAGCGCCCCGGcca                                                                                                               >  6:352163/1‑147 (MQ=255)                                          ccGTGGCCTGCGTTATCACCTGGAGGACCGCGTGCTGGTGCATGGTAACAAGACTGTGGTGTTCGACTGACACCG‑TCCGGCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCt                                                                                                           >  2:190247/1‑149 (MQ=255)                                               gCCTGCGTTATCACCTGGAGGACCGCGTGCTGGTGCATGGTAACAAGACTGTGGTGTTCGACTGACACCG‑TCCGGCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGt                                                                                                       <  1:92388/148‑1 (MQ=255)                                                   gCGTTATCACCTGGAGGACCGCGTGCTGGTGCATGGTAACAAGACTGTGGTGTTCGACTGACACCG‑TCCGGCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAAc                                                                                                  >  8:85511/1‑149 (MQ=255)                                                    ggTTTTCACCTGGAGGACCGCGTGCTGGTGCATGGAAACAAGACTGTGGTGTTCGACTGACACCG‑TCCGGCTTGCGGTGAGTTCCCGGGCTCCTGCGCAGCCCTGTCTCGACAGGGGCGCAGACCGCCCCGGCAATCCTGGGTGAAACg                                                                                                 <  2:357695/148‑1 (MQ=255)                                                         tCACCTGGAGGACCGCGTGCTGGTGCATGGTAACAAGACTGTGGTGTTCGACTGACACCG‑TCCGGCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCaa                                                                                                               >  3:138120/1‑130 (MQ=255)                                                         tCACCTGGAGGACCGCGTGCTGGTGCATGGTAACAAGACTGTGGTGTTCGACTGACACCG‑TCCGGCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCaa                                                                                                               <  4:138120/130‑1 (MQ=255)                                                                gaggaCCGCGTGCTGGTGCATGGTAACAAGACTGTGGTGTTCGACTGACACCG‑TCCGGCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCaa                                                                                      <  2:83629/148‑1 (MQ=255)                                                                  ggaCCGCGTGCTGGTGCATGGTAACAAGACTGTGGTGTTCGACTGACACCG‑TCCGGCTTGCTGTGGGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAgca                                                                                   <  8:334278/149‑1 (MQ=255)                                                                       gcgTGCTGGTGCATGGTAACAAGACTGTGGTGTTCGACTGACACCG‑TCCGGCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCgg                                                                                                                  <  4:262465/113‑1 (MQ=255)                                                                       gcgTGCTGGTGCATGGTAACAAGACTGTGGTGTTCGACTGACACCG‑TCCGGCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCgg                                                                                                                  >  3:262465/1‑113 (MQ=255)                                                                               gTGCATGGTAACAAGACTGTGGTGTTCGACTGACACCG‑TCCGGCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAAcccc                                                                      <  3:124318/149‑1 (MQ=255)                                                                                        aaCAAGACTGTGGTGTTCGACTGACACCG‑TCCGGCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAACCCCGTCGCCTCa                                                             >  4:117534/1‑149 (MQ=255)                                                                                             gACTGTGGTGTTCGACTGACACCG‑TCCGGCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAACCCCGTCGCCTCACCTc                                                         >  1:3304/1‑148 (MQ=255)                                                                                                    gtgtTCGACTGACACCG‑TCCGGCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAACCCCGTCGCCTCACCTCCACCGATc                                                 <  5:136097/149‑1 (MQ=255)                                                                                                    gtgtTCGACTGACACCG‑TCCGGCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAACCCCGTCGCCTCACCTCCACCGATc                                                 >  1:307462/1‑149 (MQ=255)                                                                                                       ttCGACTGACACCG‑TCCGGCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAACCCCGTCGCCTCACCTCCACCGATCCtt                                              >  3:220791/1‑149 (MQ=255)                                                                                                        tCGACTGACACCG‑TCCGGCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAACCCCGTCGCCTCACCTCCACCGATCCTTg                                             >  3:34799/1‑149 (MQ=255)                                                                                                         cGACTGACACCG‑TCCGGCTTGCTGTGAGGTTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAACCCCGTCGCCt                                                               <  6:250704/130‑1 (MQ=255)                                                                                                         cGACTGACACCG‑TCCGGCTTGCTGTGAGGTTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAACCCCGTCGCCt                                                               >  5:250704/1‑130 (MQ=255)                                                                                                             tGACACCG‑TCCGGCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTggg                                                                                                        <  4:151379/85‑1 (MQ=255)                                                                                                             tGACACCG‑TCCGGCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTggg                                                                                                        >  3:151379/1‑85 (MQ=255)                                                                                                             tGACACCG‑TCCGGCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGCAACGGCAATGACCAAGCAACGCAACAACCCCGTCGCCTCACCTCCACCGATCCTTGGCg                                          >  7:277815/1‑147 (MQ=255)                                                                                                                cacCG‑TCCGGCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAACCCCg                                                                     <  7:277611/117‑1 (MQ=255)                                                                                                                cacCG‑TCCGGCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAACCCCg                                                                     >  8:277611/1‑117 (MQ=255)                                                                                                                          gCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAACCCCGTCGCCTCACCTCCACCGATCCTTg                                             <  3:45744/132‑1 (MQ=255)                                                                                                                          gCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGCGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAACCCCGTCGCCTCACCTCCACCGATCCTTg                                             >  4:45744/1‑132 (MQ=255)                                                                                                                                tgtgAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAACCCCGTCGCCTCACCTCCACCGATCCTTGGCGAAGGCTGc                                  >  1:67180/1‑137 (MQ=255)                                                                                                                                tgtgAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAACCCCGTCGCCTCACCTCCACCGATCCTTGGCGAAGGCTGCCTGGCCCGCTAc                      >  3:69111/1‑149 (MQ=255)                                                                                                                                tgtgAGATTCCGGGCGGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAACCCCGTCGCCTCACCTCCACCGATCCTTGGCGAAGGCTGc                                  <  2:67180/137‑1 (MQ=255)                                                                                                                                  tgAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAACCCCGTCGCCTCACCTCCACCGATCCTTGGCGAAGGCTGCCTGGCCCGCTACGa                    <  4:220791/149‑1 (MQ=255)                                                                                                                                    agaTTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAACCCCGTCGCCTCACCTCCACCGATCCTTGGCGAAGGCTGCCTGGCCCGCTACGAcc                  >  2:228551/1‑149 (MQ=255)                                                                                                                                    agaTTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAACCCCGTCGCCTCACCTCCACCGATCCTTGGCGAAGGCTGCCTGGCCCGCTACGAcc                  <  2:307462/149‑1 (MQ=255)                                                                                                                                    agaTTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAACCCCGTCGCCTCACCTCCACCGATCCTTGGCGAAGGCTGCCTGGCCCGCTACGAcc                  <  7:85511/149‑1 (MQ=255)                                                                                                                                    agaTTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAACCCCGTCGCCTCACCTCCACCGATCCTTGGCGAAGGCTGCCTGGCCCGCTACGAc                   <  1:190247/148‑1 (MQ=255)                                                                                                                                       ttCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGCTGAAACGGCAATGACCAAGCAACGCAACAACCACGTCGCCTCACCTCCACCGATCCTTGGCGAAGGCTGCCTGGCCCGCTACGACCCTg               <  4:113813/149‑1 (MQ=255)                                                                                                                                       gacagggCTGCGCAGCAGCCCGGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAACCCCGTCGCCTCACCTCCACCGATCCTTGGCGAAGGCTGCCTGGCCCGCTACGACCCt                >  8:133082/5‑148 (MQ=255)                                                                                                                                        acagggCTGCGCAGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAACCAACGCAACACCCCCGTCGCCTCACCTCCCCCGATCCTTGGCGAAGGCTGCCTGGCCCGCTACGACCCTg               >  2:221058/4‑148 (MQ=255)                                                                                                                                          cGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGCGTGAAACGGCAATGACCAAGCAACGCAACAACCCCGTCGCCTCACCTCCACCGATCCTTGGCGAAGGCTGCCTGGCCCGCTACGACCCTGaa             <  5:352163/148‑1 (MQ=255)                                                                                                                                                    cgcAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAACCCCGTCGCCTCACCTCCACCGATCCTTGGCGAAGGCTGCCTGGCCCGCTACGACCCTGAAGCCCTGGAt    <  2:3304/147‑1 (MQ=255)                                                                                                                                                    cgcAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAACCCCGTCGCCTCACCTCCACCGATCCTTGGCGAAGGCTGCCTGGCCCGCTACGACCCTGAAGCCCTGGATGa  >  2:193476/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                     ATGGTCCGCTTCGGCCGTGATGTGGAGAAGATGGTGCTGGCCCGTGGCCTGCGTTATCACCTGGAGGACCGCGTGCTGGTGCATGGTAACAAGACTGTGGTGTTCGACTGACACCG‑TCCGGCTTGCTGTGAGATTCCGGGCTGCTGCGCAGCCCTGTCACGACAGGGGCGCAGAGCGCCCCGGCAATCCTGGGTGAAACGGCAATGACCAAGCAACGCAACAACCCCGTCGCCTCACCTCCACCGATCCTTGGCGAAGGCTGCCTGGCCCGCTACGACCCTGAAGCCCTGGATGA  >  NC_002947/1556838‑1557132

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