Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,177,663 T→A 6.0% I108F (ATC→TTC)  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,6630TA 6.0% 202.3 / 2.9 67I108F (ATC→TTC) PP_4553ECF family RNA polymerase sigma‑70 factor
Reads supporting (aligned to +/- strand):  ref base T (32/31);  new base A (2/2);  total (34/33)
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.41e-01

GGCGCGCACGATACGCTGGGCGATGGTGGCCGGGCTTTGCAGAAAGGCGCGGGCGATCTGCTCGGTGGTCAGGTCGCAGACTTCGCGAAGCGTCAACGGCACCTGAGCGTCGGCCGACAGGGCCGGATGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGGCGCAGCGCGTCAATGGCCTTGAAGCGCCCGACAG  >  NC_002947/5177520‑5177806
                                                                                                                                               |                                                                                                                                               
ggcGCGCACGATACGCTGGGCGATGGTGGCCGGGCTTTGCAGAAAGGCGCGGGCGATCTGCTCGGTGGTCAGGTCGCAGACTTCGCGAAGCGTCAACGGCACCTGAGCGTCGGCCGACAGGGCCGGATGGCATCAGGTGAAGATCAg                                                                                                                                              <  5:177466/147‑1 (MQ=255)
     gcACGATACGCTGGGCGATAGTGGCCGGGCTTTGCAGAAAGGCGCGGGCGATCTGCTCGGTGGTCAGGTCGCAGACTTCGCGAAGCGTCAACGGCACCTGAGCGTCGGCCGACAGGGCCGGATGGCAACAGGTGAAGATCAGCCGCAgg                                                                                                                                       <  5:59514/149‑1 (MQ=255)
        cGATACGCTGGGCGATGGTGGCCGGGCTTTGCAGAAAGGCGCGGGCGATCTGCTCGGTGGTCAGGTCGCAGACTTCGCGAAGCGTCAACGGCACCTGAGCGTCGGCCGACAGGGCCGGATGGCAACAGGTGAAGATCAGCCGCAGGCgg                                                                                                                                    <  6:257607/149‑1 (MQ=255)
          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                                                                                                     >  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                                                                           <  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                                     <  7:290587/132‑1 (MQ=255)
                                                                                                                        ggCCGGATGGCAACGGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCg                                     >  8:290587/1‑132 (MQ=255)
                                                                                                                         gCCGGATGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGGCGCAGCGCGTCAATGGc                   <  5:97742/149‑1 (MQ=255)
                                                                                                                           cGGATGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGGCGCAGCGCGTCAAt                      >  3:177278/1‑144 (MQ=255)
                                                                                                                           cGGATGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGGCGCAGCGCGTCAAt                      <  4:177278/144‑1 (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           <  1:181208/148‑1 (MQ=255)
                                                                                                                                  cAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGGCGCAGCGCGTCAATGGCCTTGAAgc           >  2:181208/1‑148 (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     >  1:274610/1‑148 (MQ=255)
                                                                                                                                        gTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGGCGCAGCGCGTCAATGGCCTTGAAGCGCCCGa     >  8:317195/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)
                                                                                                                                               |                                                                                                                                               
GGCGCGCACGATACGCTGGGCGATGGTGGCCGGGCTTTGCAGAAAGGCGCGGGCGATCTGCTCGGTGGTCAGGTCGCAGACTTCGCGAAGCGTCAACGGCACCTGAGCGTCGGCCGACAGGGCCGGATGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGGCGCAGCGCGTCAATGGCCTTGAAGCGCCCGACAG  >  NC_002947/5177520‑5177806

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