Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,177,671 A→T 5.7% L105Q (CTG→CAG)  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,6710AT 5.7% 208.3 / 2.6 70L105Q (CTG→CAG) PP_4553ECF family RNA polymerase sigma‑70 factor
Reads supporting (aligned to +/- strand):  ref base A (36/30);  new base T (2/2);  total (38/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 = 2.33e-01

GCACGATACGCTGGGCGATGGTGGCCGGGCTTTGCAGAAAGGCGCGGGCGATCTGCTCGGTGGTCAGGTCGCAGACTTCGCGAAGCGTCAACGGCACCTGAGCGTCGGCCGACAGGGCCGGATGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGGCGCAGCGCGTCAATGGCCTTGAAGCGCCCGACAGAAACCAGCCAC  >  NC_002947/5177525‑5177817
                                                                                                                                                  |                                                                                                                                                  
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                      >  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                                            <  3:292647/120‑1 (MQ=25)
                                                                                                                               cccagTGAAGAAGAGCTGCTGGCGGATGACCGCCTGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGgcgc                                            >  4:292647/5‑124 (MQ=25)
                                                                                                                                   gTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGGCGCAGCGCGTCAATGGCCTTGAAGCGCCCGa                >  1:274610/1‑148 (MQ=255)
                                                                                                                                   gTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGGCGCAGCGCGTCAATGGCCTTGAAGCGCCCGa                >  8:317195/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)
                                                                                                                                                  |                                                                                                                                                  
GCACGATACGCTGGGCGATGGTGGCCGGGCTTTGCAGAAAGGCGCGGGCGATCTGCTCGGTGGTCAGGTCGCAGACTTCGCGAAGCGTCAACGGCACCTGAGCGTCGGCCGACAGGGCCGGATGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGGCGCAGCGCGTCAATGGCCTTGAAGCGCCCGACAGAAACCAGCCAC  >  NC_002947/5177525‑5177817

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