Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,177,679 A→G 5.6% D102D (GAT→GAC)  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,6790AG 5.6% 210.9 / 2.1 71D102D (GAT→GAC) PP_4553ECF family RNA polymerase sigma‑70 factor
Reads supporting (aligned to +/- strand):  ref base A (37/30);  new base G (2/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.21e-01

GCTGGGCGATGGTGGCCGGGCTTTGCAGAAAGGCGCGGGCGATCTGCTCGGTGGTCAGGTCGCAGACTTCGCGAAGCGTCAACGGCACCTGAGCGTCGGCCGACAGGGCCGGATGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGGCGCAGCGCGTCAATGGCCTTGAAGCGCCCGACAGAAACCAGCCACGCTCGCGGG  >  NC_002947/5177534‑5177826
                                                                                                                                                 |                                                                                                                                                   
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                                                                                                                         >  6:153730/1‑141 (MQ=21)
                           gAAAGGCGCGGGCGATCTGCTCGGTGGTCAGGTCGCAGACTTCGCGAAGCGTCAACGGCACCTGAGCGTCGGCCGACAGGGCCGGATGGCAACAGGTGAAGAAGAGCTGCTGGCGGATGACCGCCTGCGGCTGATCTTCAcctgttg                                                                                                                         <  5:153730/147‑7 (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                                                                                               <  4:261701/149‑1 (MQ=255)
                                                   tggtCAGGTCGCAGACTTCGCGAAGCGTCAACGGCACCTGAGCGTCGGCCGACAGGGCCGGATGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATg                                                                                               <  8:156540/149‑1 (MQ=255)
                                                          gTCGCAGACTTCGCGAAGCGTCAACGGCACCTGAGCGTCGGCCGACAGGGCCGGATGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTg                                                                                                                             >  4:271036/1‑112 (MQ=255)
                                                          gTCGCAGACTTCGCGAAGCGTCAACGGCACCTGAGCGTCGGCCGACAGGGCCGGATGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTg                                                                                                                             <  3:271036/112‑1 (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                                         <  8:212644/149‑1 (MQ=255)
                                                                                                         gggCCGGATGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGGCGCAGCGCGTCAATg                                         <  2:165985/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                                    <  5:97221/149‑1 (MQ=255)
                                                                                                              ggATGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGGCGCAGCGCGTCAATGGCCtt                                    <  3:212377/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                                                                >  2:103642/1‑106 (MQ=255)
                                                                                                                             aaGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGa                                                                <  1:103642/106‑1 (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)
                                                                                                                                                 aTCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGGCGCAGCGCGTCAATGGCCTTGAAgcgc                             >  4:223840/1‑121 (MQ=255)
                                                                                                                                                 aTCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGGCGCAGCGCGTCAATGGCCTTGAAgcgc                             <  3:223840/121‑1 (MQ=255)
                                                                                                                                                 aTCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGGCGCAGCGCGTCAATGGCCTTGAAGCGCCCGACAGAAACCAGCCACGCTCGCggg  <  1:359716/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GCTGGGCGATGGTGGCCGGGCTTTGCAGAAAGGCGCGGGCGATCTGCTCGGTGGTCAGGTCGCAGACTTCGCGAAGCGTCAACGGCACCTGAGCGTCGGCCGACAGGGCCGGATGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGGCGCAGCGCGTCAATGGCCTTGAAGCGCCCGACAGAAACCAGCCACGCTCGCGGG  >  NC_002947/5177534‑5177826

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