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

CGATGGTGGCCGGGCTTTGCAGAAAGGCGCGGGCGATCTGCTCGGTGGTCAGGTCGCAGACTTCGCGAAGCGTCAACGGCACCTGAGCGTCGGCCGACAGGGCCGGATGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGGCGCAGCGCGTCAATGGCCTTGAAGCGCCCGACAGAAACCAGCCACGCTCGCGGGTTGGCGG  >  NC_002947/5177540‑5177833
                                                                                                                                                  |                                                                                                                                                   
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                                                <  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                                                 >  3:177278/1‑144 (MQ=255)
                                                                                                       cGGATGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGGCGCAGCGCGTCAAt                                                 <  4:177278/144‑1 (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                                >  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)
                                                                                                                                    ggCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCAGGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGGCGCAGCGCGTCAATGGCCTTGCAGCGCCCGACAGAAACCAGCCACgat               >  5:90605/1‑147 (MQ=255)
                                                                                                                                       ggTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGGCGCAGCGCGTGAATGGCCTTGAAGCGCCCGACAGAAAcc                        >  3:316148/1‑137 (MQ=255)
                                                                                                                                       ggTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGGCGCAGCGCGTGAATGGCCTTGAAGCGCCCGACAGAAAcc                        <  4:316148/137‑1 (MQ=255)
                                                                                                                                       ggTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGGCGCAGCGCGTCAATGGCCTTGAAGCGCCCGACAGAAACCAGCCACGCTcgc            <  4:144973/149‑1 (MQ=255)
                                                                                                                                           aTCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGGCGCAGCGCGTCAATGGCCTTGAAgcgc                                    <  3:223840/121‑1 (MQ=255)
                                                                                                                                           aTCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGGCGCAGCGCGTCAATGGCCTTGAAgcgc                                    >  4:223840/1‑121 (MQ=255)
                                                                                                                                           aTCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGGCGCAGCGCGTCAATGGCCTTGAAGCGCCCGACAGAAACCAGCCACGCTCGCggg         <  1:359716/148‑1 (MQ=255)
                                                                                                                                               gCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGGCGCAGCGCGTCAATGGCCTTGAAGCGCCCGACAGAAACCAGCCACGCTCGCGGGTTGGc    >  8:212399/1‑149 (MQ=255)
                                                                                                                                                ccAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGGCGCAGCGCGTCAATGGCCTTGAAGCGCCCGACAGAAACCAGCCACGCTCGCGGGTTGGCg   >  5:279494/1‑149 (MQ=255)
                                                                                                                                                 cagcatCTCTTCTTCAGTGGGGTACTGCGACTGACCCTCCAGCAGCATGCTCAAGTCATCCTGGGAGCGGTAGAAGCGGACGCGCCGGCGCAGCGCGTCAATGGCCTTGAAGCGCCCGACAGAAACCAGCCACGCTCGCGGGTTGGCgg  <  6:90605/149‑1 (MQ=255)
                                                                                                                                                 cagcagCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGGCGCAGCGCGTCAATGGCCTTGAAGCGCCCGACAGAAACCAGCCACGCTCGCGGGTTGGCgg  <  4:85019/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CGATGGTGGCCGGGCTTTGCAGAAAGGCGCGGGCGATCTGCTCGGTGGTCAGGTCGCAGACTTCGCGAAGCGTCAACGGCACCTGAGCGTCGGCCGACAGGGCCGGATGGCAACAGGTGAAGATCAGCCGCAGGCGGTCATCCGCCAGCAGCTCTTCTTCACTGGGGTCCTGCGCCTGACCCTCCAGCAGCATGCTCAAGTCAGCCTGGGAGCGGTCGAAGCGGACGCGCCGGCGCAGCGCGTCAATGGCCTTGAAGCGCCCGACAGAAACCAGCCACGCTCGCGGGTTGGCGG  >  NC_002947/5177540‑5177833

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