Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,864,765 G→A 6.4% A462A (GCC→GCT)  cadA‑III ← cadmium translocating P‑type ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,864,7650GA 6.4% 187.3 / 2.1 62A462A (GCC→GCT) cadA‑IIIcadmium translocating P‑type ATPase
Reads supporting (aligned to +/- strand):  ref base G (28/30);  new base A (2/2);  total (30/32)
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.23e-01

AACGAAGTCACTGACCTCGCTCAAGGCCAGGCCTTGCTGCTTG‑CCGAACTCGGCAATGGCGCGGGACACCGGGTGGTCCGAACGCTCGCCCAGGCTGGCCGCCAAGGCTTGCGCACGGTCTTCGAACAGCGGCGCCAGCACCTTGGCGTCGGTCTGCACCGGCTTGCCATGGGTGATGGTGCCGGTCTTGTCCAGGGCCAGGAAATCAAGGTGGCGGCCACCCTCCAGGTACACGCCACCCTTGATCAGGATGCCTTTGCGCGCGGCGGCGGCCAGGCCACTGACGATGG  >  NC_002947/5864621‑5864910
                                                                                                                                                 |                                                                                                                                                 
aaCGAAGTCACTGACCTCGCTCAAGGCCAGGCCTTGCTGCTTG‑CCGAACTCGGCAATGGCGCGGGACACCGGGTGGTCCGAACGCTCGCCCAGGCTGGCCGCCAAGGCTTGCGCACGGTCTTCGAACAGCGGCGCCAGCACCTTGGCGt                                                                                                                                               >  3:36913/1‑149 (MQ=255)
     aGTCACTGACCTCGCTCAAGGCCAGGCCTTGCTGCTTG‑CCGAACTCGGCAATGGCGCGGGACACCGGGTGGTCCGAACGCTCGCCCAGGCTGGCCGCCAAGGCTTGCGCACGGTCTTCGAACAGCGGCGCCAGCACCTTGGCGTCGGTc                                                                                                                                          <  5:175647/149‑1 (MQ=255)
         aCTGACCTCGCTCAAGGCCAGGCCTTGCTGCTTGTCCGAACTCGGCAATGGCGCGGGACACCGGGTGGTCCGAACGCTCGCCCAGGCTGGCCGCCAAGGCTTGCGCACGGTCTTCGAACAGCGGCGCCAGCACCTTGGCGTCGGTCTg                                                                                                                                        >  1:211768/1‑148 (MQ=255)
         aCTGACCTCGCTCAAGGCCAGGCCTTGCTGCTTG‑CCGAACTCGGCAATGGCGCGGGACACCGGGTGGTCCGAACGCTCGCCCAGGCTGGCCGCCAAGGCTTGCGCACGGTCTTCGAACAGCGGCGCCAGCACCTTGGCGTCGGt                                                                                                                                           >  4:209430/1‑144 (MQ=255)
         aCTGACCTCGCTCAAGGCCAGGCCTTGCTGCTTG‑CCGAACTCGGCAATGGCGCGGGACACCGGGTGGTCCGAACGCTCGCCCAGGCTGGCCGCCAAGGCTTGCGCACGGTCTTCGAACAGCGGCGCCAGCACCTTGGCGTCGGt                                                                                                                                           <  3:209430/144‑1 (MQ=255)
           tGACCTCGCTCAAGGCCAGGCCTTGCTGCTTG‑CCGAACTCGGCAATGGCGCGGGACACCGGGTGGTCCGAACGCTCGCCCAGGCTGGCCGCCAAGGCTTGTGCACGGTCTTCGAACAGCGGCGCCAGCACCTTGGCGTCGGTCTGCAcc                                                                                                                                    >  7:25388/1‑149 (MQ=255)
           tGACCTCGCTCAAGGCCAGGCCTTGCTGCTTG‑CCGAACTCGGCAATGGCGCGGGACACCGGGTGGTCCGAACGCTCGCCCAGGCTGGCCGCCAAGGCTTGCGCACGGTCTTCGAACAGCGGCGCCAGCACCTTGGCGTCGGTCTGCAcc                                                                                                                                    >  2:210472/1‑149 (MQ=255)
             aCCTCGCTCAAGGCCAGGCCTTGCTGCTTG‑CCGAACTCGGCAATGGCGCGGGACACCGGGTGGTCCGAACGCTCGCCCAGGCTGGCCGCCAAGGCTTGCGCACGGTCTTCGAACAGCGGCGCCAGCACCTTGGCGTCGGTCTGCACCg                                                                                                                                   <  8:1607/148‑1 (MQ=255)
                        ggccaggccTTGCTGCTTG‑CCGAACTCGGCAATGGCGCGGGACACCGGGTGGTCCGAACGCTCGCCCAGGCTGGCCGCCAAGGCTTGCGCACGGTCTTCGAACAGCGGCGCCAGCACCTTGGCGTCGGTCTGCACCGGCTTGCCAt                                                                                                                          >  2:73348/1‑146 (MQ=255)
                        ggccaggccTTGCTGCTTG‑CCGAACTCGGCAATGGCGCGGGACACCGGGTGGTCCGAACGCTCGCCCAGGCTGGCCGCCAAGGCTTGCGCACGGTCTTCGAACAGCGGCGCCAGCACCTTGGCGTCGGTCTGCACCGGCTTGCCAt                                                                                                                          <  1:73348/146‑1 (MQ=255)
                        ggccaggccTTGCTGCTTG‑CCGAACTCGGCAATGGCGC‑GGACACCGGGTGGTCCGAACGCTCGCCCAGGCTGGCCGCCAAGGCTTGCGCACGGTCTTCGAACAGCGGCGCCAGCACCTTGGCGTCGGTCTGCACCGGCTTGCCATGGGt                                                                                                                      >  6:263590/1‑149 (MQ=255)
                              gccTTGCTGCTTG‑CCGAACTCGGCAATGGCGCGGGACACCGGGTGGTCCGAACGCTCGCCCAGGCTGGCCGCCAAGGCTTGCGCACGGTCTTCGAACAGCGGCGCCAGCACCTTGGCGTCGGTCTGCACCGGCTTGCCATGGGTGATgg                                                                                                                 <  4:36913/149‑1 (MQ=255)
                                   gctgctTG‑CCGAACTCGGCAATGGCGCGGGACACCGGGTGGTCCGAACGCTCGCCCAGGCTGGCCGCCAAGGCTTGTGCACGGTCTTCGAACAGCGGCGCCAGCACCTTGGCGTCGGTCTGCACCGGCTTGCCATGGGTGATGGTGc                                                                                                              <  8:25388/147‑1 (MQ=255)
                                       ctTG‑CCGAACTCGGCAATGGCGCGGGACACCGGGTGGTCCGAACGCTCGCCCAGGCTGGCCGCCAAGGCTTGCGCACGGTCTTCGAACAGCGGCGCCAGCACCTTGGCGTCGGTCTGCACCGGCTTGCCATTGGTGATGGTGCCGGTc                                                                                                         <  8:146365/148‑1 (MQ=255)
                                       ctTG‑CCGAACTCGGCAATGGCGCGGGACACCGGGTGGTCCGAACGCTCGCCCAGGCTGGCCGCCAAGGCTTGCGCACGGTCTTCGAACAGCGGCGCCAGCACCTTGGCGTCGGTCTGCACCGGCTTGCCATGGGTGATGGTGCCGGTCt                                                                                                        >  5:52422/1‑149 (MQ=255)
                                       ctTG‑CCGAACTCGGCAATGGCGC‑GGACACCGGGTGGTCCGAACGCTCGCCCAGGCTGGCCGCCAAGGCTTGCGCACGGTCTTCGAACAGCGGCGCCAGCACCTTGGCGTCGGTCTGCACCGGCTTGCCATGGGTGAAGGTGCCGGTc                                                                                                         <  5:263590/147‑1 (MQ=255)
                                          g‑CCGAACTCGGCAATGGCGCGGGACACCGGGTGGTCCGAACGCTCGCCCAGGCTGGCCGCCAAGGCTTGCGCACGGTCTTCGAACAGCGGCGCCAGCACCTTGGCGTCGGTCTGCACCGGCTTGCCATGGGTGATGGTGCCGGTCTTg                                                                                                      <  6:128577/148‑1 (MQ=255)
                                                aCTCGGCAATGGCGCGGGACACCGGGTGGTCCGAACGCTCGCCCAGGCTGGCCGCCAAGGCTTGCGCACGGTCTTCGAACAGCGGCGCCAGCACCTTGGCGTCGGTCTGCACCGGCTTGCCATGGGTGATGGTGCCGGTCTTGTCCAgg                                                                                                >  2:47713/1‑149 (MQ=255)
                                                  tCGGCAATGGCGCGGGACACCGGGTGGTCCGAACGCTCGCCCAGGCTGGCCGCCAAGGCTTGCGCACGGTCTTCGAACAGCGGCGCCAGCACCTTGGCGTCGGTCTGCACCGGCTTGCCATGGGTGATGGTGCCGGTc                                                                                                         <  4:135508/138‑1 (MQ=255)
                                                  tCGGCAATGGCGCGGGACACCGGGTGGTCCGAACGCTCGCCCAGGCTGGCCGCCAAGGCTTGCGCACGGTCTTCGAACAGCGGCGCCAGCACCTTGGCGTCGGTCTGCACCGGCTTGCCATGGGTGATGGTGCCGGTc                                                                                                         >  3:135508/1‑138 (MQ=255)
                                                  tCGGCAATGGCGCGGGACACCGGGTGGTCCGAACGCTCGCCCAGGCTGGCCGCCAAGGCTTGCGCACGGTCTTCGAACAGCGGCGCCAGCACCTTGGCGTCGGTCTGCACCGGCTTGCCATGGGTGATGGTGCCGGTCTTGTCCAGGGc                                                                                              >  7:185052/1‑149 (MQ=255)
                                                   cGGCAATGGCGCGGGACACCGGGTGGTCCGAACGCTCGCCCAGGCTGGCCGCCAAGGCTTGCGCACGGTCTTCGAACAGCGGCGCCAGCACCTTGGCGTCGGTCTGCACCGGCTTGCCATGGGTGATGGTGCCGGTCTTGTCCAGGGc                                                                                              >  5:334390/1‑148 (MQ=255)
                                                     gCAATGGCGCGGGACACCGGGTGGTCCGAACGCTCGCCCAGGCTGGCCGCCAAGGCTTGCGCACGGTCTTCGAACAGCGGCGCCAGCACCTTGGCGTCGGTCTGCACCGGCTTGCCATGGGTGATGGTGCCGGTCTTGTCCAGGGCCAg                                                                                           <  6:52422/149‑1 (MQ=255)
                                                      cAATGGCGCGGGACACCGGGTGGTCCGAACGCTCGCCCAGGCTGGCCGCCAAGGCTTGCGCACGGTCTTCGAACAGCGGCGCCAGCACCTTGGCGTCGGTCTGCACCGGCTTGCCATGGGTGATGGTGCCGGTCTTGTCCAGGGCCAg                                                                                           >  6:167777/1‑148 (MQ=255)
                                                       aaTGGCGCGGGCCACCGGGTGGTCCGAACGCTCGCCCAGGCTGGCCGCCAAGGCTTGCGCACGGTCTTCGAACAGCGGCGCCAGCACCTTGGCGTCGGTCTGCACCGGCTTGCCATGGGTGATGGTGCCGGTCtt                                                                                                       <  7:232887/135‑1 (MQ=255)
                                                        aTGGCGCGGGACACCGGGTGGTCCGAACGCTCGCCCAGGCTGGCCGCCAAGGCTTGCGCACGGTCTTCGAACAGCGGCGCCAGCACCTTGGCGTCGGTCTGCACCGGCTTGCCATGGGTGATGGTGCCGGTCTTGTCCAGGGCCAGGaa                                                                                        <  6:334390/149‑1 (MQ=255)
                                                         tGGCGCGGGACACCGGGTGGTCCGAACGCTCGCCCAGGCTGGCCGCCAAGGCTTGCGCACGGTCTTCGAACAGCGGCGCCAGCACCTTGGCGTCGGTCTGCACCGGCTTGCCATGGGTGATGGTGCCGGTCTTGTCCAGGGCCAGGaaa                                                                                       <  8:185052/149‑1 (MQ=255)
                                                           gcgcGGGACACCGGGTGGTCCGAACGCTCGCCCAGGCTGGCCGCCAAGGCTTGCGCACGGTCTTCGAACAGCGGCGCCAGCACCTTGGCGTCGGTCTGCACCGGCTTGCCATGGGTGATGGTGCCGGTCTTGTCCAGGGCCAGGAAATc                                                                                     <  1:210472/149‑1 (MQ=255)
                                                                ggACACCGGGTGGTCCGAACGCTCGCCCAGGCTGGCCGCCAAGGCTTGCGCACGGTCTTCGAACAGCGGCGCCAGCACCTTGGCGTCGGTCTGCACCGGCTTGCCATGGGTGATGGTGCCGGTCTTGTCCAGGGCCAGGAAATCACGGt                                                                                >  2:163270/1‑149 (MQ=255)
                                                                    acCGGGTGGTCCGAACGCTCGCCCAGGCTGGCCGCCAAGGCTTGCGCACGGTCTTCGAACAGCGGCGCCAGCACCTTAGCGTCGGTCTGCACCGGCTTGCCATGGGTGATGGTGCCGGTCTTGTCCAGGGCCAGGAAATCAAGGTggcg                                                                            >  3:80866/1‑149 (MQ=255)
                                                                    acCGGGTGGTCCGAACGCTCGCCCAGGCTGGCCGCCAAGGCTTGCGCACGGTCTTCGAACAGCGGCGCCAGCACCTTAGCGTCGGTCTGCACCGGCTTGCCATGGGTGATGGTGCCGGTCTTGTCCAGGGCCAGGAAATCAAGGTggcg                                                                            >  3:76466/1‑149 (MQ=255)
                                                                       gggTGGTCCGAACGCTCGCCCAGGCTGGCCGCCAAGGCTTGCGCACGGTCTTCGAACAGCGGCGCCAGCACCTTGGCGTCGGTCTGCACCGGCTTGCCATGGGTGATGGTGCCGGTCTTGTCCAGGGCCAGGAAATCAAGGTGGCGGcc                                                                         <  4:51170/149‑1 (MQ=255)
                                                                           ggtCCGAACGCTCGCCCAGGCTGGCCGCCAAGGCTTGCGCACGGTCTTCGAACAGCGGCGCCAGCACCTTGGCGTCGGTCTGCACCGGCTTGCCATGGGTGATGGTGCCGGTCTTGTCCAGGGCCAGGAAATCAAGGTGGCGGCCAccc                                                                     >  5:223520/1‑149 (MQ=255)
                                                                               cGAACGCTCGCCCAGGCTGGCCGCCAAGGCTTGCGCACGGTCTTCGAACAGCGGCGCCAGCACCTTGGCGTCGGTCTGCACCGGCTTGCCATGGGTGATGGTGCCGGTCTTGTCCAGGGCCAGGAAATCAAGGTGGCGGCCACCCTccc                                                                 >  5:257969/1‑148 (MQ=255)
                                                                               cGAACGCTCGCCCAGGCTGGCCGCCAAGGCTTGCGCACGGTCTTCGAACAGCGGCGCCAGCACCTTGGCGTCGGTCTGCACCGGCTTGCCATGGGTGATGGTGCCGGTCTTGTCCAGGGCCAGGAAATCAAGGTGGCGGCCACCCTcc                                                                  <  3:89765/148‑1 (MQ=255)
                                                                                gAACGCTCGCCCAGGCTGGCCGCCAAGGCTTGCGCACGGTCTTCGAACAGCGGCGCCAGCACCTTGGCGTCGGTCTGCACCGGCTTGCCATGGGTGATGGTGCCGGTCTTGTCCAGGGCCAGGAAATCAAGGTGGCGGCCACCCTCCAg                                                                <  3:309016/149‑1 (MQ=255)
                                                                                gAACGCTCGCCCAGGCTGGCCGCCAAGGCTTGCGCACGGTCTTCGAACAGCGGCGCCAGCACCTTGGCGTCGGTCTGCACCGGCTTGCCATGGGTGATGGTGCCGGTCTTGTCCAGGGCCAGGAAATCAAGGTGGCGGCCACCCTCCAg                                                                <  3:313986/149‑1 (MQ=255)
                                                                                 aaCGCTCGCCCAGGCTGGCCGCCAAGGCTTGCGCACGGTCTTCGAACAGCGGCGCCAGCACCTTGGCGTCGGTCTGCACCGGCTTGCCATGGGTGATGGTGCCGGTCTTGTCCAGGGCCAGGAAATCAAGGTGGCGGCCACCCTCCa                                                                 >  4:28396/1‑147 (MQ=255)
                                                                                 aaCGCTCGCCCAGGCTGGCCGCCAAGGCTTGCGCACGGTCTTCGAACAGCGGCGCCAGCACCTTGGCGTCGGTCTGCACCGGCTTGCCATGGGTGATGGTGCCGGTCTTGTCCAGGGCCAGGAAATCAAGGTGGCGGCCACCCTCCa                                                                 <  3:28396/147‑1 (MQ=255)
                                                                                 aaCGCTCGCCCAGGCTGGCCGCCAAGGCTTGCGCACGGTCTTCGAACAGCGGCGCCAGCACCTTGGCGTCGGTCTGCACCGGCTTGCCATGGGTGATGGTGCCGGTCTTGTCCAGGGCCAGGAAATCAAGGTGGCGGCCACCCTCCAgg                                                               >  2:188031/1‑149 (MQ=255)
                                                                                         cccAGGCTGGCCGCCAAGGCTTGCGCACGGTCTTCGAACAGCGGCGCCAGCACCTTGGCGTCGGTCTGCACCGGCTTGCCATGGGTGATGGTGCCGGTCTTGTCCAGGGCCAGGAAATCAAGGTGGCGGCCACCCTCCAGGTACACGcc                                                       >  3:63974/1‑149 (MQ=255)
                                                                                                  gccgccAAGGCTTGCGCACGGTCTTCGAACAGCGGCGCCAGCACCTTAGCGTCGGTCTGCACCGGCTTGCCATGGGTGATGGTGCCGGTCTTGTCCAGGGCCAGGAAATCAAGGTGGCGGCCACCCTCCAGGTACACGCCACCCTTg                                                <  4:80866/147‑1 (MQ=255)
                                                                                                  gccgccAAGGCTTGCGCACGGTCTTCGAACAGCGGCGCCAGCACCTTAGCGTCGGTCTGCACCGGCTTGCCATGGGTGATGGTGCCGGTCTTGTCCAGGGCCAGGAAATCAAGGTGGCGGCCACCCTCCAGGTACACGCCACCCTTg                                                <  4:76466/147‑1 (MQ=255)
                                                                                                    cgccAAGGCTTGCGCACGGTCTTCGAACAGCGGCGCCAGCACCTTGGCGTCGGTCTGCACCGGCTTGCCATGGGTGATGGTGCCGGTCTTGTCCAGGGCCAGGAAATCAAGGTGGCGGCCACCCTCCAGGTACACGCCACCCTTGATCa                                            <  1:163270/149‑1 (MQ=255)
                                                                                                       cAAGGCTTGCGCACGGTCTTCGAACAGCGGCGCCAGCACCTTGGCGTCGGTCTGCACCGGCTTGCCATGGGTGATGGTGCCGGTCTTGTCCAGGGCCAGGAAATCAAGGTGGCGGCCACCCTCCAGGTACACGCCACCCTTGATCAGGa                                         <  6:54445/149‑1 (MQ=255)
                                                                                                        aaGGCTTGCGCACGGTCTTCGAACAGCGGCGCCAGCACCTTGGCGTCGGTCTGCACCGGCTTGCCATGGGTGATGGTGCCGGTCTTGTCCAGGGCCAGGAAATCAAGGTGGCGGCCACCCTCCAGGTACACGCCACCCTTGATCAGGAt                                        <  2:211768/149‑1 (MQ=255)
                                                                                                                cgcACGGTCTTCGAACAGCGGCGCCAGCACCTTGGCGTCGGTCTGCACCGGCTTGCCATGGGTGATGGTGCCGGTCTTGTCCAGGGCCAGGAAATCAAGGTGGCGGCCACCCTCCAGGTACACGCCACCCTTGATCAGGATGCCTTTgc                                <  1:188031/149‑1 (MQ=255)
                                                                                                                  cACGGTCTTCGAACAGCGGCGCCAGCACCTTGGCGTCGGTCTGCACCGGCTTGCCATGGGTGATGGTGCCGGTCTTGTCCAGGGCCAGGAAATCAAGGTGGCGGCCACCCTCCAGGTACACGCCACCCTTGATCAGGATGCCTTTgcgc                              >  1:329480/1‑149 (MQ=255)
                                                                                                                       tCTTCGAACAGCGGCGCCAGCACCTTGGCGTCGGTCTGCACCGGCTTGCCATGGGTGATGGTGCCGGTCTTGTCCAGGGCCAGGAAATCAAGGTGGCGGCCACCCTCCAGGTACACGCCACCCTTGATCAGGATGCCTTTGCGCgcggc                         <  4:63974/149‑1 (MQ=255)
                                                                                                                           cGAACAGCGGCGCCAGCACCTTGGCGTCGGTCTGCACCGGCTTGCCATGGGTGATGGTGCCGGTCTTGTCCAGGGCCAGGAAATCAAGGTGGCGGCCACCCTCCAGGTACACGCCACCCTTGATCAGGa                                         >  4:276008/1‑129 (MQ=255)
                                                                                                                           cGAACAGCGGCGCCAGCACCTTGGCGTCGGTCTGCACCGGCTTGCCATGGGTGATGGTGCCGGTCTTGTCCAGGGCCAGGAAATCAAGGTGGCGGCCACCCTCCAGGTACACGCCACCCTTGATCAGGa                                         <  3:276008/129‑1 (MQ=255)
                                                                                                                           cGAACAGCGGCGCCAGCACCTTGGCGTCGGTCTGCACCGGCTTGCCATGGGTGATGGTGCCGGTCTTGTCCAGGGCCAGGAAATCAAGGTGGCGGCCACCCTCCAGGTACACGCCACCCTTGATCAGGATGCCTTTGCGCgcggcggcg                     >  3:3259/1‑149 (MQ=255)
                                                                                                                              aCAGCGGCGCCAGCACCTTGGCGTCGGTCTGCACCGGCTTGCCATGGGTGATGGTGCCGGTCTTGTCCAGGGCCAGGAAATCAAGGTGGCGGCCACCCTCCAGGTACACGCCACCCTTGATCAGGATGCCTTTGCGCGCggcggcggc                   >  1:170570/1‑148 (MQ=255)
                                                                                                                                 gcggcgCCAGCACCTTGGCGTCGGTCTGCACCGGCTTGCCATGGGTGATGGTGCCGGTCTTGTCCAGGGCCAGGAAATCAAGGTGGCGGCCACCCTCCAGGTACACGCCACCCTTGATCAGGATGCCTTTGCGCGCGGCGGCggccagg               >  1:202611/1‑149 (MQ=255)
                                                                                                                                    gcgcCAGCACCTTGGCGTCGGTCTGCACCGGCTTGCCATGGGTGATGGTGCCGGTCTTGTCCAGGGCCAGGAAATCAAGGTGGCGGCCACCCTCCAGGTACACGCCACCCTTGATCAGGATGCCTTTGCGCGCGGCGGCggccaggcca            >  1:194911/1‑149 (MQ=255)
                                                                                                                                     cgcCAGCACCTTGGCGTCGGTCTGCACCGGCTTGCCATGGGTGATGGTGCCGGTCTTGTCCAGGGCCAGGAAATCAAGGTGGCGGCCACCCTCCAGGTACACGCCACCCTTGATCAGGATGCCTTTGCGCGCGGCGGCGGCCAGGCCAc           <  6:223520/149‑1 (MQ=255)
                                                                                                                                       ccAGCACCTTGGCGTCGGTCTGCACCGGCTTGCCATGGGTGATGGTGCCGGTCTTGTCCAGGGCCAGGAAATCAAGGTGGCGGCCACCCTCCAGGTACACGCCACCCTTGATCAGGATGCCTTTGCGCgcggcggcg                     >  7:145172/1‑137 (MQ=255)
                                                                                                                                       ccAGCACCTTGGCGTCGGTCTGCACCGGCTTGCCATGGGTGATGGTGCCGGTCTTGTCCAGGGCCAGGAAATCAAGGTGGCGGCCACCCTCCAGGTACACGCCACCCTTGATCAGGATGCCTTTGCGCgcggcggcg                     <  8:145172/137‑1 (MQ=255)
                                                                                                                                       ccAGCACCTTGGCGGCGGTCTGCACCGGCTTGCCATGGGTGATGGTGCCGGTCTTGTCCAGGGCCAGGAAATCAAGGTGGCGGCCACCCTCCAGGTACACGCCACCCTTGATCAGGATGCCTTTGCGCGCGGCGGCGGCCAGGCCACTg         <  8:306662/149‑1 (MQ=255)
                                                                                                                                        cAGCACCTTGGCGTCGGTCTGCACCGGCTTGCCATGGGTGATGGTGCCGGTCTTGTCCAGGGCCAGGAAATCAAGGTGGCGGCCACCCTCCAGGTACACGCCACCCTTGATCAGGATGCCTTTGCGCGCGGCGGCGGCCAGGCCACTGa        >  7:127312/1‑149 (MQ=255)
                                                                                                                                             ccTTGGCGTCGGTCTGCACCGGCTTGCCATGGGTGATGGTGCCGGTCTTGTCCAGGGCCAGGAAATCAAGGTGGCGGCCACCCTCCAGGTACACGCCACCCTTGATCAGGATGCCTTTGCGCGCGGCGGCGGCCAGGCCACTGACGATg   <  1:47713/149‑1 (MQ=255)
                                                                                                                                              cTTGGCGTCGGTCTGCACCGGCTTGCCATGGGTGATGGTGCCGGTCTTGTCCAGGGCCAGGAAATCAAGGTGGCGGCCACCCTCCAGGTACACGCCACCCTTGATCAGGATGCCTTTGCGCGCGGCGGCGGCCAGGCCACTGACGATgg  >  5:295694/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
AACGAAGTCACTGACCTCGCTCAAGGCCAGGCCTTGCTGCTTG‑CCGAACTCGGCAATGGCGCGGGACACCGGGTGGTCCGAACGCTCGCCCAGGCTGGCCGCCAAGGCTTGCGCACGGTCTTCGAACAGCGGCGCCAGCACCTTGGCGTCGGTCTGCACCGGCTTGCCATGGGTGATGGTGCCGGTCTTGTCCAGGGCCAGGAAATCAAGGTGGCGGCCACCCTCCAGGTACACGCCACCCTTGATCAGGATGCCTTTGCGCGCGGCGGCGGCCAGGCCACTGACGATGG  >  NC_002947/5864621‑5864910

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