Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,762,284 A→T 7.1% 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,2840AT 7.1% 1836.0 / 87.9 562L87L (CTA→CTTyfjQCP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (343/179);  new base T (26/13);  total (369/193)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.10e-01

Only 100 of 579 total aligned reads displayed.

TCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCA  >  CP009273/2762217‑2762351
                                                                   |                                                                   
tCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTa                                                                     <  1:18373887/68‑1 (MQ=255)
tCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTa                                                                     <  1:5124325/68‑1 (MQ=255)
 ctatGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTAc                                                                    <  1:21668779/65‑1 (MQ=255)
 cTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTAc                                                                    <  1:10032809/68‑1 (MQ=255)
 cTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTAc                                                                    <  1:19073638/68‑1 (MQ=255)
   ttGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGt                                                                  <  1:22664346/68‑1 (MQ=255)
   ttGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGt                                                                  <  1:31878040/68‑1 (MQ=255)
      ccTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCgg                                                               >  1:15720812/1‑68 (MQ=255)
         gTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAg                                                            <  1:16852016/68‑1 (MQ=255)
           cAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAggg                                                          <  1:297422/68‑1 (MQ=255)
               cacGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAgggg                                                         <  1:12510921/65‑1 (MQ=255)
               cacGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGc                                                        <  1:23961986/66‑1 (MQ=255)
                acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                     <  1:21903299/68‑1 (MQ=255)
                acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                     <  1:8936032/68‑1 (MQ=255)
                acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                     >  1:29004030/1‑68 (MQ=255)
                      aCGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTAcc                                                <  1:12303324/67‑1 (MQ=255)
                      aCGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCg                                               >  1:1602480/1‑68 (MQ=255)
                       cGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCgg                                              <  1:19163438/68‑1 (MQ=255)
                       cGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCgg                                              <  1:27236874/68‑1 (MQ=255)
                       cGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCgg                                              <  1:26318369/68‑1 (MQ=255)
                         tGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTa                                            >  1:28249535/1‑68 (MQ=255)
                          gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTTCGTCGGGAGGGGCAAATTACCGGTaa                                           >  1:26082806/1‑68 (MQ=255)
                          gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                           >  1:33057792/1‑68 (MQ=255)
                          gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                           >  1:20177580/1‑68 (MQ=255)
                            cccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaaa                                          <  1:22291424/67‑1 (MQ=255)
                             ccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAAATACCGGTAAACa                                        >  1:19495155/1‑68 (MQ=255)
                               ccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAgg                                      <  1:3565559/68‑1 (MQ=255)
                               ccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAgg                                      <  1:28579574/68‑1 (MQ=255)
                                      cgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGCCAAATTACCGGTAAACAGGTTCCGGa                               >  1:15414509/1‑68 (MQ=255)
                                       gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                              <  1:15817757/68‑1 (MQ=255)
                                           gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAAtt                           <  1:12892075/67‑1 (MQ=255)
                                              tACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                         >  1:16528053/1‑66 (MQ=255)
                                                cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                      >  1:11371619/1‑67 (MQ=255)
                                                cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                      >  1:17935102/1‑67 (MQ=255)
                                                 tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTATACAGGTTCCGGAAATTATCCTGc                    <  1:5884287/68‑1 (MQ=255)
                                                 tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                      >  1:33749268/1‑66 (MQ=255)
                                                 tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                    <  1:10464644/68‑1 (MQ=255)
                                                 tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                    <  1:29947137/68‑1 (MQ=255)
                                                 tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                    <  1:31701875/68‑1 (MQ=255)
                                                 tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                    <  1:2920553/68‑1 (MQ=255)
                                                 tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                    <  1:31219621/68‑1 (MQ=255)
                                                 tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                    <  1:2949069/68‑1 (MQ=255)
                                                  aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                    <  1:17470413/64‑1 (MQ=255)
                                                  aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCt                   >  1:4805744/4‑68 (MQ=255)
                                                  aaGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCt                   >  1:8150180/1‑68 (MQ=255)
                                                   aGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCtt                  >  1:3718783/1‑68 (MQ=255)
                                                    gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                 >  1:9493130/1‑68 (MQ=255)
                                                     cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                >  1:1513669/1‑68 (MQ=255)
                                                     cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                >  1:21078742/1‑68 (MQ=255)
                                                     cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                <  1:7816833/68‑1 (MQ=255)
                                                     cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                >  1:13967956/1‑68 (MQ=255)
                                                     cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTaa                >  1:14295274/1‑68 (MQ=255)
                                                     cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                 >  1:13167253/1‑67 (MQ=255)
                                                     cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                 >  1:34372039/1‑67 (MQ=255)
                                                     cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                 >  1:1638561/1‑67 (MQ=255)
                                                     cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                 >  1:27495277/1‑67 (MQ=255)
                                                     cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                 >  1:4403429/1‑67 (MQ=255)
                                                     cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                 >  1:6982336/1‑67 (MQ=255)
                                                     cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                 >  1:25765953/1‑67 (MQ=255)
                                                      atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTaa                <  1:32700738/67‑1 (MQ=255)
                                                      atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc               <  1:19866412/68‑1 (MQ=255)
                                                      atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc               <  1:13354688/68‑1 (MQ=255)
                                                      atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc               <  1:24761432/68‑1 (MQ=255)
                                                       tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAct              <  1:27994078/68‑1 (MQ=255)
                                                       tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc               >  1:25179442/1‑67 (MQ=255)
                                                        atGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc             <  1:3863056/68‑1 (MQ=255)
                                                         tgcttcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATACTGCTTAActc             >  1:11738291/1‑67 (MQ=255)
                                                         tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc             >  1:4253471/1‑67 (MQ=255)
                                                         tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc             >  1:32403103/1‑67 (MQ=255)
                                                         tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc             >  1:7097306/1‑67 (MQ=255)
                                                         tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc             >  1:9818991/1‑67 (MQ=255)
                                                         tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc             >  1:24387524/1‑67 (MQ=255)
                                                          gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATc                        >  1:30875966/1‑55 (MQ=255)
                                                          gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc           >  1:18708903/1‑68 (MQ=255)
                                                          gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc           >  1:4191224/1‑68 (MQ=255)
                                                          gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc           >  1:20320524/1‑68 (MQ=255)
                                                          gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc           >  1:23001747/1‑68 (MQ=255)
                                                          gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc           >  1:5679865/1‑68 (MQ=255)
                                                          gctgcGCCTACGGCGGGAGGGGCAAATTACCGGTAAACAGGTGCCGGAAATTATCCTGCTTAACtctc           >  1:33377422/1‑68 (MQ=18)
                                                           ctgcGCCTACGTCGGGAGGGGCAAATTCCCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc           >  1:14581571/1‑67 (MQ=255)
                                                           ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc           >  1:12697126/1‑67 (MQ=255)
                                                           ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc           >  1:12088966/1‑67 (MQ=255)
                                                           ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc           >  1:7428873/1‑67 (MQ=255)
                                                           ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc           >  1:26739706/1‑67 (MQ=255)
                                                           ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc           >  1:6212445/1‑67 (MQ=255)
                                                           ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc           >  1:11121851/1‑67 (MQ=255)
                                                           ctgcGCCTACGTCGGGAGGGGCAAATTACCCGTAAACAGGTTCCGGAAAATATCCTGCTTCACtctc           >  1:26921056/1‑67 (MQ=255)
                                                            tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCa          <  1:30365560/67‑1 (MQ=255)
                                                            tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCAc         >  1:30775073/1‑68 (MQ=255)
                                                            tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCAc         >  1:21516314/1‑68 (MQ=255)
                                                               gcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGa       <  1:33511931/67‑1 (MQ=255)
                                                               gcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGAt      >  1:26462730/1‑68 (MQ=255)
                                                               gcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGAt      >  1:23683015/1‑68 (MQ=255)
                                                                 cTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATgg    <  1:34016511/68‑1 (MQ=255)
                                                                  tACTTCGGGAGTGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc   <  1:1162213/68‑1 (MQ=255)
                                                                  tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc   <  1:4988719/68‑1 (MQ=255)
                                                                  tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc   <  1:6568937/68‑1 (MQ=255)
                                                                   aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa  <  1:27800958/68‑1 (MQ=255)
                                                                   aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa  <  1:20708331/68‑1 (MQ=255)
                                                                   aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa  <  1:10772910/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: