Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,657,965 G→A 12.9% A34T (GCC→ACC)  ynfH → oxidoreductase, membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,657,9650GA12.9% 80.2 / 5.4 31A34T (GCC→ACC) ynfHoxidoreductase, membrane subunit
Reads supporting (aligned to +/- strand):  ref base G (15/12);  new base A (2/2);  total (17/14)
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.87e-01

GAAGAGCTGCGCCAGAAGCACGGCACTCTTGCTGCCGTCGCGCCACTGCCGCGCGCGCATTTCACAAAACCCAATATCGTTATCAAACCTAACGCCAACAGCCGCCCGACCGGCGATACCACTGGTTATCTGGCTAATCCGGAGGAAGTGTAAAATGGGAAATGGATGGCATGAGTGGCCACTGGTTATCTTTACCGTATTGGGTCAATGCGTAGTAGGCGCGCTGATTGTTAGCGGAATTGGCTGGTTCGCCGCCAAAAATGACGCAGACCGTCAGCGTATTGTCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGT  >  CP009273/1657712‑1658217
                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                            
gggcagcggncccAGAAGCACGGCACCCTTGCTGCCGTCGCGCCACTGCCGCGCGCGCATTTCACAAAACCCAATATCGTTATCAAACCTAACGCCAACAGCCGCCCGACCGGCGATACCACTGGTTATCTGGCTAATCCGGAGGAAGTGTAAAATGGGAAATGGATGGCATGAGTGGCCACTGGTTATCTTTACCGTATTGGGTCAATGCGTAGTAGGCGCGCTGATTGTTAGCGGAATTGGCTGGTTCGCCGCCAAAAATGACGCAGACCGTc                                                                                                                                                                                                                                         <  2:275153/264‑1 (MQ=255)
                                tgcCCTCACCCCACCGCCGCGCCGGCACTGGCCAAAACCCAATACTCGTATCAAACCTAACCCCAACAGCCGCCCGACTGGCGATACCACTGGTTATCGGGCTAATCCGGAGGAAGTGTAAAATGGGAAATGGATGGCATGAGTGGCCACTGGTTATCTTTACCGTATTGGGTCAATGCGTAGTAGGCGCGCTGATTGTTAGCGGAATTGGCTGGTTCGCCGCCAAAAATGTCGCAGACCGTCAGCGTATTGTCCGTGGCATGTTTTTCCTCTgg                                                                                                                                                                                                         <  2:190110/275‑1 (MQ=255)
                                                 cgcgcgcgCATTTCACAAAACCCAATATCGTTATCAAACCTAACGCCAACAGCCGCCCGACCGGCGATACCACTGGTTCTCTGGCTAATCCGGAGGAAGTGTAAAATGGGAAATGGATGGCATGAGTGGCCACTGGTTATCTTTACCGTATTGGGTCAATGCGTAGTAGGCGCGCTGATTGTTAGCGGAATTGGCTGGTTCGCCGCCAAAAATGACGCAGACCGTCAGCGTATTGTCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCg                                                                                                                                                                                         <  1:719217/274‑1 (MQ=255)
                                                                   cgaacacacctCGTTATAAAACCCAACGCCCACAGCCGGCCGACAGGCGATGCCACTGGTTATTGGGCTAAGCCGGAGGAAGTGTAAAATGGGAAATGGATGGCATGAGTGGCCACTGGTTATCTTTACCGTATTGGGTCAATGCGTAGTAGGCGCGCTGATTGTTAGCGGAATTGGCTGGTTCGCCGCCAAAAATGACGCAGACCGTCAGCGTATTGTCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCa                                                                                                                                                                       <  2:650845/265‑1 (MQ=255)
                                                                                       aCCTAACGCCAACAGCCGCCCGACCGGCGATACCACCGGTTATCTGGCTAATCCGGAGGAAGTGTAAAATGGGAAATGGATGGCATGAGTGGCCACTGGTTATCTTTACCGTATTGGGTCAATGCGTAGTAGGCGCGCTGATTGTTAGCGGAATTGGCTGGTTCGCCGCCAAAAATGTCGCAGACCGTCAGCGTATTGTCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCACt                                                                                                                                                    >  1:399994/1‑274 (MQ=255)
                                                                                                          ccGACCGGCGATACCACTGGTTATCTGGCTAATCCGGAGGAAGTGTAAAATGGGAAATGGATGGCATGAGTGGCCACTGGTTATCTTTACCGTATTGGGTCAATGCGTAGTAGGCGCGCTGATTGTTAGCGGAATTGGCTGGTTCGCCGCCAAAAATGACGCAGACCGTCAGCGTATTGTCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTg                                                                                                                                                                          >  1:682459/1‑233 (MQ=255)
                                                                                                              ccGGCGATACCACTGGTTATCTGGCTAATCCGGAGGAAGTGTAAAATGGGAAATGGATGGCATGAGTGGCCACTGGTTATCTTTACCGTATTGGGTCAATGCGTAGTAGGCGCGCTGATTGTTAGCGGAATTGGCTGGTTCGCCGCCAAAAATGACGCAGACCGTCAGCGTATTGTCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCAt                                                                                                                                                                      <  2:682450/233‑1 (MQ=255)
                                                                                                                                       tAATCCGGAGGAAGTGTAAAATGGGAAATGGATGGCATGAGTGGCCACTGGTTATCTTTACCGTATTGGGTCAATGCGTAGTAGGCGCGCTGATTGTTAGCGGAATTGGCTGGTTCGCCGCCAAAAATGACGCAGACCGTCAGCGTATTGTCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATTCACCTGTGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCGGGTCTGAGTCAc                                                                                                   >  2:232380/1‑275 (MQ=255)
                                                                                                                                                                          cATGAGTGGCCACTGGTTATCTTTACCGTATTGGGTCAATGCGTAGTAGGCGCGCTGATTGTTAGCGGAATTGGCTGGTTCGCCGCCAAAAATGACGCAGACCGTCAGCGTATTGTCCGTGGCATGTTTTTCctct                                                                                                                                                                                                           >  1:602303/1‑136 (MQ=255)
                                                                                                                                                                                  gCCACTGGTTATCTTTACCGTATTGGGTGAATGCGTAGTAGGCGCGCTGATTGTTAGCGGAATTGGCTGGTTCGCCGCCAAAAATGACGCAGACCGTCAGCGTATTGTCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCGGTCATGCACCTGGGATCGCCACTNCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTAt                                                        >  1:524460/1‑275 (MQ=255)
                                                                                                                                                                                         gTTATCTTTACCGTATTGGGTCAATGCGTAGTAGGCGCGCTGATTGTTAGCGGAATTGGCTGGTTCGCCGCCAAAAATGACGCAg                                                                                                                                                                                                                                               >  1:472893/1‑85 (MQ=255)
                                                                                                                                                                                             tCTTTACCGTATTGGGTCAATGCGTAGTAGGCGCGCTGATTGTTAGCGGAATTGGCTGGTTCGCCGCCAAAAATGACGCAGACCg                                                                                                                                                                                                                                           <  2:472886/85‑1 (MQ=255)
                                                                                                                                                                                                 tACCGTATTGGGTCAATGCGTAGTAGGCGCGCTGATTGTTAGCGGAATTGGCTGGTTCGCCGCCAAAAATGACGCAGACCGTCAGCGTATTGTCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGACACTGCGGGCGTTTAACTCCCTgaatcg                                                                                                                            >  2:387043/1‑192 (MQ=255)
                                                                                                                                                                                                   ccGTATTGGGTCAATGCGTAGTAGGCGCGCTGATTGTTAGCGGAATTGGCTGGTTCGCCGCCAAAAATGACGCAGACCGTCAGCGTATTGTCCGTGGCCTGtttt                                                                                                                                                                                                                 >  2:619815/1‑105 (MQ=255)
                                                                                                                                                                                                     gTATTGGGTCAATGCGTAGTAGGCGCGCTGATTGTTAGCGGAATTGGCTGGTTCGCCGCCAAAAATGACGCAGACCGTCAGCGTATTGTCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTgaatcgaatc                                                                                                                        <  1:387049/192‑1 (MQ=255)
                                                                                                                                                                                                       aTTGGGTCAATGCGTAGTAGGCGCGCTGATTGTTAGCGGAATTGGCTGGTTCGCCGCCAAAAATGACGCAGACCGTCAGCGTATTGTCCGTGGCATGTTTTTCct                                                                                                                                                                                                             <  1:619824/105‑1 (MQ=255)
                                                                                                                                                                                                                aTGCGTAGTAGGCGCGCTGATTGTTAGCGGAATTGGCTGGTTCGCCGCCAAAAATGACGCAGACCGTCAGCGTATTGTCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGc                          >  1:79532/1‑275 (MQ=255)
                                                                                                                                                                                                                aTGCGTAGTAGGCGCGCTGATTGTTAGCGGAATTGGCTGGTTCGCCGCCAAAAATGACGCAGACCGTCAGCGTATTGTCCGTGGCATGCTTTTACTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGTCACTGCGGGCGGTTAACTCCCGGAATCGAATCGGTGCGGCGGGTCTGAGTAACGGGATTGCGTCCGGTTACCTATTTTTGGCCGGTGGCGGATTGTGGGTGGGGGGGGCAGTAATCGGTAAAATgg                          >  2:632020/1‑274 (MQ=255)
                                                                                                                                                                                                                 tGCGTAGTAGGCGCGCTGATTGTTAGCGGAATTGGCTGGTTCGCCGCCAAAAATGACGCAGACCGTCAGCGTATTGTCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCg                                                                                                                                             >  2:601550/1‑159 (MQ=255)
                                                                                                                                                                                                                 tGCGTAGTAGGCGCGCTGATTGTTAGCGGAATTGGCTGGTTCGCCACCAAAAATGACGCAg                                                                                                                                                                                                                                               >  1:537079/1‑61 (MQ=255)
                                                                                                                                                                                                                     tagtagGCGCGCTGATTGTTAGCGGAATTGGCTGGTTCGCCGCCAAAAATGACGCAGACCGTCAGCGTATTGTCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTa                                                                                                                                         <  1:601559/159‑1 (MQ=255)
                                                                                                                                                                                                                     tagtagGCGCGCTGATTGTTAGCGGAATTGGCTGGTTCGCCACCAAAAATGACGCAGACCg                                                                                                                                                                                                                                           <  2:537071/61‑1 (MQ=255)
                                                                                                                                                                                                                     tagtagGCGCGCTGATTGTTAGCGGAATTGGCTGGTTCGCCACCAAAAATGACGCAGACCGTCAGCGTATTGTCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTg                                                                                                              >  2:262128/1‑186 (MQ=255)
                                                                                                                                                                                                                         agGCGCGCTGATTGTTAGCGGAATTGGCTGGTTCGCCACCAAAAATGACGCAGACCGTCAGCGTATTGTCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCt                                                                                                          <  1:262133/186‑1 (MQ=255)
                                                                                                                                                                                                                              cgcTGATTGTTAGCGGAATTGGCTGGTTCGCCGCCAAAAATGACGCAGACCGTCAGCGTATTGTCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCgg                                                                                                                                                                                        >  2:281197/1‑103 (MQ=255)
                                                                                                                                                                                                                                cTGATTGTTAGCGGAATTGGCTGGTTCGCCGCCAAAAATGACGCAGACCGTCAGCGTATTGTCCGTGGCATGtttt                                                                                                                                                                                                                 >  1:579801/1‑76 (MQ=255)
                                                                                                                                                                                                                                  gATTGTTAGCGGAATTGGCTGGTTCGCCGCCAAAAATGACGCAGACCGTCAGCGTATTGTCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGttt                                                                                                                                                                                    <  1:281202/103‑1 (MQ=255)
                                                                                                                                                                                                                                    ttGTTAGCGGAATTGGCTGGTTCGCCGCCAAAAATGACGCAGACCGTCAGCGTATTGTCCGTGGCATGTTTTTCct                                                                                                                                                                                                             <  2:579792/76‑1 (MQ=255)
                                                                                                                                                                                                                                         aGCGGAATTGCCTGGTTCGCCGCCTAAAATTACGCAACCCGTCAGCGTATTGTCCGCGTCATGTTTTTCCTCTGTCTGGTGATGGGCGTCGGGTTTATTGCCTCTGTCCTCCACCTGGGCTCGCCACTGCGGGCGTTTAACTACCTGACTCGAATCGGATCTTCTGGTCTGACTAACGATACTGCTGCCGGCTCCCTATTTTTTGCCTTTGGCGGATTCTTGTTGTGGGGGTCTGTCATACGTTCAATGCCGCCGTCGGTGGGGAAAGTCTggg  >  2:209198/1‑273 (MQ=255)
                                                                                                                                                                                                                                               aTTGGCTGGTTCGCCGCCAAAAATGACGCAGACCGTCAGCGTATTGTCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTAtt                                                                                                                                                                                 >  1:392824/1‑93 (MQ=255)
                                                                                                                                                                                                                                                   gCTGGTTCGCCGCCAAAAATGACGCAGACCGTCAGCGTATTGTCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCct                                                                                                                                                                             <  2:392818/93‑1 (MQ=255)
                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                            
GAAGAGCTGCGCCAGAAGCACGGCACTCTTGCTGCCGTCGCGCCACTGCCGCGCGCGCATTTCACAAAACCCAATATCGTTATCAAACCTAACGCCAACAGCCGCCCGACCGGCGATACCACTGGTTATCTGGCTAATCCGGAGGAAGTGTAAAATGGGAAATGGATGGCATGAGTGGCCACTGGTTATCTTTACCGTATTGGGTCAATGCGTAGTAGGCGCGCTGATTGTTAGCGGAATTGGCTGGTTCGCCGCCAAAAATGACGCAGACCGTCAGCGTATTGTCCGTGGCATGTTTTTCCTCTGGCTGTTGATGGGCGTCGGGTTTATTGCCTCTGTCATGCACCTGGGATCGCCACTGCGGGCGTTTAACTCCCTGAATCGAATCGGAGCTTCTGGTCTGAGTAACGAGATTGCTGCCGGTTCCATATTTTTTGCCGTTGGCGGATTATGGTGGCTGGTGGCAGTCATCGGTAAAATGCCGCAAGCGTTGGGGAAACTCTGGT  >  CP009273/1657712‑1658217

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: