Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,763,506	C→T	9.2%	intergenic (+26/+1649)	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,506	0	C	T	9.2%	182.5 / 6.6	67	intergenic (+26/+1649)	PP_3326/PP_3328	hypothetical protein/ring‑cleaving dioxygenase
Reads supporting (aligned to +/- strand):  ref base C (33/27);  new base T (3/3);  total (37/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 = 3.75e-01

CGCAATGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGG‑TATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCATGTCGATTACCTGCATGATGGGCACTT  >  NC_002947/3763358‑3763647                                                                                                                                                     |                                                                                                                                               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)                                                                                                                                                     |                                                                                                                                               CGCAATGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGG‑TATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCATGTCGATTACCTGCATGATGGGCACTT  >  NC_002947/3763358‑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: ‑