Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,628,638 T→G 7.0% Y50D (TAC→GAC)  mucB → sigma factor AlgU regulator MucB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,628,6380TG 7.0% 216.0 / 5.6 71Y50D (TAC→GAC) mucBsigma factor AlgU regulator MucB
Reads supporting (aligned to +/- strand):  ref base T (37/29);  new base G (2/3);  total (39/32)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

CGCTACCTCTCTTGTCGCTGCTGATTGGCAGCTGCATGACGGTGCCAGCATTGGCGGCCAATTCTTCGCCTGAGGCGAGTGAGTGGCTGAACAAGCTGGCACAGGCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTCGACGACGGCAAGGTCAGCGAGCGGCTGTTGCAGCTTGATGGCGCCGCCCAGGAAATCGTGCGCGTGGATGGCAAGGTGGAATGTGTCAGCGGGGCCTTG  >  NC_002947/1628498‑1628784
                                                                                                                                             |                                                                                                                                                  
cgcTACCTCTCTTGTCGCTGCTGATTGGCAGCTGCATGACGGTGCCAGCATTGGCGGCCAATTCTTCGCCTGAGGCGAGTGAGTGGCTGAACAAGCTGGCACAGGCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTAGACGAACga                                                                                                                                            <  6:152400/149‑2 (MQ=255)
  cTACCTCTCTTGTCGCTGCTGATTGGCAGCTGCATGACGGTGCCAGCATTGGCGGCCAATTCTTCGCCTGAGGCGAGCGAGTGGCTGAACAAGCTGGCACAGGCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTCTACGAACGCaa                                                                                                                                          >  3:190237/1‑149 (MQ=255)
   tACCTCTCTTGTCGCTGCTGATTGGCAGCTGCATGACGGTGCCAGCATTGGCGGCCAATTCTTCGCGTGACGCGAGTCAGTGGCTGCACAAGCTGGCACAGGCCGAGCAAAAACAGAGCGATCAAGGGTCGTTCGTCTACGACCGCAAc                                                                                                                                         >  7:263091/1‑149 (MQ=255)
     ccTCTCTTGTCGCTGCTGATTGGCAGCTGCATGACGGTGCCAGCATTGGCGGCCAATTCTTCGCCTGAGGCGAGTGAGTGGCTGAACAAGCTGGCACAGGCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACgg                                                                                                                                       >  5:156111/1‑149 (MQ=255)
     ccTCTCTTGTCGCTGCTGATTGGCAGCTGCATGACGGTGCCAGCATTGGCGGCCAATTCTTCGCCTGAGGCGAGTGAGTGGCTGAACAAGCTGGCACAGGCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTAGACGAACgaccc                                                                                                                                         <  1:231699/147‑5 (MQ=255)
           ttGTCGCTGCTGATTGGCAGCTGCATGACGGTGCCAGCATTGGCGGCCAATTCTTCGCCTGAGGCGAGTGAGTGGCTGAACAAGCTGGCACAGGCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGttt                                                                                                                                 >  4:253386/1‑149 (MQ=255)
              tCGCTGCTGATTGGCAGCTGCATGACGGTGCCAGCATTGGCGGCCAATTCTTCGCCTGAGGCGAGTGAGTGGCTGAACAAGCTGGCACAGGCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTtctc                                                                                                                              <  1:5507/149‑1 (MQ=255)
                   gctgATTGGCAGCTGCATGACGGTGCCAGCATTGGCGGCCAATTCTTCGCCTGAGGCGAGTGAGTGGCTGAACAAGCTGGCACAGGCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTccc                                                                                                                         >  8:145650/1‑149 (MQ=255)
                   gctgATTGGCAGCTGCATGACGGTGCCAGCATTGGCGGCCAATTCTTCGCCTGAGGCGAGTGAGTGGCTGAACAAGCTGGCACAGGCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTAGACGAACgac                                                                                                                                           >  8:277537/1‑129 (MQ=255)
                   gctgATTGGCAGCTGCATGACGGTGCCAGCATTGGCGGCCAATTCTTCGCCTGAGGCGAGTGAGTGGCTGAACAAGCTGGCACAGGCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTAGACGAACgac                                                                                                                                           <  7:277537/131‑3 (MQ=255)
                      gATTGGCAGCTGCATGACGGTGCCAGCATTGGCGGCCAATTCTTCGCCTGAGGCGAGTGAGTGGCTGAACAAGCTGGCACAGGCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCAc                                                                                                                       >  4:230329/1‑148 (MQ=255)
                       aTTGGCAGCTGCATGACGGTGCCAGCATTGGCGGCCAATTCTTCGCCTGAGGCGAGTGAGTGGCTGAACAAGCTGGCACAGGCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGa                                                                                                                     >  5:210112/1‑149 (MQ=255)
                           gCAGCTGCATGACGGTGCCAGCATTGGCGGCCAATTCTTCGCCTGAGGCGAGTGAGTGGCTGAACAAGCTGGCACAGGCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATc                                                                                                                 >  6:315269/1‑149 (MQ=255)
                           gCAGCTGCATGACGGTGCCAGCATTGGCGGCCAATTCTTCGCCTGAGGCGAGTGAGTGGCTGAACAAGCTGGCACAGGCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATc                                                                                                                 >  4:349155/1‑149 (MQ=255)
                              gCTGCATGACGGTGCCAGCATTGGCGGCCAATTCTTCGCCTGAGGCGAGTGAGTGGCTGAACAAGCTGGCACAGGCCTAGCAAAAACAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTgg                                                                                                              <  4:61042/149‑1 (MQ=255)
                              gCTGCATGACGGTGCCAGCATTGGCGGCCAATTCTTCGCCTGAGGCGAGTGAGTGGCTGAACAAGCTGGCACAGGCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTgg                                                                                                              <  2:203816/149‑1 (MQ=255)
                                 gCATGACGGTGCCAGCATTGGCGGCCAATTCTTCGCCTGAGGCGAGTGAGTGGCTGAACAAGCTGGCACAGGCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGAGATCTGGccc                                                                                                           >  4:109599/1‑147 (MQ=255)
                                  cATGACGGTGCCAGCATTGGCGGCCAATTCTTCGCCTGAGGCGAGTGAGTGGCTGAACAAGCTGGCACAGGCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGca                                                                                                            >  6:83922/1‑147 (MQ=255)
                                  cATGACGGTGCCAGCATTGGCGGCCAATTCTTCGCCTGAGGCGAGTGAGTGGCTGAACAAGCTGGCACAGGCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGca                                                                                                            <  5:165254/147‑1 (MQ=255)
                                          tGCCAGCATTGGCGGCCAATTCTTCGCCTGAGGCGAGTGAGTGGCTGAACAAGCTGGCACAGGCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTcgac                                                                                                  <  4:167741/149‑1 (MQ=255)
                                                cATTGGCGGCCAATTCTTCGCCTGAGGCGAGTGAGTGGCTGAACAAGCTGGCACAGGCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTCGACGACgg                                                                                             >  7:78211/1‑148 (MQ=255)
                                                    ggcggcCAATTCTTCGCCTGAGGCGAGTGAGTGGCTGAACAAGCTGGCACAGGCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGGCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTCGACGACGGCAAg                                                                                         <  7:5795/148‑1 (MQ=255)
                                                           aaTTCTTCGCCTGAGGCGAGTGAGTGGCTGAACAAGCTGGCACAGGCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTCGACGACGGCAAGGTCagc                                                                                   >  3:34171/1‑147 (MQ=255)
                                                               cttcGCCTGAGGCGAGTGAGTGGCTGAACAAGCTGGCACAGGCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTCGACGACGGCAAGGTCagcgagcg                                                                              <  5:315269/148‑1 (MQ=255)
                                                               cttcGCCTGAGGCGAGTGAGTGGCTGAACAAGCTGGCACAGGCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTCGACGACGGCAAGGTCagcgagcg                                                                              <  8:189525/148‑1 (MQ=255)
                                                                ttcGCCTGAGGCGAGTGAGTGGCTGAACAAGCTGGCACAGGCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTCGACGACGGCAAg                                                                                         <  3:4305/136‑1 (MQ=255)
                                                                ttcGCCTGAGGCGAGTGAGTGGCTGAACAAGCTGGCACAGGCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTCGACGACGGCAAg                                                                                         >  4:4305/1‑136 (MQ=255)
                                                                     cTGAGGCGAGTGAGTGGCTGAACAAGCTGGCACAGGCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTCGACGACGGCAAGGTCa                                                                                     >  1:307979/1‑135 (MQ=255)
                                                                     cTGAGGCGAGTGAGTGGCTGAACAAGCTGGCACAGGCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTCGACGACGGCAAGGTCa                                                                                     <  2:307979/135‑1 (MQ=255)
                                                                     cTGAGGCGAGTGAGTGGCTGAACAAGCTGGCACAGGCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTCGACGACGGCAAGGTCAGCGAGCGGCTGTTg                                                                       >  7:230960/1‑149 (MQ=255)
                                                                      tGAGGCGAGTGAGTGGCTGAACAAGCTGGCACAGGCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTCGACGACGGCAAGGTCAGCGAGCGGCTGTTg                                                                       >  1:282216/1‑148 (MQ=255)
                                                                       gAGGCGAGTGAGTGGCTGAACAAGCTGGCACAGGCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTAGCCTCCCACGACATCTGGCACAGGGTCGACGACGGAAAGGTCAGCGAGCGGGTTTTg                                                                       >  8:295761/1‑147 (MQ=255)
                                                                        aGGCGAGTGAGTGGCTGAACAAGCTGGCACAGGCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTCGACGACGGCAAGGTCAGCGAGCGGCTGTTGCAg                                                                    <  8:78211/149‑1 (MQ=255)
                                                                           cGAGTGAGTGGCTGAACAAGCTGGCACAGGCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTCGACGACGGCAAGGTCAGCGAGCGGCTGTTGCAGCtt                                                                 >  7:360503/1‑149 (MQ=255)
                                                                                gagtgGCTGAACAAGCTGGCACAGGCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTCGACGACGGCAAGGTCAGCGAGCGGCTGTTGCAGCTTGATgg                                                            >  7:229487/1‑149 (MQ=255)
                                                                                  gtgGCTGAACAAGCTGGCACAGGCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTCGACGACGGCAAGGTCAGCGAGCGGCTGTTGCAGCTTGATGgcg                                                          >  7:102479/1‑149 (MQ=255)
                                                                                  gtgGCTGAACAAGCTGGCACAGGCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTCGACGACGGCAAGGTCAGCGAGCGGCTGTTGCAGCTTGATGgcg                                                          >  8:212385/1‑149 (MQ=255)
                                                                                   tgGCTGAACAAGCTGGCACAGGCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTCGACGACGGCAAGGTCAGCGAGCGGCTGTTGCAGCTTGATGgcgc                                                         <  3:230329/149‑1 (MQ=255)
                                                                                      cTGAACAAGCTGGCACAGGCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTCGACGACGGCAAGGTCAGCGAGCGGCTGTTGCAGCTTGATGGcgccgc                                                      <  6:156111/149‑1 (MQ=255)
                                                                                         aaCAAGCTGGCACAGGCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTCGACGACGGCAAGGTCAGCGAGCGGCTGTTGCAGCTTGATGGCgccgcc                                                     >  7:4485/1‑147 (MQ=255)
                                                                                          aCAAGCTGGCACAGGCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTCGACGACGGCAAGGTCAGCGAGCGGCTGTTGCAGCTTGATGGCGCCGCCCAg                                                  >  7:103245/1‑149 (MQ=255)
                                                                                            aaGCTGGCACAGGCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTCGACGACGGCAAGGTCAGCGAGCGGCTGTTGCAGCTTGATGgcgc                                                         >  1:25786/1‑140 (MQ=255)
                                                                                            aaGCTGGCACAGGCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTCGACGACGGCAAGGTCAGCGAGCGGCTGTTGCAGCTTGATGgcgc                                                         <  2:25786/140‑1 (MQ=255)
                                                                                              gCTGGCACAGGCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCtcctcc                                                                                                                          <  6:102489/73‑1 (MQ=255)
                                                                                              gCTGGCACAGGCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCtcctcc                                                                                                                          >  5:102489/1‑73 (MQ=255)
                                                                                                tGGCACAGGCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTCGACGACGGCAAGGTCAGCGAGCGGCTGTTGCAGCTTGATGGCGCCGCCCAGGAAATc                                            >  5:306446/1‑149 (MQ=255)
                                                                                                  gCACAGGCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTCGACGACGGCAAGGTCAGCGAGCGGCTGTTGCAGCTTGATGGCGCCGCCCAGGAAATc                                            <  6:210112/147‑1 (MQ=255)
                                                                                                  gCACAGGCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTCGACGACGGCAAGGTCAGCGAGCGGCTGTTGCAGCTTGATGGCGCCGCCCAGGAAATc                                            <  8:229487/147‑1 (MQ=255)
                                                                                                    acaGGCCGAGCAAAAACAGAGCTATCAAGGGATCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTCGACGACGGCAAGGTCAGCGAGCGGCTGTTGCAGCTTGATGGCGCCGCCCAGGAAATCGTg                                         >  1:42148/1‑149 (MQ=255)
                                                                                                        gCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTCGACGACGGCAAGGTCAGCGAGCgg                                                                             >  6:346281/1‑108 (MQ=255)
                                                                                                        gCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTCGACGACGGCAAGGTCAGCGAGCgg                                                                             <  5:346281/108‑1 (MQ=255)
                                                                                                        gCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTCGACGACGGCAAGGTCAGCGAGCGGCTGTTGCAGCTTGATGGCGCCGCCCAGGAAATCGTGCGCGt                                    <  4:34171/149‑1 (MQ=255)
                                                                                                          cGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTCGACGACGGCAAGGTCAGCGAGCGGCTGTTGCAGCTTGATGGCGCCGCCCAGGAAATCGTGCGCGTg                                   >  2:261852/1‑148 (MQ=255)
                                                                                                             gCAAAAACAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTCGACGACGGCAAGGTCAGCGAGCGGCTGTTGCAGCTTGATGGCGCCGCCCAGGAAATCGTGCGCGTGGATg                               <  3:253386/149‑1 (MQ=255)
                                                                                                                 aaaCAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTCGACGACGGCAAGGTCAGCGAGCGGCTGTTGCAGCTTGATGGCGCCGCCCAGGAAATCGTGCGCGTGGATgga                             >  4:291147/1‑146 (MQ=255)
                                                                                                                 aaaCAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTCGACGACGGCAAGGTCAGCGAGCGGCTGTTGCAGCTTGATGGCGCCGCCCAGGAAATCGTGCGCGTGGATGGCaa                           <  4:269352/149‑1 (MQ=255)
                                                                                                                     agagCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTCGACGACGGAAAGGTCAGCGAGCGGCTGTTGCAGCGTGATGGCGCCGCCCAGGAACTCGTGCGCGGGGATGGCAAGGTg                       >  8:271686/1‑149 (MQ=255)
                                                                                                                        gCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTCGACGACGGCAAGGTCAGCGAGCGGCTGTTGCAGCTTGATGGCGCCGCCCAGGAAATCGTGCGCGTGGATGGCAAGGTGGaa                    <  7:295761/149‑1 (MQ=255)
                                                                                                                         cTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTCGACGAGGGCAAGGTCAGCGAGCGGCt                                                                           >  7:183144/1‑93 (MQ=255)
                                                                                                                         cTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTCGACGACGGCAAGGTCAGCGAGCGGCt                                                                           <  8:183144/93‑1 (MQ=255)
                                                                                                                         cTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTCGACGACGGCAAGGTCAGCGAGCGGCTGTTGCAGCTTGATGGCGCCGCCCAGGAAATCGTGCGCGTGGATGGCAAGGTGGAAt                   <  2:42148/149‑1 (MQ=255)
                                                                                                                          tATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTCGACGACGGCAAGGTCAGCGAGCGGCTGTTGCAGCTTGATGGCGCCGCCCAGGAAATCGTGCGCGTGGATGGCAAGGTGGAAtg                  >  5:364366/1‑149 (MQ=255)
                                                                                                                            tCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTCGACGACGGCAAGGTCAGCGAGCGGCTGTTGCAGCTTGATGTCGCCGCCCAGGAAATCGTGCGCGTTGATGGCAAGGTTGAAtgtg                <  8:103245/149‑1 (MQ=255)
                                                                                                                            tCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTCGACGACGGCAAGGTCAGCGAGCGGCTGTTGCAGCTTGATGGCGCCGCCCAGGAAATCGTGCGCGTGGATGGCAAGGTGGAAtgtg                >  7:348066/1‑149 (MQ=255)
                                                                                                                              aaGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTCGACGACGGCAAGGTCAGCGAGCGGCTGTTGCAGCTTGATGGCGCCGCCCAGGAAATCGTGCGCGTGGATGGCAAGGTGGAATgtgt               <  1:261852/148‑1 (MQ=255)
                                                                                                                               aGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTCGACGACGGCAAGGTCAGCGAGCGGCTGTTGCAGCTTGATGGCGCCGCCCAGGAAATCGTGCGCGTGGATGGCAAGGTGGAATGTGTCa             <  7:212385/149‑1 (MQ=255)
                                                                                                                                 gttgcgttcgtAGACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTCGACGACGGCAAGGTCAGCGAGCGGCTGTTGCAGCTTGATGGCGCCGCCCAGGAAATCGTGCGCGTGGATGGCAAGGTGGAATGTGTCAg            >  2:392512/5‑149 (MQ=255)
                                                                                                                                 ggTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTCGACGACGGCAAGGTCAGCGAGCGGCTGTTGCAGCTTGATGGCGCCGCCCAGGAAATCGTGCGCGTGGATGGCAAGGTGGAATGTGTCAGc           <  8:360503/149‑1 (MQ=255)
                                                                                                                                       ttcgtCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTCGACGACGGCAAGGTCAGCGAGCGGCTGTTGCAGCTTGATGGCGCCGCCCAGGAAATCGTGCGCGTGGATGGCAAGGTGGAATGTGTCAGCGGGGc      >  6:278823/1‑149 (MQ=255)
                                                                                                                                       ttcgtCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTCGACGACGGCAAGGTCAGCGAGCGGCTGTTGCAGCTTGATGGCGCCGCCCAGGAAATCGTGCGCGTGGATGGCAAGGTGGAATGTGTCAGCGGGGc      >  6:125102/1‑149 (MQ=255)
                                                                                                                                           tCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTCGACGACGGCAAGGTCAGCGAGCGGCTGTTGCAGCTTGATGGCGCCGCCCAGGAAATCGTGCGCGTGGATGGCAAGGTGGAATGTGTCAGCGGGGCCTTg  <  7:145650/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                  
CGCTACCTCTCTTGTCGCTGCTGATTGGCAGCTGCATGACGGTGCCAGCATTGGCGGCCAATTCTTCGCCTGAGGCGAGTGAGTGGCTGAACAAGCTGGCACAGGCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTCGACGACGGCAAGGTCAGCGAGCGGCTGTTGCAGCTTGATGGCGCCGCCCAGGAAATCGTGCGCGTGGATGGCAAGGTGGAATGTGTCAGCGGGGCCTTG  >  NC_002947/1628498‑1628784

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