Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,762,304 C→A 8.3% T94N (ACC→AAC)  yfjQ → CP4‑57 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,762,3040CA 8.3% 2242.1 / 134.2 642T94N (ACC→AAC) yfjQCP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (405/182);  new base A (38/16);  total (444/198)
Fisher's exact test for biased strand distribution p-value = 8.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

Only 100 of 648 total aligned reads displayed.

ACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTG  >  CP009273/2762239‑2762371
                                                                 |                                                                   
aCGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTAcc                                                                    <  1:9134562/67‑1 (MQ=255)
aCGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTAcc                                                                    <  1:38211686/67‑1 (MQ=255)
 cGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCgg                                                                  <  1:17680986/68‑1 (MQ=255)
  gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                                                 <  1:30255016/68‑1 (MQ=255)
  gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                                                 >  1:34115584/1‑68 (MQ=255)
   tGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTa                                                                >  1:15468172/1‑68 (MQ=255)
    gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                                               >  1:13820386/1‑68 (MQ=255)
     aCCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                                               >  1:26642169/1‑67 (MQ=255)
     aCCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                                               >  1:38454370/1‑67 (MQ=255)
      cccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaaa                                                              <  1:38708419/67‑1 (MQ=255)
      cccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAAc                                                             >  1:29534834/1‑68 (MQ=255)
      cccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAAc                                                             >  1:23915625/1‑68 (MQ=255)
            tcgtcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGtt                                                        >  1:28936068/1‑67 (MQ=255)
                cgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGa                                                   >  1:32967020/1‑68 (MQ=255)
                cgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGa                                                   >  1:34998042/1‑68 (MQ=255)
                 gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                                  <  1:31130885/68‑1 (MQ=255)
                 gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                                  <  1:16665553/68‑1 (MQ=255)
                 gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                                  <  1:34714029/68‑1 (MQ=255)
                 gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                                  <  1:20368943/68‑1 (MQ=255)
                  tGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaaa                                                 >  1:22638424/1‑68 (MQ=255)
                     gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTa                                              >  1:13188643/1‑68 (MQ=255)
                         aCTAAGCTTATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                                          <  1:424228/68‑1 (MQ=255)
                         aCTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                                             >  1:5935841/1‑65 (MQ=255)
                           tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                                        <  1:31866954/68‑1 (MQ=255)
                           tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                                        <  1:11504983/68‑1 (MQ=255)
                            aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                                        <  1:10872406/64‑1 (MQ=255)
                            aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                                        <  1:12469876/64‑1 (MQ=255)
                            aaGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCAGCt                                       >  1:35414272/1‑68 (MQ=255)
                              gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                     >  1:1279666/1‑68 (MQ=255)
                              gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                     >  1:2795211/1‑68 (MQ=255)
                              gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                     >  1:25459953/1‑68 (MQ=255)
                              gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                     >  1:18352508/1‑68 (MQ=255)
                              acATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTa                                     >  1:26324433/2‑68 (MQ=255)
                               cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                    >  1:26048485/1‑68 (MQ=255)
                               cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                    >  1:4017360/1‑68 (MQ=255)
                               cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                    >  1:27032757/1‑68 (MQ=255)
                               cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                    >  1:2468997/1‑68 (MQ=255)
                               cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                    >  1:12154951/1‑68 (MQ=255)
                               cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                    >  1:21753957/1‑68 (MQ=255)
                               cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                    >  1:21290694/1‑68 (MQ=255)
                               cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                     >  1:27237331/1‑67 (MQ=255)
                               cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                     >  1:19393385/1‑67 (MQ=255)
                               cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                     >  1:3892408/1‑67 (MQ=255)
                               cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                     >  1:28409686/1‑67 (MQ=255)
                               cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                     >  1:27741480/1‑67 (MQ=255)
                               cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                     >  1:11758459/1‑67 (MQ=255)
                               cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                     >  1:17902881/1‑67 (MQ=255)
                               cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                     >  1:18100352/1‑67 (MQ=255)
                               cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                     >  1:15092909/1‑67 (MQ=255)
                                atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                     >  1:5690446/1‑66 (MQ=255)
                                atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc                                   <  1:4748301/68‑1 (MQ=255)
                                   tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActct                                <  1:8536071/68‑1 (MQ=255)
                                   tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                 >  1:25271490/1‑67 (MQ=255)
                                   tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                 >  1:36767365/1‑67 (MQ=255)
                                   tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                 >  1:18696107/1‑67 (MQ=255)
                                   tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                 >  1:21579813/1‑67 (MQ=255)
                                   tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                 >  1:17304808/1‑67 (MQ=255)
                                   tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                 >  1:33856206/1‑67 (MQ=255)
                                   tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                 >  1:33393002/1‑67 (MQ=255)
                                   tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                 >  1:17000313/1‑67 (MQ=255)
                                   tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                 >  1:30647426/1‑67 (MQ=255)
                                    gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                               >  1:9627290/1‑68 (MQ=255)
                                    gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                               >  1:3744360/1‑68 (MQ=255)
                                    gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                               >  1:14801249/1‑68 (MQ=255)
                                    gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                               >  1:3213419/1‑68 (MQ=255)
                                    gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                               >  1:20033362/1‑68 (MQ=255)
                                    gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                               >  1:10161998/1‑68 (MQ=255)
                                    gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAAAAGGTTCCGGAAATTATCCTGCTTAACtctc                               >  1:837377/1‑68 (MQ=255)
                                     ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAATCAGGTTCCGGAAATTATCCTGCTGAACtctc                               >  1:16227905/1‑67 (MQ=255)
                                     ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTACCtctc                               >  1:35216674/1‑67 (MQ=255)
                                     ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                               >  1:3261895/1‑67 (MQ=255)
                                     ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                               >  1:14309690/1‑67 (MQ=255)
                                     ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                               >  1:21163820/1‑67 (MQ=255)
                                     ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                               >  1:7774288/1‑67 (MQ=255)
                                     ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                               >  1:34342693/1‑67 (MQ=255)
                                     ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                               >  1:24542710/1‑67 (MQ=255)
                                     ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                               >  1:29298502/1‑67 (MQ=255)
                                     ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAAAAGGTTCCGGAAATTATCCTGCTTAACtctc                               >  1:5031035/1‑67 (MQ=255)
                                       gcgcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACg                            <  1:37937187/68‑1 (MQ=255)
                                       gcgcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACg                            <  1:19101686/68‑1 (MQ=255)
                                       gcgcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACg                            <  1:32377299/68‑1 (MQ=255)
                                          ccTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATg                         >  1:23595783/1‑68 (MQ=255)
                                           cTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATg                         >  1:15695252/1‑67 (MQ=255)
                                            tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc                       <  1:35885234/68‑1 (MQ=255)
                                            tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc                       <  1:11171146/68‑1 (MQ=255)
                                            tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc                       <  1:6975211/68‑1 (MQ=255)
                                            tACGTCGGGAGGGGCAAATTAACGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc                       <  1:6493028/68‑1 (MQ=255)
                                              cGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGcc                      >  1:10014583/1‑66 (MQ=255)
                                                tCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCa                   >  1:24262915/1‑68 (MQ=255)
                                                 cGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAg                  >  1:36498219/1‑68 (MQ=255)
                                                 cGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAg                  >  1:36316041/1‑68 (MQ=255)
                                                  gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                 <  1:13061926/68‑1 (MQ=255)
                                                  gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                 <  1:19627252/68‑1 (MQ=255)
                                                  gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                 <  1:7428150/68‑1 (MQ=255)
                                                        ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt             <  1:8228079/66‑1 (MQ=255)
                                                           aaaTTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAg        <  1:8691163/68‑1 (MQ=255)
                                                             aTTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAt      <  1:15903884/68‑1 (MQ=255)
                                                             aTTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAt      <  1:20820617/68‑1 (MQ=255)
                                                             aTTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAt      <  1:22301445/68‑1 (MQ=255)
                                                                 ccGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAtgctg  <  1:29927868/68‑1 (MQ=255)
                                                                 |                                                                   
ACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTG  >  CP009273/2762239‑2762371

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: