Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,846,846 C→A 8.6% A532A (GCC→GCA)  PP_4261 → cation‑transporting P‑type ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,846,8460CA 8.6% 205.0 / 9.0 70A532A (GCC→GCA) PP_4261cation‑transporting P‑type ATPase
Reads supporting (aligned to +/- strand):  ref base C (32/32);  new base A (3/3);  total (35/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.47e-01

ACCCTGCACAAGCTTGGCCTGCTGGTGACCCGCGGCCACGTGCTGGAAGGCCTGAACCAGATCGACACGGTGATCTTCGACAAGACCGGCACACTGACCGAGGGCCGCCTGACCTTGCGCAGCATCCGCCCGCTCGGCAGCCAGGCCGCCGACCGCTGCCTGGCCCTGGCCGCGGCGCTGGAAAACCGCTCTGAACACCCCATCGCTCG‑‑CGCCTTTGGTCGTACCGCCACACCTGCCGATGACGTACAGAGCGTACCCGGCCTGGGCCTGGAGGGGGTGGTCGACGGCCAGC  >  NC_002947/4846700‑4846991
                                                                                                                                                  |                                                                                                                                                   
aCCCTGCACAAGCTTGGCCTGCTGGTGACCCGCGGCCACGTGCTGGAAGGCCTGAACCAGATCGACACGGTGATCTTCGACAAGACCGGCACACTGACCGAGGGCCGCCTGACCTTGCGCAGCATCCGCCCGCTCGGCAGCCAGgccgc                                                                                                                                                   <  2:142658/149‑1 (MQ=255)
   cTGCACAAGCTTGGCCTGCTGGTGACCCGCGGCCACGTGCTGGAAGGCCTGAACCAGATCGACACGGTGATCTTCGACAAGACCGGCACACTGACCGAGGGCCGCCTGACCTTGCGCAGCATCCGCCCGCTCGGCAGCCAGGCCGCCGa                                                                                                                                                >  5:198502/1‑149 (MQ=255)
      cacaAGCTTGGCCTGCTGGTGACCCGCGGCCACGTGCTGGAAGGCCTGAACCAGATCGACACGGTGATCTTCGACAAGACCGGCACACTGACCGAGGGCCGCCTGACCTTGCGCAGCATCCGCCCGCTCGGCAGCCAGGCCGCCGACCg                                                                                                                                             >  5:248229/1‑149 (MQ=255)
       acaAGCTTGGCCTGCTGGTGACCCGCGGCCACGTGCTGGAAGGCCTGAACCAGATCGACACGGTGATCTTCGACAAGACCGGCACACTGACCGAGGGCCGCCTGACCTTGCGCAGCATCCGCCCGCTCGGCAGCCAGGCCGCCGACCGc                                                                                                                                            >  1:47991/1‑149 (MQ=255)
       acaAGCTTGGCCTGCTGGTGACCCGCGGCCACGTGCTGGAAGGCCTGAACCAGATCGACACGGTGATCTTCGACAAGACCGGCACACTGACCGAGGGCCGCCTGACCTTGCGCAGCATCCGCCCGCTCGGCAGCCAGGCCGCCGACCGc                                                                                                                                            <  6:198502/149‑1 (MQ=255)
         aaGCTTGGCCTGCTGGTGACCCGCGGCCACGTGCTGGAAGGCCTGAACCAGATCGACACGGTGATCTTCGACAAGACCGGCACACTGACCGAGGGCCGCCTGACCTTGCGCAGCATCCGCCCGCTCGGCAGCCAGGCCGCCGACCGCTg                                                                                                                                          >  5:359330/1‑149 (MQ=255)
          aGCTTGGCCTGCTGGTGACCCGCGGCCACGTGCTGGAAGGCCTGAACCAGATCGACACGGTGATCTTCGACAAGACCGGCACACTGACCGAGGGCCGCCTGACCTTGCGCAGCATCCGCCCGCTCGGCAGCCAGGCCGCCGACCGCTGc                                                                                                                                         <  4:136042/149‑1 (MQ=255)
           gCTTGGCCTGCTGGTGACCCGCGGCCACGTGCTGGAAGGCCTGAACCAGATCGACACGGTGATCTTCGACAAGACCGGCACACTGACCGAGGGCCGCCTGACCTTGCGCAGCATCCGCCCGCTCGGCAGCCAGGCCGCCGACCGCTGcc                                                                                                                                        >  7:192281/1‑149 (MQ=255)
               ggCCTGCTGGTGACCCGCGGCCACGTGCTGGAAGGCCTGAACCAGATCGACACGGTGATCTTCGACAAGACCGGCACACTGACCGAGGGCCGCCTGACCTTGCGCAGCATCCGCCCGCTCGGCAGCCAGGCCGCCGACCGCTGcctgg                                                                                                                                     <  3:54689/148‑1 (MQ=255)
                gCCTGCTGGTGACCCGCGGCCACGTGCTGGAAGGCCTGAACCAGATCGACACGGTGATCTTCGACAAGACCGGCACACTGACCGAGGGCCGCCTGACCTTGCGCAGCATCCGCCCGCTCGGCAGCCAGGCAGCGGTCGGCGGcctggc                                                                                                                                    <  6:72611/148‑1 (MQ=255)
                 ccTGCTGGTGACCCGCGGCCACGTGCTGGAAGGCCTGAACCAGATCGACACGGTGATCTTCGACAAGACCGGCACACTGACCGAGGGCCGCCTGACCTTGCGCAGCATCCGCCCGCTCGGCAGCCAGGCCGCCGACCGCTGcctggccc                                                                                                                                  >  1:154202/1‑149 (MQ=255)
                 ccTGCTGGTGACCCGCGGCCACGTGCTGGAAGGCCTGAACCAGATCGACACGGTGATCTTCGACAAGACCGGCACACTGACCGAGGGCCGCCTGACCTTGCGCAGCATCCGCCCGCTCGGCAGCCAGGCCGCCGACCGCTGcctggccc                                                                                                                                  >  7:294787/1‑149 (MQ=255)
                        gTGACCCGCGGCCACGTGCTAGAAGGCCTGAACCAGATCGACACGGTGATCTTCGACAAGACCGGCACACTGACCGAGGGCCGCCTGACCTTGCGCAGCATCCGCCCGCTCGGCAGCCAGGCCGCCGACCGCTGCCTGGCCCTGGCcgc                                                                                                                           >  7:53128/1‑149 (MQ=255)
                           aCCCGCGGCCACGTGCTGGAAGGCCTGAACCAGATCGACACGGTGATCTTCGACAAGACCGGCACACTGACCGAGGGCCGCCTGACCTTGCGCAGCATCCGCCCGCTCGGCAGCCAGGCCGCCGACCGCTGCCTGGCCCTGGCCgcggc                                                                                                                        >  2:344354/1‑149 (MQ=255)
                           aCCCGCGGCCACGTGCTGGAAGGCCTGAACCAGATCGACACGGTGATCTTCGACAAGACCGGCACACTGACCGAGGGCCGCCTGACCTTGCGCAGCATCCGCCCGCTCGGCAGCCAGGCCGCCGACCGCTGCCTGGCCCTGGCCgcgg                                                                                                                         <  6:248229/148‑1 (MQ=255)
                             ccGCGGCCACGTGCTGGAAGGCCTGAACCAGATCGACACGGTGATCTTCGACAAGACCGTCACACTGACCGAGGGCCGCCTGACCTTGCGCAGCATCCGCCCGCTCGGCAGCCAGGCCGCCGACCGCTGCCTGGCCCTGTCCGCGgcgc                                                                                                                      <  2:47991/149‑1 (MQ=255)
                              cgcgGCCACGTGCTGGAAGGCCTGAACCAGATCGACACGGTGATCTTCGACAAGACCGGCACACTGACCGAGGGCCGCCTGACCTTGCGCAGCATCCGCCCGCTCGGCAGCCAGGCCGCCGACCGCTGCCTGGCCCTGGCCGCGGCGCt                                                                                                                     <  2:306784/149‑1 (MQ=255)
                                 ggCCACGTGCTGGAAGGCCTGAACCAGATCGACACGGTGATCTTCGACAAGACCGGCACACTGACCGAGGGCCGCCTGACCTTGCGCAGCATCCGCCCGCTCGGCAGCCAGGCAGCGGTCGGCGGcctggct                                                                                                                                   <  6:173100/132‑2 (MQ=255)
                                 ggCCACGTGCTGGAAGGCCTGAACCAGATCGACACGGTGATCTTCGACAAGACCGGCACACTGACCGAGGGCCGCCTGACCTTGCGCAGCATCCGCCCGCTCGGCAGCCAGGCAGCGGTCGGCGGcctggct                                                                                                                                   >  5:173100/1‑131 (MQ=255)
                                     aCGTGCTGGAAGGCCTGAACCAGATCGACACGGTGATCTTCGACAAGACCGGCACACTGACCGAGGGCCGCCTGACCTTGCGCAGCATCCGCCCGCTCGGCAGCCAGGCCGCCGACCGCTGCCTGGCCCTGGCcgc                                                                                                                           <  2:10217/136‑1 (MQ=255)
                                     aCGTGCTGGAAGGCCTGAACCAGAACGACACGGTGATCTTCGACAAGACCGGCACACTGACCGAGGGCCGCCTGACCTTGCGCAGCATCCGCCCGCTCGGCAGCCAGGCCGCCGACCGCTGCCTGGCCCTGGCcgc                                                                                                                           >  1:10217/1‑136 (MQ=255)
                                             gAAGGCCTGAACCAGATCGACACGGTGATCTTCGACAAGACCGGCACACTGACCGAGGGCCGCCTGACCTTGCGCAGCATCCGCCCGCTCGGCAGCCAGGCCGCCGACCGCTGCCTGGCCCTGGCCGCGGCGCTGGAAAACCGCTCTg                                                                                                       >  5:335666/1‑148 (MQ=255)
                                              aaGGCCTGAACCAGATCGACACGGTGATCTTCGACAAGACCGGCACACTGACCGAGGGCCGCCTGACCTTGCGCAGCATCCGCCCGCTCGGCAGCCAGGCCGCCGACCGCTGCCTGGCCCTGGCCGCGGCGCTGGAAAACCGCTCTGaa                                                                                                     <  1:344354/149‑1 (MQ=255)
                                                   cTGAACCAGATCGACACGGTGATCTTCGACAAGACCGGCACACTGACCGAGGGCCGCCTGACCTTGCGCAGCATCCGCCCGCTCGGCAGCCAGGCCGCCGACCGCTGCCTGGCCCTGGCCGCGGCGCTGGAAAACCGCTCTGAACAccc                                                                                                >  4:75852/1‑149 (MQ=255)
                                                       aCCAGATCGACACGGTGATCTTCGACAAGACCGGCACACTGACCGAGGGCCGCCTGACCTTGCGCAGCATCCGCCCGCTCGGCAGCCAGGCCGCCGACCGCTGcct                                                                                                                                       >  1:375444/1‑106 (MQ=255)
                                                       aCCAGATCGACACGGTGATCTTCGACAAGACCGGCACACTGACCGAGGGCCGCCTGACCTTGCGCAGCATCCGCCCGCTCGGCAGCCAGGCCGCCGACCGCTGcct                                                                                                                                       <  2:375444/106‑1 (MQ=255)
                                                         cAGATCGACACGGTGATCTTCGACAAGACCGGCACACTGACCGAGGGCCGCCTGACCTTGCGCAGCATCCGCCCGCTCGGCAGCCAGGCCGCCGACCGCTGCCTGGCCCTGGCCGCGGCGCTGGAAAACCGCTCTGAACACCCCAtcgc                                                                                          >  8:42557/1‑149 (MQ=255)
                                                          aGATCGACACGGTGATCTTCGACAAGACCGGCACACTGACCGAGGGCCGCCTGACCTTGCGCAGCATCCGCCCGCTCGGCAGCCAGGCCGCCGACCGCTGCCTGGCCCTGGCCGCGGCGCTGGAAAACCGCTCTGAACACCCCAtcgct                                                                                         >  1:160690/1‑149 (MQ=255)
                                                           gATCGACACGGTGATCTTCGACAAGACCGGCACACTGACCGAGGGCCGCCTGACCTTGCGCAGCATCCGCCCGCTCGGCAGCCAGGCCGCCGACCGCTGCCTGGCCCTGGCCGCGGCGCTGGAAAACCGCTCTGAACACCCCAtcgctc                                                                                        <  6:246485/149‑1 (MQ=255)
                                                           gATCGACACGGTGATCTTCGACAAGACCGGCACACTGACCGAGGGCCGCCTGACCTTGCGCAGCATCCGCCCGCTCGGCAGCCAGGCCGCCGACCGCTGCCTGGCCCTGGCCGCGGCGCTGGAAAACCGCTCTGAACACCCCAtcgctc                                                                                        <  1:186973/149‑1 (MQ=255)
                                                               gACACGGTGATCTTCGACAAGACCGGCACACTGACCGAGGGCCGCCTGACCTTGCGCAGCATCCGCCCGCTCGGCAGCCAGGGCGCCGACCGCTGCCTGGCCCTGGCCGCGGCGCTTGAAAACCTCTCTGAACACCCCATCGCTCg‑‑cgc                                                                                  <  8:53128/149‑1 (MQ=255)
                                                               gACACGGTGATCTTCGACAAGACCGGCACACTGACCGAGGGCCGCCTGACCTTGCGCAGCATCCGCCCGCTCGGCAGCCAGGCCGCCGACCGCTGCCTGGCCCTGGCCGCGGCGCTGGAAAACCGCTCTGAACACCCCATCGCTCg‑‑cgc                                                                                  <  8:133991/149‑1 (MQ=255)
                                                                 cacGGTGATCTTCGACAAGACCGGCACACTGACCGCGGGCCGCCTGACCTTGCGCAGCATCCGCCCGCTCGGCAGCCAGGCCGCCGACCGCTGCCTGGCCCTGGCCGCGGCGCTGGAAAACCGCGCTGAACACCCCATCGCACG‑‑CGCCt                                                                                >  7:141158/1‑149 (MQ=255)
                                                                   cGGTGATCTTCGACAAGACCGGCACACTGACCGAGGGCCGCCTGACCTTGCGCAGCATCCGCCCGCTCGGCAGCCAGGCCGCCGACCGCTGCCTGGCCCTGGCCGCGGCGCTGGAAAACCGCTCTGAACACCCCATCGCTCG‑‑CGCCttt                                                                              <  4:304679/149‑1 (MQ=255)
                                                                    ggTGATCTTCGACAAGACCGGCACACTGACCGAGGGCCGCCTGACCTTGCGCAGCATCCGCCCGCTCGGCAGCCAGGCCGCCGACCGCTGCCTGGCCCTGGCCGCGGCGCTGGAAAACCGCTCTGAACACCCCATCGCTCG‑‑CGCCTTTg                                                                             <  2:271875/149‑1 (MQ=255)
                                                                     gTGATCTTCGACAAGACCGGCACACTGACCGAGGGCCGCCTGACCTTGCGCAGCATCCGCCCGCTCGGCAGCCAGGCCGCCGACCGCTGCCTGGCCCTGGCCGCGGCGCTGGAAAACCGCTCTGAACACCCCATCGCTCG‑‑CGCCTTTgg                                                                            >  2:310079/1‑149 (MQ=255)
                                                                      tGATCTTCGACAAGACCGGCACACTGACCGAGGGCCGCCTGACCTTGCGCAGCATCCGCCCGCTCGGCAGCCAGGCCGCCGACCGCTGCCTGGCCCTGGCCGCGGCGCTGGAAAACCGCTCTGAACACCCCATCGCTCG‑‑CGCCTTTGGt                                                                           <  5:344333/149‑1 (MQ=255)
                                                                             cGACAAGACCGGCACACTGACCGAGGGCCGCCTGACCTTGCGCAGCATCCGCCCGCTCGGCAGCCAGGCCGCCGACCGCTGCCTGGCCCTGGCCGCGGCGCTGGAAAACCGCTCTGAACACCCCATCGCTCG‑‑CGCCTTTGGTCGt                                                                        <  7:42557/145‑1 (MQ=255)
                                                                             cGACAAGACCGGCACACTGACCGAGGGCCGCCTGACCTTGCGCAGCATCCGCCCGCTCGGCAGCCAGGCCGCCGACCGCTGCCTGGCCCTGGCCGCGGCGCTGGAAAACCGCTCTGAACACCCCATCGCTCG‑‑CGCCTTTGGTCGTAc                                                                      <  4:198984/147‑1 (MQ=255)
                                                                              gACAAGACCGGCACACTGACCGAGGGCCGCCTGACCTTGCGCAGCATCCGCCCGCTCGGCAGCCAGGCCGCCGACCGCTGCCTGGCCCTGGCCGCGGCGCTGGAAAACCGCTCTGAACACCCCATCGCTCG‑‑CGCCTTTGGTCGTACCGc                                                                   >  5:168857/1‑149 (MQ=255)
                                                                                  aGACCGGCACACTGACCGAGGGCCGCCTGACCTTGCGCAGCATCCGCCCGCTCGGCAGCCAGGCCGCCGACCGCTGCCTGGCCCTGGCCGCGGCGCTGGAAAACCGCTCTGAACACCCCATCGCTCG‑‑CGCCTTTGGTCGTACCGCcaca                                                               >  5:21519/1‑149 (MQ=255)
                                                                                    aCCGGCACACTGACCGAGGGCCGCCTGACCTTGCGCAGCATCCGCCCGCTCGGCAGCCAGGCCGCCGACCGCTGCCTGGCCCTGGCCGCGGCGCTGGAAAACCGct                                                                                                          <  8:184512/106‑1 (MQ=255)
                                                                                      cGGCACACTGACCGAGGGCCGCCTGACCTTGCGCAGCATCCGCCCGCTCGGCAGCCAGGCCGCCGACCGCTGCCTGGCCCTGGCCGCGGCGCTGGAAAACCGCTCTGAACACCCCATCGCTCG‑‑CGCCTTTGGTCGTACCGCCACACCTg                                                           <  3:75852/149‑1 (MQ=255)
                                                                                       ggCACACTGACCGAGGGCCGCCTGACCTTGCGCAGCATCCGCCCGCTCGGCAGCCAGGCCGCCGACCGCTGCCTGGCCCTGGCCGCGGCGCTGGAAAACCGCTCTGAACACCCCATCGCTCG‑‑CGCCTTTGGTCGTACCGc                                                                   >  5:336935/1‑140 (MQ=255)
                                                                                       ggCACACTGACCGAGGGCCGCCTGACCTTGCGCAGCATCCGCCCGCTCGGCAGCCAGGCCGCCGACCGCTGCCTGGCCCTGGCCGCGGCGCTGGAAAACCGCTCTGAACACCCCATCGCTCG‑‑CGCCTTTGGTCGTACCGc                                                                   <  6:336935/140‑1 (MQ=255)
                                                                                           cacTGACCGAGGGCCGCCTGACCTTGCGCAGCATCCGCCCGCTCGGCAGCCAGGCCGCCGACCGCTGCCTGGCCCTGGCCGCGGCGCTGGAAAACCGCTCTGAACACCCCATCGCTCG‑‑CGCCTTTGGTCGTACCGCCACACCTGCCGAt                                                      <  8:192281/149‑1 (MQ=255)
                                                                                            acTGACCGAGGGCCGCCTGACCTTGCGCAGCATCCGCCCGCTCGGCAGCCAGGCCGCCGACCGCTGCCTGGCCCTGGCCGCGGCGCTGGAAAACCGCTCTGAACACCCCATCGCTCG‑‑CGCCTTTGGTCGTACCGCCACACCTGCCGAt                                                      >  4:179842/1‑148 (MQ=255)
                                                                                              tGACCGAGGGCCGCCTGACCTTGCGCAGCATCCGCCCGCTCGGCAGCCAGGCCGCCGACCGCTGCCTGGCCCTGGCCGCGGCGCTGGAAAACCGCTCTGAACACCCCATCGCTCG‑‑CGCCTTTGGTCGTACCGCCACACCTGCCGATGAc                                                   >  8:30167/1‑149 (MQ=255)
                                                                                                aCCGAGGGCCGCCTGACCTTGCGCAGCATCCGCCCGCTCGGCAGCCAGGCCGCCGACCGCTGCCTGGCCCTGGCCGCGGCGCTGGAAAACCGCTCTGAACACCCCATCGCTCG‑‑CGCCTTTGGTCGTACCGCCACACCTGCCGATGACg                                                  <  3:234191/148‑1 (MQ=255)
                                                                                                  cGAGGGCCGCCTGACCTTGCGCAGCATCCGCCCGCTCGGCAGCCAGGCCGCCGACCGCTGCCTGGCCCTGGCcgc                                                                                                                           >  1:182494/1‑75 (MQ=255)
                                                                                                  cGAGGGCCGCCTGACCTTGCGCAGCATCCGCCCGCTCGGCAGCCAGGCCGCCGACCGCTGCCTGGCCCTGGCcgc                                                                                                                           <  2:182494/75‑1 (MQ=255)
                                                                                                      ggCCGCCTGACCTTGCACAGCATCCGCCCGCTCGGCAGCCAGGCCGCCGACCGCTGCCTGGCCCTGGCCGCGGCGCTGGAAAACCGCTCTGAACACCCCATCGCTCG‑‑CGCCTTTGGTCGTACCGCCACACCTGCCGATGACGTACagag                                           >  2:109685/1‑149 (MQ=255)
                                                                                                       gccgccTGACCTTGCGCAGCATCCGCCCGCTCGGCAGCCAGGCCGCCGACCGCTGCCTGGCCCTGGCCGCGGCGCTGGAAAACCGCTCTGAACACCCCATCGCTCG‑‑CGCCTTTGGTCGTACCGCCACACCTGCCGATGACGTACaga                                            <  2:213983/147‑1 (MQ=255)
                                                                                                       gccgccTGACCTTGCGCAGCATCCGCCCGCTCGGCAGCCAGGCCGCCGACCGCTGCCTGGCCCTGGCCGCGGCGCTGGAAAACCGCTCTGAACACCCCATCGCTCG‑‑CGCCTTTGGTCGTACCGCCACACCTGCCGATGACGTACaga                                            >  1:213983/1‑147 (MQ=255)
                                                                                                               aCCTTGCGCAGCATCCGCCCGCTCGGCAGCCAGGCCGCCGACCGCTGCCTGGCCCTGGCCGCGGCGCTGGAAAACCGCTCTGAACACCCCATCGCTCG‑‑CGCCTTTGGTCGTACCGCCACACCTGCCGATGACGTACAGAGCGTACCCgg                                  >  5:129562/1‑149 (MQ=255)
                                                                                                                  ttGCGCAGCATCCGCCCGCTCGGCAGCCAGGCCGCCGACCGCTGCCTGGCCCTGGCCGCGGCGCTGGAAAACCGCTCTGAACACCCCATCGCTCG‑‑CGCCTTTGGTCGTACCGCCACACCTGCCGATGACGTACAGAGCGTACCCggcct                               >  2:134055/1‑149 (MQ=255)
                                                                                                                  ttGCGCAGCATCCGCCCGCTCGGCAGCCAGGCCGCCGACCGCTGCCTGGCCCTGGCCGCGGCGCTGGAAAACCGCTCTGAACACCCCATCGCTCG‑‑CGCCTTTGGACGTACCGCCACACCTGCCGATGACGTACAGAGCGTACCCggcct                               >  6:219613/1‑149 (MQ=255)
                                                                                                                     cgcAGCATCCGCCCGCTCGGCAGCCAGGCCGCCGACCGCTGCCTGGCCCTGGCCGCGGCGCTGGAAAACCGCTCTGAACACCCCATCGCTCG‑‑CGCCTTTGGTCGTACCGCCACACCTGCCGATGACGTACAGAGCGTACCCggcctgg                             <  8:277747/148‑1 (MQ=255)
                                                                                                                       cagcaTCCGCCCGCTCGGCAGCCAGGCCGCCGACCGCTGCCTGGCCCTGGCCGCGGCGCTGGAAAACCGCTCTGAACACCCCATCGCTCG‑‑CGCCTTTGGTCGTACCGCCACACCTGCCGATGACGTACAGAGCGTACCCggcctgggcc                          <  2:160690/149‑1 (MQ=255)
                                                                                                                          caTCCGCCCGCTCGGCAGCCAGGCCGCCGACCGCTGCCTGGCCCTGGCCGCGGCGCTGGAAAACCGCTCTGAACACCCCATCGCTCG‑‑CGCCTTTGGTCGTACCGCCACACCTGCCGATGACGTACAGAGCGTACCCGgcctgggcctgg                       >  7:6104/1‑149 (MQ=255)
                                                                                                                            tCCGCCCGCTCGGCAGCCAGGCCGCCGACCGCTGCCTGGCCCTGGCCGCGGCGCTGGAAAACCGCTCTGAACACCCCATCGCTCG‑‑CGCCTTTGGTCGTACCGCCACACCTGCCGATGACGTACAGAGCGTACCCGGCCTGGGCCTGGAg                     <  6:264786/149‑1 (MQ=255)
                                                                                                                             ccgcccgcTCGGCAGCCAGGCCGCCGACCGCTGCCTGGCCCTGGCCGCGGCGCTGGAAAACCGCTCTGAACACCCCATCGCTCG‑‑CGCCTTTGGTCGTACCGCCACACCTGCCGATGACGTACAGAGCGTACCCGGCCTGGGCCTGGAgg                    >  3:73187/1‑149 (MQ=255)
                                                                                                                              ctcccgcGCGGCAGCCAGGCCGCCGACCGCGGCCTGGCCCTGGCCGCGGCGCTGGAAAACCGCTCTGAACACCCCATCGCTCG‑‑CGCCTTTGGTCGTACCGCCACCCCTGCCGATGACGTACAGAGCGTACCCGGCCTGGGCCTGGAggg                   <  8:294787/147‑1 (MQ=255)
                                                                                                                                  cgcTCCGGACCCAGGCCCCCGACCGCTGCCCGGCCCTGGCCCCGGCGCTGGCAAACCCCTCTGAACACCCCCTCGCTCGCCCTCCTTTGGTCGTACCGCCACACCTGCCGATGACGTACAGAGCGTACCCGGCCTGGGCCTGGAggggg                 <  1:310079/149‑1 (MQ=255)
                                                                                                                                      ccagggCCAGGCAGCGGTCGGCGGCCTGGCCCTGGCCGCGGCGCTGGAAAACCGCTCTGAACACCCCATCGCTCG‑‑CGCCTTTGGTCGTACCGCCACACCTGCCGATGACGTACAGAGCGTACCCGGCCTGGGCCTGGAGGGGGTGGTc            >  3:179458/6‑148 (MQ=255)
                                                                                                                                        gCAGCCAGGCCGCCGACCGCTGCCTGGCCCTGGCCGCGGCGCTGGAAAACCGCTCTGAACACCCCATCGCTCG‑‑CGCCTTTGGTCGTACCGCCACACCTGCCGATGACGTACAGAGCGTACCCGGCCTGGGCCTGGAGGGGGTGGTCGAc         >  8:183479/1‑149 (MQ=255)
                                                                                                                                        agggCCAGGCAGCGGTCGGCGGCCTGGCCCTGGCCGCGGCGCTGGAAAACCGCTCTGAACACCCCATCGCTCG‑‑CGCCTTTGGTCGTACCGCCACACCTGCCGATGACGTACAGAGCGTACCCGGCCTGGGCCTGGAGGGGGTGGTCGAc         <  4:179458/146‑1 (MQ=255)
                                                                                                                                         cAGCCAGGCCGCCGACCGCTGCCTGGCCCTGGCCGCGGCGCTGGAAAACCGCTCTGAACACCCCATCGCTCG‑‑CGCCTTTGGTCGTACCGCCACACCTGCCGATGACGTACAGAGCGTACCCGGCCTGGGCCTGGAGGGGGTGGTCGACg        <  6:359330/149‑1 (MQ=255)
                                                                                                                                           gCCAGGCAGCGGTCGGCGGCCTGGCCCTGGCCGCGGCGCTGGAAAACCGCTCTGAACACCCCATCGCTCG‑‑CGCCATAGGTCGTACCGCCACACCTGCCGATGACGAACAGCGCGTACCCGGCCTGGGCCTGGAGGGCGTGGTCGACGGc      >  8:287365/1‑149 (MQ=255)
                                                                                                                                              aGGCCGCCGACCGCTGCCTGGCCCTGGCCGCGGCGCTGGAAAACCGCTCTGACCACCCCCTCGCTCG‑‑CGCCTTTGGTCGTACCGCCACACCTGCCGATGACGTACAGCGCGTACCCGGCCTGGGCCTGGAGGGGGTGGTCGACGGCCAg   >  5:320464/1‑149 (MQ=255)
                                                                                                                                               ggCCGCCGACCGCTGCCTGGCCCTGGCCGCGGCGCTGGAAAACCGCTCTGAACACCCCATCGCTCG‑‑CGCCTTTGGTCGTACCGCCACACCTGCCGATGACGTACAGAGCGTACCCGGCCTGGGCCTGGAGGGGGTGGTCGACGGCCAgc  <  6:168857/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
ACCCTGCACAAGCTTGGCCTGCTGGTGACCCGCGGCCACGTGCTGGAAGGCCTGAACCAGATCGACACGGTGATCTTCGACAAGACCGGCACACTGACCGAGGGCCGCCTGACCTTGCGCAGCATCCGCCCGCTCGGCAGCCAGGCCGCCGACCGCTGCCTGGCCCTGGCCGCGGCGCTGGAAAACCGCTCTGAACACCCCATCGCTCG‑‑CGCCTTTGGTCGTACCGCCACACCTGCCGATGACGTACAGAGCGTACCCGGCCTGGGCCTGGAGGGGGTGGTCGACGGCCAGC  >  NC_002947/4846700‑4846991

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