Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,762,284 A→T 10.9% L87L (CTA→CTT yfjQ → CP4‑57 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,762,2840AT10.9% 1674.5 / 160.2 569L87L (CTA→CTTyfjQCP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (341/164);  new base T (55/7);  total (396/173)
Fisher's exact test for biased strand distribution p-value = 3.63e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.75e-01

Only 100 of 585 total aligned reads displayed.

TCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCA  >  CP009273/2762217‑2762351
                                                                   |                                                                   
  tttGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACg                                                                   <  1:30248620/68‑1 (MQ=255)
  tttGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACg                                                                   <  1:15376274/68‑1 (MQ=255)
   ttGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGt                                                                  <  1:27686083/68‑1 (MQ=255)
     gCCTGTCAGACACGGGTACCTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCg                                                                >  1:5757109/1‑68 (MQ=255)
      ccTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCgg                                                               >  1:22359485/1‑68 (MQ=255)
        tGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGa                                                             <  1:2660911/68‑1 (MQ=255)
        tGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGa                                                             >  1:8810666/1‑68 (MQ=255)
          tCAGACCCGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAgg                                                           >  1:4911190/1‑68 (MQ=255)
           cAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAggg                                                          <  1:16760626/68‑1 (MQ=255)
           cAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAggg                                                          <  1:10951462/68‑1 (MQ=255)
           cAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAggg                                                          <  1:1337079/68‑1 (MQ=255)
               cacGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGc                                                        <  1:1305292/66‑1 (MQ=255)
                acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                     <  1:26956558/68‑1 (MQ=255)
                acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                     <  1:4405646/68‑1 (MQ=255)
                      aCGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTAcc                                                <  1:28203190/67‑1 (MQ=255)
                      aCGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTAcc                                                <  1:28335381/67‑1 (MQ=255)
                      aCGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTAcc                                                <  1:6528016/67‑1 (MQ=255)
                      aCGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCg                                               <  1:18576514/68‑1 (MQ=255)
                        gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                             >  1:17880873/1‑68 (MQ=255)
                        gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                             >  1:17447716/1‑68 (MQ=255)
                        gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                             >  1:11450201/1‑68 (MQ=255)
                        gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                             >  1:4045404/1‑68 (MQ=255)
                         tGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTa                                            >  1:25265406/1‑68 (MQ=255)
                          gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                           >  1:20262975/1‑68 (MQ=255)
                          gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                           >  1:29638909/1‑68 (MQ=255)
                            cccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaatc                                         >  1:34873354/1‑66 (MQ=255)
                            cccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaaa                                          <  1:3242139/67‑1 (MQ=255)
                            cccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAAc                                         >  1:5439384/1‑68 (MQ=255)
                            cccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAAc                                         >  1:12372924/1‑68 (MQ=255)
                             ccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACa                                        >  1:259095/1‑68 (MQ=255)
                                      cgtGAGCATACTAAGCATATGCTGCGCCTACGGCGGGAGGGGCAAATTACCGGTAAACATCTACCGGa                               >  1:26709626/1‑68 (MQ=18)
                                       gtgtGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                              <  1:15178376/68‑1 (MQ=255)
                                       gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                              <  1:28575796/68‑1 (MQ=255)
                                       gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                              <  1:26396327/68‑1 (MQ=255)
                                       gtGAGCATACTAAGCATATCCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                              <  1:34182949/68‑1 (MQ=255)
                                               aCTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCTGTAAACAGGTTCCGGAAATTATCCt                      <  1:30015463/68‑1 (MQ=255)
                                                cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTg                     <  1:937904/68‑1 (MQ=255)
                                                 tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                    <  1:27466868/68‑1 (MQ=255)
                                                 tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                    <  1:31424205/68‑1 (MQ=255)
                                                 tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                    <  1:14772896/68‑1 (MQ=255)
                                                 tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                    <  1:16098467/68‑1 (MQ=255)
                                                  aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                    <  1:2183009/64‑1 (MQ=255)
                                                  aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                    <  1:20907736/64‑1 (MQ=255)
                                                  aaGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCt                   >  1:1971406/1‑68 (MQ=255)
                                                    gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                 >  1:10034764/1‑68 (MQ=255)
                                                    gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                 >  1:7819397/1‑68 (MQ=255)
                                                    acATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTa                 >  1:71282/2‑68 (MQ=255)
                                                    acATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTa                 >  1:27887648/2‑68 (MQ=255)
                                                     cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                >  1:25073811/1‑68 (MQ=255)
                                                     cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                >  1:11915352/1‑68 (MQ=255)
                                                     cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                >  1:14179275/1‑68 (MQ=255)
                                                     cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                >  1:30905236/1‑68 (MQ=255)
                                                     cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                >  1:17193476/1‑68 (MQ=255)
                                                     cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                >  1:29343515/1‑68 (MQ=255)
                                                     cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                >  1:19175032/1‑68 (MQ=255)
                                                     cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                >  1:20514560/1‑68 (MQ=255)
                                                     cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                >  1:8126553/1‑68 (MQ=255)
                                                     cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                >  1:26820894/1‑68 (MQ=255)
                                                     cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGATCCTGAAATTATCCTGCTTaa                >  1:15632952/1‑68 (MQ=255)
                                                     cATATGCTGCGCCTAGGTCGGGAGGGGCAATTTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                 >  1:11712770/1‑67 (MQ=255)
                                                     cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTaa                >  1:33949800/1‑68 (MQ=255)
                                                     cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                 >  1:12836046/1‑67 (MQ=255)
                                                     cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                 >  1:22740350/1‑67 (MQ=255)
                                                     cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                 >  1:3351933/1‑67 (MQ=255)
                                                     cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                 >  1:16921040/1‑67 (MQ=255)
                                                     cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                 >  1:6737259/1‑67 (MQ=255)
                                                     cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTACGGAAATTATCCTGCTTa                 >  1:32257214/1‑67 (MQ=255)
                                                      atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                 >  1:19898946/1‑66 (MQ=255)
                                                      atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                 >  1:29114777/1‑66 (MQ=255)
                                                      atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc               <  1:28884586/68‑1 (MQ=255)
                                                      atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc               <  1:23284781/68‑1 (MQ=255)
                                                       tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAct              <  1:11133528/68‑1 (MQ=255)
                                                       tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc               >  1:24728653/1‑67 (MQ=255)
                                                        atGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc             <  1:13700081/68‑1 (MQ=255)
                                                         tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc             >  1:2120389/1‑67 (MQ=255)
                                                         tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc             >  1:1655718/1‑67 (MQ=255)
                                                         tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc             >  1:31237878/1‑67 (MQ=255)
                                                          gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActct            <  1:24024384/67‑1 (MQ=255)
                                                          gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActcac           >  1:17707702/1‑66 (MQ=255)
                                                          gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc           >  1:32609843/1‑68 (MQ=255)
                                                          gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc           >  1:30667107/1‑68 (MQ=255)
                                                          gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc           >  1:26284667/1‑68 (MQ=255)
                                                          gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc           >  1:8940982/1‑68 (MQ=255)
                                                          gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc           >  1:34571840/1‑68 (MQ=255)
                                                          gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc           >  1:18457913/1‑68 (MQ=255)
                                                          gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc           >  1:2213418/1‑68 (MQ=255)
                                                           ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc           >  1:9800992/1‑67 (MQ=255)
                                                           ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc           >  1:7503084/1‑67 (MQ=255)
                                                           ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc           >  1:5291392/1‑67 (MQ=255)
                                                           ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc           >  1:25591834/1‑67 (MQ=255)
                                                           ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc           >  1:21441790/1‑67 (MQ=255)
                                                           ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc           >  1:1812363/1‑67 (MQ=255)
                                                            tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCa          <  1:14391303/67‑1 (MQ=255)
                                                             gcgcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACg        <  1:6094453/68‑1 (MQ=255)
                                                               gcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGAt      >  1:8420731/1‑68 (MQ=255)
                                                                 cTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATg     >  1:33335086/1‑67 (MQ=255)
                                                                   aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa  <  1:3185012/68‑1 (MQ=255)
                                                                   aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa  <  1:23057966/68‑1 (MQ=255)
                                                                   aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa  <  1:10777752/68‑1 (MQ=255)
                                                                   aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa  <  1:10311697/68‑1 (MQ=255)
                                                                   |                                                                   
TCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCA  >  CP009273/2762217‑2762351

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