Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,177,677 T→A 6.5% D103V (GAC→GTC)  PP_4553 ← ECF family RNA polymerase sigma‑70 factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,177,6770TA 6.5% 217.0 / 3.2 71D103V (GAC→GTC) PP_4553ECF family RNA polymerase sigma‑70 factor
Reads supporting (aligned to +/- strand):  ref base T (36/30);  new base A (3/2);  total (39/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

ATACGCTGGGCGATGGTGGCCGGGCTTTGCAGAAAGGCGCGGGCGATCTGCTCGGTGGTCAGGTCGCAGACTTCGCGAAGCGTCAACGGCACCTGAGCGTCGGCCGACAGGGCCGGATGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGGCGCAGCGCGTCAATGGCCTTGAAGCGCCCGACAGAAACCAGCCACGCTCGC  >  NC_002947/5177530‑5177823
                                                                                                                                                   |                                                                                                                                                  
aTACGCTGGGCGATGGTGGCCGGGCTTTGCAGAAAGGCGCGGGCGATCTGCTCGGTGGTCAGGTCGCAGACTTCGCGAAGCGTCAACGGCACCTGAGCGTCGGCCGACAGGGCCGGATGGCAACAGGTGAAGATCAGCCGCAGGCGGt                                                                                                                                                    >  7:212644/1‑148 (MQ=255)
    gCTGGGCGATGGTGGCCGGGCTTTGCAGAAAGGCGCGGGCGATCTGCTCGGTGGTCAGGTCGCAGACTTCGCGAAGCGTCAACGGCACCTGAGCGTCGGCCGACAGGGCCGGATGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATcc                                                                                                                                               >  2:326176/1‑149 (MQ=255)
          cGATGGTGGCCGGGCTTTGCAGAAAGGCGCGGGCGATCTGCTCGGTGGTCAGGTCGCAGACTTCGCGAAGCGTCAACGGCACCTGAGCGTCGGCCGACATGGCCGGATGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCcagc                                                                                                                                         <  1:270782/149‑1 (MQ=255)
           gATGGTGGCCGGGCTTTGCAGAAAGGCGCGGGCGATCTGCTCGGTGGTCAGGTCGCAGACTTCGCGAAGCGTCAACGGCACCTGAGCGTCGGCCGACAGGGCCGGATGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCcagca                                                                                                                                        <  7:279117/149‑1 (MQ=255)
             tggtggCCGGGCTTTGCAGAAAGGCGCGGGCGATCTGCTCGGTGGTCAGGTCGCAGACTTCGCGAAGCGTCAACGGCACCTGAGCGTCGGCCGACAGGGCCGGATGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCagcagc                                                                                                                                      >  3:306487/1‑149 (MQ=255)
                tggCCGGGCTTTGCAGAAAGGCGCGGGCGATCTGCTCGGTGGTCAGGTCGCAGACTTCGCGAAGCGTCAACGGCACCTGAGCGTCGGCCGACAGGGCCGGATGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGctct                                                                                                                                   <  3:94326/149‑1 (MQ=255)
                   ccGGGCTTTGCAGAAAGGCGCGGGCGATCTGCTCGGTGGTCAGGTCGCAGACTTCGCGAAGCGTCAACGGCACCTGAGCGTCGGCCGACAGGGCCGGATGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCtcttct                                                                                                                                >  8:183826/1‑149 (MQ=255)
                    cGGGCTTTGCAGAAAGGCGCGGGCGATCTGCTCGGTGGTCAGGTCGCAGACTTCGCGAAGCGTCAACGGCACCTGAGCGTCGGCCGACAGGGCCGGATGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTcttctt                                                                                                                               >  3:192719/1‑149 (MQ=255)
                      ggCTTTGCAGAAAGGCGCGGGCGATCTGCTCGGTGGTCAGGTCGCAGACTTCGCGAAGCGTCAACGGCACCTGAGCGTCGGCCGACAGGGCCGGATGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCttcttc                                                                                                                              >  3:261701/1‑148 (MQ=255)
                      ggCTTTGCAGAAAGGCGCGGGCGATCTGCTCGGTGGTCAGGTCGCAGACTTCGCGAAGCGTCAACGGCACCTGAGCGTCGGCCGACAGGGCCGGATGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCa                                                                                                                             <  4:306487/149‑1 (MQ=255)
                         tttGCAGAAAGGCGCGGGCGATCTGCTCGGTGGTCAGGTCGCAGACTTCGCGAAGCGTCAACGGCACCTGAGCGTCGGCCGACAGGGCCGGATGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTg                                                                                                                          >  5:217572/1‑149 (MQ=255)
                           tGCAGAAAGGCGCGGGCGATCTGCTCGGTGGTCAGGTCGCAGACTTCGCGAAGCGTCAACGGCACCTGAGCGTCGGCCGACAGGGCCGGATGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTg                                                                                                                          <  5:276846/147‑1 (MQ=255)
                               gAAAGGCGCGGGCGATCTGCTCGGTGGTCAGGTCGCAGACTTCGCGAAGCGTCAACGGCACCTGAGCGTCGGCCGACAGGGCCGGATGGCAACAGGTGAAGAAGAGCTGCTGGCGGATGACCGCCTGCGGCTGATCTTCAcctgttg                                                                                                                      <  5:153730/147‑7 (MQ=21)
                               gAAAGGCGCGGGCGATCTGCTCGGTGGTCAGGTCGCAGACTTCGCGAAGCGTCAACGGCACCTGAGCGTCGGCCGACAGGGCCGGATGGCAACAGGTGAAGAAGAGCTGCTGGCGGATGACCGCCTGCGGCTGATCTTCAcctgttg                                                                                                                      >  6:153730/1‑141 (MQ=21)
                                    gcgcGGGCGATCTGCTCGGTGGTCAGGTCGCAGACTTCGCGAAGCGTCAACGGCACCTGAGCGTCGGCCGACAGGGCCGGATGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTgcgc                                                                                                               <  8:201358/149‑1 (MQ=255)
                                            gATCTGCTCGGTGGTCAGGTCGCAGACTTCGCGAAGCGTCAACGGCACCTGAGCGTCGGCCGACAGGGCCGGATGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCt                                                                                                       >  3:134037/1‑149 (MQ=255)
                                                 gCTCGGTGGTCAGGTCGCAGACTTCGCGAAGCGTCAACGGCACCTGAGCGTCGGCCGACAGGGCCGGATGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCcagc                                                                                                  >  7:156540/1‑149 (MQ=255)
                                                 gCTCGGTGGTCAGGTCGCAGACTTCGCGAAGCGTCAACGGCACCTGAGCGTCGGCCGACAGGGCCGGATGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCcag                                                                                                   >  1:165985/1‑148 (MQ=255)
                                                  cTCGGTGGTCAGGTCGCAGACTTCGCGAAGCGTCAACGGCACCTGAGCGTCGGCCGACAGGGCCGGATGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCcagca                                                                                                 <  5:85770/149‑1 (MQ=255)
                                                  cTCGGTGGTCAGGTCGCAGACTTCGCGAAGCGTCAACGGCACCTGAGCGTCGGCCGACAGGGCCGGATGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCca                                                                                                    <  1:326176/146‑1 (MQ=255)
                                                       tggtCAGGTCGCAGACTTCGCGAAGCGTCAACGGCACCTGAGCGTCGGCCGACAGGGCCGGATGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATg                                                                                            <  8:156540/149‑1 (MQ=255)
                                                       tggtCAGGTCGCAGACTTCGCGAAGCGTCAACGGCACCTGAGCGTCGGCCGACAGGGCCGGATGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATg                                                                                            <  4:261701/149‑1 (MQ=255)
                                                              gTCGCAGACTTCGCGAAGCGTCAACGGCACCTGAGCGTCGGCCGACAGGGCCGGATGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTg                                                                                                                          <  3:271036/112‑1 (MQ=255)
                                                              gTCGCAGACTTCGCGAAGCGTCAACGGCACCTGAGCGTCGGCCGACAGGGCCGGATGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTg                                                                                                                          >  4:271036/1‑112 (MQ=255)
                                                              gTCGCAGACTTCGCGAAGCGTCAACGGCACATGAGCGTCGCCCGACAGGGCCGGATGGCAAAAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCATGCGCCTGACCCTCCAGCAGCATGATCAAg                                                                                      >  3:30278/1‑148 (MQ=255)
                                                                 gCAGACTTCGCGAAGCGTCAACGGCACCTGAGCGTCGGCCGACAGGGCCGGATGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGc                                                                                           >  7:324838/1‑140 (MQ=255)
                                                                 gCAGACTTCGCGAAGCGTCAACGGCACCTGAGCGTCGGCCGACAGGGCCGGATGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGc                                                                                           <  8:324838/140‑1 (MQ=255)
                                                                  cAGACTTCGCGAAGCGTCAACGGCACCTGAGCGTCGGCCGACAGGGCCGGATGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGc                                                                                 >  6:123588/1‑149 (MQ=255)
                                                                     aCTTCGCGAAGCGTCAACGGCACCTGAGCGTCGGCCGACAGGGCCGGATGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCCAGGCCGCCTg                                                                              >  4:73270/1‑149 (MQ=255)
                                                                     aCTTCGCGAAGCGTCAACGGCACCTGAGCGTCGGCCGACAGGGCCGGATGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCACGTCAGCCTg                                                                              >  6:97742/1‑149 (MQ=255)
                                                                     aCTTCGCGAAGCGTCAACGGCACCTGAGCGTCGGCCGACAGGGCCGGATGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTg                                                                              >  4:333602/1‑149 (MQ=255)
                                                                                  tCAACTGCACCTGAGCGTCGGCCGACAGGGCCGGATGGCAACAGGTGAAGATCAGCCGCATGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGaa                                                                  <  7:168564/148‑1 (MQ=255)
                                                                                  tCAACGGCACCTGAGCGTCGGCCGACAGGGCCGGATGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCATGCCCCTCCAGCAGCATCATCAAGTCAGACTGGGATCGGtag                                                                    >  8:168564/1‑144 (MQ=255)
                                                                                    aaCGGCACCTGAGCGTCGGCCGACAGGGCCGGATGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCg                                                               >  3:85019/1‑149 (MQ=255)
                                                                                       ggCACCTGAGCGTCGGCCGACAGGGCCGGATGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGAc                                                            <  4:134037/149‑1 (MQ=255)
                                                                                              gAGCGTCGGCCGACAGGGCCGGATGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCgg                                                     >  6:339291/1‑149 (MQ=255)
                                                                                                      gCCGACAGGGCCGGATGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGGCGCAgcgc                                             >  4:212377/1‑149 (MQ=255)
                                                                                                        cGACAGGGCCGGATGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGGCGCAGCGCGt                                           <  3:333602/149‑1 (MQ=255)
                                                                                                            aGGGCCGGATGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGGCGCAGCGCGTCAAt                                       >  3:144973/1‑149 (MQ=255)
                                                                                                             gggCCGGATGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGGCGCAGCGCGTCAATg                                      <  2:165985/149‑1 (MQ=255)
                                                                                                             gggCCGGATGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGGCGCAGCGCGTCAATg                                      <  8:212644/149‑1 (MQ=255)
                                                                                                              ggCCGGATGGCAACGGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCg                                                      >  8:290587/1‑132 (MQ=255)
                                                                                                              ggCCGGATGGCAACGGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCg                                                      <  7:290587/132‑1 (MQ=255)
                                                                                                               gCCGGATGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGGCGCAGCGCGTCAATGGc                                    <  5:97742/149‑1 (MQ=255)
                                                                                                                 cGGATGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGGCGCAGCGCGTCAAt                                       <  4:177278/144‑1 (MQ=255)
                                                                                                                 cGGATGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGGCGCAGCGCGTCAAt                                       >  3:177278/1‑144 (MQ=255)
                                                                                                                  ggATGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGGCGCAGCGCGTCAATGGCCtt                                 <  3:212377/149‑1 (MQ=255)
                                                                                                                  ggATGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGGCGCAGCGCGTCAATGGCCtt                                 <  5:97221/149‑1 (MQ=255)
                                                                                                                   gATGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGGCGCAGCGCGTCAATGGCCtt                                 >  4:311592/1‑148 (MQ=255)
                                                                                                                    aTGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGGCGCAGCGCGTCAATGGCCTTg                                >  2:150368/1‑148 (MQ=255)
                                                                                                                     tGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGGCGCAGCGCGTCAATGGCCTTGaa                              <  4:192719/149‑1 (MQ=255)
                                                                                                                       gCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGGCGCAGCGCGTCAATGGCCTTGAAgc                            >  4:351513/1‑149 (MQ=255)
                                                                                                                        cAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGGCGCAGCGCGTCAATGGCCTTGAAgc                            >  2:181208/1‑148 (MQ=255)
                                                                                                                        cAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGGCGCAGCGCGTCAATGGCCTTGAAgc                            <  1:181208/148‑1 (MQ=255)
                                                                                                                         aaCAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGGCGCAGCGCGTCAATGGCCTTGAAgcgc                          <  7:183826/149‑1 (MQ=255)
                                                                                                                          tcAGGTGAATATCAGCCGCAGGCGGTCATCGGCCAGCAGCTCTTCTTCAGTGGTGTCCTGCTCCTGACCCTCCATCAGCATGCTCAAGTCAGCCTTGGAGCGGTCGAAGCGGACGCGCCGGCGCAGCGCGTCAATGGCCTTGAAgcgc                          <  4:30278/147‑1 (MQ=255)
                                                                                                                          cccagTGAAGAAGAGCTGCTGGCGGATGACCGCCTGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGgcgc                                                  >  4:292647/5‑124 (MQ=25)
                                                                                                                          cccagTGAAGAAGAGCTGCTGGCGGATGACCGCCTGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGgcgc                                                  <  3:292647/120‑1 (MQ=25)
                                                                                                                              gTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGGCGCAGCGCGTCAATGGCCTTGAAGCGCCCGa                      >  8:317195/1‑148 (MQ=255)
                                                                                                                              gTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGGCGCAGCGCGTCAATGGCCTTGAAGCGCCCGa                      >  1:274610/1‑148 (MQ=255)
                                                                                                                               tGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGGCGCAGCGCGTCAATGGCCTTGAAGCGCCCGACa                    <  3:311592/149‑1 (MQ=255)
                                                                                                                                 aaGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGa                                                             <  1:103642/106‑1 (MQ=255)
                                                                                                                                 aaGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGa                                                             >  2:103642/1‑106 (MQ=255)
                                                                                                                                 aaGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGAACGCGCCGGCGCAGCGCGTCAATGGCCTTGAAGCGCCCGACAg                   >  6:245549/1‑148 (MQ=255)
                                                                                                                                      cAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGGCGCAGCGCGTCAATGGCCTTGAAGCGCCCGACAGAAAc               >  3:189393/1‑147 (MQ=255)
                                                                                                                                        gCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGGCGCAGCGCGTCAATGGCCTTGAAGCGCCCGACAGAAACCAGc           >  2:359716/1‑149 (MQ=255)
                                                                                                                                          cGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGGCGCAGCGCGTCAATGGCCTTGAAGCGCCCGACAGAAACCAGCCa         >  8:45362/1‑149 (MQ=255)
                                                                                                                                           tcAGGCGTACATCAGCCAGCAGCTATTCTTCACTGGGGTCCTGCCCCTGACCCTCCAGCAGCATGCTCAAGTAAGCCTGGGAGAGGTCGAAGCGGACGCGCCGGCGCAGCGCGACAATGGCCTTGAACCGCCCGACAGAACCCAgaccc        >  8:266256/2‑145 (MQ=255)
                                                                                                                                              ggCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCAGGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGGCGCAGCGCGTCAATGGCCTTGCAGCGCCCGACAGAAACCAGCCACgat     >  5:90605/1‑147 (MQ=255)
                                                                                                                                                 ggTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGGCGCAGCGCGTGAATGGCCTTGAAGCGCCCGACAGAAAcc              <  4:316148/137‑1 (MQ=255)
                                                                                                                                                 ggTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGGCGCAGCGCGTGAATGGCCTTGAAGCGCCCGACAGAAAcc              >  3:316148/1‑137 (MQ=255)
                                                                                                                                                 ggTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGGCGCAGCGCGTCAATGGCCTTGAAGCGCCCGACAGAAACCAGCCACGCTcgc  <  4:144973/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ATACGCTGGGCGATGGTGGCCGGGCTTTGCAGAAAGGCGCGGGCGATCTGCTCGGTGGTCAGGTCGCAGACTTCGCGAAGCGTCAACGGCACCTGAGCGTCGGCCGACAGGGCCGGATGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGGCGCAGCGCGTCAATGGCCTTGAAGCGCCCGACAGAAACCAGCCACGCTCGC  >  NC_002947/5177530‑5177823

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑