Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,763,502	A→G	8.7%	intergenic (+22/+1653)	PP_3326 → / ← PP_3328	hypothetical protein/ring‑cleaving dioxygenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,763,502	0	A	G	8.7%	187.4 / 6.7	69	intergenic (+22/+1653)	PP_3326/PP_3328	hypothetical protein/ring‑cleaving dioxygenase
Reads supporting (aligned to +/- strand):  ref base A (36/27);  new base G (3/3);  total (39/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.94e-01

TACGCAATGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGG‑TATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCATGTCGATTACCTGCATGATGGGCACTT  >  NC_002947/3763356‑3763647                                                                                                                                                   |                                                                                                                                                   tACGCAATGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCAt                                                                                                                                                   >  5:329112/1‑148 (MQ=255)   cGCAATGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGAc                                                                                                                                                >  8:157423/1‑149 (MQ=255)    gCAATGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACg                                                                                                                                               >  7:146730/1‑149 (MQ=255)     cAATGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACgg                                                                                                                                              >  5:165050/1‑149 (MQ=255)      aaTGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCGTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGc                                                                                                                                             >  6:132559/1‑149 (MQ=255)      aaTGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGc                                                                                                                                             >  6:137505/1‑149 (MQ=255)       aTGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGcc                                                                                                                                            <  7:157423/149‑1 (MQ=255)        tGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGcc                                                                                                                                            >  2:300031/1‑148 (MQ=255)        tGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCg                                                                                                                                           <  7:61449/149‑1 (MQ=255)            ggACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCGTCATGGCCGAGG‑t                                                                                                                                      <  6:230127/149‑1 (MQ=255)             gACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCGTCATGGCCGAGG‑Ta                                                                                                                                     <  2:296933/149‑1 (MQ=255)             gACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGG‑Ta                                                                                                                                     >  1:278811/1‑149 (MQ=255)                tCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGG‑TATgg                                                                                                                                  >  1:272893/1‑149 (MQ=255)                  aaCATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAgg                                                                                                                                        >  4:159691/1‑142 (MQ=255)                  aaCATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAgg                                                                                                                                        <  3:159691/142‑1 (MQ=255)                  aaCATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGG‑TATgg                                                                                                                                  >  6:155407/1‑147 (MQ=255)                      tcatcGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGG‑TATGGCAGTTg                                                                                                                            >  7:3070/1‑149 (MQ=255)                         tcGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGG‑TATGGCAGTTGGAg                                                                                                                         <  5:155407/149‑1 (MQ=255)                         tcGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGG‑TATGGCAGTTGGAg                                                                                                                         <  2:272893/149‑1 (MQ=255)                            cgGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGG‑TATGGCAGTTGGAGGCg                                                                                                                      <  3:90591/149‑1 (MQ=255)                              gCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGG‑TATGGCAGTTGGAGGCGtt                                                                                                                    >  3:147338/1‑149 (MQ=255)                               cccAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGG‑TATGGCAGTTGGAGGCGTTc                                                                                                                   >  1:25557/1‑149 (MQ=255)                               cccAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGG‑TATGGCAGTTGGAGGCGTTc                                                                                                                   >  2:251906/1‑149 (MQ=255)                                   ggCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGG‑TATGGCAGTTGGAGGCGTTCGGaa                                                                                                               >  1:286500/1‑149 (MQ=255)                                   ggCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGG‑TATGGCAGTTGGAGGCGTTCGGaa                                                                                                               <  6:329112/149‑1 (MQ=255)                                        cccTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGG‑TATGGCAGTTGGAGGCGTTCGGAATGa                                                                                                            <  8:146730/147‑1 (MQ=255)                                          cTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGG‑TATGGCAGTTGGAg                                                                                                                         >  1:222841/1‑132 (MQ=255)                                          cTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGG‑TATGGCAGTTGGAg                                                                                                                         <  2:222841/132‑1 (MQ=255)                                           tGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAAGTAAGCGACATACCTCGGCCATGAGGGCCGCGG‑TATGGCAGTTGGAGGCGTTCGGAATGAGTGcc                                                                                                       >  7:212697/1‑149 (MQ=255)                                             ccAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGG‑TATGGCAGTTGGAGGCGTTCGGAATGAGTGCCt                                                                                                      <  1:300031/148‑1 (MQ=255)                                             ccAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGG‑TATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTg                                                                                                     >  3:104674/1‑149 (MQ=255)                                              cAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCCGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGG‑TATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTg                                                                                                     >  2:121276/1‑148 (MQ=255)                                                   gCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGG‑TATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAAc                                                                                                 >  5:157376/1‑147 (MQ=255)                                                   gCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGG‑TATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAAc                                                                                                 <  6:157376/147‑1 (MQ=255)                                                       aCGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGG‑TATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCAt                                                                                           >  5:110823/1‑149 (MQ=255)                                                             gCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGG‑TATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCAc                                                                                     <  2:8200/149‑1 (MQ=255)                                                                 cGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGG‑TATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATg                                                                                 >  4:204799/1‑149 (MQ=255)                                                                     ggccAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGAc                                                                               <  1:251906/146‑1 (MQ=255)                                                                            agagCCCGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGG‑TATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGcac                                                                      <  1:121276/149‑1 (MQ=255)                                                                               gCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGG‑TATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACaa                                                                    <  2:278811/148‑1 (MQ=255)                                                                                ccTGTTCAGTTCACTGTCGTTCCTGGGTAAGTCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGG‑TATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGc                                                                  >  3:7616/1‑149 (MQ=255)                                                                                ccTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGG‑TAt                                                                                                                                    <  7:12516/83‑1 (MQ=255)                                                                                ccTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGG‑TAt                                                                                                                                    >  8:12516/1‑83 (MQ=255)                                                                                      cagttcaCTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAGGCGACATACCTCGGCCATGACGGGCGAGG‑TATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCt                                                            >  1:67550/1‑149 (MQ=255)                                                                                      cagttcaCTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGG‑TATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCg                                                                                        <  6:181269/121‑1 (MQ=255)                                                                                        gttcaCTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGG‑TATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGc                                                          >  2:246356/1‑149 (MQ=255)                                                                                           caCTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGG‑TATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGc                                                       >  7:7686/1‑149 (MQ=255)                                                                                                  gTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGG‑TATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAACCCCGGCTGCGGCCACACgg                                                >  2:275590/1‑149 (MQ=255)                                                                                                   ttCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGG‑TATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGc                                               <  8:7686/149‑1 (MQ=255)                                                                                                   ttCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGG‑TATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGc                                               <  3:204799/149‑1 (MQ=255)                                                                                                          gTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGG‑TATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGaa                                        <  4:147338/149‑1 (MQ=255)                                                                                                               ccGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGG‑TATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGa                                   >  3:161315/1‑149 (MQ=255)                                                                                                                 gCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGG‑TATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATg                                 <  2:218566/149‑1 (MQ=255)                                                                                                                   accacGAGTATAGTAAGCGACATACCTCGGCCATGACGGCCGAGG‑TATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCATCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGa                                   <  4:147917/145‑1 (MQ=255)                                                                                                                   accacGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGG‑TATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGAtag                                >  3:147917/1‑146 (MQ=255)                                                                                                                        gAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGG‑TATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCATGt                          >  4:154566/1‑149 (MQ=255)                                                                                                                           tAAACTAAGCGACATACCTCGGCCATGACGGCCGAGG‑TATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCATGTCGa                       <  3:154566/149‑1 (MQ=255)                                                                                                                              aCTAAGCGACATACCTCGGCCATGACGGCCGAGG‑TATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCATGTCGATTa                    <  6:165050/149‑1 (MQ=255)                                                                                                                                     gACATACCTCGGCCATGACGGCCGAGG‑TATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCATGTCGATTACCTGCat             <  2:67550/149‑1 (MQ=255)                                                                                                                                     gACATACCTCGGCCATGACGGCCGAGG‑TATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCATGTCGATTACCTGCat             >  4:51450/1‑149 (MQ=255)                                                                                                                                      aCATACCTCGGCCATGACGGGCGAGG‑TTTGGCAGTTGGAGGCGTTGGGAATGAGTGCCTGCACCCAGCATCCGCACCATGAGCTTTGGCACAAGCCCGGGTGCGGGCACAGGGCTATCGACGACGATGGGGAGGACGATTACCtcaagg            >  6:51102/1‑144 (MQ=255)                                                                                                                                       cATACCTCGGCCATGACGGCCGAGG‑TATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGCCAATCGAACACGATGGGCATGTCGATTACCTGCatga           >  7:257282/1‑149 (MQ=255)                                                                                                                                        aTACCTCGGCCGTCATGGCCGAGG‑TATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCATGTCGATTACCTGCatgat          <  7:217938/149‑1 (MQ=255)                                                                                                                                        aTACCTCGGCCGTCATGGCCGAGG‑TATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCATGTCGATTACCTGCatgat          >  8:217938/1‑149 (MQ=255)                                                                                                                                          aCCTCGGCCGTCATGGCCGAGG‑TATGGAAGTTGGCGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCATGTCGATTACCTGCATGATgg        >  2:231531/1‑149 (MQ=255)                                                                                                                                           ccTCGGCCGTCATGGCCGAGG‑TATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCATGTCGATTACCTGCATGATggg       >  3:69655/1‑149 (MQ=255)                                                                                                                                            cTCGGCCATGACGGCCGAGG‑TATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCATGTCGATTACCTGCATGATggg       >  3:144160/1‑148 (MQ=255)                                                                                                                                               ggCCATGACGGCCGAGG‑TATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCATGTCGATTACCTGCATGATGGGCAc    <  4:144160/148‑1 (MQ=255)                                                                                                                                                gCCATGACGGCCGAGG‑TATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCATGTCGATTACCTGCATGATGGGCACtt  <  4:161315/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                   TACGCAATGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGG‑TATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCATGTCGATTACCTGCATGATGGGCACTT  >  NC_002947/3763356‑3763647

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