Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 652,767 G→C 5.3% intergenic (‑101/+17) dsbD‑I ← / ← ctpL thiol/disulfide interchange protein 1/methyl‑accepting chemotaxis protein CtpL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947652,7670GC 5.3% 279.0 / 4.9 88intergenic (‑101/+17)dsbD‑I/ctpLthiol/disulfide interchange protein 1/methyl‑accepting chemotaxis protein CtpL
Reads supporting (aligned to +/- strand):  ref base G (45/38);  new base C (3/2);  total (48/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

GCCAGCAGTAGGAACAGCATGAAAAACAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCACCCAGTTCGTGAATCCGCTCGCTGTGGTCA  >  NC_002947/652629‑652916
                                                                                                                                          |                                                                                                                                                     
gCCAGCAGTAGGAACAGCATGAAAAACAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGC                                                                                                                                                      <  1:376000‑M1/140‑1 (MQ=255)
gCCAGCAGTAGGAACAGCATGAAAAACAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGC                                                                                                                                                      >  2:376000‑M1/1‑140 (MQ=255)
 ccAGCAGTAGGAACAGCATGAAAAACAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTC                                                                                                                                           <  2:1302195‑M1/150‑1 (MQ=255)
  cagcagTAGGAACAGCATGAAAAACAGGCGGAGCATGGCGGACTCGCGAGGTCGAACGTGCCCGGCATGATAACGGATGCGCAGGCATAATTGGCCCGAAACCAACGTACGACCGGGGGTACCCGCGCAGCGCCGCCCGAGCGGCGATCG                                                                                                                                          >  2:889208‑M1/1‑150 (MQ=255)
    gcagTAGGAACAGCATGAAAAACAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGAT                                                                                                                                        <  1:719895‑M1/150‑1 (MQ=255)
     cagTAGGAACAGCATGAAAAACAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATC                                                                                                                                       <  2:788874‑M1/150‑1 (MQ=255)
          ggAACAGCATGAAAAACAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGa                                                                                                                                  >  2:535045/1‑150 (MQ=255)
             aCAGCATGAAAAACAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGAccc                                                                                                                               >  1:217778/1‑150 (MQ=255)
                gCATGAAAAACAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGC                                                                                                                                             <  1:239089‑M1/133‑1 (MQ=255)
                gCATGAAAAACAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGC                                                                                                                                             >  2:239089‑M1/1‑133 (MQ=255)
                 cATGAAAAACAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCGGGAATAGGTGGCCGGAAACAAACGTAGGAGCGGGTTTCCCGACGAAGCGCCGCGCGGGCGGCGCTCGAGCTCAGACCCTGaa                                                                                                                           >  2:187291/1‑150 (MQ=37)
                 cATGAAAAACAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGaa                                                                                                                           >  1:1453185/1‑150 (MQ=255)
                  aTGAAAAACAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAAc                                                                                                                          >  2:1242026/1‑150 (MQ=255)
                         aCAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCa                                                                                                                   <  2:1413278/150‑1 (MQ=255)
                             gcggcgCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACAGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGc                                                                                                               >  2:22318/1‑150 (MQ=255)
                               ggcgCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGccc                                                                                                             >  1:486420/1‑150 (MQ=255)
                                 cgcATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCgg                                                                                                           <  1:1856536/150‑1 (MQ=255)
                                    aTGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCCCGCGCGGCGCTtcgcgg                                                                                                                                      >  2:1864195‑M1/1‑114 (MQ=255)
                                    aTGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCCCGCGCGGCGCTtcgcgg                                                                                                                                      <  1:1864195‑M1/120‑7 (MQ=255)
                                      ggCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAg                                                                                                      <  2:486420/150‑1 (MQ=255)
                                             tCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCg                                                                                               >  2:1908480/1‑150 (MQ=255)
                                                       gAACGGGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTca                                                                                     <  1:1154224/150‑1 (MQ=255)
                                                         aCGTGCCCGGCATGATAACGGATGCGCCGGAATAATTGGCCCGAAACAAACGTAGCAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTcagc                                                                                   >  1:31763/1‑150 (MQ=255)
                                                          cGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAACCAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTcagca                                                                                  <  2:1066450/150‑1 (MQ=255)
                                                           gTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGc                                                                                                                    >  1:623072/1‑115 (MQ=255)
                                                           gTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTcagcag                                                                                 >  1:1150754/1‑150 (MQ=255)
                                                           gTGACCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGc                                                                                                                    <  2:623072/115‑1 (MQ=255)
                                                                     tGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCa                                                                                                                                    >  2:854458/1‑89 (MQ=255)
                                                                     tGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCa                                                                                                                                    <  1:854458/89‑1 (MQ=255)
                                                                       aTAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTg                                                                     <  2:972760/150‑1 (MQ=255)
                                                                         aaCGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGtt                                                                   <  1:1242026/150‑1 (MQ=255)
                                                                           aaatcccgcaggaatAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCccaggcgctgaaactctcccggcgaacacctcccaggcgctgaaaatcccccggcgaacacctcccag                                                                 <  1:406768‑M1/144‑69 (MQ=255)
                                                                           aaatcccgcaggaatAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCccaggcgctgaaactctcccggcgaacacctcccaggcgctgaaaatcccccggcgaacacctcccag                                                                 <  1:1461929‑M1/144‑69 (MQ=255)
                                                                              tcccgcaggAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCccaggcgctgaaactctcccggcgaacacctcccaggcgctgaaaatccctgggcgaacccctcccaggcg                                                              >  2:803397‑M1/4‑79 (MQ=255)
                                                                              aTGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGc                                                              <  1:1467554/150‑1 (MQ=255)
                                                                               cccgcagGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCccaggcgctgaaactctcccggcgaacacctcccaggcgctgaaaatcccccggcgaacacctcccaggcgc                                                             <  1:551933‑M1/148‑73 (MQ=255)
                                                                                  gcAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGGAGCGCCGCCCGCGCGGCGCTt                                                                                                                                           <  2:1706420‑M1/69‑2 (MQ=255)
                                                                                  gcAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGGAGCGCCGCCCGCGCGGCGCTt                                                                                                                                           >  1:1706420‑M1/1‑68 (MQ=255)
                                                                                  gcAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCt                                                          >  1:1657371/1‑150 (MQ=255)
                                                                                    aGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCccaggcgctgaaactctcccggcgaacacctcccaggcgctgaaaatcccccggcgaacacctcccaggcgctgaaa                                                        <  1:803397‑M1/150‑78 (MQ=255)
                                                                                     ggAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCccaggcgctgaaactctcccggcgaacacctcccaggcgctgaaaatcccccggcgaacacctcccaggcgctgaaaa                                                       <  1:920874‑M1/150‑79 (MQ=255)
                                                                                          aatTGGCCCGAAACAAACGTAGGAGCGGGTCTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCccaggcgctgaaactctcccggcgaacacctcccaggcgctgaaaatcccccggcgaacacctcccaggcgctgaaaacctcc                                                  <  2:342898‑M1/150‑84 (MQ=255)
                                                                                           atTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCccaggcgctgaaactctcccggcgaacacctcccaggcgctgaaaatcccccggcgaacacctcccaggcgctgaaaacctccc                                                 <  2:1515961‑M1/150‑85 (MQ=255)
                                                                                           atTGGCCCGAAACAAACGTAGCAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCccaggcgctgaaactctcccggcgaacacctcccaggcgctgaaaatcccccggcgaacacctcccaggcgctgaaaacctccc                                                 >  1:1515961‑M1/1‑66 (MQ=255)
                                                                                            ttGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCccaggcgctgaaactctcccggcgaacacctcccaggcgctgaaaatcccccggcgaacacctcccaggcgctgaaaacctcccg                                                >  1:1858180‑M1/1‑65 (MQ=255)
                                                                                            ttGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCccaggcgctgaaactctcccggcgaacacctcccaggcgctgaaaatcccccggcgaacacctcccaggcgctgaaaacctcccg                                                >  2:1051569‑M1/1‑65 (MQ=255)
                                                                                               gCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCccaggcgctgaaactctcccggcgaacacctcccaggcgctgaaaatcccccggcgaacacctcccaggcgctgaaaacctcccgcca                                             >  2:369221‑M1/1‑62 (MQ=255)
                                                                                               gCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCAcc                                                                                                >  2:172464/1‑99 (MQ=255)
                                                                                               gCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCAcc                                                                                                <  1:172464/99‑1 (MQ=255)
                                                                                               gCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACAGACCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTg                                             <  1:535045/150‑1 (MQ=255)
                                                                                                 ccGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCccaggcgctgaaactctcccggcgaacacctcccaggcgctgaaa                                                                                        <  1:1782693‑M1/105‑46 (MQ=255)
                                                                                                 ccGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCccaggcgctgaaactctcccggcgaacacctcccaggcgctgaaa                                                                                        >  2:1782693‑M1/1‑60 (MQ=255)
                                                                                                 ccGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCa                                           >  2:1061690/1‑150 (MQ=255)
                                                                                                  cGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAg                                          >  1:1317943/1‑150 (MQ=255)
                                                                                                     atcaatcGTTGGAGCGGGTTTACCCGCGAACGGCCGCGCGGGCGGCGCTCGATATCAGACCCTGCACGCTCGCGCCGCCCGGGTCAGCTCACCGCCCATGCTCAGCAGCGGTTCACCCTGTTTACACCCCTCGCCAATGCGCTGCAGGtt                                       <  2:33086/148‑1 (MQ=37)
                                                                                                      acaaacGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCccaggcgctgaaactctcccggcgaacacctcccaggcgctgaaaatcccccggcgaacacctcccaggcgctgaaaacctcccgccaaacacct                                      >  2:1758920‑M1/1‑55 (MQ=255)
                                                                                                       caaacGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCa                                                                                                  <  2:142690/89‑1 (MQ=255)
                                                                                                       caaacGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCa                                                                                                  >  1:142690/1‑89 (MQ=255)
                                                                                                          acGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCtt                                  >  1:714268/1‑150 (MQ=255)
                                                                                                            gTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCa                                >  1:1796571/1‑150 (MQ=255)
                                                                                                               ggAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCCCCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCAccc                             >  1:121935/1‑150 (MQ=255)
                                                                                                                gAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTca                                                                                     >  2:1358460/1‑93 (MQ=255)
                                                                                                                gAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTca                                                                                     <  1:1358460/93‑1 (MQ=255)
                                                                                                                gAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACGCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCACCCa                            >  2:1398670/1‑150 (MQ=255)
                                                                                                                 aGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCccaggcgctgaaactctcccggcgaacacctcccaggcgctgaaaatcccccggcgaacacctcccaggcgctgaaaacctcccgccaaacacctacccgccaaag                           >  2:698247‑M1/1‑44 (MQ=255)
                                                                                                                  gCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTACGTTGCAAAGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCACCCCGt                          >  1:336375/1‑150 (MQ=37)
                                                                                                                    gggTTTACCCGCGAATCGCCGCGCGGGGGGCGCTCGATCTCccaggcgctgaaactctcccggcgaccacctcccaggcgctgaaaatcccccggcgaacacctcccaggcgctgaaaacctcccgccaaacacctacccgcccaagcgt                        <  2:485507‑M1/150‑110 (MQ=255)
                                                                                                                    gggTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCccaggcgctgaaactctcccggcgaacacctcccaggcgctgaaaatcccccggcgaacacctcccaggcgctgaaaacctcccgccaaacacctacccgccaaagcgc                        <  1:698247‑M1/150‑110 (MQ=255)
                                                                                                                    gggTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCACCCAGTTc                        <  2:1453185/150‑1 (MQ=255)
                                                                                                                     ggTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCa                                                                                                       >  2:1130682/1‑70 (MQ=255)
                                                                                                                     ggTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCa                                                                                                       <  1:1130682/70‑1 (MQ=255)
                                                                                                                      gTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGCCCATCAAACCCCCCACCGCCCAGATCATCTCACCGCCGAGGCTCAGCAGCTGTTCACCCCGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCACCCAGTTCGt                      <  2:217778/150‑1 (MQ=37)
                                                                                                                      gTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCACCCAGTTCGt                      >  1:914189/1‑150 (MQ=255)
                                                                                                                           cccGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCccaggcgctgaaactctcccggcgaacacctcccaggcgctgaaaatcccccggcgaacacctcccaggcgctgaaaacctcccgccaaacacctacccgccaaagcgcgaattcc                 <  1:1758920‑M1/150‑117 (MQ=255)
                                                                                                                            ccGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCccaggcgctgaaactctcccggcgaacacctcccaggcgctgaaaatcccccggcgaacacctcccaggcgctgaaaacctcccgccaaacacctacccgccaaagcgcgaattacg                >  1:969917‑M1/1‑33 (MQ=255)
                                                                                                                            ccGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCACCCAGTTCGTGAATcc                >  2:961300/1‑150 (MQ=255)
                                                                                                                             cgcgAAGCGCCGCGCGGGCGGCGCTCGATCTCccaggcgctgaaactctcccggcgaacacctcccaggcgctgaaaatcccccggcgaacacct                                                                      >  2:650302‑M1/1‑32 (MQ=255)
                                                                                                                             cgcgAAGCGCCGCGCGGGCGGCGCTCGATCTCccaggcgctgaaactctcccggcgaacacctcccaggcgctgaaaatcccccggcgaacacct                                                                      <  1:650302‑M1/95‑64 (MQ=255)
                                                                                                                             cgcgAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCACCCAGTTCGTGAATCcg               >  1:1346211/1‑150 (MQ=255)
                                                                                                                              gcgAAGCGCCGCGCGGGCGGCGCTCGATCTCccaggcgctgaaactctcccggcgaacacctcccaggcgctgaaaatcccccggcgaacacctcccaggcgctgaaaacctcccgccaaacacctacccgccaaagcgcgaattccgcc              >  2:1527949‑M1/1‑31 (MQ=255)
                                                                                                                              gcgAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCCCACCGCCCAGGc                                                                                        <  1:987290/76‑1 (MQ=255)
                                                                                                                              gcgAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCCCACCGCCCAGGc                                                                                        >  2:987290/1‑76 (MQ=255)
                                                                                                                                 aaGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCACCCAGTTCGTGAATCcgctcg           <  1:22318/150‑1 (MQ=255)
                                                                                                                                   gcgcCGCGCGGGCGGCGCTCGATCTCccaggcgctgaaactctcccggcgaacacctcccaggcgctgaaaatcccccggcgaacacctcccaggcgctgaaaacctcccgccaaacacctacccgccaaagcgcgaattccgcccattg         >  2:1219649‑M1/1‑26 (MQ=255)
                                                                                                                                     gccgcGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCACCCAGTTCGTGAATCCGCTCGCtgt       >  1:419603/1‑150 (MQ=255)
                                                                                                                                      ccgcGCGGGCGGCGCTCGATCTCccaggcgctgaaactctcccggcgaacacctcccaggcgctgaaaatcccccggcgaacacctcccaggcgctgaaaacctcccgccaaacacctacccgccaaagcgcgaattccgcccattgccc      <  2:1858180‑M1/150‑128 (MQ=255)
                                                                                                                                      ccgcGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCACCCAGTTCGTGAATCCGCTCGCtgtg      >  1:652903/1‑150 (MQ=255)
                                                                                                                                      ccgcGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCACCCAGTTCGTGAATCCGCTCGCtgtg      >  1:363238/1‑150 (MQ=255)
                                                                                                                                         cgcgGGCGGCGCTCGATCTCccaggcgctgaaactctcccggcgaacacctcccaggcgctgaaaatcccccggcgaacacctcccaggcgctgaaaccctcccgccaaacacctacccgccaaagcgcgaattccgcccattgccctat   <  1:1805436‑M1/150‑131 (MQ=255)
                                                                                                                                          gcggtcggCGCTCGATCTCccagccgatatacgcctcccgtctaatgcctcccaggagctgaaaatcccccggcgaacacctcccaggcgctgaaaacctcccgccaaacacctacccgccaaagcgcgaattccgcccattgccctata  <  1:492490‑M1/150‑132 (MQ=255)
                                                                                                                                          gcgggcggCGCTCGATCTCccaggcgctgaaactctcccggcgaacacctcccaggcgctgaaaatcccccggcgaacacctcccaggcgctgaaaacctcccgccaaacacctacccgccaaagcgcgaattccgcccattgccctata  <  1:412994‑M1/150‑132 (MQ=255)
                                                                                                                                          |                                                                                                                                                     
GCCAGCAGTAGGAACAGCATGAAAAACAGGCGGCGCATGGCGGACTCGCGAGGTGGAACGTGCCCGGCATGATAACGGATGCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCAGACCCTGAACGCTCGCACCGCCCGGTTCAGCTCACCGCCCAGGCTCAGCAGCTGTTCACCCTGTTTACGCCCCTCGCCAATGCGCTGCAGGTTGTCTTCACCCAGTTCGTGAATCCGCTCGCTGTGGTCA  >  NC_002947/652629‑652916

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