Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,658,664 A→T 7.6% A208A (GCT→GCA)  yfkO ← FMN reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,658,6640AT 7.6% 123.0 / 2.8 48A208A (GCT→GCA) yfkOFMN reductase
Reads supporting (aligned to +/- strand):  ref base A (16/27);  new base T (2/2);  total (19/29)
Fisher's exact test for biased strand distribution p-value = 6.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.60e-01

TTTTCACGCAAAAGCCTGCGATTCATTTATGTAATCGCAGGCTTTTTTATTGCTAACTTTTTAATTTTATATATTTAGACATACGGTATTCTATTTATTCTTTTGGTCCAACCCATTCAATAACATCTTCATAAGCTTGGCGTGTTTTATTTTTCGGTGGCTCTTGTTGTCTGTAGCCAAATGCGACCATCACTGATAAACCAAATTGCTCAGTATCTAAGATCCCTTTATTTGCTAAAATGTCTGTCACTGTATCCAGACTAAAACCTTCC  >  CP000730/2658531‑2658802
                                                                                                                                     |                                                                                                                                          
ttttttttttttttccTGCGATTCATTTATGTAATCGCAGGCTTTTTTATTGCTAACTTTTTAATTTTATATATTTAGACATACGGTATTCTATTTATTCTTTTGTTCCAACCCATTCAATAACATCTTCATAAGCTTTGc                                                                                                                                     >  2:53555/15‑141 (MQ=255)
ttttCACGCAAAAGCCTGCGATTCATTTACGTAATCGCAGGCTTTTTTATTGCTAACTTTTTAATTTTATATATTTAGAAATACGGTATTCTATTTATTCTTTTGGTCCAACCCATTCAATAACATCTTCATAAGCTTGGc                                                                                                                                     <  2:421863/141‑1 (MQ=255)
     aCGCAAAAGCCTGCGATTCATTTATGTAATCGCAGGCTTTTTTATTGCTAACTTTTTAATTTTATATATTTAGACATACGGTATTCTATTTATTCTTTTGGTCCAACCCATTCAATAACATCTTCATAAGCTTGGCGTGtt                                                                                                                                >  2:57677/1‑141 (MQ=255)
     aCGCAAAAGCCTGCGATTCATTTATGTAATCGCAGGCTTTTTTATTGCTAACTTTTTAATTTTATATATTTAGACATACGGTATTCTATTTATTCTTTTGGTCCAACCCATTCAATAACATCTTCATAAGCTTGGCGTGtt                                                                                                                                <  1:230182/141‑1 (MQ=255)
          aaaGCGTGCGATTCATTTATGTAATCGCAGGCTTTTTTATTGCTAACTTTTTAATTTTATATATTTAGACATACGGTATTCTATTTATTCTTTTGGTCCAACCCATTCAATAACATCTTCATAAGCTTGGCGTGTTTTAtt                                                                                                                           <  2:395076/141‑1 (MQ=255)
          aaaGCCTGCGATTCATTTATGTAATCGCAGGCTTTTTTATTGCTAACTTTTTAATTTTATATATTTAGACATACGGTATTCTATTTATTCTTTTGGTCCAACCCATTCAATAACATCTTCATAAGCTTGGCGTGTTTTAtt                                                                                                                           <  2:166713/141‑1 (MQ=255)
          aaaGCCTGCGATTCATTTATGTAATCGCAGGCTTTTTTATTGCTAACTTTTTAATTTTATATATTTAGACATACGGTATTCTATTTATTCTTTTGGTCCAACCCATTCAATAACATCTTCATAAGCTTGGCGTGTTTTAtt                                                                                                                           <  1:43467/141‑1 (MQ=255)
                tGCGATTCATTTATGTAATCGCAGGCTTTTTTATTGCTAACTTTTTAATTTTATATATTTAGACATACGGTATTCTATTTATTCTTTTGGTCCAACCCATTCAATAACATCTTCATAAGCTTGGCGTGTTTTATTTTTCgg                                                                                                                     >  1:254907/1‑141 (MQ=255)
                cgCGATTCATTTATGTAATCGCGGGCTTCTTCATTGCTAACTTTTTAATTTTATCTATTTAGCCATACGGTATTTTATTTATTCTTTTGTTCCCACCCCCTCAATAACATCTTCATACGCTTAGCGCGTCTTATTTTTctt                                                                                                                     >  2:448949/2‑139 (MQ=255)
                      tCATTTATGTAATCGCAGGCTTTTTTATTGCTAACTTTTTAATTTTATATATTTAGACATACGGTATTCTATTTATTCTTTTGGTCCAACCCATTCAATAACATCTTCTTAAGCTTGGCGTGTTTTATTTTTCGGTGGctc                                                                                                               <  2:371252/141‑1 (MQ=255)
                        cataaatgaaTCGCAGGCTTTTTTATTGCTAACTTTTTAATTTTATATATTTAGACATACGGTATTCTATTTATTCTTTTGGTCCAACCCATTCAATAACATCTTCATAAGCTTGGCGTGTTTTATTTTTCGGTGGCTCtt                                                                                                             >  1:264639/9‑141 (MQ=255)
                          ttATGTAATCGCAGGCTTTTTTCTTGCTAACTTTTTAATTTTATATATTTAGACATACGGTATTCTATTTATTCTTTTGGTCCAACCCATTCAATATCATCTTCATAAGCTTGGCGTGTTTTATTTTTCGGTGGCTCttgt                                                                                                           <  2:213894/141‑1 (MQ=255)
                          ttATGTAATCGCAGGCTTTTTTATTGCTAACTTTTTAATTTTATATATTTAGACATACGGTATTCTATTTATTCTTTTGGTCCAACCCATTCAATAACATCTTTATAAGCTTGGCGTGTTTTTTTTTTTGGTGGCTCttgt                                                                                                           >  1:300720/1‑141 (MQ=255)
                                 aacgatgttttttttttttCTAACTTTTTAATTTTATATATTTAGACATACGGTATTCTATTTATTCTTTTGGTCCAACCCATTCAATAACATCTTCATAAGCTTGGCGTGTTTTATTTTTCGGTGGCTCTTGTTGTCTGt                                                                                                    >  2:245725/10‑141 (MQ=255)
                                 aTCGCAGGCTTTTTTATTGCTAACTTTTTAATTTTATATATTTAGACATACGGTATTCTATTTATTCTTTTGGTCCAACCCATTCAATAACATCTTCATAAGCTTGGCGTGTATTATTTTTCGGTGGCTCTTGTTGTCTGt                                                                                                    >  1:95114/1‑141 (MQ=255)
                                           tttttATTGCTAACTTTTTAATTTTATATATTTAGACATACGGTATTCTATTTATTCTTTTGGTCCAACCCATTCAATAACATCTTCATAAGCTTGGCGTGTTTTATTTTTCGGTGGCTCTTGTTGTCTGTAGCCAAATGc                                                                                          <  1:298537/141‑1 (MQ=255)
                                           tttttATTGCTAACTTTTTAATTTTATATATTTAGACATACGGTATTCTATTTATTCTTTTGGTCCAACCCATTCAATAACATCTTCATAAGCTTGGCGTGTTTTATTTTTCGGTGGCTCTTGTTGTCTGTAGCCAAATGc                                                                                          <  1:280980/141‑1 (MQ=255)
                                                  tGCTAACTTTTTAATTTTATATATTTAGACATACGGTATTCTATTTATTCTTTTGGTCCAACCCATTCAATAACATCTTCATAAGCTTGGCGTGTTTTATTTTTCGGTGGCTCTTGTTGTCTGTAGCCAAATGCGACCATc                                                                                   <  1:412990/141‑1 (MQ=255)
                                                   gCTAACTTTTTAATTTTATATATTTAGACATACGGTATTCTATTTATTCTTTTGGTCCAACCCATTCAATAACATCTTCATAAGCTTGGCGTGTTTTATTTTTCGGTGGCTCTTGTTGTCTGTAGCCAAATGCGACCATCa                                                                                  >  1:175138/1‑141 (MQ=255)
                                                    cTAACTTTTTAATTTTATATATTTAGACATACGGTATTCTATTTATTCTTTTGGTCCAACCCATTCAATAACATCTTCATAAGCTTGGCGTGTTTTATTTTTCGGTGGCTCTTGTTGTCTGTAGCCAAATGCGACCATCAc                                                                                 <  2:110997/141‑1 (MQ=255)
                                                       aCTTTTTAATTTTATATATTTAGACATACTGTATTCTATTTATTCTTTTTTTCCAAGCCAGTTAATAACATCTTCATATGCTTGGCGTGTTTTATTTTTCGGTGGATCTTGTTGTCTGTAGCCAAATGCAACCATCACTgg                                                                              <  1:268383/141‑2 (MQ=255)
                                                         tttttAATTTTATATATTTAGACATACGGTATTCTATTTATTCTTTTGGTCCAACCCATTCAATAACATCTTCATAAGCTTGGCGTGTTTTATTTTTCGGTGGCTCTTGTTGTCTGTAGCCAAATGCGACCATCACTGATa                                                                            <  2:268343/141‑1 (MQ=255)
                                                          ttttAATTTTATATATTTAGACATACGGTATTCTATTTATTCTTTTGGTCCAACCCATTCAATAACATCTTCATAAGCTTGGCGTGTTTTATTTTTCGGTGGCTCTTGTTGTCTGTAGCCAAATGCGACCATCACTGATaa                                                                           >  1:96204/1‑141 (MQ=255)
                                                            ttAATTTTATATATTTAGACATACGGTATTCTATTTATTCTTTTGGTCCAACCCATTCAATAACATCTTCATAAGCTTGGCGTGTTTTATTTTTCGGTGGCTCTTGTTGTCTGTAGCCAAATGCGACCATCACTGATAAAc                                                                         <  1:179183/141‑1 (MQ=255)
                                                            ttAATTTTATATATTTAGACATACGGTATTCTATTTATTCTTTTGGTCCAACCCATTCAATAACATCTTCATAAGCTTGGCGTGTTTTATTTTTCGGTGGCTCTTGTTGTCTGTAGCCAAATGCGACCATCACTGATAAAc                                                                         <  1:337250/141‑1 (MQ=255)
                                                             tAATTTTATATATTTAGACATACGGTATTCTATTTATTCTTTTGGTCCAACCCATTCAATAACATCTTCATAAGCTTGGCGTGTTTTATTTTTCGGTGGCTCTTGTTGTCTGTAGCCAAATGCGACCATCACTGATAAccc                                                                        >  2:323337/1‑141 (MQ=255)
                                                              aaTTTTATATATTTAGACATACGGTATTCTATTTATTCTTTTGGTCCAACCCATTCAATAACATCTTCATAAGCTTGGCGTGTTTTATTTTTCGGTGGCTCTTGTTGTCTGTAGCCAAATGCGACCATCACTGATAAACCa                                                                       <  1:156142/141‑1 (MQ=255)
                                                                     tatatTTAGACATACGGTATTCTATTTATTCTTTTGGTCCAACCCATTCAATAACATCTTCATAAGCTTGGCGTGTTTTATTTTTCGGTGGCTCTTGTTGTCTGTAGCCAAATGCGACCATCACTGATAAACCAAATTGCt                                                                >  2:14725/1‑141 (MQ=255)
                                                                         tttAGACATACGGTATTCTATTTATTCTTTTGGTCCATCCCATTCAATAACATCTTCATAAGCTTGGCGTGTTTTATTTTTCGGTGGCTCTTGTTGTCTGTAGCCAAATGCGACCATCACTGATAAACCAAATTGCTCAGt                                                            <  2:246164/141‑1 (MQ=255)
                                                                          ttAGACATACGGTATTCTATTTATTCTTTTGGTCCAACCCATTCAATAACATCTTCATAAGCTTGGCGTGTTTTATTTTTCGGTGGCTATTGTTGTCTGTAGCCAAATGCGACCATCACTGATAAACCAAATTGCTCAGTa                                                           >  1:153487/1‑141 (MQ=255)
                                                                                 tACGGTATTCTATTTATTCTTTTGGTCCAACCCATTCAATAACATCTTCATAAGCTTGGCGTGTTTTATTTTTCGGTGGCTCTTGTTGTCTGTAGCCAAATGCGACCATCACTGATAAACCAAATTGCTCAGTATCTAAGa                                                    <  2:257334/141‑1 (MQ=255)
                                                                                   cGGTATTCTCTTTATTCTTTTGGTCCAACCCATTCAATAACATCTTCATATGCTTGGCGTGTTTTATTTTTCGGTGGCTCTTGTTGTCTGTAGCCAAATGCGACCATTACTGATAAACCAAATTGCTCAGTATCTAAGATc                                                  >  2:159949/1‑141 (MQ=255)
                                                                                         tCTATTTATTCTTTTGGTCCAACCCATTCAATAACATCTTCATAAGCTTGGCGTGTTTTATTTTTCGGTGGCTCTTGTTGTCTGTAGCCAAATGCGACCATCACTGATAAACCAAATTGCTCAGTATCTAAGATCCCTTTa                                            >  2:187839/1‑141 (MQ=255)
                                                                                               tattCTTTTGGTCCAACCCATTCAATAACATCTTCATAAGCTTGGCGTGTTTTATTTTTCGGTGGCTCTTGTTGTCTGTAGCCAAATGCGACCATCACTGATAAACCAAATTGCTCAGTATCTAAGATCCCTTTATTTGCt                                      <  1:147393/141‑1 (MQ=255)
                                                                                                    ttttGGTCCAACCCATTCAATAACATCTTCATAAGCTTGGCGTGTTTTATTTTTCGGTGGCTCTTGTTGTCTGTAGCCAAATGCGACCATCACTGATAAACCAAATTGCTCAGTATCTAAGATCCCTTTATTTGCTAAAAt                                 <  2:148017/141‑1 (MQ=255)
                                                                                                       tGGTCCAACCCATTCAATAACATCTTCATAAGCTTGGCGTGTTTTATTTTTCGGTGGCTCTTGTTGTCTGTAGCCAAATGCGACCATCACTGATAAACCAAATTGCTCAGTATCTAAGATCCCTTTATTTGCtctaaatat                              >  2:216409/1‑133 (MQ=255)
                                                                                                       tGGTCCAACCCATTCAATAACATCTTCATAAGCTTGGCGTGTTTTATTTTTCGGTGGCTCTTGTTGTCTGTAGCCAAATGCGACCATCACTGATAAACCAAATTGCTCAGTATCTAAGATCCCTTTATTTGCTAAAAtgtc                              >  1:258634/1‑141 (MQ=255)
                                                                                                           ccAACCCATTCAATAACATCTTCATAAGCTTGGCGTGTTTTATTTGTCGGTGGCTCTTGTTGTCTGTAGCCAAATGCGACCATCACTGATAAACCAAATTGCTCAGTATCTAAGATCCCTTTATTTGCTAAAAtgtctgtc                          <  2:390347/141‑1 (MQ=255)
                                                                                                              aCCCATTCAATAACATCTTCATAAGCTTGGCGTGTTTTATTTTTCGGTGGCTCTTGTTGTCTGTAGCCAAATGCGACCATCACTGATAAACCAAATTGCTCAGTATCTAAGATCCCTTTATTTGCTAAAATGTCTGTCACt                       <  2:384231/141‑1 (MQ=255)
                                                                                                                  aTTCAATAACATCTTCATAAGCTTGGCGTGTTTTATTTTTCGGTGGCTCTTGTTGTCTGTAGCCAAATGCGACCATCACTGATAAACCAAATTGCTCAGTATCTAAGATCCCTTTATTTGCTAAAATGTCTGTCACTGTAt                   >  2:344358/1‑141 (MQ=255)
                                                                                                                   ttCAATAACATCTTCATATGCTTGGCGTGTTTTATTTTTCGGTGGATCTTGTTGTCTGTAGCCAAATGCAACCATCACTGATAAACCAAATTGCTCAGTATCTAAGATACCTTTATTTGCTAAAATGTCTGTCACTGTATc                  >  1:62304/1‑141 (MQ=255)
                                                                                                                        tAACATCTTCATAAGCTTGGCGTGTTTTATTTTTCGGTGGCTCTTGTTGTCTGTAGCCAAATGCGACCATCACTGATAAACCAAATTGCTCAGTATCTAAGATCCCTTTATTTGCTAAATTGTCTGTCACTGTATCCAGAc             <  1:449682/141‑1 (MQ=255)
                                                                                                                        tAACATCTTCATAAGCTTGGCGTGTTTTATTTTTCGGTGGCTCTTGTTGTCTGTAGCCAAATGCGACCATCACTGATAAACCAAATTGCTCAGTATCTAAGATCCCTTTATTTGCTAAAATGTCTGTCACTGTATCCAGAc             <  1:364779/141‑1 (MQ=255)
                                                                                                                            aTCTTCATATGCTTGGCGTGTTTTATTTTTCGGTGGATCTTGTTGTCTGTAGCCAAATGCAACCATCACTGATAAACCAAATTGCTCAGTATCTAAGATACCTTTATTTGCTAACATGTCTGTCACTGTATCCAGACTaaa         <  2:219842/141‑1 (MQ=255)
                                                                                                                              cTTCATAAGCTTGGCGTGTTTTATTTTTCGGTGGCTCTTGTTGTCTGTAGCCAAATGCGACCATCACTGATAAACCAAATTGCTCAGTATCTAAGATCCCTTTATTTGCTAAAATGTCTGTCACTGTATCCAGACTAAAAc       <  2:215688/141‑1 (MQ=255)
                                                                                                                              cTTCATAAGCTTGGCGTGTTTTATTTTTCGGTGGCTCTTGTTGTCTGTAGCCAAATGCGACCATCACTGATAAACCAAATTGCTCAGTATCTAAGATCCCTTTATTTGCTAAAATGTCTGTCACTGTATCCAGACTAAAAc       <  2:90053/141‑1 (MQ=255)
                                                                                                                                tCATAAGCTTGGCGTGTTTTATTTTTCGGTGGCTCTTGTTGTCTGTAGCCAAATGCGACCATCACTGATAAACCAAATTGCTCAGTATCTAAGATCCCTTTATTTGCTAAAATGTCTGTCACTGTATCCAGACTAAAACCt     <  1:89660/141‑1 (MQ=255)
                                                                                                                                   tAAGCTTGGCGTGTTTTATTTTTCGGTGGCTCTTGTTGTCTGTAGCCAAATGCGACCATCACTGATAAACCAAATTGCTCAGTATCTAAGATCCCTTTATTTGCNAAAATGTCTGTCACTGTATCCAGACTAAAACCTTcc  <  1:86083/141‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                          
TTTTCACGCAAAAGCCTGCGATTCATTTATGTAATCGCAGGCTTTTTTATTGCTAACTTTTTAATTTTATATATTTAGACATACGGTATTCTATTTATTCTTTTGGTCCAACCCATTCAATAACATCTTCATAAGCTTGGCGTGTTTTATTTTTCGGTGGCTCTTGTTGTCTGTAGCCAAATGCGACCATCACTGATAAACCAAATTGCTCAGTATCTAAGATCCCTTTATTTGCTAAAATGTCTGTCACTGTATCCAGACTAAAACCTTCC  >  CP000730/2658531‑2658802

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑