Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 2,871,022 C→A 22.4% E61* (GAA→TAA)  U712_14805 ← Putative tRNA‑binding protein ytpR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068812,871,0220CA22.4% 48.1 / 2.5 33E61* (GAA→TAA) U712_14805Putative tRNA‑binding protein ytpR
Reads supporting (aligned to +/- strand):  ref base C (4/18);  new base A (6/2);  total (13/20)
Fisher's exact test for biased strand distribution p-value = 7.23e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTTCTTTTGATTCTACATAACCCACAACAAACTTCGGAGATAAATCAACGACTAATGTTTCTTCAACGCCGTTTCTGTTTAAAATCTCATTCACATCTTGCACGAATGTTTCTGAGAGCGCAACAGGGCCGTTTTCATCAATGGTCAAGTAAGAAG‑ACGCATTGAAAATGTTGAAGCCCGTTGTTTCTTTTGTTTCGTTATTGAAAATTTTGACGACGTCGCCGTGTTTTTCATAGCCTAATTGTTCGCGTGTCACATCTTGAAGA  >  CP006881/2870886‑2871152
                                                                                                                                        |                                                                                                                                   
ttttCTTTTGATTCTACATAACCCACAACAAACTTCGGAGATAAATCAACGACTAATGTTTCTTCAACGCCGTTTCTGTTTAAAAACTCATTCACATCTTGCACGAATGTTTCTGAGAGCGCAACAGGGCCGTTTtca                                                                                                                                    >  5:85590/1‑138 (MQ=255)
      tttGATTCTACATAACCCACAACAAACTTCGGAGATAAATCAACGACTAATGTTTCTTCAACGCCGTTTCTGTTTAAAATCTCATTCACATCTTGCACGAATGTTTCTGAGAGCGCAACAGGGCCGTTTTCCTGAATgg                                                                                                                             >  6:262990/1‑139 (MQ=255)
        tGATTATACATAACCCACAACAAACTTCGGAGATAAATCAACGACTAATGTTTCTTCAACGCCGTTTCTGTTTAAAATCTCATTCACATCTTGCAAGAATGTTTCTTAGAGCGCAACAGGGCCGTTTTCATCAATGGt                                                                                                                            <  4:110005/138‑1 (MQ=255)
         gATTCTACATAACCCACAACAAACTTCGGAGATAAATCAACGACTAATGTTTCTTCAACGCCGTTTCTGTTTAAAATCTCATTCACATCTTGCACGAATGTTTCTGAGAGCGCAACAGGGCCGTTTTCATCAATGGTc                                                                                                                           <  8:210469/138‑1 (MQ=255)
         gATTCTACATAACCCACAACAAACTTCGGAGATAAATCAACGACTAATGTTTCTTCAACGCCGTTTCTGTTTAAAATCTCATTCACATCTTGCACGAATGTTTCTGAGAGCGCAACAGGGCCGTTTTACTCAATTGtg                                                                                                                           >  7:135410/1‑137 (MQ=255)
              tACATAACCCACAACAAACTTCGGAGATAAATCAACGACTAATGTTTCTTCAACGCCGTTTCTGTTTAAAATCTCATTCACATCTTGCACGAATGTTTCTGAGAGCGCAACAGGGCCGTTTTCATCAATGGTCaatgaa                                                                                                                     >  8:307820/1‑135 (MQ=255)
               aCATAACCCACAACAAACTTCGGAGATAAATCAACGACTAATGTTTCTTCAACGCCGTTTCTGTTTAAAATCTCATTCACATCTTGCACGAATGTTTCTGAGAGCGCAACAGGGCCGTTTTCATCAATGGTCAAGTaag                                                                                                                    <  5:170594/139‑1 (MQ=255)
                  tAACCCACAACAAACTTCGGAGATAAATCAACGACTAATGTTTCTTCAACGCCGTTTCTGTTTAAAATCTCATTCACATCTTGCACGAATGTTTCTGAGAGCGCAACAGGGCCGTTTTCATCAATGGTCAAGTaagaag                                                                                                                 <  7:101427/139‑1 (MQ=255)
                    aCCCACAACAAACTTCGGAGATAAATCAACGACTAATGTTTCTTCAACGCCGTTTCTGTTTAAAATCTAATTCACATCTTGCACGAATGTTTCTGAGAGCGCAACAGGGCCGTTTTCATCAATGGTCAAGTAagaag‑a                                                                                                               <  5:291591/138‑1 (MQ=255)
                     cccACCACAAACTTCGGAGAGAAATCAACGACTTCAGTTTCTTCAACGCCGTTTCTGGTTAAAATCTCATTCACATCTTGCACGAATGTTGCTGAGAGCGCAACAGGGCCGTTTACATCAATGGGCAAGTAAGAAG‑ACg                                                                                                             <  4:125181/139‑1 (MQ=255)
                          aacaaACTTCGGAGATAAATCAACGACTAATGTTTCTTCAACGCCGTTTCTGTTTAAAATCTCATTCACATCTTGCACGAATGTTTCTGAGAGCGCAACAGGGCCGTTTTAATAAATGGTGAAGTAAAAAG‑ACGCAtt                                                                                                         >  6:291591/1‑138 (MQ=255)
                                   cGGAGATAAATCAACGACTAATGTTTCTTCAACGCCGTTTCTGTTTAAAATCTCATTCACATCTTGCACGAATGTTTCTGAGAGCGCAACAGGGCCGTGTTGATCAATTGTCAAGTAACAAT‑CCGCAATTAAAAAGTTg                                                                                               >  3:310011/1‑139 (MQ=255)
                                           aaTCAACGACTAATGTTTCTTCAACGCCGTTTCTGTTTAAAATCTCATTCACATCTTGCACGAATGTTTCTGAGAGCGCAACAGGGCCGTGTTAATAAAAGTGCAAGTAAAAAA‑ACGCCTTTAAAAATTTGAATccccg                                                                                       >  6:116161/1‑137 (MQ=255)
                                             acagcGACTTATGTTTCTTCATCGCAGTTTCTGTTTAAAATCTCATTCACTTCTTGCACTAATGTTTCTGAGAGCGCAACAGGTCCTTTTTAATCAATTGTCAAGTAAGAAG‑ACGCATTGAAAATGTTGAAGCCCgtt                                                                                      <  5:95052/134‑1 (MQ=255)
                                               aaCGACTAATGTTTCTTCAACGCCGTTTCTGTTTAAAATCTCATTCACATCTTGCACGAATGTTTCTGAGAGCGCAACAGGGCCGTTTTCATCAATGGTCAAGTAAGAAG‑ACGCATTGAAAATGTTGAAGCCCgttgt                                                                                    <  6:348878/138‑1 (MQ=255)
                                                                  aCGCCTTTTCTCTTTAAAAACTCATTCACATCTTGCACGAATGTTTCTGAGAGCGCAACAGGGCCGTTTTAAAAAAATGGGAATTAATAAC‑ACGCAATGAAAAACTTGAACCCCGTTGTTTTTTTTTTTTTATTATTg                                                                 >  2:272145/1‑138 (MQ=255)
                                                                   cGCCGTTTCTGTTTAAAATCTCATTCACATCTTGAACGAATGTTTCTGAGAGCGCAAAAGGGCCGTTTTAAAAAATTGTAAATAAAGAAC‑ACCCATTGAAAAATTATAAGCCCGGTGTTTATTTTTTTtatttattta                                                                >  2:128306/1‑136 (MQ=255)
                                                                        tttCTGTTTAAAATCTCATTCACATCTTGCACGAATGTTTCTGAGAGCGCAACAGGGCCGTGTTGATCAAATTGGAAGTAATAAGAACCCATTAAAAATTTTAACGCCC‑TTGTTTTGTTTGTTTTGGTTTTTTAAAttt                                                          >  1:70840/1‑139 (MQ=255)
                                                                              tttAAAATCTCATTCACATCTTGCACGAATGTTTCTGAGAGCGCAACAGGGCCGTTTTAATCAATTGTCAAGTAAGAAG‑AAGCATTGAAAATTTTGAAGCCCGGTGTTTGTTTTTTTTCGTTATTTAAAAATTTtacga                                                    >  5:145583/1‑139 (MQ=255)
                                                                                               aTCTTGCACGAATGTTTCTGAGAGCGCAACAGGGCCGTTTTCATCAATGGTCAAGTAAGAAG‑ACGCATTGAAAATGTTGAAGCCCGTTGTTTCTTTTGTTTCGTTATTTAAAATTTTGACGACGTCGCCGTGTTTTTCa                                   <  5:116161/139‑1 (MQ=255)
                                                                                               aTCTTGCACCAATGGTTCGTCGCCCGCACCAGGGCCGTTTTCATCATTGGTCAAGTATGAAG‑ACGCATTGAAAATGTTGAAGCCCGTTGTTTCTTTTGTTTCGTTATTTAAAATTTTGACGACGTCGCCGTGTTTTTCa                                   <  2:197825/139‑1 (MQ=255)
                                                                                                     cACGAATGTTTCTGAGAGCGCAACCGGGCCGTTTTCATCAATGGTCAAGTAAGAAG‑ACGCATTGAAAATGTTGAAGCCCGTTGTTTCTTTTGTTTCGTTATTGAAAATTTTGACGACGTCGCCGTGTTTTTCATAGCCt                             <  7:307820/139‑1 (MQ=255)
                                                                                                      aCGAATGTTTCTGAGAGCGCAACAGGGCCGTTTTCATCAATGGTCAAGTAAGAAG‑ACGCATTGAAAATGTT‑AAGCCCGTTGTTTCTTTTGTTTCGTTATTGAAAATTTTGACGACGTCGCCGTGTTTTTCATAGCCTaa                           <  1:6810/139‑1 (MQ=255)
                                                                                                         aaTGTTTCTGAGAGCGCAACAGGGCCGTTTTCATAAATGGGCAAGTAAGAAG‑TAGCCTTTAAAATTTTGAAGCCCGTTGTTTGTTTTGTTTCGGTTTTGAATATTTTGACGACGTCGCCGTGTTGTTCTGAGCGTAAtt                         >  6:353422/1‑139 (MQ=255)
                                                                                                         aaTGTTTCTGAGAGCGCAACAGGGCCGTGTTGATCAATGGTCAAGTTAGAAG‑ACGCCTTGAAAATTTTGAAGCCCGGTGGTGGTTTTGTTTGGTTTTTGAAAATTTTTTCGCCGTCGGCGGGGTTTTTATAGACGAAtt                         >  2:214044/1‑139 (MQ=255)
                                                                                                               tCTGAGAGCGCAACAGGGCCGTTTTCATCAATGGTCAAGTAAGAAG‑ACGCATTGAAAATGTTGAAGCCCGTTGTTTCTTTTGTTTCGTTTTTGAAAATTTTGACGACGTCGCCGTGGTTTTGTTAGGCGAATTTTTggc                   >  3:257773/1‑139 (MQ=255)
                                                                                                                cTGAGAGCGCAACAGGGCCGTTTTCATCAATGGTCAAGTAAGAAG‑ACGCATTGAAAATGTTGAAGCCCGTTGTTTCTTTTGTTTCGTTATTGAAAATTTTGACGACGTCGCCGTGTTTTTCATAGCCTAATTGTTcgc                   <  6:145583/138‑1 (MQ=255)
                                                                                                                    gcgcgcAACAGGGCCGTTTTCATCAATGGTCAAGAAAGAAG‑ACGCATTGAAACTGTTGAAGCCCGTTGTTTCTTTTGTTTCGTTATTGAAAATTTTGACGACGTCGCCGTGTTTTTCATAGCCTAATTGTTCGCGTGTc              <  8:135410/137‑1 (MQ=255)
                                                                                                                      gcgcAACATGGCATTTTTCATCATTGGTCAAGTAAGAAG‑ACGCATTGAAAATGTTGAAGCCCGTTGTTTCTTTTGTTTCTTTTTTTAAAATTTTGACGACGTCGCCGTGTTTTTCATAGCCTAATTGTTCGCGTGTcac            <  1:128306/139‑1 (MQ=255)
                                                                                                                      gcgcAACAGGGCCGTTTTCATCAATGGTCAAGTAAGAAG‑ACGCATTGAAAATGTTGAAGCCCGTTGTTTCTTTTGTTTCGTTATTGAAAATTTTGACGACGTCGCCGTGTTTTTCATAGCCTAATTGTTCGCGTGTcac            <  6:246493/139‑1 (MQ=255)
                                                                                                                                gCCGTTTTCATCAATGGTCAAGTAAGAAG‑ACGCATTGAAAATGTTGAAGCCCGTTGTTTCTTTTGTTTCGTTATTGAAAATTTTGACGACGTCGCCGTGTTTTTCATAGCCTAATTGTTCGCGTGTCACATCTTGAaga  <  5:70276/139‑1 (MQ=255)
                                                                                                                                gCCGTTTTCATCAATGGTCAAGTAAGAAG‑ACGCATTGAAAATGTTGAAGCCCGTTGTTTCTTTTGTTTCGTTATTGAAAATTTTGACGACGTCGCCGTGTTTTTCATAGCCTAATTGTTCGCGTGTCACATCTTGAaga  <  2:129423/139‑1 (MQ=255)
                                                                                                                                gCCGTTTTCATCAATGGTCAAGTAAGAAC‑ACGCATTGAAAATGTGGAAGCCCGTTGTTTCTTTTGTTTCGTTATTGAAAATTTTGACGACGTCGCCGTGTTTTTCATAGCCTAATTGTTCGCGTGTCACATCTTGAaga  <  1:93043/139‑1 (MQ=255)
                                                                                                                                gCCGTTTTAATCAATGGTAAAGTAAGAAC‑ACGCATTTAAAATTTTGAAGCCCGTTGTTTCTTTTGTTTCGTTATTGAAAATTTTGACGACGTCGCCGTGTTTTTCATAGCCTAATTGTTCGCGTGTCACATCTTGAaga  <  6:85590/139‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                   
TTTTCTTTTGATTCTACATAACCCACAACAAACTTCGGAGATAAATCAACGACTAATGTTTCTTCAACGCCGTTTCTGTTTAAAATCTCATTCACATCTTGCACGAATGTTTCTGAGAGCGCAACAGGGCCGTTTTCATCAATGGTCAAGTAAGAAG‑ACGCATTGAAAATGTTGAAGCCCGTTGTTTCTTTTGTTTCGTTATTGAAAATTTTGACGACGTCGCCGTGTTTTTCATAGCCTAATTGTTCGCGTGTCACATCTTGAAGA  >  CP006881/2870886‑2871152

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: ‑