Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,467,260 T→C 5.3% D401D (GAT→GAC) ‡ 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,2600TC 5.3% 858.5 / 21.9 264D401D (GAT→GAC) ‡PP_4798membrane‑bound lytic murein transglycosylase
Reads supporting (aligned to +/- strand):  ref base T (106/144);  new base C (10/4);  total (116/148)
Fisher's exact test for biased strand distribution p-value = 5.01e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.34e-01

Only 100 of 265 total aligned reads displayed.

CTGGCGGTGAGCCTGCTCGGGGATCGCTTCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCT  >  NC_002947/5467111‑5467408
                                                                                                                                                     |                                                                                                                                                    
cTGGCGGTGAGCCTGCTCGGGGATCGCTTCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGAt                                                                                                                                                      <  1:4105670/150‑1 (MQ=255)
 tGGCGGTGAGCCTGCTCGGGGATCGCTTCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATg                                                                                                                                                     <  2:2852986/150‑1 (MQ=255)
 tGGCGGTGAGCCTGCTCGGGGATCGCTTCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATg                                                                                                                                                     <  2:3044198/150‑1 (MQ=255)
 tGGCGGTGAGCCTGCTCGGGGATCGCTTCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATg                                                                                                                                                     >  1:21077/1‑150 (MQ=255)
     ggTGAGCCTGCTCGGGGATCGCTTCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCAt                                                                                                                                                 >  2:2177661/1‑150 (MQ=255)
     ggTGAGCCTGCTCGGGGATCGCTTCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGAACTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCAt                                                                                                                                                 >  2:2240053/1‑150 (MQ=255)
            cTGCTCGGGGATCGCTTCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCggg                                                                                                                                          <  1:2160403/150‑1 (MQ=255)
                  ggggATCGCTTCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAAt                                                                                                                                    >  2:3396391/1‑150 (MQ=255)
                     gATCGCTTCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATAccc                                                                                                                                <  2:3088684/150‑1 (MQ=255)
                      aTCGCTTCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATAccc                                                                                                                                >  2:1826549/1‑150 (MQ=255)
                         gCTTCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCa                                                                                                                             >  1:1283619/1‑150 (MQ=255)
                          cgtCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCaa                                                                                                                            <  2:5099399/148‑1 (MQ=255)
                          cTTCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCaa                                                                                                                            <  1:2988534/150‑1 (MQ=255)
                           ttCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAg                                                                                                                           >  2:3642447/1‑150 (MQ=255)
                           ttCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAg                                                                                                                           >  1:2484387/1‑150 (MQ=255)
                           ttCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAg                                                                                                                           >  1:905308/1‑150 (MQ=255)
                             cAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCATCGGCATTac                                                                                                                                              >  2:2117431/1‑128 (MQ=255)
                               gCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCaa                                                                                                                       <  1:747397/150‑1 (MQ=255)
                                      gggCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCaa                                                                                                                >  1:4187440/1‑150 (MQ=255)
                                         cTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGc                                                                                                             <  1:2625872/150‑1 (MQ=255)
                                          ttCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGcc                                                                                                            >  2:453385/1‑150 (MQ=255)
                                          ttCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGcc                                                                                                            >  1:518428/1‑150 (MQ=255)
                                           tCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGccc                                                                                                           <  1:1555344/150‑1 (MQ=255)
                                           tCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGccc                                                                                                           <  1:2808724/150‑1 (MQ=255)
                                            cATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCa                                                                                                          <  2:5550756/150‑1 (MQ=255)
                                            cATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCa                                                                                                          <  2:1151756/150‑1 (MQ=255)
                                             aTTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAg                                                                                                         >  1:3234816/1‑150 (MQ=255)
                                             aTTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAg                                                                                                         >  1:627210/1‑150 (MQ=255)
                                              ttGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGc                                                                                                        >  2:2602852/1‑150 (MQ=255)
                                                    gaaGCTGGCCGAAAGAGGACCTGCAGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGACGCGGGCATGATGCGCGTAATGCCGATGGCATTATCGGGGCCAATACCCGAAAGGCAATCCGCAATGCCCAGCAAgcgc                                                                                                  >  2:3966032/3‑150 (MQ=255)
                                                       gCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTgg                                                                                               <  2:1737128/150‑1 (MQ=255)
                                                         tGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGa                                                                                             >  2:3081671/1‑150 (MQ=255)
                                                         tGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGa                                                                                             >  2:1806228/1‑150 (MQ=255)
                                                           gCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATg                                                                                           <  1:4269633/150‑1 (MQ=255)
                                                           gCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATg                                                                                           <  2:2466001/150‑1 (MQ=255)
                                                           gCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATg                                                                                           <  2:1027993/150‑1 (MQ=255)
                                                              gAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAGGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGcc                                                                                        >  1:2842914/1‑150 (MQ=255)
                                                              gAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGcc                                                                                        >  2:2839785/1‑150 (MQ=255)
                                                                aaGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCcgg                                                                                      <  1:3396391/150‑1 (MQ=255)
                                                                  gagGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCcggcg                                                                                    <  1:1501077/150‑1 (MQ=255)
                                                                   agGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCaa                                                                                                                <  2:919039/121‑1 (MQ=255)
                                                                     gACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGAt                                                                                 <  1:5624600/150‑1 (MQ=255)
                                                                     gACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGAt                                                                                 <  2:4152170/150‑1 (MQ=255)
                                                                     gACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGAt                                                                                 <  1:4614607/150‑1 (MQ=255)
                                                                     gACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGAt                                                                                 >  1:5076908/1‑150 (MQ=255)
                                                                     gACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGAt                                                                                 <  1:1105351/150‑1 (MQ=255)
                                                                     gACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGAt                                                                                 <  2:4685700/150‑1 (MQ=255)
                                                                      aCCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATg                                                                                <  2:1573800/150‑1 (MQ=255)
                                                                       ccTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATgg                                                                               <  2:2756681/150‑1 (MQ=255)
                                                                       ccTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATgg                                                                               <  2:256643/150‑1 (MQ=255)
                                                                         tGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCt                                                                             <  2:3710998/150‑1 (MQ=255)
                                                                          gCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTACTGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTa                                                                            <  2:5286139/150‑1 (MQ=255)
                                                                          gCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTa                                                                            >  1:3065729/1‑150 (MQ=255)
                                                                           ccGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTAc                                                                           >  2:3305517/1‑150 (MQ=255)
                                                                                cAGCCGCACCGAGTGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGAc                                                                      <  2:1793609/150‑1 (MQ=255)
                                                                                  gCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGAccc                                                                    >  2:48679/1‑150 (MQ=255)
                                                                                      cACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACaa                                                                <  1:703384/150‑1 (MQ=255)
                                                                                         cGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAgct                                                             <  1:541625/150‑1 (MQ=255)
                                                                                         cGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAgct                                                             <  1:3354700/150‑1 (MQ=255)
                                                                                             cgcATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCtt                                                         <  2:3682090/150‑1 (MQ=255)
                                                                                             cgcATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCtt                                                         <  1:4314111/150‑1 (MQ=255)
                                                                                             cgcATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCtt                                                         <  1:881712/150‑1 (MQ=255)
                                                                                               cATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGAAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTga                                                       <  1:395180/150‑1 (MQ=255)
                                                                                                 tGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTgaga                                                     <  2:1962710/150‑1 (MQ=255)
                                                                                                 tGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTgaga                                                     <  1:2005504/150‑1 (MQ=255)
                                                                                                  ggAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGagag                                                    >  1:2117013/1‑150 (MQ=255)
                                                                                                     gTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTgaga                                                     >  2:1313941/1‑146 (MQ=255)
                                                                                                       tGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTgc                                               <  1:5348318/150‑1 (MQ=255)
                                                                                                        gCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTgcg                                              >  2:605350/1‑150 (MQ=255)
                                                                                                         cAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTgcgc                                             >  2:7682/1‑150 (MQ=255)
                                                                                                            aaCTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCcag                                          >  1:4419507/1‑150 (MQ=255)
                                                                                                              cTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGGGCCTGCGCcagca                                        <  2:4476842/150‑1 (MQ=255)
                                                                                                              cTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCcagca                                        <  1:3621600/150‑1 (MQ=255)
                                                                                                              cTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCcagca                                        <  1:1360200/150‑1 (MQ=255)
                                                                                                               ttGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCcagcag                                       >  2:4930343/1‑150 (MQ=255)
                                                                                                                  cTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTaa                                    <  2:3427520/150‑1 (MQ=255)
                                                                                                                  cTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAAGACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTaa                                    <  2:4702579/150‑1 (MQ=255)
                                                                                                                   tGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAg                                   >  2:5017530/1‑150 (MQ=255)
                                                                                                                   tGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAg                                   >  1:5590448/1‑150 (MQ=255)
                                                                                                                   tGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAg                                   >  2:873806/1‑150 (MQ=255)
                                                                                                                    gAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGc                                  >  1:4557123/1‑150 (MQ=255)
                                                                                                                     aaCAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCa                                 <  1:2378384/150‑1 (MQ=255)
                                                                                                                        aGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTc                              <  2:4331938/150‑1 (MQ=255)
                                                                                                                         gCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTcc                             <  1:2358139/150‑1 (MQ=255)
                                                                                                                           cgcgGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCgg                           >  1:4792068/1‑150 (MQ=255)
                                                                                                                            gcgGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCggg                          <  1:2602852/150‑1 (MQ=255)
                                                                                                                            gcgGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCggg                          <  2:1202896/150‑1 (MQ=255)
                                                                                                                                gCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCg                      <  2:5403593/150‑1 (MQ=255)
                                                                                                                                    gatgCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTgc                  <  1:1826549/150‑1 (MQ=255)
                                                                                                                                       cccgatAATGCCATCGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGcgcg               >  1:4954628/6‑150 (MQ=255)
                                                                                                                                        ccgatAATGCCATCGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCgg              >  2:2335708/5‑150 (MQ=255)
                                                                                                                                        ccgatAATGCCATCGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCgg              >  1:4393680/5‑150 (MQ=255)
                                                                                                                                        cGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGAt                                                                                 >  2:4321884/1‑83 (MQ=255)
                                                                                                                                        cGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCgg              <  2:4407767/150‑1 (MQ=255)
                                                                                                                                           gTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGccc           <  2:1445191/150‑1 (MQ=255)
                                                                                                                                              aTGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCgg        <  1:3830947/150‑1 (MQ=255)
                                                                                                                                                 ccGATGGCATTATCGGGGCCAAGACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCGCAGCGGCCCCggaat     <  1:5015401/150‑4 (MQ=255)
                                                                                                                                                  cGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTatg    <  1:4927862/150‑3 (MQ=255)
                                                                                                                                                  cGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAAt    >  2:413752/1‑150 (MQ=255)
                                                                                                                                                  cGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCAGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTaac    <  1:1776297/150‑2 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CTGGCGGTGAGCCTGCTCGGGGATCGCTTCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGTTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCT  >  NC_002947/5467111‑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: