Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,467,258 G→A 6.2% D401N (GAT→AAT) ‡ 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,2580GA 6.2% 1074.8 / 37.7 340D401N (GAT→AAT) ‡PP_4798membrane‑bound lytic murein transglycosylase
Reads supporting (aligned to +/- strand):  ref base G (124/195);  new base A (14/7);  total (138/202)
Fisher's exact test for biased strand distribution p-value = 1.97e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 343 total aligned reads displayed.

CCCTGGCGGTGAGCCTGCTCGGGGATCGCTTCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATC  >  NC_002947/5467109‑5467407
                                                                                                                                                     |                                                                                                                                                     
 ccTGGCGGTGAGCCTGCTCGGGGATCGCTTCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGa                                                                                                                                                      <  2:8754788/150‑1 (MQ=255)
 ccTGGCGGTGAGCCTGCTCGGGGATCGCTTCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGa                                                                                                                                                      <  2:3617709/150‑1 (MQ=255)
  cTGGCGGTGAGCCTGCTCGGGGATCGCTTCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGAt                                                                                                                                                     >  1:3824928/1‑150 (MQ=255)
  cTGGCGGTGAGCCTGCTCGGGGATCGCTTCAGCGGCTGGGGCTTCATTGCCGGCAGCGGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGAt                                                                                                                                                     <  2:8199131/150‑1 (MQ=255)
     gCGGTGAGCCTGCTCGGGGATCGCTTCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGc                                                                                                                                                  <  2:6951041/150‑1 (MQ=255)
           aGCCTGCTCGGGGATCGCTTCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTacc                                                                                                                                            <  2:4840008/150‑3 (MQ=255)
           aGCCTGCTCGGGGATCGCTTCAGCGGCAGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCATCGGCATTacc                                                                                                                                            <  2:141349/150‑3 (MQ=255)
             ccTGCTCGGGGATCGCTTCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCgg                                                                                                                                          <  1:3736389/150‑1 (MQ=255)
             ccTGCTCGGGGATCGCTTCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCgg                                                                                                                                          >  2:60078/1‑150 (MQ=255)
             ccTGCTCGGGGATCGCTTCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCgg                                                                                                                                          >  2:4494262/1‑150 (MQ=255)
             ccTGCTCGGGGATCGCTTCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCgg                                                                                                                                          >  2:2722131/1‑150 (MQ=255)
               tGCTCGGGGATCGCTTCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCgggg                                                                                                                                        >  2:4119517/1‑150 (MQ=255)
                gCTCGGGGATCGCTTCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGc                                                                                                                                       <  2:565319/150‑1 (MQ=255)
                gCTCGGGGATCGCTTCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGATAATGCCATCGGCATTacccg                                                                                                                                          >  1:8107358/1‑143 (MQ=255)
                      ggATCGCTTCAGCGGCTGGGGCTTCATTGCCGGCAGCTGACCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATAc                                                                                                                                 <  1:7380823/150‑1 (MQ=255)
                        aTCGCTTCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATAccc                                                                                                                               >  1:396346/1‑150 (MQ=255)
                        aTCGCTTCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATAccc                                                                                                                               >  1:1150460/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                                                                                                                    <  2:1150460/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)
                                            ttCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGcc                                                                                                           >  2:367662/1‑150 (MQ=255)
                                             tCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGccc                                                                                                          <  1:4900647/150‑1 (MQ=255)
                                              cATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCa                                                                                                         <  1:7951955/150‑1 (MQ=255)
                                                ttGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGc                                                                                                       <  1:1597715/150‑1 (MQ=255)
                                                 tGCCGGCAGCTGGCCGACAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCa                                                                                                      <  1:6943898/150‑1 (MQ=255)
                                                  gCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCaa                                                                                                     >  1:7426171/1‑150 (MQ=255)
                                                    cGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAgg                                                                                                   >  2:1959599/1‑150 (MQ=255)
                                                    cGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACCTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAgg                                                                                                   >  1:4802219/1‑150 (MQ=255)
                                                     ggCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAggg                                                                                                  >  2:5513266/1‑150 (MQ=255)
                                                     ggCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCCGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAggg                                                                                                  >  2:2499123/1‑150 (MQ=255)
                                                        aGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTg                                                                                               <  2:3890305/150‑1 (MQ=255)
                                                         gCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTgg                                                                                              <  2:5208251/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)
                                                                 aaaGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCcg                                                                                      >  2:1001414/1‑150 (MQ=255)
                                                                      ggACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGa                                                                                 <  1:4434298/150‑1 (MQ=255)
                                                                       gACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGAt                                                                                <  1:802540/150‑1 (MQ=255)
                                                                          cTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGc                                                                             >  2:6566752/1‑150 (MQ=255)
                                                                             ccGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTAc                                                                          <  1:433343/150‑1 (MQ=255)
                                                                             ccGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTAc                                                                          <  1:2241039/150‑1 (MQ=255)
                                                                               gCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTAccc                                                                        <  1:2056647/150‑1 (MQ=255)
                                                                               gCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTAccc                                                                        >  1:6885702/1‑150 (MQ=255)
                                                                               gCTCAGCCGCACCGAGCGCATGGAGTAGCAGACCTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAACACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCCAccc                                                                        <  1:3296219/150‑1 (MQ=255)
                                                                                  cAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGAc                                                                     <  1:3545874/150‑1 (MQ=255)
                                                                                     ccGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTgg                                                                                              <  2:7934202/122‑1 (MQ=255)
                                                                                     ccGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCca                                                                  >  2:4442480/1‑150 (MQ=255)
                                                                                      cGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCcac                                                                 >  2:8297998/1‑150 (MQ=255)
                                                                                        cACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACaa                                                               <  1:5839771/150‑1 (MQ=255)
                                                                                        cACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACaa                                                               <  2:5037198/150‑1 (MQ=255)
                                                                                        cACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACaa                                                               >  1:1888436/1‑150 (MQ=255)
                                                                                        cACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACaa                                                               <  1:5142611/150‑1 (MQ=255)
                                                                                             aGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAgctgc                                                          <  2:4567796/150‑1 (MQ=255)
                                                                                             aGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAgctgc                                                          <  1:8494248/150‑1 (MQ=255)
                                                                                               cgcATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCtt                                                        <  1:4481369/150‑1 (MQ=255)
                                                                                               cgcATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCtt                                                        <  1:7215844/150‑1 (MQ=255)
                                                                                                  aTGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTgag                                                     <  2:6393895/150‑1 (MQ=255)
                                                                                                   tGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTgaga                                                    <  1:4935172/150‑1 (MQ=255)
                                                                                                   tGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTgaga                                                    <  1:8299775/150‑1 (MQ=255)
                                                                                                   tGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTgaga                                                    <  2:1098690/150‑1 (MQ=255)
                                                                                                   tGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTgaga                                                    <  1:6538986/150‑1 (MQ=255)
                                                                                                   tGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTgaga                                                    <  2:1247277/150‑1 (MQ=255)
                                                                                                   tGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTgaga                                                    <  2:6852364/150‑1 (MQ=255)
                                                                                                   tGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTgaga                                                    <  2:5769335/150‑1 (MQ=255)
                                                                                                     gAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGc                                                  <  2:3417794/150‑1 (MQ=255)
                                                                                                      aGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGcc                                                 <  1:1382533/150‑1 (MQ=255)
                                                                                                       gTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCt                                                <  2:2730582/150‑1 (MQ=255)
                                                                                                           cAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTgcgc                                            >  2:2803946/1‑150 (MQ=255)
                                                                                                           cAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTgcgc                                            <  1:6644011/150‑1 (MQ=255)
                                                                                                            aGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGcc                                           >  2:1736911/1‑150 (MQ=255)
                                                                                                            aGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGcc                                           >  2:8026574/1‑150 (MQ=255)
                                                                                                              aaCTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCcag                                         <  1:4709183/150‑1 (MQ=255)
                                                                                                              aaCTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCcag                                         >  2:6298985/1‑150 (MQ=255)
                                                                                                                cTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCcagca                                       >  1:6309163/1‑150 (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)
                                                                                                                        aCAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCAt                               <  1:2781653/150‑1 (MQ=255)
                                                                                                                         cAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCAtt                              <  2:8098717/150‑1 (MQ=255)
                                                                                                                               cgGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCgggg                        <  2:2201496/150‑1 (MQ=255)
                                                                                                                                        tgCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTgcgc               <  2:3271637/150‑1 (MQ=255)
                                                                                                                                         cccgatAATGCCATCGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGcgcg              >  1:595869/6‑150 (MQ=255)
                                                                                                                                         cccgatAATGCCATCGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGcgcg              >  1:7048371/6‑150 (MQ=255)
                                                                                                                                          ccgatAATGCCATCGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCgg             >  2:3995297/5‑150 (MQ=255)
                                                                                                                                          ccgatAATGCCATCGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCgg             >  1:2745039/5‑150 (MQ=255)
                                                                                                                                           gggTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGc            >  1:1844347/1‑150 (MQ=255)
                                                                                                                                           cgatAATGCCATCGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCTCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGc            >  2:1643515/4‑150 (MQ=255)
                                                                                                                                              tAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCGCAGCGGcccc         <  1:3590598/150‑1 (MQ=255)
                                                                                                                                                aTGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCgg       >  2:2571924/1‑150 (MQ=255)
                                                                                                                                                  gCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTa     <  1:2942263/150‑1 (MQ=255)
                                                                                                                                                   ccGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTaa    <  1:6506616/150‑1 (MQ=255)
                                                                                                                                                   ccGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTaa    <  1:8340330/150‑1 (MQ=255)
                                                                                                                                                   ccGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTaa    <  2:8527131/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCCTGGCGGTGAGCCTGCTCGGGGATCGCTTCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATC  >  NC_002947/5467109‑5467407

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: