Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,467,259 A→T 6.3% D401V (GAT→GTT) ‡ PP_4798 → membrane‑bound lytic murein transglycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,467,2590AT 6.3% 1119.2 / 42.7 335D401V (GAT→GTT) ‡PP_4798membrane‑bound lytic murein transglycosylase
Reads supporting (aligned to +/- strand):  ref base A (123/191);  new base T (14/7);  total (137/198)
Fisher's exact test for biased strand distribution p-value = 2.00e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

Only 100 of 342 total aligned reads displayed.

CCTGGCGGTGAGCCTGCTCGGGGATCGCTTCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCT  >  NC_002947/5467110‑5467408
                                                                                                                                                     |                                                                                                                                                     
ccTGGCGGTGAGCCTGCTCGGGGATCGCTTCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGa                                                                                                                                                       <  2:8754788/150‑1 (MQ=255)
ccTGGCGGTGAGCCTGCTCGGGGATCGCTTCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGa                                                                                                                                                       <  2:3617709/150‑1 (MQ=255)
 cTGGCGGTGAGCCTGCTCGGGGATCGCTTCAGCGGCTGGGGCTTCATTGCCGGCAGCGGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGAt                                                                                                                                                      <  2:8199131/150‑1 (MQ=255)
    gCGGTGAGCCTGCTCGGGGATCGCTTCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGc                                                                                                                                                   <  2:6951041/150‑1 (MQ=255)
         gAGCCTGCTCGGGGATCGCTTCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCATCGGCATTac                                                                                                                                              <  1:3806358/150‑2 (MQ=255)
          aGCCTGCTCGGGGATCGCTTCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTacc                                                                                                                                             <  2:4840008/150‑3 (MQ=255)
          aGCCTGCTCGGGGATCGCTTCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATc                                                                                                                                             >  2:2986536/1‑150 (MQ=255)
          aGCCTGCTCGGGGATCGCTTCAGCGGCAGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCATCGGCATTacc                                                                                                                                             <  2:141087/150‑3 (MQ=255)
            ccTGCTCGGGGATCGCTTCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCgg                                                                                                                                           <  1:3736389/150‑1 (MQ=255)
            ccTGCTCGGGGATCGCTTCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCgg                                                                                                                                           >  2:60078/1‑150 (MQ=255)
            ccTGCTCGGGGATCGCTTCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCgg                                                                                                                                           >  2:2722131/1‑150 (MQ=255)
            ccTGCTCGGGGATCGCTTCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCgg                                                                                                                                           >  2:4494262/1‑150 (MQ=255)
               gCTCGGGGATCGCTTCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGc                                                                                                                                        <  2:5269075/150‑1 (MQ=255)
               gCTCGGGGATCGCTTCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGc                                                                                                                                        <  2:565319/150‑1 (MQ=255)
                 tCGGGGATCGCTTCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATAGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCa                                                                                                                                      <  1:989602/150‑1 (MQ=255)
                     ggATCGCTTCAGCGGCTGGGGCTTCATTGCCGGCAGCTGACCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATAc                                                                                                                                  <  1:7380823/150‑1 (MQ=255)
                      gATCGCTTCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATAcc                                                                                                                                 >  1:6536966/1‑150 (MQ=255)
                       aTCGCTTCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATAccc                                                                                                                                >  1:396346/1‑150 (MQ=255)
                       aTCGCTTCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATAccc                                                                                                                                >  1:1150460/1‑150 (MQ=255)
                       aTCGCTTCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATAccc                                                                                                                                <  2:2230265/150‑1 (MQ=255)
                        tCGCTTCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCg                                                                                                                               >  2:3498631/1‑150 (MQ=255)
                           cTTCAGCGGCTGGGGCTTCAGTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCaa                                                                                                                            <  2:6667243/150‑1 (MQ=255)
                            ttCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAg                                                                                                                           <  1:1443614/150‑1 (MQ=255)
                            ttCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAg                                                                                                                           <  1:5432116/150‑1 (MQ=255)
                                  ggCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATc                                                                                                                     <  1:1276550/150‑1 (MQ=255)
                                   gCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATcc                                                                                                                    <  2:4396676/150‑1 (MQ=255)
                                    cTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCg                                                                                                                   <  2:772963/150‑1 (MQ=255)
                                       gggCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCaa                                                                                                                >  1:7793060/1‑150 (MQ=255)
                                        ggCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAAt                                                                                                               >  1:2545719/1‑150 (MQ=255)
                                        ggCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAAt                                                                                                               <  1:2888925/150‑1 (MQ=255)
                                           ttCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGcc                                                                                                            >  2:367662/1‑150 (MQ=255)
                                             cATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCTTTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCa                                                                                                          <  2:640090/150‑1 (MQ=255)
                                             cATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCa                                                                                                          <  2:114403/150‑1 (MQ=255)
                                              aTTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAg                                                                                                         >  1:4883513/1‑150 (MQ=255)
                                              aTTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAg                                                                                                         >  1:4150089/1‑150 (MQ=255)
                                               ttGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGc                                                                                                        <  1:1597715/150‑1 (MQ=255)
                                                   cGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGAGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAgg                                                                                                    >  1:5121128/1‑150 (MQ=255)
                                                   cGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACCTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAgg                                                                                                    >  1:4802219/1‑150 (MQ=255)
                                                    ggCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCCGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAggg                                                                                                   >  2:2499123/1‑150 (MQ=255)
                                                       aGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTg                                                                                                <  2:3890305/150‑1 (MQ=255)
                                                        gCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTgg                                                                                               >  1:5302337/1‑150 (MQ=255)
                                                        gCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTgg                                                                                               <  2:7349967/150‑1 (MQ=255)
                                                         cTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTggg                                                                                              <  1:8603976/150‑1 (MQ=255)
                                                                     ggACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGa                                                                                  <  1:4434298/150‑1 (MQ=255)
                                                                      gACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGAt                                                                                 <  1:802540/150‑1 (MQ=255)
                                                                         cTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGc                                                                              >  2:6566752/1‑150 (MQ=255)
                                                                           gCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTa                                                                            >  1:6912171/1‑150 (MQ=255)
                                                                            ccGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTAc                                                                           >  1:8041289/1‑150 (MQ=255)
                                                                            ccGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTAc                                                                           <  1:2241039/150‑1 (MQ=255)
                                                                              gCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTAccc                                                                         <  1:2056647/150‑1 (MQ=255)
                                                                              gCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTAccc                                                                         >  1:6885702/1‑150 (MQ=255)
                                                                                 cAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGAc                                                                      <  1:3545874/150‑1 (MQ=255)
                                                                                    ccGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTgg                                                                                               >  1:7934202/1‑122 (MQ=255)
                                                                                    ccGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCca                                                                   >  1:7198254/1‑150 (MQ=255)
                                                                                     cGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCcac                                                                  >  2:8297998/1‑150 (MQ=255)
                                                                                      gCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAACACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCcaca                                                                 >  2:1486545/1‑150 (MQ=255)
                                                                                       cACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACaa                                                                <  2:5037198/150‑1 (MQ=255)
                                                                                       cACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACaa                                                                <  1:5839771/150‑1 (MQ=255)
                                                                                            aGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAgctgc                                                           <  1:8494248/150‑1 (MQ=255)
                                                                                            aGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAgctgc                                                           <  2:4567796/150‑1 (MQ=255)
                                                                                            aGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAgctgc                                                           <  1:2750823/150‑1 (MQ=255)
                                                                                              cgcATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCtt                                                         <  1:4481369/150‑1 (MQ=255)
                                                                                                  tGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTgaga                                                     <  2:6852364/150‑1 (MQ=255)
                                                                                                  tGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTgaga                                                     <  2:1098690/150‑1 (MQ=255)
                                                                                                  tGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTgaga                                                     <  1:4935172/150‑1 (MQ=255)
                                                                                                  tGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTgaga                                                     <  2:5769335/150‑1 (MQ=255)
                                                                                                  tGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTgaga                                                     <  1:8299775/150‑1 (MQ=255)
                                                                                                  tGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTgaga                                                     <  2:1247277/150‑1 (MQ=255)
                                                                                                  tGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTgaga                                                     <  2:7963815/150‑1 (MQ=255)
                                                                                                  tGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTgaga                                                     <  2:556040/150‑1 (MQ=255)
                                                                                                     aGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGcc                                                  <  1:1382533/150‑1 (MQ=255)
                                                                                                      gTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCt                                                 <  2:2730582/150‑1 (MQ=255)
                                                                                                          cAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTgcgc                                             <  1:6644011/150‑1 (MQ=255)
                                                                                                           aGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGcc                                            >  2:8026574/1‑150 (MQ=255)
                                                                                                           aGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGcc                                            >  2:1736911/1‑150 (MQ=255)
                                                                                                             aaCTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCcag                                          >  2:6298985/1‑150 (MQ=255)
                                                                                                             aaCTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCcag                                          <  1:4709183/150‑1 (MQ=255)
                                                                                                                   cTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTaa                                    <  2:5361910/150‑1 (MQ=255)
                                                                                                                   cTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTaa                                    <  1:1299062/150‑1 (MQ=255)
                                                                                                                    tGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAg                                   <  1:3373329/150‑1 (MQ=255)
                                                                                                                       aCAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCAt                                >  1:7597781/1‑150 (MQ=255)
                                                                                                                        cAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCAtt                               <  2:8098717/150‑1 (MQ=255)
                                                                                                                                       tgCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTgcgc                <  2:3271637/150‑1 (MQ=255)
                                                                                                                                       tgCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTgcgc                <  2:2647357/150‑1 (MQ=255)
                                                                                                                                        cccgatAATGCCATCGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGcgcg               >  2:1923315/6‑150 (MQ=255)
                                                                                                                                        cccgatAATGCCATCGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGcgcg               >  1:595869/6‑150 (MQ=255)
                                                                                                                                        cccgatAATGCCATCGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGcgcg               >  1:7048371/6‑150 (MQ=255)
                                                                                                                                        cccgatAATGCCATCGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGcgcg               >  2:4152495/6‑150 (MQ=255)
                                                                                                                                         ccgatAATGCCATCGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCgg              >  1:7409046/5‑150 (MQ=255)
                                                                                                                                         ccgatAATGCCATCGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCgg              >  2:3995297/5‑150 (MQ=255)
                                                                                                                                          gggTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGc             >  1:1844347/1‑150 (MQ=255)
                                                                                                                                            gTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGccc           <  1:623635/150‑1 (MQ=255)
                                                                                                                                             tAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCGCAGCGGcccc          <  1:3590598/150‑1 (MQ=255)
                                                                                                                                                tGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGt       <  1:3218350/150‑1 (MQ=255)
                                                                                                                                                tGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGt       >  2:4449666/1‑150 (MQ=255)
                                                                                                                                                 gCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTa      >  1:1866932/1‑150 (MQ=255)
                                                                                                                                                  ccGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTaa     <  1:8340267/150‑1 (MQ=255)
                                                                                                                                                  ccGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTaa     <  1:6506616/150‑1 (MQ=255)
                                                                                                                                                  ccGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTaa     <  2:8527131/150‑1 (MQ=255)
                                                                                                                                                     aTGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAAtct  <  1:2986536/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCTGGCGGTGAGCCTGCTCGGGGATCGCTTCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCT  >  NC_002947/5467110‑5467408

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