Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,203,380 G→C 5.3% A142P (GCA→CCA)  stfP → e14 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,203,3800GC 5.3% 1574.8 / 9.8 478A142P (GCA→CCA) stfPe14 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (296/153);  new base C (2/27);  total (298/180)
Fisher's exact test for biased strand distribution p-value = 1.76e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 479 total aligned reads displayed.

ATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTG  >  CP009273/1203314‑1203446
                                                                  |                                                                  
  aCCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGc                                                                   >  1:23870684/1‑66 (MQ=255)
  aCCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCa                                                                  >  1:19246346/1‑67 (MQ=255)
  aCCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCa                                                                  >  1:26612206/1‑67 (MQ=255)
  aCCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCa                                                                  >  1:33768725/1‑67 (MQ=255)
  aCCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCa                                                                  >  1:5988694/1‑67 (MQ=255)
  aCCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCa                                                                  >  1:12757249/1‑67 (MQ=255)
  aCCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCa                                                                  >  1:12989319/1‑67 (MQ=255)
  aCCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCa                                                                  >  1:13173453/1‑67 (MQ=255)
  aCCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCa                                                                  >  1:13396051/1‑67 (MQ=255)
  aCCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCa                                                                  >  1:23338566/1‑67 (MQ=255)
    cGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCa                                                                  >  1:17660831/1‑65 (MQ=255)
    cGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCa                                                                  >  1:17408200/1‑65 (MQ=255)
    cGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCa                                                                  >  1:25507837/1‑65 (MQ=255)
    cGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCa                                                                  >  1:15423055/1‑65 (MQ=255)
       ttATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGTGTGGAATGCCGCAAAAGCATCTGc                                                             <  1:27161141/67‑1 (MQ=255)
       ttATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTTGAATGCCGCAAAAGCATCTGc                                                             <  1:14996150/67‑1 (MQ=255)
        tATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGcc                                                            <  1:1950459/67‑1 (MQ=255)
        tATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGcc                                                            <  1:1844201/67‑1 (MQ=255)
        tATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGcc                                                            <  1:4048362/67‑1 (MQ=255)
        tATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGcc                                                            <  1:11144227/67‑1 (MQ=255)
          tGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATATGcc                                                            >  1:20437343/1‑65 (MQ=255)
          tGATTGGGGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCctct                                                         <  1:12359360/68‑1 (MQ=255)
           gATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTg                                                        >  1:36440220/1‑68 (MQ=255)
            aTTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTgg                                                       >  1:11530281/1‑68 (MQ=255)
            aTTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTgg                                                       >  1:8398730/1‑68 (MQ=255)
            aTTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTgg                                                       >  1:26949501/1‑68 (MQ=255)
            aTTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTgg                                                       >  1:31539527/1‑68 (MQ=255)
            aTTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTgg                                                       >  1:24135469/1‑68 (MQ=255)
            aTTGGTGGTATTGAGCAATAAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTgg                                                       <  1:32212807/68‑1 (MQ=255)
            aTTGGTGGGATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTgg                                                       >  1:33226070/1‑68 (MQ=255)
             ttGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTggg                                                      >  1:28956458/1‑68 (MQ=255)
             ttGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTggg                                                      >  1:24819734/1‑68 (MQ=255)
             ttGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTggg                                                      >  1:16574685/1‑68 (MQ=255)
                gtggtATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTgggggg                                                   >  1:18638602/1‑68 (MQ=255)
                 tggtATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTggggggg                                                  >  1:11728403/1‑68 (MQ=255)
                 tggtATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTggggggg                                                  >  1:28617371/1‑68 (MQ=255)
                  ggtATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGa                                                 >  1:6472580/1‑68 (MQ=255)
                  ggtATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGa                                                 >  1:8200996/1‑68 (MQ=255)
                  ggtATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGa                                                 >  1:12076767/1‑68 (MQ=255)
                    tATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATg                                               >  1:21551152/1‑68 (MQ=255)
                     aTTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATgg                                              >  1:7160348/1‑68 (MQ=255)
                     aTTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATg                                               <  1:33529078/67‑1 (MQ=255)
                     aTTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATg                                               <  1:30915154/67‑1 (MQ=255)
                      ttGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGa                                             <  1:31269025/68‑1 (MQ=255)
                      ttGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGa                                             <  1:34341981/68‑1 (MQ=255)
                       tGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGaa                                            >  1:28529282/1‑68 (MQ=255)
                       tGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGaa                                            >  1:28333793/1‑68 (MQ=255)
                       tGACCAATCAGTAATCATTCAGGGGGGGGGGCATGCCGCAAAACCATCTGCCTCTGGGGGGGATGGaa                                            <  1:6810162/68‑1 (MQ=255)
                          gCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAAt                                          >  1:22817837/1‑67 (MQ=255)
                          gCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAAt                                          >  1:7574355/1‑67 (MQ=255)
                          gCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAAt                                          >  1:32996722/1‑67 (MQ=255)
                          gCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAAt                                          >  1:35827317/1‑67 (MQ=255)
                          gCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAAt                                          >  1:37000585/1‑67 (MQ=255)
                          gCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAAt                                          <  1:5210633/67‑1 (MQ=255)
                          gCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGACTCTGGGGGGGATGGAAAt                                          >  1:4504731/1‑67 (MQ=255)
                          ccAAGCAGTAATCAATCAGTGGGCGTGGAATGCCGCAAAAGCTTCTGCCTCTGGGGGGGATGGAAAt                                          <  1:27590110/66‑1 (MQ=255)
                            aaTCAGTAATCATTCCGTGGGGGGGGAATGCCGCAAAAGCATCTGCCTCGGGGGGGGATGGAAATACa                                       <  1:5532839/68‑1 (MQ=255)
                            aaTCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACa                                       <  1:15265020/68‑1 (MQ=255)
                            aaTCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACa                                       <  1:14232330/68‑1 (MQ=255)
                            aaTCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACa                                       >  1:32028704/1‑68 (MQ=255)
                              tCAGTATTCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGt                                     <  1:22239821/68‑1 (MQ=255)
                              tCAGTACTCATTCAGGGGGGGTGCAATCCCGCAAAAGCAGCTGCCTCTGGGGGGGATGGAAATACAGt                                     <  1:9698873/68‑1 (MQ=255)
                              tCAGTAAGCATTCAGGGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGt                                     <  1:20938899/68‑1 (MQ=255)
                               cAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGtt                                    >  1:8518173/1‑68 (MQ=255)
                               cAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGtt                                    >  1:19847201/1‑68 (MQ=255)
                                aGTAATCATTCAGTGGGGGCTGAATGCCGCAAAAGCATCTGCCGCCGGGGGGGATGGAAATACAGt                                     <  1:10027410/66‑1 (MQ=255)
                                   aaTCATTCCGTGGGGGTGGACTGCCGCAAAAGCATCTGCCCCTGGGGGGGATGGAAATACAGTTg                                   <  1:36708875/65‑1 (MQ=255)
                                    ttCATTCCGTGGGGCTGGAATGCCGCAAAAGCAGCTGCCTCTGGGGGGGATGGAAATACAGTTGTAt                                <  1:24543858/66‑1 (MQ=255)
                                    aTCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGTAAATACAGTTGTAt                                >  1:8917967/1‑67 (MQ=255)
                                    aTCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTAt                                >  1:202523/1‑67 (MQ=255)
                                    aTCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTAt                                >  1:2183428/1‑67 (MQ=255)
                                    aTCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTAt                                >  1:29215503/1‑67 (MQ=255)
                                    aTCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTAt                                >  1:34929557/1‑67 (MQ=255)
                                      cATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTc                              >  1:17253780/1‑67 (MQ=255)
                                      cATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTc                              >  1:18807588/1‑67 (MQ=255)
                                        ttCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCg                           >  1:15873278/1‑68 (MQ=255)
                                        ttCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCg                           >  1:5325678/1‑68 (MQ=255)
                                            gTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGtt                        >  1:32611070/1‑67 (MQ=255)
                                             tGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTAtt                               >  1:26301272/1‑59 (MQ=255)
                                             tGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGc                      >  1:10731334/1‑68 (MQ=255)
                                              ggggggGCAATGCCGCAAAAGCATCTGCCTCTGGGGGGGCTGGAAATACAGTTGTATTCCCGGTTGc                      <  1:10352259/67‑1 (MQ=255)
                                              gggggTGGACTGCCGCAAAACCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGc                      <  1:9465344/67‑1 (MQ=255)
                                              gggggTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCg                     >  1:17893714/1‑68 (MQ=255)
                                              gggggTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCg                     >  1:4854466/1‑68 (MQ=255)
                                               ggggTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGt                    >  1:30490538/1‑68 (MQ=255)
                                                ggggggACTGCCGCAAACCCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGt                    <  1:29645688/67‑1 (MQ=255)
                                                ggggggAACGCCGCAAAAGCACCTGCCGCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGt                    <  1:1377658/67‑1 (MQ=255)
                                                ggggggAAAGCCGCAAAAGCCACTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGt                    <  1:10096560/67‑1 (MQ=255)
                                                gggTGGAATGCCGCACACCCATCAGCCTCGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGt                    <  1:1820872/67‑1 (MQ=255)
                                                gggTGGAATGCCGCAAACCCAGCTGCCTCTGGGGGGGATGGAAATACAGTCGTATTCCCGGTTGCGt                    <  1:26243805/67‑1 (MQ=255)
                                                 ggTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGttt                  >  1:21226841/1‑68 (MQ=255)
                                                  gTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTta                 >  1:10951367/1‑68 (MQ=255)
                                                    ggAATCCCGCAAAACCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTtaa                <  1:29969208/67‑1 (MQ=255)
                                                     caaTGCCCCAACAGCCCCTGCCCCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTtaat               <  1:35606315/66‑1 (MQ=255)
                                                      aaTGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTtaataa             >  1:16415860/1‑68 (MQ=255)
                                                       aTGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTaataat            >  1:6314170/1‑68 (MQ=255)
                                                           cGCAAACCCCTCTGCCGCTGGGGGCGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGcc         <  1:1617594/67‑1 (MQ=255)
                                                                 aGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGtt   >  1:34610992/1‑67 (MQ=255)
                                                                 aGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGtt   >  1:2575381/1‑67 (MQ=255)
                                                                  gCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTg  >  1:251426/1‑67 (MQ=255)
                                                                  |                                                                  
ATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTG  >  CP009273/1203314‑1203446

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