Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,628,637 C→A 7.1% V49V (GTC→GTA 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,6370CA 7.1% 218.3 / 6.4 71V49V (GTC→GTAmucBsigma factor AlgU regulator MucB
Reads supporting (aligned to +/- strand):  ref base C (37/29);  new base A (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 = 1.47e-01

TCATGCGCGCGCTACCTCTCTTGTCGCTGCTGATTGGCAGCTGCATGACGGTGCCAGCATTGGCGGCCAATTCTTCGCCTGAGGCGAGTGAGTGGCTGAACAAGCTGGCACAGGCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTCGACGACGGCAAGGTCAGCGAGCGGCTGTTGCAGCTTGATGGCGCCGCCCAGGAAATCGTGCGCGTGGATGGCAAGGTGGAATGTGTCAGCGGGGCCTTG  >  NC_002947/1628489‑1628784
                                                                                                                                                     |                                                                                                                                                   
tCATGCGCGCGCTACCTCTCTTGTCGCTGCTGATTGGCAGCTGCATGACGGTGCCAGCATTGGCGGCCAATTCTTCGGCTGAGGCGAGTGAGTGGCTGAACAAGCTGGCACAGGCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTc                                                                                                                                                     >  1:124375/1‑149 (MQ=255)
tCATGCGCGCGCTACCTCTCTTGTCGCTGCTGATTGGCAGCTGCATGACGGTGCCAGCATTGGCGGCCAATTCTTCGCCTGAGGCGAGTGAGTGGCTGAACAAGCTGGCACAGGCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTc                                                                                                                                                     <  4:19205/149‑1 (MQ=255)
tCATGCGCGCGCTACCTCTCTTGTCGCTGCTGATTGGCAGCTGCATGACGGTGCCAGCATTGGCGGCCAATTCTTCGCCTGAGGCGAGTGAGTGGCTGAACAAGCTGGCACAGGCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTc                                                                                                                                                     <  2:41072/149‑1 (MQ=255)
         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                                                                                                                 >  4:349155/1‑149 (MQ=255)
                                    gCAGCTGCATGACGGTGCCAGCATTGGCGGCCAATTCTTCGCCTGAGGCGAGTGAGTGGCTGAACAAGCTGGCACAGGCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATc                                                                                                                 >  6:315269/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                                                                                                            <  5:165254/147‑1 (MQ=255)
                                           cATGACGGTGCCAGCATTGGCGGCCAATTCTTCGCCTGAGGCGAGTGAGTGGCTGAACAAGCTGGCACAGGCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGca                                                                                                            >  6:83922/1‑147 (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                                                                              <  8:189525/148‑1 (MQ=255)
                                                                        cttcGCCTGAGGCGAGTGAGTGGCTGAACAAGCTGGCACAGGCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTCGACGACGGCAAGGTCagcgagcg                                                                              <  5:315269/148‑1 (MQ=255)
                                                                         ttcGCCTGAGGCGAGTGAGTGGCTGAACAAGCTGGCACAGGCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTCGACGACGGCAAg                                                                                         <  3:4305/136‑1 (MQ=255)
                                                                         ttcGCCTGAGGCGAGTGAGTGGCTGAACAAGCTGGCACAGGCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTCGACGACGGCAAg                                                                                         >  4:4305/1‑136 (MQ=255)
                                                                              cTGAGGCGAGTGAGTGGCTGAACAAGCTGGCACAGGCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTCGACGACGGCAAGGTCa                                                                                     <  2:307979/135‑1 (MQ=255)
                                                                              cTGAGGCGAGTGAGTGGCTGAACAAGCTGGCACAGGCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTCGACGACGGCAAGGTCa                                                                                     >  1:307979/1‑135 (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                                                          >  8:212385/1‑149 (MQ=255)
                                                                                           gtgGCTGAACAAGCTGGCACAGGCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTCGACGACGGCAAGGTCAGCGAGCGGCTGTTGCAGCTTGATGgcg                                                          >  7:102479/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                                            <  8:229487/147‑1 (MQ=255)
                                                                                                           gCACAGGCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTCGACGACGGCAAGGTCAGCGAGCGGCTGTTGCAGCTTGATGGCGCCGCCCAGGAAATc                                            <  6:210112/147‑1 (MQ=255)
                                                                                                             acaGGCCGAGCAAAAACAGAGCTATCAAGGGATCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTCGACGACGGCAAGGTCAGCGAGCGGCTGTTGCAGCTTGATGGCGCCGCCCAGGAAATCGTg                                         >  1:42148/1‑149 (MQ=255)
                                                                                                                 gCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTCGACGACGGCAAGGTCAGCGAGCgg                                                                             <  5:346281/108‑1 (MQ=255)
                                                                                                                 gCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTCGACGACGGCAAGGTCAGCGAGCgg                                                                             >  6:346281/1‑108 (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)
                                                                                                                                                     |                                                                                                                                                   
TCATGCGCGCGCTACCTCTCTTGTCGCTGCTGATTGGCAGCTGCATGACGGTGCCAGCATTGGCGGCCAATTCTTCGCCTGAGGCGAGTGAGTGGCTGAACAAGCTGGCACAGGCCGAGCAAAAACAGAGCTATCAAGGGTCGTTCGTCTACGAACGCAACGGCAGTTTCTCCTCCCACGACATCTGGCACAGGGTCGACGACGGCAAGGTCAGCGAGCGGCTGTTGCAGCTTGATGGCGCCGCCCAGGAAATCGTGCGCGTGGATGGCAAGGTGGAATGTGTCAGCGGGGCCTTG  >  NC_002947/1628489‑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: