Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 684,846:1 +A 10.9% coding (1401/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,8461.A10.9% 225.7 / 24.9 64coding (1401/2400 nt)cadA‑IIcadmium translocating P‑type ATPase
Reads supporting (aligned to +/- strand):  ref base . (27/30);  new base A (4/3);  total (31/33)
Fisher's exact test for biased strand distribution p-value = 7.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCAGGCGGTAAAGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTACCCCAGCCAGCCACCAGCC  >  NC_002947/684703‑684992
                                                                                                                                                |                                                                                                                                                       
gCCAGGCGGTAAAGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGAt                                                                                                                                                    >  1:242931/1‑149 (MQ=255)
  cAGGCGGTAAAGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGc                                                                                                                                                  >  3:109496/1‑149 (MQ=255)
  cAGGCGGTAAAGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTGATCAAGGATGc                                                                                                                                                  <  5:289600/149‑1 (MQ=255)
  cAGGCGGTAAAGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTGATCAAGGATGc                                                                                                                                                  >  6:289600/1‑149 (MQ=255)
      cGGTAAAGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCG‑‑Tg                                                                                                                                            >  2:87488/1‑149 (MQ=255)
          aaaGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCG‑‑Tggcgg                                                                                                                                        >  2:290629/1‑149 (MQ=255)
               tCGGCGCTGTTGCCGTCCAGCGCCTGGCTATGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCG‑‑TGGCGggcggc                                                                                                                                   >  6:204899/1‑149 (MQ=255)
                 ggCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCG‑‑TGGCGGGCGGCG‑‑a                                                                                                                               >  3:289659/1‑149 (MQ=255)
                     cTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTGGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGTATGCCG‑‑TGGCGGGCGGCG‑‑Acccc                                                                                                                          <  8:260441/149‑1 (MQ=255)
                        ttGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCa                                                                                                                                                        <  1:331172/121‑1 (MQ=255)
                        ttGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCa                                                                                                                                                        >  2:331172/1‑121 (MQ=255)
                         tGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTa                                                                                                                      >  7:26292/1‑149 (MQ=255)
                            cGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCg                                                                                                                   >  5:105949/1‑149 (MQ=255)
                             gTCCAGCGCCTGGCTATGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCgg                                                                                                                  <  5:204899/149‑1 (MQ=255)
                                cAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCGGcc                                                                                                                <  6:105949/148‑1 (MQ=255)
                                         gCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCGGCCATGATCgcgg                                                                                                      <  1:87488/149‑1 (MQ=255)
                                           tgtgAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCG‑‑TGGCGGGCGGCG‑‑Acc                                                                                                                             >  7:240485/1‑125 (MQ=255)
                                           tgtgAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCG‑‑TGGCGGGCGGCG‑‑Acc                                                                                                                             <  8:240485/125‑1 (MQ=255)
                                           tgtgAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCGGCCATGATCgcggcg                                                                                                    >  3:304395/1‑149 (MQ=255)
                                                accaccTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCGGCCATGATCGCGGCGGGCGt                                                                                               <  3:150953/149‑1 (MQ=255)
                                                 ccaccTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCGGCCATGATCgcggcg                                                                                                    <  8:165571/143‑1 (MQ=255)
                                                 ccaccTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCGGCCATGATCgcggcg                                                                                                    >  7:165571/1‑143 (MQ=255)
                                                    ccTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGc                                                                                            <  5:79364/148‑1 (MQ=255)
                                                    ccTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGc                                                                                            >  6:79364/1‑148 (MQ=255)
                                                     cTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCa                                                                                          <  2:242931/149‑1 (MQ=255)
                                                     cTGCGGGCTGCCGGAGGTGAGGGTGCCGGTATTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGc                                                                                            >  1:30047/1‑147 (MQ=255)
                                                         gggCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGcc                                                                                      <  4:109496/149‑1 (MQ=255)
                                                           gCTGCCGGAGGGGAGGGTGCCGGGTTTGTCGACCCCTCCGCGGTTGACCGCATGGGCACGCGCCCGGGCTTCGGCAGCCTTGAGCAGGATGCCG‑‑GGGCGGGCGGCG‑‑ACCCCGGTACCGGCCATGATCGCGGCGGGCGGGGCCAGGCCGa                                                                                    <  5:218693/149‑1 (MQ=255)
                                                             tGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAgg                                                                                  <  2:30047/149‑1 (MQ=255)
                                                             gcccGCAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAgg                                                                                  <  2:257864/147‑1 (MQ=255)
                                                               ccGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGgc                                                                                >  2:121127/1‑149 (MQ=255)
                                                                  gAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCGGCCATGATCGcggcgg                                                                                                   <  3:303528/127‑1 (MQ=255)
                                                                  gAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCGGCCATGATCGcggcgg                                                                                                   >  4:303528/1‑127 (MQ=255)
                                                                      tGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAgggg                                                                         <  1:121127/149‑1 (MQ=255)
                                                                           gTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCG‑‑TGGCGggcggc                                                                                                                                   >  5:339175/1‑89 (MQ=255)
                                                                           gTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCG‑‑TGGCGggcggc                                                                                                                                   <  6:339175/89‑1 (MQ=255)
                                                                             gtcGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGa                                                                  >  7:308866/3‑149 (MQ=255)
                                                                             gCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCG‑‑TGGCGGGCCGCC‑‑ACCCCGGTACCCGCCACGATCCCGGCGGGCGTGGCCAGGCCGAGGGCGGAGGGGCAGGCg                                                                   >  8:225647/1‑148 (MQ=255)
                                                                                               tACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCacggcgacgg                                                >  7:165892/1‑149 (MQ=255)
                                                                                               tACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCacggcgacgg                                                >  7:189852/1‑149 (MQ=255)
                                                                                                 cgcgGTTGCCCGCATGGGCACGCTCCACGGCTTAGGCATCCTTGTTCATGAAGGCGCCTGGCGGGCGGCGGCTCCACCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCacggc                                                     <  1:290629/147‑1 (MQ=255)
                                                                                                 cgcgGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCacggcgacggcg                                              >  1:202686/1‑149 (MQ=255)
                                                                                                   cgGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGtt                                            <  8:308866/149‑1 (MQ=255)
                                                                                                   cgGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGtt                                            <  2:36306/149‑1 (MQ=255)
                                                                                                   cgGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGtt                                            <  4:289659/149‑1 (MQ=255)
                                                                                                    ggTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTg                                           <  8:26292/149‑1 (MQ=255)
                                                                                                    ggTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTg                                           <  5:3742/149‑1 (MQ=255)
                                                                                                            gCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCGGCCATGATCGcggcgg                                                                                                   <  6:23807/85‑1 (MQ=255)
                                                                                                            gCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCGGCCATGATCGcggcgg                                                                                                   >  5:23807/1‑85 (MQ=255)
                                                                                                              aTGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCgg                                 <  2:202686/149‑1 (MQ=255)
                                                                                                                 ggCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGttt                              <  8:165892/149‑1 (MQ=255)
                                                                                                                          cAGGGCTTCGGCATCCTTGATCAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTa                     >  8:323679/1‑149 (MQ=255)
                                                                                                                           aGGGCTTCGGCATCCTTGATCAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCacggcgac                                                  <  4:1847/119‑1 (MQ=255)
                                                                                                                           aGGGCTTCGGCATCCTTGATCAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCacggcgac                                                  >  3:1847/1‑119 (MQ=255)
                                                                                                                              gCTTCGGCATCCTTGATCAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTAccc                  >  7:101284/1‑148 (MQ=255)
                                                                                                                               ccacGGCATCCTGATCAAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTa                     <  5:303908/141‑1 (MQ=255)
                                                                                                                               ccacGGCATCCTGATCAAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTa                     >  6:303908/4‑144 (MQ=255)
                                                                                                                                ttCGGCATCCTTGATCAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTAccc                  >  5:247522/1‑146 (MQ=255)
                                                                                                                                ttCGGCATCCTTGATCAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTAccc                  <  6:247522/146‑1 (MQ=255)
                                                                                                                                 acGGCATCCTGATCAAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGc                                                                                       >  6:98024/2‑76 (MQ=38)
                                                                                                                                 acGGCATCCTGATCAAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGc                                                                                       <  5:98024/75‑1 (MQ=38)
                                                                                                                                   ggCATCCTGATCAAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTACCccagc              >  6:105883/1‑147 (MQ=255)
                                                                                                                                    gCATCCTTGATCAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTACCccagccag           <  8:101284/149‑1 (MQ=255)
                                                                                                                                     cATCCTTGATCAGGCTGGCG‑‑TGGGGGGCGGCG‑‑ACCCCGGTACCGGCCATGATCGCGGCGGGCGGGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTACCCcagccagc          <  4:90671/149‑1 (MQ=255)
                                                                                                                                       tCCTTGATCAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTACCCCAgccagcca        <  7:323679/149‑1 (MQ=255)
                                                                                                                                              tCAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTACCCCAGccagccaccagcc  >  2:108309/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                       
GCCAGGCGGTAAAGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTACCCCAGCCAGCCACCAGCC  >  NC_002947/684703‑684992

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: