Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,957,995 G→A 6.9% P812S (CCG→TCG)  mnxG ← manganese‑oxidizing multicopper oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,957,9950GA 6.9% 165.0 / 2.8 58P812S (CCG→TCG) mnxGmanganese‑oxidizing multicopper oxidase
Reads supporting (aligned to +/- strand):  ref base G (28/26);  new base A (2/2);  total (30/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.12e-01

GCGGTCGGCAGTGAAGGCCGAGGCGCCACGGAAGTTCACGCCCGTCAGGCTCTCTCCACTGTTGAACTCCCCAACCCCGGAAGACTGGGTCAGGCGTTTTCCACGGTCATCCATGCAGGGCTCGTAGAACGGCGCACCCGCAGTCGGCAAGGCGCCGTTGGTGCGGAACGCGCTGGCCACTGGCTGGCTGCCGGGGATCAATGCATAGCTGGGGTGTTCACCCATGGCGTGGAACTGCATGGCTGCCTGCTCGACATCGGTGCCCTCTTCCGGCAGGTAGATGGGCTTGGCCTT  >  NC_002947/3957849‑3958142
                                                                                                                                                  |                                                                                                                                                   
gCGGTCGGCAGTGAAGGCCGAGGCGCCACGGAAGTTCACGCCCGTCAGGCTCTCTCCACTGTTGAACTCCCCAACCCCGGAAGACTGGGTCAGGCGTTTTCCACGGTCATCCATGCAGGGCTCGTAGAACGGCGCACCCGCAGTCGGCaa                                                                                                                                                  <  1:1237311/150‑1 (MQ=255)
gCGGTCGGCAGTGAAGGCCGAGGCGCCACGGAAGTTCACGCCCGTCAGGCTCTCTCCACTGTTGAACTCCCCAACCCCGGAAGACTGGGTCAGGCGTTTTCCACGGTCATCCATGCAGGGCTCGTAGAACGGCGCACCCGCAGTCGGCaa                                                                                                                                                  >  1:149327/1‑150 (MQ=255)
    tCGGCAGTGAAGGCCGAGGCGCCACGGAAGTTCACGCCCGTCAGGCTCTCTCCACTGTTGAACTCCCCAACCCCGGAAGACTGGGTCAGGCGTTTTCCACGGTCATCCATGCAGGGCTCGTAGAACGGCGCACCCGCAGTCGGCAAGgcg                                                                                                                                              >  2:876217/1‑150 (MQ=255)
    tCGGCAGTGAAGGCCGAGGCGCCACGGAAGTTCACGCCCGTCAGGCTCTCTCCACTGTTGAACTCCCCAACCCCGGAAGACTGGGTCAGGCGTTTTCCACGGTCATCCATGCAGGGCTCGTAGAACGGCGCACCCGCAGTCGGCAAGgcg                                                                                                                                              >  2:875952/1‑150 (MQ=255)
           tGAAGGCCGAGGCGCCACGGAAGTTCACGCCCGTCAGGCTCTCTCCACTGTTGAACTCCCCAACCCCGGAAGACTGGGTCAGGCGTTTTCCACGGTCATCCATGCAGGGCTCGTAGAACGGCGCACCCGCAGTCGGCAAGGCGCCGTTgg                                                                                                                                       >  2:917139/1‑150 (MQ=255)
           tGAAGGCCGAGGCGCCACGGAAGTTCACGCCCGTCAGGCTCTCTCCACTGTTGAACTCCCCAACCCCGGAAGACTGGGTCAGGCGTTTTCCACGGTCATCCATGCAGGGCTCGTAGAACGGCGCACCCGCAGTCGGCAAGGCGCCGTTgg                                                                                                                                       >  2:480258/1‑150 (MQ=255)
             aaGGCCGAGGCGCCACGGAAGTTCACGCCCGTCAGGCTCTCTCCACTGTTGAACTCCCCAACCCCGGAAGACTGGGTCAGGCGTTTTCCACGGTCATCCATGCAGGGCTCGTAGAACGGCGCACCCGCAGTCGGCAAGGCGCCGTTGGTg                                                                                                                                     <  2:251130/150‑1 (MQ=255)
               ggCCGAGGCGCCACGGAAGTTCACGCCCGTCAGGCTCTCTCCACTGTTGAACTCCCCAACCCCGGAAGACTGGGTCAGGCGTTTTCCACGGTCATCCATGCAGGGCTCGTAGAACGGCGCACCCGCAGTCGGCAAGGCGCCGTTGGTGCg                                                                                                                                   >  2:884454/1‑150 (MQ=255)
                      gcgcCACGGAAGTTCACGCCCGTCAGGCTCTCTCCACTGTTGAACTCCCCAACCCCGGAAGACTGGGTCAGGCGTTTTCCACGGTCATCCATGCAGGGCTCGTAGAACGGCGCACCCGCAGTCGGCAAGGCGCCGTTGGTGCGGAAcgcg                                                                                                                            >  2:318899/1‑150 (MQ=255)
                                     ccgcccgccAGACTCTCTCCTCTGTTGAACTCCCCAACCCCGGAAGACTGGGTCAGGCGTTTTCCCCGCTCATCCATGCAGGGCTCGTAGACCGGCACCATCGAAGGCGGCAAGGCGCCGTGGGTGCTGAACGCGCTGGCCActggctgg                                                                                                             <  1:292338/149‑1 (MQ=255)
                                     aCGCCCGTCAGGCTCTCTCCACTGTTGAACTCCCCAACCCCGGAAGACTGGGTCAGGCGTTTTCCACGGTCATCCATGCAGGGCTCGTAGAACGGCGCACCCGCAGTCGGCAAGGCGCCGTTGGTGCGGAACGCGCTGGCCActggctgg                                                                                                             <  1:318899/150‑1 (MQ=255)
                                      cGCCCGTCAGGCTCTCTCCACTGTTGAACTCCCCAACCCCGGAAGACTGGGTCAGGCGTTTTCCACGGTCATCCATGCAGGGCTCGTAGAACGGCGCACCCGCAGTCGGCAAGGCGCCGTTGGTGCGGAACGCGCTGGCCACtggctggc                                                                                                            >  1:1236970/1‑150 (MQ=255)
                                                 ctctctCCACTGTTGAACTCCCCAACCCCGGAAGACTGGGTCAGGCGTTTTCCACGGTCATCCATGCAGGGCTCGTAGAACGGCGCACCCGCAGTCGGCAAGGCGCCGTTGGt                                                                                                                                      <  1:722228/113‑1 (MQ=255)
                                                 ctctctCCACTGTTGAACTCCCCAACCCCGGAAGACTGGGTCAGGCGTTTTCCACGGTCATCCATGCAGGGCTCGTAGAACGGCGCACCCGCAGTCGGCAAGGCGCCGTTGGt                                                                                                                                      >  2:722228/1‑113 (MQ=255)
                                                        cACTGTTGAACTCCCCAACCCCGGAAGACTGGGTCAGGCGTTTTCCACGGTCATCCATGCAGGGCTCGTAGAACGGCGCACCCGCAGTCGGCAAGGCGCCGTTGGTGCGGAACGCGCTGGCCACTGGCTGGCTGCCGGGGATCAATGCAt                                                                                          >  1:1001730/1‑150 (MQ=255)
                                                        cACTGTTGAACTCCCCAACCCCGGAAGACTGGGTCAGGCGTTTTCCACGGTCATCCATGCAGGGCTCGTAGAACGGCGCACCCGCAGTCGACAAGGCGCCGTTGGTGCGGAACGCGCTGGCCACTGGCTGGCTGCCGGGGATCAATGCAt                                                                                          >  2:899590/1‑150 (MQ=255)
                                                        cACTGTTGAACTCCCCAACCCCGGAAGACTGGGTCAGGCGTTTTCCACGGTCATCCATGCAGGGCTCGTAGAACGGCGCACCCGCAGTCGACAAGGCGCCGTTGGTGCGGAACGCGCTGGCCACTGGCTGGCTGCCGGGGATCAATGCAt                                                                                          >  2:899700/1‑150 (MQ=255)
                                                           tGTTGAACTCCCCAACCCCGGAAGACTGGGTCAGGCGTTTTCCACGGTCATCCATGCAGGGCTCGTAGAACGGCGCACCCGCAGTCGGCAAGGCGCCGTTGGTGCGGAACGCGCTGGCCACTGGCTGGCTGCCGGGGATCAATGCATAGc                                                                                       >  2:220156/1‑150 (MQ=255)
                                                               gAACTCCCCAACCCCGGAAGACTGGGTCAGGCGTTTTCCACGGTCATCCATGCAGGGCTCGTAGAACGGCGCACCCGCAGTCGGCAAGGCGCCGTTGGTGCGGAACGCGCTGGCCACTGGCTGGCTGCCGGGGATCAATGCATAGCTggg                                                                                   <  2:1001730/150‑1 (MQ=255)
                                                                aaCTCCCCAACCCCGGAAGACTGGGTCAGGCGTTTTCCACGGTCATCCATGCAGGGCTCGTAGAACGGCGCACCCGCAGTCGGCAAGGCGCCGTTGGTGCGGAACGCGCTGGCCACTGGCTGGCTGCCGGGGATCAATGCATAGCTg                                                                                     <  1:571573/147‑1 (MQ=255)
                                                                aaCTCCCCAACCCCGGAAGACTGGGTCAGGCGTTTTCCACGGTCATCCATGCAGGGCTCGTAGAACGGCGCACCCGCAGTCGGCAAGGCGCCGTTGGTGCGGAACGCGCTGGCCACTGGCTGGCTGCCGGGGATCAATGCATAGCTg                                                                                     >  2:571573/1‑147 (MQ=255)
                                                                aaCTCCCCAACCCCGGAAGACTGGGTAAGGCGTTTTCCACGGTCATCCATACAGGGCTCGTAGTACGGCGCACCCGCAGTCGGCAAGGCGTCGTTGGTGCGGAACGCGCTGGCCGCTGGCTGGGTTCGGGGGGTCAATGAATAGCGgggg                                                                                  >  2:610772/1‑150 (MQ=255)
                                                                 aCTCCCCAACCCCGGAAGACTGGGTCAGGCGTTTTCCACGGTCATCCATGCAGGGCTCGTAGAACGGCGCACCCGCAGTCGGCAAGGCGCCGTTGGTGCGGAACGCGCTGGCCACTGGCTGGCTGCCGGGGATCAATGCATAGCTGGGgt                                                                                 <  1:220156/150‑1 (MQ=255)
                                                                 aCTCCCCAACCCCGGAAGACTGGGTCAGGCGTTTTCCACGGTCATCCATGCAGGGCTCGTAGAACGGCGCACCCGCAGTCGACAAGGCGCCGTTGGTGCGGAACGCGCTGGCCACTGGCTGGCTGCCGGGGATCAATGCATAGCTGGGgt                                                                                 <  1:899590/150‑1 (MQ=255)
                                                                 aCTCCCCAACCCCGGAAGACTGGGTCAGGCGTTTTCCACGGTCATCCATGCAGGGCTCGTAGAACGGCGCACCCGCAGTCGACAAGGCGCCGTTGGTGCGGAACGCGCTGGCCACTGGCTGGCTGCCGGGGATCAATGCATAGCTGGGgt                                                                                 <  1:899700/150‑1 (MQ=255)
                                                                   tCCCCAACCCCGGAAGACTGGGTCAGGCGTTTTCCACGGTCATCCATGCAGGGCTCGTAGAACGGCGCACCCGCAGTCGGCAAGGCGCCGTTGGTGCGGAACGCGCTGGCCACTGGCTGGCTGCCGGGGATCAATGCATAGCTGGGgtgt                                                                               <  1:917139/150‑1 (MQ=255)
                                                                     cccAACCCCGGAAGACTGGGTCAGGCGTTTTCCACGGTCATCCATGCAGGGCTCGTAGAACGGCGCACCCGCAGTCGGCAAGGCGCCGTTGGTGCGGAACGCGCTGGCCACTGGCTGGCTGCCGGGGATCAATGCATAGCTGGGGTGTTc                                                                             <  2:929569/150‑1 (MQ=255)
                                                                         aCCCCGGAAGACTGGGTCAGGCGTTTTCCACGGTCATCCATGCAGGGCTCGTAGAACGGCGCACCCGCAGTCGGCAAGGCGCCGt                                                                                                                                          <  2:380406/85‑1 (MQ=255)
                                                                         aCCCCGGAAGACTGGGTCAGGCGTTTTCCACGGTCATCCATGCAGGGCTCGTAGAACGGCGCACCCGCAGTCGGCAAGGCGCCGt                                                                                                                                          >  1:380406/1‑85 (MQ=255)
                                                                            ccGGAAGACTGGGTCAGGCGTTTTCCACGGTCATCCTTGCAGGGCTCGTAGAACGGCGCACCCGCAGTCGGCAAGGCGCCGTTGGTGCGGAACGCGCTGGCCACTGGCTGGCTGCCGGGGATCAATGCATAGCTGGGGTGTTCACCCATg                                                                      <  2:580088/150‑1 (MQ=255)
                                                                            ccGGAAGACTGGGTCAGGCGTTTTCCACGGTCATCCATGCAGGGCTCGTAGAACGGCGCACCCGCAGTCGGCAAGGCGCCGTTGGTGCGGAACGCGCTGGCCACTGGCTGGCTGCCGGGGATCAATGCATAGCTGGGGTGTTCACCCATg                                                                      >  2:1033584/1‑150 (MQ=255)
                                                                                   aCTGGGTCAGGCGTTTTCCAAGGTCATCCATGCAGGGCTCGTAGAACGGCGCACCCGCAGTCGGCAAGGCGCCGTTGGTGCGGAACGCGCTGGCCACTGGCTGGCTGCCGGGGATCAATGCATAGCTGGGGTGTTCACCCATGGCGTGGa                                                               >  1:996202/1‑150 (MQ=255)
                                                                                     tGGGTCAGGCGTTTTCCACGGTCATCCATGCAGGGCTCGTAGAACGGCGCACCCGCAGTCGGCAAGGCGCCGTTGGTGCGGAACGCGCTGGCCACTGGCTGGCTGCCGGGGATCAATGCATAGCTGGGGTGTTCACCCATGGCGTTGAAc                                                             >  2:620639/1‑150 (MQ=255)
                                                                                     tGGGTCAGGCGTTTTCCACGGTCATCCATGCAGGGCTCGTAGAACGGCGCACCCGCAGTCGGCAAGGCGCCGTTGGTGCGGAACGCGCTGGCCACTGGCTGGCTGCCGGGGATCAATGCATAGCTGGGGTGTTCACCCATGGCGTGGAAc                                                             <  1:1033584/150‑1 (MQ=255)
                                                                                            ggCGTTTTCCACGGTCATCCATGCAGGGCTCGTAGAACGGCGCACCCGCAGTCGGCAAGGCGCCGTTGGTGCGGAACGCGCTGGCCACTGGCTGGCTGCCGGGGATCAATGCATAGCTGGGGTGTTCACCCATGGCGTGGAACTGc                                                          >  2:1206742/1‑146 (MQ=255)
                                                                                            ggCGTTTTCCACGGTCATCCATGCAGGGCTCGTAGAACGGCGCACCCGCAGTCGGCAAGGCGCCGTTGGTGCGGAACGCGCTGGCCACTGGCTGGCTGCCGGGGATCAATGCATAGCTGGGGTGTTCACCCATGGCGTGGAACTGc                                                          <  1:1206742/146‑1 (MQ=255)
                                                                                                    ccACGGTCATCCATGCAGGGCTCGTAGAACGGCGCACCCGCAGTCGGCAAGGCGCCGTTGGTGCGGAACGCGCTGGCCACTGGCTGGCTGCCGGGGATCAATGCATAGCTGGGGTGTTCACCCATGGCGTGGAACTGCATGGctgcctgc                                              <  2:489031/150‑1 (MQ=255)
                                                                                                      aCGGTCATCCATGCAGGGCTCGTAGAACGGCGCACCCGCAGTCGGCAAGGCGCCGTTGGTGCGGAACGCGCTGGCCACTGGCTGGCTGCCGGGGATCAATGCATAGCTGGGGTGTTCACCCATGGCGTGGAACTGCATGGCTGCCTGCTc                                            <  1:480258/150‑1 (MQ=255)
                                                                                                         gTCATCCATGCAGGGCTCGTAGAACGGCGCACCCGCAGTCGGCAAGGCGCCGTTGGTGCGGAACGCGCTGGCCACTGGCTGGCTGCCGGGGATCAATGCATAGCTGGGGTGTTCACCCATGGCGTTGAACTGCATGGCTGCCTGCTCGAc                                         <  1:620639/150‑1 (MQ=255)
                                                                                                            aTCCATGCAGGGCTCGTAGAACGGCGCACCCGCAGTCGGCAAGGCGCCGTTGGTGCGGAACGCGCTGGCCACTGGCTGGCTGCCGGGGATCAATGCATAGCTGGGGTGTTCACCCATGGCGTGGAACTGCATGGCTGCCTGCTCGACATc                                      >  2:907729/1‑150 (MQ=255)
                                                                                                             tCCATGCAGGGCTCGTAGAACGGCGCACCCGCAGTCGGCAAGGCGCCGTTGGTGCGGAACGCGCTGGCCACTGGCTGGCTGCCGGGGATCAATGCATAGCTGGGGTGTTCACCCATGGCGTGTTATTGCATGGCTGCCTGCTCGACATCg                                     >  1:295/1‑150 (MQ=255)
                                                                                                              ccATGCAGGGCTCGTAGAACGGCGCACCCGCAGTCGGCAAGGCGCCGTTGGTGCGGAACGCGCTGGCCACTGGCTGGCTGCCGGGGATCAATGCATAGCTGGGGTGTTCACCCATGGCGTGGAACTGCATGGCTGCCTGCTCGACATCgg                                    >  2:971016/1‑150 (MQ=255)
                                                                                                               cATGCAGGGCTCGTAGAACGGCGCACCCGCAGTCGGCAAGGCGCCGTTGGTGCGGAACGCGCTGGCCACTGGCTGGCTGCCGGGGATCAATGCATAGCTGGGGTGTTCACCCATGGCGTGGAACTGCATGGCTGCCTGCTCGACATCGGt                                   >  2:142424/1‑150 (MQ=255)
                                                                                                                 tGCAGGGCTCGTAGAACGGCGCACCCGCAGTCGGCAAGGCGCCGTTGGTGCGGAACGCGCTGGCCACTGGCTGGCTGCCGGGGATCAATGCATAGCTGGGGTGTTCACCCATGGCGTGGAACTGCATGGCTGCCTGCTCGACATCGGTGc                                 <  2:996202/150‑1 (MQ=255)
                                                                                                                        cTCGTAGAACGGCGCACCCGCAGTCGGCAAGGCGCCGTTGGTGCGGAACGCGCTGGCCACTGGCTGGCTGCCGGGGATCAATGCATAGCTGGGGTGTTCACCCATGGCGTGGAACTGCATGGCTGCCTGCTCGACATCGGTGCCCTCTTc                          <  1:875952/150‑1 (MQ=255)
                                                                                                                        cTCGTAGAACGGCGCACCCGCAGTCGGCAAGGCGCCGTTGGTGCGGAACGCGCTGGCCACTGGCTGGCTGCCGGGGATCAATGCATAGCTGGGGTGTTCACCCATGGCGTGGAACTGCATGGCTGCCTGCTCGACATCGGTGCCCTCTTc                          <  1:876217/150‑1 (MQ=255)
                                                                                                                        cTCGTAGAACGGCGCACCCGCAGTCGGCAAGGCGCCGTTGGTGCGGAACGCGCTGGCCACTGGCTGGCTGCCGGGGATCAATGCATAGCTGGGGTGTTCACCCATGGCGTGGAACTGCATGGCTGCCTGCTCGACATCGGTGCCCTCTTc                          <  1:907729/150‑1 (MQ=255)
                                                                                                                         tCGTAGAACGGCGCACCCGCAGTCGGCAAGGCGCCGTTGGTGCGGAACGCGCTGGCCACTGGCTGGCTGCCGGGGATCAATGCATAGCTGGGGTGTTCACCCATGGCGTGGAACTGCATGGCTGCCTGCTCGACATCGGTGCCCTCTTcc                         <  1:884454/150‑1 (MQ=255)
                                                                                                                           gTAGAACGGCGCACCCGCAGTCGGCAAGGCGCCGTTGGTGCGGAACGCGCTGGCCACTGGCTGGCTGCCGGGGATCAATGCATAGCTGGGGTGTTCACCCATGGCGTGGAACTGCATGGctgcc                                                 >  1:364621/1‑124 (MQ=255)
                                                                                                                           gTAGAACGGCGCACCCGCAGTCGGCAAGGCGCCGTTGGTGCGGAACGCGCTGGCCACTGGCTGGCTGCCGGGGATCAATGCATAGCTGGGGTGTTCACCCATGGCGTGGAACTGCATGGctgcc                                                 <  2:364621/124‑1 (MQ=255)
                                                                                                                           gTAGAACGGCGCACCCGCAGTCGGCAAGGCGCCGTTGGTGCGGAACGCGCTGGCCACTGGCTGGCTGCCGGGGATCAATGCATAGCTGGGGTGTTCACCCATGGCGTGGAACTGCATGGCTGCCTGCTCGACATCGGTGCCCTCTTCCgg                       <  2:721913/150‑1 (MQ=255)
                                                                                                                           gTAGAACGGCGCACCCGCAGTCGGCAAGGCGCCGTTGGTGCGGAACGCGCTGGCCACTGGCTGGCTGCCGGGGATCAATGCATAGCTGGGGTGTTCACCCATGGCGTGGAACTGCATGGCTGCCTGCTCGACATCGGTGCCCTCTTCCgg                       >  1:1085469/1‑150 (MQ=255)
                                                                                                                           gTAGAACGGCGCACCCGCAGTCGGCAAGGCGCCGTTGGTGCGGAACGCGCTGGCAACTGGCTGGCTGCCGGGGATCAATGCATAGCTGGGGTGTTCACCCATGGCGTGGAACTGCATGGCTGCCTGCTCGACATCGGTGCCCTCTTCCgg                       >  1:869482/1‑150 (MQ=255)
                                                                                                                            tAGAACGGCGCACCCGCAGTCGGCAAGGCGCCGTTGGTGCGGAACGCGCTGGCCACTGGCTGGCTGCCGGGGATCAATGCATAGCTGGGGTGTTCACCCATGGCGTGGAACTGCATGGCTGCCTGCTCGACATCGGTGCCCTCTTCCGGc                      <  1:142424/150‑1 (MQ=255)
                                                                                                                             aGAACGGCGCACCCGCAGTCGGCAAGGCGCCGTTGGTGCGGAACGCGCTGGCCACTGGCTGGCTGCCGGGGATCAATGCATAGCTGGGGTGTTCACCCATGGCGTGGAACTGCATGGCTGCCTGCTCGACATCGGTGCCCTCTTCCggga                     >  1:150240/1‑148 (MQ=255)
                                                                                                                                   gcgcACCCGCAGTCGGCAAGGCGCCGTTGGTGCGGAACGCGCTGGCCACTGGCTGGCTGCCGGGGATCAATGCATAGCTGGGGTGTTCACCCATGGCGTGGAACTGCATGGCTGCCTGCTCGACATCGGTGCCCTCTTCCGGCAGGTAGa               >  2:78470/1‑150 (MQ=255)
                                                                                                                                           gCAGTCGGCAAGGCGCCGTTGGTGCGGAACGCGCTGGCCACTGGCTGGCTGCCGGGGATCAATGCATAGCTGGGGTGTTCACCCATGGCGTGGAACTGCATGGCTGCCTGCTCGACATCGGTGCCCTCTTCCGGCAGGTAGATGGGCTTg       >  2:751294/1‑150 (MQ=255)
                                                                                                                                                cGGCAAGGCGCCGTTGGTGCGGAACGCGCTGGCCACTGGCTGGCTGCCGGGGATCAATGCATAGCTGGGGTGTTCACCCATGGCGTGGAACTGCATGGCTGCCTGCTCGACATCGGTGCCCTCTTCCGGCAGGTAGATGGGCTTGGCCtt  <  1:78470/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GCGGTCGGCAGTGAAGGCCGAGGCGCCACGGAAGTTCACGCCCGTCAGGCTCTCTCCACTGTTGAACTCCCCAACCCCGGAAGACTGGGTCAGGCGTTTTCCACGGTCATCCATGCAGGGCTCGTAGAACGGCGCACCCGCAGTCGGCAAGGCGCCGTTGGTGCGGAACGCGCTGGCCACTGGCTGGCTGCCGGGGATCAATGCATAGCTGGGGTGTTCACCCATGGCGTGGAACTGCATGGCTGCCTGCTCGACATCGGTGCCCTCTTCCGGCAGGTAGATGGGCTTGGCCTT  >  NC_002947/3957849‑3958142

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑