Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 684,844 Δ1 bp 6.3% coding (1406/2400 nt) cadA‑II ← cadmium translocating P‑type ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947684,8440T. 6.3% 191.6 / 7.1 63coding (1406/2400 nt)cadA‑IIcadmium translocating P‑type ATPase
Reads supporting (aligned to +/- strand):  ref base T (24/34);  new base . (2/2);  total (27/36)
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.47e-01

GGGCGCCGGCCAGGCGGTAAAGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTACCCCAGCCAGCCACCA  >  NC_002947/684698‑684992
                                                                                                                                                  |                                                                                                                                                     
gggCGCCGGCCAGGCGGTAAAGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGa                                                                                                                                                    <  2:706833/150‑1 (MQ=255)
gggCGCCGGCCAGGCGGTAAAGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGGGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGa                                                                                                                                                    >  2:933781/1‑150 (MQ=255)
 ggCGCCGGCCAGGCGGTAAAGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGAt                                                                                                                                                   <  2:57128/150‑1 (MQ=255)
            ggCGGTAAAGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGt                                                                                                                                       <  1:639390/150‑1 (MQ=255)
             gCGGTAAAGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTg                                                                                                                                      >  1:686106/1‑150 (MQ=255)
                      gTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGggcggc                                                                                                                             <  1:904625/150‑1 (MQ=255)
                        cGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGa                                                                                                                           <  1:253281/150‑1 (MQ=255)
                                        cAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGTCCGCATGGGCACGCTCCAGGGCTTCGGATTCCTTGATCAGGAAGCCGTGGCGGGCGGCGAACCCGGTACCGGGCAt                                                                                                           <  2:496531/150‑1 (MQ=255)
                                        cAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCAt                                                                                                           <  1:378206/150‑1 (MQ=255)
                                            gcCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATc                                                                                                       >  2:487277/1‑150 (MQ=255)
                                                 gCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACGGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCgcggc                                                                                                  <  2:208670/150‑1 (MQ=255)
                                                 gCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCgcggc                                                                                                  <  2:208931/150‑1 (MQ=255)
                                                   tgtgAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAgg                                                                                                                                              >  1:638013/1‑104 (MQ=255)
                                                   tgtgAACCACCTGCGGGCTGCCGGAGGTGAGGGTGACGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAgg                                                                                                                                              <  2:638013/104‑1 (MQ=255)
                                                     tgAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCggcgggc                                                                                              >  1:459525/1‑150 (MQ=255)
                                                       aaCCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCGCGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGt                                                                                            >  1:207861/1‑150 (MQ=255)
                                                        cccaccTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTg                                                                                           <  2:1013323/149‑1 (MQ=255)
                                                            ccTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCa                                                                                       <  2:213454/150‑1 (MQ=255)
                                                             cTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAg                                                                                      <  1:352846/150‑1 (MQ=255)
                                                             cTGCGGGCTGCAGGCGGTGGGGGTGCCGGGTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCGTGATCAGGATGCCGTGGCGGGCGGCGACCACGGTACCCGCCATGATCGCGGCGGGCGTGGCCac                                                                                      >  2:272060/1‑149 (MQ=255)
                                                              tGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGAccc                                                                                                                        <  2:1075234/115‑1 (MQ=255)
                                                              tGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGAccc                                                                                                                        >  1:1075234/1‑115 (MQ=255)
                                                                    cTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAgc                                                                               <  1:487277/150‑2 (MQ=255)
                                                                     tGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAggg                                                                              >  1:268046/1‑150 (MQ=255)
                                                                         ggAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTTGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGgcgaa                                                                          >  1:538077/1‑148 (MQ=255)
                                                                          gAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAg                                                                         <  2:246964/150‑1 (MQ=255)
                                                                                aGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAg                                                                   <  2:908043/150‑1 (MQ=255)
                                                                                 gggTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAgg                                                                  <  2:330944/150‑1 (MQ=255)
                                                                                 gggTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAgg                                                                  <  1:1084144/150‑1 (MQ=255)
                                                                                  ggTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGc                                                                 <  1:73234/150‑1 (MQ=255)
                                                                                     gCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGAt                                                              >  2:168455/1‑150 (MQ=255)
                                                                                      ccGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGAcccc                                                                                                                       <  2:887148/92‑1 (MQ=255)
                                                                                      ccGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGAcccc                                                                                                                       >  1:887148/1‑92 (MQ=255)
                                                                                         gTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGAcc                                                          <  2:107859/150‑1 (MQ=255)
                                                                                           tttGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAg                                                        <  1:933781/150‑1 (MQ=255)
                                                                                               tCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCacg                                                    >  1:361373/1‑150 (MQ=255)
                                                                                                         cgcgGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCCGCACGGCGACGGCGt                                          >  1:763921/1‑150 (MQ=255)
                                                                                                          gcgGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCTGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGtt                                         <  2:686106/150‑1 (MQ=255)
                                                                                                            ggTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGGGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGa                                       >  1:541810/1‑150 (MQ=255)
                                                                                                             gTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGTCCATGACCAGCACGGCGACGGCGTTGAt                                      <  1:869572/150‑1 (MQ=255)
                                                                                                              ttGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCgg                                                                                                               <  2:885263/76‑1 (MQ=255)
                                                                                                              ttGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCgg                                                                                                               >  1:885263/1‑76 (MQ=255)
                                                                                                                  ccGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGGTCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAgcg                                 >  2:1103340/1‑150 (MQ=255)
                                                                                                                     cATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGTATGCCGTGGCGGGCGGCGACCCCGGTACCCGCCCTGATGGCGGCGGGCGTGGGCAGGCCGAGGGAGCAGGGGCAGGCGATGACCAGCACGGCGCCGGCCTTGATCAGCGCgg                              >  1:719918/1‑150 (MQ=255)
                                                                                                                     cATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCgg                              >  2:255146/1‑150 (MQ=255)
                                                                                                                       gggggACGCGTCAGGGCGTCGGCAGCCTTGAGCAGGATGCCGTGGGGGGCGGCGCGCCCGGTACCGGGCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCCGCGCGGtt                            <  2:507611/149‑1 (MQ=255)
                                                                                                                         ggCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTc                          <  1:972707/150‑1 (MQ=255)
                                                                                                                          gCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTcc                         <  2:252978/150‑1 (MQ=255)
                                                                                                                              gCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCgg                              <  2:707877/141‑1 (MQ=255)
                                                                                                                              gCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCgg                              >  1:707877/1‑141 (MQ=255)
                                                                                                                                  cAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTAc                 <  1:594821/150‑1 (MQ=255)
                                                                                                                                   aGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTAcc                >  2:537898/1‑150 (MQ=255)
                                                                                                                                      gCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGc                                                                 >  2:1026493/1‑98 (MQ=255)
                                                                                                                                      gCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGc                                                                 <  1:1026493/98‑1 (MQ=255)
                                                                                                                                       ccacGGCATCCTGATCAAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCacggcga                                                <  2:559893/111‑1 (MQ=255)
                                                                                                                                       ccacGGCATCCTGATCAAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCacggcga                                                >  1:559893/4‑114 (MQ=255)
                                                                                                                                       cTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTACCccag            >  1:835251/1‑150 (MQ=255)
                                                                                                                                        ttCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTACCccagc           >  2:816251/1‑150 (MQ=255)
                                                                                                                                           ggCATCCTGATCAAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGACCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGtt                            >  1:940725/1‑130 (MQ=255)
                                                                                                                                           ggCATCCTGATCAAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGACCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGtt                            <  2:940725/130‑1 (MQ=255)
                                                                                                                                              aTCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTACCCCAgccagcca     <  1:537898/150‑1 (MQ=255)
                                                                                                                                               tCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTACCCCAGccagccac    <  1:816251/150‑1 (MQ=255)
                                                                                                                                                 cTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTATCCAGCGGTATCCCAGtcagactcca  >  2:1085994/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GGGCGCCGGCCAGGCGGTAAAGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTACCCCAGCCAGCCACCA  >  NC_002947/684698‑684992

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