Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 3,837,361 A→C 26.2% F34C (TTT→TGT)  U712_19795 ← Uncharacterized protein yxiK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068813,837,3610AC26.2% 79.8 / 13.5 48F34C (TTT→TGT) U712_19795Uncharacterized protein yxiK
Reads supporting (aligned to +/- strand):  ref base A (33/2);  new base C (2/11);  total (35/13)
Fisher's exact test for biased strand distribution p-value = 2.43e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTCCCTTTTGGCCAGCTTCGTCACCTTGTAAGCATGACCTTCAAAATACTCGGACAGGTCTTCTATAAATGAACCAGTCGTTATTACCTGTATCATGATTCGATCGTAGGCAAGCCCTGCCCAAAAACGATAAAACGACTCTGCATCTTCAAGTTCGCATGCCATAACGGTTTGGTCAGGCTTAAGCTCGATCATGTCCAAATTGCCAATACTGCTTATAGAGGTTATTTCCATAAAATCTTCCCCTTATCAGTTGGATAGCGCAACG  >  CP006881/3837227‑3837495
                                                                                                                                        |                                                                                                                                      
aTTCCCTTTTGGCCAGCTTCGCCACCTTGAAAGCATGCCATACAAAAAACGCGCACAGGTCTTATAAAAATCACCCCGTCTTTATACCCTGTACCATGATGCGATCGTCGGCACGCCCTCCCAAAAAACGATAAAACGa                                                                                                                                    <  8:216854/139‑1 (MQ=255)
      ttttGGCCAGCTTCGTCACCTTGTAAGCATGACCTTCAAAATACTCGGACAGGTCTTCTATAAATGAACCAGTCGTTATTACCTGTATCATGATTAGATCGTAGGCAAGCCCTGCCCAAAAACGATAAAAAGACTATg                                                                                                                               >  7:209946/1‑138 (MQ=255)
                cTTCGTCACCTTGTAAGCATGACCTTCAAAATACTCGGACAGGTCTTCTATAAATGAACCAGTCGTTATTACCTGTATCATGATTCGATCGTAGGCAAGCCCTGCCCAAAAACGATAAAACGACTCTGCATCTTCAAGt                                                                                                                    >  3:45475/1‑139 (MQ=255)
                     tCACGTTGTACCCATGACATTAAAAAAACTCCGACAGGTTTTATAAAAATCACCCAGTGTTTATTACGTTTAAAATGTTCCTTCCGACGCACCCCCCTCCCAAAAACCTAAAACACGACTCTGCATCTTCAAGTTCGCa                                                                                                               <  8:156220/139‑1 (MQ=255)
                                   tGACCTTCAAAATACTCGGACAGGTCTTCTATAAATGAACCAGTCGTTTTTACCTGTTTAATGATTGGATCGTATGGAAGCCCAGCCCAAAAACGATAAAACGACTCTGGATCTTCAACTTCGCATGCCAAAACGGttt                                                                                                 >  3:179667/1‑139 (MQ=255)
                                   tGACCTTCAAAATACTCGGACAGGTCTTCTATAAATGAACCAGTCGTTATTACCTGTATCATGATTCGATCGTAGGCAAGCCCTGCCCAAAAACGATAAAACGACTCTGCATCTTCAAGTTCGCATGCCATAACGGttt                                                                                                 >  7:223260/1‑139 (MQ=255)
                                   tGACCTTCAAAATACTCGGACAGGTCTTCTATAAATGAACCAGTCGTTATTACCTGTAACATTATTCGATCGTAGGCAAGCCAGCCCCAAAAACAAAAACACCACTATGCAAATTAAAGTTCGAATGCCATAACCGttt                                                                                                 >  4:264557/1‑139 (MQ=255)
                                    gACCTTCAAAATACTCGGACAGGTCTTCTATAAATGAACCAGTCGTTATTACCTGTATCATGATTAGATCGTAGGCAAGCCCTGCCCAAAAACGATAAAACGACTCTGCATCTTCAAGTTCGCATGCCATAACGGTTTg                                                                                                >  8:196243/1‑139 (MQ=255)
                                                 cTCGGACAGGTCTTCTATAAATGAACCAGTCGTGATTAACTGTATAATGATTAGATCGTAGGCAAGCCCGGACCCAAAACGCTAAAAAGACAGTGCATCGTCAATTTCGCATGCCAATACAGTTTGGTCAGGCTTAAGc                                                                                   >  6:296323/1‑139 (MQ=255)
                                                    ggACAGGTCTTCTATAAATGAACCAGTCGTTATTTCCTGTATCATGATTCGATCGTAGGCAAGCCCTGCCCAAAAACGATAAAACGACTCTGCATCTTCAACTTCGCATGCCATAACGGTTTGGTCAGGCTGAAGcgc                                                                                 >  8:95722/1‑136 (MQ=255)
                                                    ggACAGGTCTTCTATAAATGAACCAGTCGTTATTACCTGTATCATGATTCGATCGTAGGCAAGCCCTGCCCAAAAACGATAAAACGACTCAGCATCTTCAACTTCGCATGCCATAACGGTTTGGTCAGGCTTAAGCTCg                                                                                >  5:237076/1‑139 (MQ=255)
                                                      aCAGGTCTTATAAAAATCAACCCGTCTTTATTACCTTTATCATTATGCGATCTTCGGCACGCCCTGCCAAAAAACTAAAACACGACTCTGCATCTTCAAGTTCGCATGCCATAACGGTTTGGTCAGGCTTAAGCTCGAt                                                                              <  7:293086/139‑1 (MQ=255)
                                                                aataAATGAACCAGTCGTTATTACCTTTATCATTATGCGATCGTAGGCAACCCCTGCCCAAAAACTAAAACACGACTCTGCATCTTCAAGTTCGCATGCCATAACGGTTTGGTCAGGCTTAAGCTCGATCATGTCCaa                                                                     <  8:223260/137‑1 (MQ=255)
                                                                  aaaaTCAACCAGTCTTAATAACCTGTATCATGATTCGATCGTAGCCAACCCCTCCCAAAAAACCAAAACACGACTCTGCATCTTCAAGTTCGCATGCCATAACGGTTTGGTCAGGCTTAAGCTCGATCATGTCCaa                                                                     <  3:124645/135‑1 (MQ=255)
                                                                    aaTGAACCAGTCGTTATTTCCTGTATCATTAATCGATCGTAGGCAAGACAAGCCCAAAAAAGATAAAACGACTCTGCATCTTCAAGTTCGCTATCCAGAACCGTTTTGTAAGGCTTTAGCTCGATCATTTACAAATTGc                                                                >  2:203268/1‑139 (MQ=255)
                                                                    aaTGAACCAGTCGTTATTACCTGTATCATGATTCGATCGTAGGCAAGCCCTGCCCAAAAACGATAAAACGACTCTGCATCTTCAAGTTCGCATGCCATAACGGTTTGGTCAGGCTTAAGCTCGATCATGTCCAAATTGc                                                                >  5:211832/1‑139 (MQ=255)
                                                                    aaTGAACCAGTCGTTATTACCTGTATCATGATTCGATCGTAGGCAAGCCCTGCCCAAAAACGATAAAACGACTCTGCATCTTCAAGTTCGCATGCCATAACGGTTTGGTCAGGCTTAAGCTCGATCATGTCCAAATTGc                                                                >  5:127628/1‑139 (MQ=255)
                                                                     aTGAACCAGTCGTTATTACCTGTATCATGATTCGATCGTAGGCAAGCCCTGCCCAAAAACGATAACACGACTCTGCATCTTCAAGTTCGCATGCCATAACGGTTTGGTCAGGCTTAAGCTCGATCATGTCCAAATTGc                                                                >  4:169920/1‑138 (MQ=255)
                                                                        gacatgTCGTTATTACCTGTATCATGATTCGATCGTAGGCAAGCCCTGCCCAAAAACGATAAAACGACTCTGCATCTTCAAGTTCGCATGCCATAACGGTTTGGTCAGGCTTAAGCTCGATCATGTCCAAATTGCCaa                                                             >  5:15202/6‑138 (MQ=255)
                                                                         aCCAGTCGTTATTACCTGTATCATGATTCGATCGTAGGCAAGCCCTGCCCAAAAACGATAAAACGACTCTGCATCTTCAAGTTCGCATGCCATAACGGTTTGGTCAGGCTTAAGCTCGATCATGTCCAAATTGCCAATa                                                           >  6:109338/1‑139 (MQ=255)
                                                                            cggcGTTATTACCTGTATCATGTTGCGAGCGTCGGCAAGCCCGGCCAAAAAACTAAAACACGACTCTGCATCTTCAAGTTCGCATGCCATAACGGTTTGGTCAGGCTTAAGCTCGATCATGTCCAAATTGCCAATACTg                                                        <  7:98161/136‑1 (MQ=255)
                                                                             gTCGTTATTTCCTGTATCATGATTCGATCGTAGGCAAGCCCTGCCCAAAAACGATAAAACGACTCTGCCTCTTCAAGTTCGCATGCCATAACGGTTTGGTCAGGCTTAAGCTCGATCATGTCCAAATTGCCAATACTGc                                                       >  7:252540/1‑139 (MQ=255)
                                                                                gTTATTACCTGTATCATGATTCGATCGTAGGCAAGCCCTGCCCAAAAACGATAAAACGACACGGCATCTTCAAGTTCGCAAGCCATAACGGTTTGGTCAGGCTTAAGCTCGATCAAGTCCAAATTGCCAATACTGCTta                                                    >  7:203682/1‑139 (MQ=255)
                                                                                          gTATCATGATTCGATCGTAGGCAAGCCCTGCCCAAAAACGATAAAACGACTCTGCATCTTCAAGTTCGCATGCCATAACGGTTTGGTCAGGCTTAAGCTCGATCATGTCCAAATTGCCAATACTGCTTATAGAGGTTAt                                          >  7:109695/1‑139 (MQ=255)
                                                                                           tATCATGATTCGATCGTAGGCAAGCCCTGCCCAAAAACGATAAAACGACTCTGCATCTTCACGTGCGCATGCCATATCGGGTGGGTCAGGCTTAAGCTCGCTCATGTCCAACTTGCCAATCCTGCTTATCGAGGTTAtt                                         >  6:49694/1‑139 (MQ=255)
                                                                                           tATCATGATTCGATCGTAGGCAAGCCCTGCCCAAAAACGATAAAACGACTCTGCATCTTCAAGTTCGCATGCCATAACGGTTTGGTCAGGCTTAAGCTCGATCATGTCCAAATTGCCAATACTGCTTATATAGGTTAtt                                         >  7:252692/1‑139 (MQ=255)
                                                                                           tATCATGATTCGATCGTAGGCAAGCCCTGCCCAAAAACGATAAAACGACTCTGCATCTTCAAGTTCGCATGCCATAACGGTTTGGTCAGGCTTAAGCTCGATCATGTCCAAATTGCCAATACTGCTTATAGAGGTTAtt                                         >  3:275013/1‑139 (MQ=255)
                                                                                            aTCTTGATTCGAGCGTAGCCACCCCCGCCCAAAAAACCAAAACACGACTCTGCATCTTCAAGTTCGCATGCCATAACGGTTTGGTCAGGCTTAAGCTCGATCATGTCCAAATTGCCAATACTGCTTATAGAGGTTAtt                                         <  4:275013/138‑1 (MQ=255)
                                                                                                 gATTCGATCGTAGGCAAGCCCTGCCCAAAAACGATAAAACGACTCTGCATCTTCAAGTTCGCATGCCATAACGGTTTGGTCAGGCTTAAGCTCGATCATGTCCAAATTACCAATACTGCTTATAGAGGTTATTTCCAt                                    >  3:209693/1‑138 (MQ=255)
                                                                                                  aTTCGTTCTTAGCCACGCGCTCCCAAAAAAGAATAACACGACTCTGCATCTTCAAGTTCGCATGCCATAACGGTTTGGTCAGGCTTAAGCTCGATCATGTCCAAATTGCCAATACTGCTTATAGAGGTTATTTCCATa                                   <  4:21646/138‑1 (MQ=255)
                                                                                                  aTTCGATCGTAGGCAAGCCCTGCCCAAAAACGATAAAACGACTCTGCATCTTCAAGTTCGCATGCCATAACGGTTTGGTCAGGCTTAAGCTCGATCATGTCCAAATTGCCAATACTGCTTATAGAGGTTATTTCCATa                                   >  2:12296/1‑138 (MQ=255)
                                                                                                  aTTCGATCGTAGGCAAGCCCTGCCCAAAAACGATAAAACGACTCTGCATCTTCAAGTTCGCATGCCATAACGGTTTGGTCAGGCTTAAGCTCGATCATGTCCAAATTGCCAATACTGCTTATAGAGGTTATTTCCATa                                   >  8:46403/1‑138 (MQ=255)
                                                                                                       aTCGTAGGCAAGCCCTGCCCAAAAACGATAAAACGACTCTGCATCTTCAAGTTCGCATGCCATAACGGTTTGGTCAGGCTTAAGCTCGATCATGTCCAAATTGCCAATACTGCTTATAGAGGTTATTTCCATAAAAtat                             >  2:250237/1‑137 (MQ=255)
                                                                                                       aTCGTAGGCAAGCCCTGCCCAAAAACGATAAAACGACTCTGCATCTTCAAGTTCGCATGCCATAACGGTTTGGTCAGGCTTAAGCTCGATCATGTCCAAATTGCCAATACTGCTTATAGAGGTTATTTCCATAAAATCt                             >  6:68285/1‑139 (MQ=255)
                                                                                                       aTCGTAGGCAAGCCCTGCCCAAAAACGATAAAACGACTCTGCATCTTCAAGTTCGCATGCCATAACGGTTTGGTCAGGCTTAAGCTCGATCATGTCCAAATTGCCAATACTGCTTATAGAGGTTATTTCCATAAAATCt                             >  4:10617/1‑139 (MQ=255)
                                                                                                       aTCGTAGGCAAGCCCTGCCCAAAAACGATAAAACGACTCTGCATCTTCAAGTTCGCATGCCATAACGGCTTGGTCAGGCTTAAGCTCGATCATGTCCAAATTGCCAATACTGCTTATAGAGGTTATTTCCATAAAATCt                             >  6:212957/1‑139 (MQ=255)
                                                                                                             ggCAAGCCCTGCCCAAAAACGATAAAACGACTCTGCATCTTCAAGTTCGCATGCCATAACGGTTTGGTCAGGCTTAAGCTCGATCATGTCCAAATTGCCAATACTGCTTATAGAGGTTATTTCCATAAAATCTTCCCCt                       >  5:267047/1‑139 (MQ=255)
                                                                                                             ggCAAGCCCTGCCCAAAAACGATAAAACGACTCTGCATCTTCAAGTTCGCATGCCATAACGGTTTGGTCAGGCTTAAGCTCGATCATGTCCAAATTGCCAATACTGCTTATAGAGGTTAGTTCCATAAAATCTTCCCCt                       >  8:83131/1‑139 (MQ=255)
                                                                                                              gCATGCCCTCCCCAAAAACTAAAACACGACTCTGCATCTTCAAGTTCGCATGCCATAACGGTTTGGTCAGGCTTAAGCTCGATCATGTCCAAATTGCCAATACTGCTTATAGAGGTTATTTCCATAAAATCTTCCCCtt                      <  8:203682/139‑1 (MQ=255)
                                                                                                                acgCCCTCCCAAAACACTAAAAAACGACTCTGCATCTTCAAGTTCGCATGCCATAACGGTTTGGTCAGGCTTAAGCTCGATCATGTCCAAATTGCCAATACTGCTTATAGAGGTTATTTCCATAAAATCTTCCCCTTAt                    <  7:265943/137‑1 (MQ=255)
                                                                                                                   cccTGCCCAAAAACTAAAACACGACTCTGCATCTTCAAGTTCGCATGCCATAACGGTTTGGTCAGGCTTAAGCTCGATCATGTCCAAATTGCCAATACTGCTTATAGATGTTATTTCCATAAAATCTTCCCCTTATCAg                 <  2:79844/139‑1 (MQ=255)
                                                                                                                   cccTGCCCAAAAACGAAAACACGACTCTGCATCTTCAAGTTCGCATGCCATAACGGTTTGGTCAGGCTTAAGCTCGATCATGTCCAAATTGCCAATACTGCTTATAGAGGTTATTTCCATAAAATCTTCCCCTTATCAg                 <  5:109338/139‑1 (MQ=255)
                                                                                                                        cccAAAAACGATAAAACGACTCTGCATCTTCAAGTTCGCATGCCATAACGGTTTGGTCAGGCTTAAGCTCGATCATGTCCAAATTGCCAATACTGCTTATAGAGGTTATTTCCATAAAATCTTCCCCTTATCAGTTGGa            >  7:1488/1‑139 (MQ=255)
                                                                                                                        cccAAAAACGATAAAACGACTCTGCATCTTCAAGTTCGCATGCCATAACGGTTTGGTCAGGCTTAAGCTCGATCATGTCCAAATTGCCAATACTGCTTATAGAGGTTATTTCCATAAAATCTTCCCCTTATCAGTTGGa            >  2:233893/1‑139 (MQ=255)
                                                                                                                        cccAAAAACGATAAAACGACTCTGCATCTTCAAGTTCGCATGCCATAACGGTTTGGTCAGGCTTAAGCTCGATCATGTCCAAATTGCCAATACTGCTTATAGAGGTTATTTCCATAAAATCTTCCCCTTATCAGTTGGa            >  1:230394/1‑139 (MQ=255)
                                                                                                                         ccAAAAACATATAACACGACTCTGCATCTTCAAGTTCGCATGCCATAACGGTTTGGTCAGGCTTAAGCTCGATCATGTCCAAATTGCCAATACTGCTTATAGAGGTTATTTCCATAAAATCTTCCCCTTATCAGTTGGa            <  6:76944/139‑1 (MQ=255)
                                                                                                                                    aTAAAACGACTCTGCATCTTCAAGTTCGCATGCCATAACGGTTTGGTCAGGCTTAAGCTCGATCATGTCCAAATTGCCAATACTGCTTATAGAGGTTATTTCCATAAAATCTTCCCCTTATCAGTTGGATAGcgcaacg  >  1:121155/1‑139 (MQ=255)
                                                                                                                                    aTAAAACGACTCTGCATCTTCAAGTTCGCATGCCATAACGGTTTGGTCAGGCTTAAGCTCGATCATGTCCAAATTGCCAATACTGCTTATAGAGGTTATTTCCATAAAATCTTCCCCTTATCAGTTGGATAGcgcaacg  >  5:130486/1‑139 (MQ=255)
                                                                                                                                        |                                                                                                                                      
ATTCCCTTTTGGCCAGCTTCGTCACCTTGTAAGCATGACCTTCAAAATACTCGGACAGGTCTTCTATAAATGAACCAGTCGTTATTACCTGTATCATGATTCGATCGTAGGCAAGCCCTGCCCAAAAACGATAAAACGACTCTGCATCTTCAAGTTCGCATGCCATAACGGTTTGGTCAGGCTTAAGCTCGATCATGTCCAAATTGCCAATACTGCTTATAGAGGTTATTTCCATAAAATCTTCCCCTTATCAGTTGGATAGCGCAACG  >  CP006881/3837227‑3837495

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