Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 132,467 G→C 6.5% intergenic (+41/‑111) cc → / → dsbA cytochrome c4/thiol:disulfide interchange protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947132,4670GC 6.5% 545.9 / 18.7 166intergenic (+41/‑111)cc/dsbAcytochrome c4/thiol:disulfide interchange protein
Reads supporting (aligned to +/- strand):  ref base G (65/89);  new base C (4/7);  total (69/97)
Fisher's exact test for biased strand distribution p-value = 7.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 169 total aligned reads displayed.

AGGCAACCGCACCAACGATGGCGATGCCATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAACTGATTCTCAGCGCTGCGCTGGTCGCCGC  >  NC_002947/132318‑132615
                                                                                                                                                     |                                                                                                                                                    
 ggCAACCGCACCAACGATGGCGATGCCATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCgg                                                                                                                                                     <  1:3824768/150‑1 (MQ=255)
    aaCCGCACCAACGATGGCGATGCCATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCt                                                                                                                                                  >  2:5322251/1‑150 (MQ=255)
     aCCGCACCAACGATGGCGATGCCATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTg                                                                                                                                                 <  1:2975949/150‑1 (MQ=255)
      ccGCACCAACGATGGCGATGCCATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGc                                                                                                                                                <  2:4948917/150‑1 (MQ=255)
         cACCAACGATGGCGATGCCATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCAGCCCGGCCGCCCt                                                                                                                                             <  2:579431/150‑1 (MQ=255)
          aCCAACGATGGCGATGCCATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCAGCCCGGCCGCCCtt                                                                                                                                            <  2:3020120/150‑1 (MQ=255)
           ccAACGATGGCGATGCCATGACCCTGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCttt                                                                                                                                           <  1:4713141/150‑1 (MQ=255)
           ccAACGATGGCGATGCCATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCttt                                                                                                                                           <  1:4713046/150‑1 (MQ=255)
           ccAACGATGGCGATGCCATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCttt                                                                                                                                           >  1:2812640/1‑150 (MQ=255)
            cAACGATGGCGATGCCATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCtttt                                                                                                                                          <  2:1580940/150‑1 (MQ=255)
             aaCGATGGCGATGCCATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCttttt                                                                                                                                         <  1:2700188/150‑1 (MQ=255)
             aaCGATGGCGATGCCATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCAGCCCGGCCGCCCtttca                                                                                                                                         <  1:2138350/150‑3 (MQ=255)
                gagggCGATGCCATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTtcc                                                                                                                                      <  1:660143/147‑1 (MQ=255)
                gATGGCGATGCCATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTtcc                                                                                                                                      <  1:919085/150‑1 (MQ=255)
                  tGGCGATGCCATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTtccgt                                                                                                                                    <  1:1240018/150‑1 (MQ=255)
                      gATGCCATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGc                                                                                                                                >  1:1763765/1‑150 (MQ=255)
                       aTGCCATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCa                                                                                                                               >  1:5370141/1‑150 (MQ=255)
                          ccATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGcc                                                                                                                            <  1:4981955/150‑1 (MQ=255)
                           cATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCg                                                                                                                           >  2:1860969/1‑150 (MQ=255)
                              gACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTa                                                                                                                        <  2:3750767/150‑1 (MQ=255)
                              gACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTa                                                                                                                        <  1:10351/150‑1 (MQ=255)
                                ccATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTaca                                                                                                                      >  1:1383762/1‑150 (MQ=255)
                                  aTGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACt                                                                                                                    >  1:4212598/1‑150 (MQ=255)
                                     cgcACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACa                                                                                                                 <  1:801105/150‑1 (MQ=255)
                                      gcgccATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAg                                                                                                                <  2:2455345/147‑1 (MQ=255)
                                        aCCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGaa                                                                                                              <  1:36655/150‑1 (MQ=255)
                                           aTCGCTGGCAAGCTGAGCAACCATGACATCGAAGTGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTc                                                                                                           <  2:342823/150‑1 (MQ=255)
                                             cGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGa                                                                                                         <  2:2332013/150‑1 (MQ=255)
                                             cGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGa                                                                                                         >  1:3698987/1‑150 (MQ=255)
                                              gCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGcc                                                                                                                            >  2:1645572/1‑130 (MQ=255)
                                              gCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGcc                                                                                                                            <  1:1645572/130‑1 (MQ=255)
                                              gCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAg                                                                                                        <  1:2871907/150‑1 (MQ=255)
                                               cTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCAGTCAGCAGACGTTACACTACAGAACTCGAGc                                                                                                       >  2:2366934/1‑150 (MQ=255)
                                                  gCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCt                                                                                                    >  2:514986/1‑150 (MQ=255)
                                                    aaGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTc                                                                                                  <  2:2552979/150‑1 (MQ=255)
                                                    aaGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTc                                                                                                  >  1:2144829/1‑150 (MQ=255)
                                                     aGCTGAGCAACCATGTCATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTcc                                                                                                 <  1:3448813/150‑1 (MQ=255)
                                                        tGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCAGCCCGGCCGCCCtt                                                                                                                                            <  2:5353997/104‑1 (MQ=255)
                                                           gCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCt                                                                                           <  2:4063058/150‑1 (MQ=255)
                                                             aaCCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGt                                                                                         >  1:2292892/1‑150 (MQ=255)
                                                              aCCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGtc                                                                                        >  2:5282557/1‑150 (MQ=255)
                                                                cATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTGTCGGTTAATGTTGAAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCAGTCCGCAGCCGATACACTACAGAACTCGAGCCCTTCCCGGCCTGtctc                                                                                      >  2:4396592/1‑150 (MQ=255)
                                                                 aTGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCa                                                                                     >  1:4683857/1‑150 (MQ=255)
                                                                 aTGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCa                                                                                     >  2:416213/1‑150 (MQ=255)
                                                                   gACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGt                                                                                   <  1:3734254/150‑1 (MQ=255)
                                                                    aCATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTg                                                                                  <  2:4694750/150‑1 (MQ=255)
                                                                    aCATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTg                                                                                  <  1:1371258/150‑1 (MQ=255)
                                                                       tCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAg                                                                               >  2:3947524/1‑150 (MQ=255)
                                                                         gAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGt                                                                             >  1:5345158/1‑150 (MQ=255)
                                                                           aGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTcg                                                                           <  1:997443/150‑1 (MQ=255)
                                                                                 ggCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTACTGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGcgtcgt                                                                     <  2:3273469/150‑1 (MQ=255)
                                                                                  gCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTg                                                                    >  1:1200343/1‑150 (MQ=255)
                                                                                  gCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTg                                                                    >  1:1200348/1‑150 (MQ=255)
                                                                                   ccAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTgg                                                                   >  1:4504601/1‑150 (MQ=255)
                                                                                       cTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGAc                                                               <  2:318646/150‑1 (MQ=255)
                                                                                            tCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGt                                                                             <  2:955450/131‑1 (MQ=255)
                                                                                            tCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATg                                                          >  1:2635278/1‑150 (MQ=255)
                                                                                            tCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATg                                                          >  1:2635195/1‑150 (MQ=255)
                                                                                            tCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCCACGACGCGAcctgca                                                          <  2:2812640/150‑5 (MQ=255)
                                                                                             ccAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGa                                                         >  1:1196296/1‑150 (MQ=255)
                                                                                             ccAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCGTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGa                                                         <  2:2799011/150‑1 (MQ=255)
                                                                                              cAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGAc                                                        <  2:1763765/150‑1 (MQ=255)
                                                                                                    cTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCa                                                  <  1:1005190/150‑1 (MQ=255)
                                                                                                    cTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCa                                                  <  2:2182623/150‑1 (MQ=255)
                                                                                                     tGCATGAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCGTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAg                                                 <  2:3680076/150‑1 (MQ=255)
                                                                                                       cATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGa                                               <  2:1122662/150‑1 (MQ=255)
                                                                                                       cATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGa                                               <  1:3330129/150‑1 (MQ=255)
                                                                                                        aTTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAg                                              >  1:3006970/1‑150 (MQ=255)
                                                                                                         ttAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGt                                             >  1:15834/1‑150 (MQ=255)
                                                                                                           aaCTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTaa                                           >  2:4171493/1‑150 (MQ=255)
                                                                                                            aCTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTaaa                                          <  2:4212598/150‑1 (MQ=255)
                                                                                                              tttCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGcgtcgt                                                                     <  2:3257827/121‑1 (MQ=255)
                                                                                                              tttCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGc                                        <  2:654328/150‑1 (MQ=255)
                                                                                                              tttCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGc                                        >  1:2446356/1‑150 (MQ=255)
                                                                                                               ttCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCa                                       <  1:4666116/150‑1 (MQ=255)
                                                                                                                  ggTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGc                                    <  2:1807242/150‑1 (MQ=255)
                                                                                                                  ggTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGc                                    <  1:5279277/150‑1 (MQ=255)
                                                                                                                  ggTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGc                                    <  1:632066/150‑1 (MQ=255)
                                                                                                                  ggTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGc                                    <  2:841367/150‑1 (MQ=255)
                                                                                                                   gTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCg                                   >  1:5172284/1‑150 (MQ=255)
                                                                                                                       aTGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTaaa                               >  2:3365723/1‑150 (MQ=255)
                                                                                                                        tgttgtAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAAc                              <  2:3376520/150‑1 (MQ=255)
                                                                                                                        tgttgtAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAAc                              <  2:1925473/150‑1 (MQ=255)
                                                                                                                           tgtAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAACTGa                           <  2:3611296/150‑1 (MQ=255)
                                                                                                                             tAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAACTGAtt                         <  2:3796316/150‑1 (MQ=255)
                                                                                                                             tAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAACTGAtt                         <  1:4452964/150‑1 (MQ=255)
                                                                                                                                gTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAACTGATtctc                      >  2:5175042/1‑150 (MQ=255)
                                                                                                                                 tCAATGTTGAATGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCTGTGCAGGTCGCGTCGTGGCGACCCATGAGTGCTCAGGAGTAATGCATGCGTAATCTGATtgtct                     >  2:3995360/1‑149 (MQ=255)
                                                                                                                                 tCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAACTGATTCTCa                     <  2:3122115/150‑1 (MQ=255)
                                                                                                                                     tgatgaAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCCTGCGTAAACTGATTCTCAgcgc                 >  2:3651240/1‑150 (MQ=255)
                                                                                                                                         aaaaGGGCAGCCCGGCCGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGc                                    >  1:1536479/2‑127 (MQ=255)
                                                                                                                                          aaaGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAACTGATTCTCAgcgctgcgc            >  1:2951921/1‑150 (MQ=255)
                                                                                                                                           aaGGGCAGCCCGGCCGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAAATCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAACTGATTCTCAgcgctgcgct           >  1:5062393/1‑150 (MQ=255)
                                                                                                                                            aGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAACTGATTCTCAGcgctgcgctg          <  2:3338664/150‑1 (MQ=255)
                                                                                                                                            aGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAACTGATTCTCAGcgctgcgctg          <  2:1383762/150‑1 (MQ=255)
                                                                                                                                             gggCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAACTGATTCTCAGCGCTGCGCTgg         <  1:766937/150‑1 (MQ=255)
                                                                                                                                                 ggcCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAACTGATTCTCAGCGCTGCGCTGGTcgc     >  1:486200/1‑150 (MQ=255)
                                                                                                                                                  gcCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAACTGATTCTCAGCGCTGCGCTGGTcgcc    <  1:52944/150‑1 (MQ=255)
                                                                                                                                                  gcCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAACTGATTCTCAGCGCTGCGCTGGTcgcc    <  2:2311381/150‑1 (MQ=255)
                                                                                                                                                   ccGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAACTGATTCTCAGCGCTGCGCTGGTcgccg   >  1:7299/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AGGCAACCGCACCAACGATGGCGATGCCATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAACTGATTCTCAGCGCTGCGCTGGTCGCCGC  >  NC_002947/132318‑132615

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: