Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 684,841 Δ1 bp 10.8% 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,8410T.10.8% 193.2 / 18.1 65coding (1406/2400 nt)cadA‑IIcadmium translocating P‑type ATPase
Reads supporting (aligned to +/- strand):  ref base T (27/31);  new base . (4/3);  total (31/34)
Fisher's exact test for biased strand distribution p-value = 7.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.19e-01

GCCAGGCGGTAAAGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTACCCCAGCCAGCCA  >  NC_002947/684703‑684986
                                                                                                                                          |                                                                                                                                                       
gCCAGGCGGTAAAGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGAt                                                                                                                                              >  1:242931/1‑149 (MQ=255)
  cAGGCGGTAAAGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGc                                                                                                                                            >  3:109496/1‑149 (MQ=255)
  cAGGCGGTAAAGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTGATCAAGGATGc                                                                                                                                            >  6:289600/1‑149 (MQ=255)
  cAGGCGGTAAAGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTGATCAAGGATGc                                                                                                                                            <  5:289600/149‑1 (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)
                          gCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCtt                                                                                                                                                        >  5:305782/1‑114 (MQ=255)
                          gCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCtt                                                                                                                                                        <  6:305782/114‑1 (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                                                                                              >  7:165571/1‑143 (MQ=255)
                                                 ccaccTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCGGCCATGATCgcggcg                                                                                              <  8:165571/143‑1 (MQ=255)
                                                    ccTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGc                                                                                      >  6:79364/1‑148 (MQ=255)
                                                    ccTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGc                                                                                      <  5:79364/148‑1 (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:189852/1‑149 (MQ=255)
                                                                                               tACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCacggcgacgg                                          >  7:165892/1‑149 (MQ=255)
                                                                                                 cgcgGTTGCCCGCATGGGCACGCTCCACGGCTTAGGCATCCTTGTTCATGAAGGCGCCTGGCGGGCGGCGGCTCCACCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCacggc                                               <  1:290629/147‑1 (MQ=255)
                                                                                                 cgcgGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCacggcgacggcg                                        >  1:202686/1‑149 (MQ=255)
                                                                                                   cgGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGtt                                      <  4:289659/149‑1 (MQ=255)
                                                                                                   cgGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGtt                                      <  8:308866/149‑1 (MQ=255)
                                                                                                   cgGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGtt                                      <  2:36306/149‑1 (MQ=255)
                                                                                                    ggTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTg                                     <  5:3742/149‑1 (MQ=255)
                                                                                                    ggTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTg                                     <  8:26292/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               >  6:303908/4‑144 (MQ=255)
                                                                                                                               ccacGGCATCCTGATCAAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTa               <  5:303908/141‑1 (MQ=255)
                                                                                                                                ttCGGCATCCTTGATCAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTAccc            <  6:247522/146‑1 (MQ=255)
                                                                                                                                ttCGGCATCCTTGATCAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTAccc            >  5:247522/1‑146 (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)
                                                                                                                                          |                                                                                                                                                       
GCCAGGCGGTAAAGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCG‑‑TGGCGGGCGGCG‑‑ACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTACCCCAGCCAGCCA  >  NC_002947/684703‑684986

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: