Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 684,844 Δ1 bp 5.5% 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. 5.5% 286.5 / 9.2 91coding (1406/2400 nt)cadA‑IIcadmium translocating P‑type ATPase
Reads supporting (aligned to +/- strand):  ref base T (44/42);  new base . (2/3);  total (46/45)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

AGGGCGCCGGCCAGGCGGTAAAGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTACCCCAGCCAGCCACCAG  >  NC_002947/684697‑684993
                                                                                                                                                   |                                                                                                                                                       
aGGGCGCCGGCCAGGCGGTAAAGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTg                                                                                                                                                       <  2:1326175/150‑1 (MQ=255)
 gggCGCCGGCCAGGCGGTAAAGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGa                                                                                                                                                      <  2:699324/150‑1 (MQ=255)
 gggCGCCGGCCAGGCGGTAAAGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGATGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGa                                                                                                                                                      >  1:1067521/1‑150 (MQ=255)
       cGGCCAGGCGGTAAAGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGa                                                                                                                                               <  2:629749/150‑1 (MQ=255)
          ccAGGCGGTAAAGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGc                                                                                                                                            <  1:1165985/150‑1 (MQ=255)
           cAGGCGGTAAAGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTGATCAAGGATGcc                                                                                                                                           <  2:349596/150‑1 (MQ=255)
           cAGGCGGTAAAGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTGATCAAGGATGcc                                                                                                                                           <  1:1089558/150‑1 (MQ=255)
            aGGCGGTAAAGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCg                                                                                                                                          >  1:1012771/1‑150 (MQ=255)
              gCGGTAAAGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTg                                                                                                                                        <  2:426484/150‑1 (MQ=255)
              gCGGTAAAGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTg                                                                                                                                        <  2:32227/150‑1 (MQ=255)
                       gTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGggcggc                                                                                                                               >  2:95517/1‑150 (MQ=255)
                       gTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGggcggc                                                                                                                               >  2:1551531/1‑150 (MQ=255)
                       gTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGggcggc                                                                                                                               <  1:1480318/150‑1 (MQ=255)
                                     cGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCgg                                                                                                                 <  1:540600/150‑1 (MQ=255)
                                       tCCAGCGCCTGGCTGTGAACCACATGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCATGATGCCGTGGCGGGCGGCGACCCCGGTACCGGcc                                                                                                               >  1:450189/1‑150 (MQ=255)
                                         cAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGACACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCAt                                                                                                             <  1:1541881/150‑1 (MQ=255)
                                          aGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGGGGCGACCCCGGTACCGGCCATg                                                                                                            >  2:26992/1‑150 (MQ=255)
                                           gcgcCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGa                                                                                                           <  2:964337/150‑1 (MQ=255)
                                              ccTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATcg                                                                                                        >  1:1582609/1‑150 (MQ=255)
                                               cTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATcgc                                                                                                       >  1:705398/1‑150 (MQ=255)
                                                  gCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGGGGCGGGCGAATGCCCCGGTACCGGCCATGATCgcggc                                                                                                    >  1:1282/1‑150 (MQ=255)
                                                    tgtgAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGCCCGCATGGGCAGCCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGcggcgg                                                                                                  <  1:1176145/150‑1 (MQ=255)
                                                    tgtgAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGcggcgg                                                                                                  <  1:1775781/150‑1 (MQ=255)
                                                      tgAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATc                                                                                                                                                    >  2:1719000/1‑99 (MQ=255)
                                                      tgAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATc                                                                                                                                                    <  1:1719000/99‑1 (MQ=255)
                                                        aaCCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCTTTGGAATCCGCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGt                                                                                              <  1:798425/150‑1 (MQ=255)
                                                        aaCCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGt                                                                                              >  1:110917/1‑150 (MQ=255)
                                                         accaccTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTg                                                                                             >  1:1866797/1‑150 (MQ=255)
                                                           caccTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGc                                                                                           <  1:882894/150‑1 (MQ=255)
                                                            accTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCTTGATCGCGGCGGGCGTGGcc                                                                                          >  1:169089/1‑150 (MQ=255)
                                                            accTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGcc                                                                                          >  2:1773382/1‑150 (MQ=255)
                                                            accTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGcc                                                                                          >  1:759011/1‑150 (MQ=255)
                                                              cTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAg                                                                                        <  1:95517/150‑1 (MQ=255)
                                                                    gCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAg                                                                                 >  1:1430980/1‑150 (MQ=255)
                                                                     cTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGAGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGCCCAgtgctggg                                                                                 <  2:1012771/150‑8 (MQ=255)
                                                                     cTGCCGGAGGAGAGGGTGCCAGTTTTGTAGAACTTTATTCGGCGTACCGCATTGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGCTCGCGGCGGGCGTGGCCAGGCCGAgg                                                                                <  2:958595/150‑1 (MQ=255)
                                                                          ggAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGTCCAGGCCGAGGGCGCa                                                                           <  1:353780/150‑1 (MQ=255)
                                                                             ggTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGACTACCATGCCGTGGCCGCAAACCATCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAggg                                                                        <  1:1673846/150‑1 (MQ=255)
                                                                                      gctgggtttGGCGAACACGACGCGGTTGACCGCATGGGCTCGCTCCAGGGCTTCGGCATCCTGGATCAGGATGCCGGGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCAAGGCGGAGGTCGCAGGGGCAGGCGAt                                                               <  2:169089/144‑1 (MQ=255)
                                                                                      gCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGATTTATTTCATCACTGACCCTTTCAAAGCGGCGACCCCGGTACCGACCATGCTCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGAt                                                               <  1:1875724/150‑1 (MQ=255)
                                                                                        cGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCTCAGGGGCAGGCGATGa                                                             <  1:1004178/150‑1 (MQ=255)
                                                                                           ttttGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTGATCAAGGATGc                                                                                                                                            >  2:150815/1‑69 (MQ=38)
                                                                                           ttttGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTGATCAAGGATGc                                                                                                                                            <  1:150815/69‑1 (MQ=38)
                                                                                                  gAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGCCCCCGGTACCGGCCATGATCGCGGCGGGCGTgg                                                                                            <  1:1077496/110‑1 (MQ=255)
                                                                                                  gAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTgg                                                                                            >  2:1077496/1‑110 (MQ=255)
                                                                                                  gAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCTGGCCGAGGGCGCAGGGGCAGGCGAt                                                               >  1:1535266/1‑138 (MQ=255)
                                                                                                  gAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGAt                                                               <  2:1535266/138‑1 (MQ=255)
                                                                                                  gAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCacggc                                                   <  2:1067521/150‑1 (MQ=255)
                                                                                                  gAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCacggc                                                   <  2:295413/150‑1 (MQ=255)
                                                                                                     cacTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTg                                                                                             <  2:468300/106‑1 (MQ=255)
                                                                                                     cacTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTg                                                                                             >  1:468300/1‑106 (MQ=255)
                                                                                                      acTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCacggcgacg                                               >  2:1830024/1‑150 (MQ=255)
                                                                                                        tACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTg                                                                                             >  2:1901227/1‑103 (MQ=255)
                                                                                                        tACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTg                                                                                             <  1:1901227/103‑1 (MQ=255)
                                                                                                        tACGCGGTTGACCGCATGGACACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCa                                                                                                              <  1:892914/86‑1 (MQ=255)
                                                                                                        tACGCGGTTGACCGCATGGACACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCa                                                                                                              >  2:892914/1‑86 (MQ=255)
                                                                                                           gcgGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGtt                                          <  2:1635422/150‑1 (MQ=255)
                                                                                                             ggTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGa                                        >  1:1478951/1‑150 (MQ=255)
                                                                                                                  aCCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAgc                                   <  2:110917/150‑1 (MQ=255)
                                                                                                                  aCCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAgc                                   >  1:11061/1‑150 (MQ=255)
                                                                                                                       aTGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAg                                                                                 <  2:128559/99‑1 (MQ=255)
                                                                                                                       aTGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAg                                                                                 >  1:128559/1‑99 (MQ=255)
                                                                                                                       aTGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGAAGGCCTTGATCAGCGCGGt                              >  2:725383/1‑150 (MQ=255)
                                                                                                                        tGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGtt                             <  2:670203/150‑1 (MQ=255)
                                                                                                                         gggCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATc                                                                                                         <  2:8338/74‑1 (MQ=255)
                                                                                                                         gggCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATc                                                                                                         >  1:8338/1‑74 (MQ=255)
                                                                                                                         gggCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATAATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGa                                        >  2:743177/1‑138 (MQ=255)
                                                                                                                         gggCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATAATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGa                                        <  1:743177/138‑1 (MQ=255)
                                                                                                                           gCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGGTGATCAGCGCGGTGTcc                          >  2:1144416/1‑150 (MQ=255)
                                                                                                                            cACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAgg                                                                                <  1:477943/95‑1 (MQ=255)
                                                                                                                            cACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAgg                                                                                >  2:477943/1‑95 (MQ=255)
                                                                                                                            cACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGAt                                       <  1:56487/136‑1 (MQ=255)
                                                                                                                            cACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGAt                                       >  2:56487/1‑136 (MQ=255)
                                                                                                                            cACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCa                         >  2:1155849/1‑150 (MQ=255)
                                                                                                                              cGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGACCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGc                       <  1:77832/150‑1 (MQ=255)
                                                                                                                               gCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCACGACCGGGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGTGCGGTTTCCAGCg                      >  2:1598096/1‑150 (MQ=255)
                                                                                                                                 tCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGGGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAACGCGGTTTCCCGCggg                    >  1:365310/1‑149 (MQ=255)
                                                                                                                                 tCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGGGGCGACCCCGGTCCCGGCCATGATCGCGGGGGGCGTGGCCAGGCCGAGGGCGCAGGGGCCGGCGCTGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGt                    >  2:820151/1‑150 (MQ=255)
                                                                                                                                 tCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGt                    >  1:977367/1‑150 (MQ=255)
                                                                                                                                   cAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTAc                  <  1:1578417/150‑1 (MQ=255)
                                                                                                                                   cAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTAc                  >  1:1306480/1‑150 (MQ=255)
                                                                                                                                       gCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTACCcca              <  1:1874275/150‑1 (MQ=255)
                                                                                                                                         cacGGCATCCTGATCAAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTACCccagc            >  1:821399/3‑150 (MQ=255)
                                                                                                                                             gCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTACCCcagccagc        >  1:1350531/1‑150 (MQ=255)
                                                                                                                                             gCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTACCCcagccagc        >  2:1428771/1‑150 (MQ=255)
                                                                                                                                               aTCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTACCCCAgccagcca      >  2:1692127/1‑150 (MQ=255)
                                                                                                                                               aTCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTACCCCAgccagcca      >  2:1691221/1‑150 (MQ=255)
                                                                                                                                                  cTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTCCCGGCCATGATCGCGGCGGGCAAGGCGCGGGCGAGGGGCGACAGCGTGGCCGAGGCGAGCGCGGCGACGGCGGGGGTCGGCGCGGGTTCCAGCGGTACCCCAGccagccacca   >  2:1160305/1‑150 (MQ=255)
                                                                                                                                                  cTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTACCCCAGccagccacca   <  1:1428771/150‑1 (MQ=255)
                                                                                                                                                  cTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTACCCCAGccagccacca   >  1:1116319/1‑150 (MQ=255)
                                                                                                                                                  cTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTACCCCAGccagccacca   >  1:932371/1‑150 (MQ=255)
                                                                                                                                                   ttGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTACCCCAGccagccacccc  >  2:1443991/1‑148 (MQ=255)
                                                                                                                                                   ttGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTACCCCAGccagccaccag  >  2:1760516/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                       
AGGGCGCCGGCCAGGCGGTAAAGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTACCCCAGCCAGCCACCAG  >  NC_002947/684697‑684993

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: