Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,368,101 T→C 100% A21A (GCT→GCC arnA → fused UDP‑4‑amino‑4‑deoxy‑L‑arabinose formyltransferase/UDP‑glucuronate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,368,1010TC96.3% 99.1 / ‑5.2 51A21A (GCT→GCCarnAfused UDP‑4‑amino‑4‑deoxy‑L‑arabinose formyltransferase/UDP‑glucuronate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (0/3);  new base C (27/19);  total (28/23)
Fisher's exact test for biased strand distribution p-value = 8.36e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TCTACACCGATGTCCGCGCCCGCCCCCGCTATTTTGTTCAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTGGCGGCGGAAAGAGGCATTCCGGTTTATGCGCCGGATAACGTTAATCATCCGCTGTGGGTG  >  NC_000913/2367963‑2368242
                                                                                                                                          |                                                                                                                                             
tCTACACCGATGTCCGCGCCCGCCCCCGCTATTTTGTTCAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTCCCACGATATGCGCTGCCTCGGTATCGAAGCCCCGCTGGccgccg                                                                                                                                           >  2:270114/1‑143 (MQ=255)
    cacCGATGTCCGCGCCCGCCCCCGCTATTTTGTTCAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAATACTGTTGTTTTTGCCTAACAAGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTa                                                                                                                                       <  1:1059226/143‑1 (MQ=255)
     acCGATGTCCGCGCCCGCCCCCGCTATTTTGTTCAGCAAGTTATCCGTCCATCCAGCACGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTag                                                                                                                                      >  2:1014907/1‑142 (MQ=255)
     acCGATGTCCGCGCCCGCCCCCGCTATTTTGTTCAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACTGTGGTTGTCCCCAACCACGATATGGGCTGCCTCGGTATTGTCGCCCTGCTGGCCGCCCGtttc                                                                                                                                      >  1:519839/1‑141 (MQ=255)
     acCGATGTCCGCGCCCGCCCCCGCTATTTTGTTCAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGACGCCCTGCTTGCCGCCGGTTAc                                                                                                                                      >  1:857921/1‑143 (MQ=255)
       cGATGTCCGCGCCCGCCCCCGCTATTTTGTTCAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCATACCACGATATGGGATGCCTCGGTATTGAAGCCCTACTGCCCGCCGCTTACgg                                                                                                                                    >  2:162711/1‑142 (MQ=255)
       cGATGTCCGCGCCCGCCCCCGCTATTTTGTTAAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATTGGCTGCTTCGGTATTCAACACATGCTGGCGCCCAGTTActa                                                                                                                                    >  2:556460/1‑141 (MQ=255)
            tCCGCGCCCGCCCCCGCTATTTTGTTCAGCAAGTTATCAGTCCATGCAGCAAGGAAAATGAATAGTGAAAACAGTCGTTTTTGAGTACCACGATAGGGGATGGCTCGGAAATGAAGCCTGGGCGGCCGCCGGTTGCGGAAtga                                                                                                                               >  2:1030865/1‑141 (MQ=255)
             ccGCGCCCGCCCCCGCTATTTTGTTCAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGATTACTGTTGTTTTTGTTTATTTTGATATGGGATGACTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAg                                                                                                                              <  1:1020163/143‑1 (MQ=255)
              cgtgcccgcctCGGCTATTGGGTGTAGGAGGTTAGCCGTCCATCCAGCAAGGAAGATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAgc                                                                                                                             <  2:771747/140‑1 (MQ=255)
                   ccgcccCCGCTATTTTGTTCAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGCGTCCCTGCTGGCCGCCGTCTGCGTACTTAGCGCCAt                                                                                                                        >  1:837036/1‑143 (MQ=255)
                      ccccccGTATTTTGTTCGGGAGGTTTGCCGGCCAGCCAGCAAGGAAAATGGATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCGGTTACGAAATTAGCGCCAtttt                                                                                                                     <  2:1078757/143‑1 (MQ=255)
                        cccGCTATTTTGTTCAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACTGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTCACGAGATTAGCGCCATTCTTa                                                                                                                   >  1:1073314/1‑143 (MQ=255)
                        cccGCTATTTTGTTCAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCACAATTTTTa                                                                                                                   >  1:557919/1‑143 (MQ=255)
                                        gCAAGTTATCCGTACATCCAGCAAGGAAAATGAATAATGAAAACCGTCGGTTTTGCCTACCACGATATGGGATGCCTAGGTATTGAAGCCCTGCCGAACATCGCGTACGAAATCAGCGCCATTTTTACGCATACTGATAATcc                                                                                                   >  2:497374/1‑143 (MQ=255)
                                           tgTTCTCCGTCCTTCCAGCAAGGAAACTGTATAATGACAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCgg                                                                                                <  2:164252/142‑1 (MQ=255)
                                           atttATCCGTCCATCCAGCAAAGAAAATGATTAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCgg                                                                                                <  2:1072268/141‑1 (MQ=255)
                                           aGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTAAGAAATTAGCGCCATTTTTACCCATACTGATAAGCGCgg                                                                                                >  2:46050/1‑143 (MQ=255)
                                           aGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGATACGAAATTAGCGCCCTTTTTACCCATACTGATAATCCCgg                                                                                                >  2:144113/1‑143 (MQ=255)
                                            gTTATCAGTCTATCCAGCAAGGAAAATGAATAATGAAATCCGTCTTTTTTGTCTACCACGATATGTGATGTTTCGGTATTTAAGCTCTGCTGGCAGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGt                                                                                               <  1:640201/143‑1 (MQ=255)
                                                  cGTCCATCCAGCAAGGAAAATGAATAATGACAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGaaaaa                                                                                         <  2:746028/143‑1 (MQ=255)
                                                       atccaGCAAGGAAAATGGATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAATCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCtt                                                                                    <  2:409388/143‑1 (MQ=255)
                                                           aGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTTGTGCCACTGTTTCTTATCCTGATAATCCCGATGAAAAAGCCATGTAt                                                                                >  1:516275/1‑143 (MQ=255)
                                                            gCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTAACACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGACCGACGGGTACGAAATTAGCGCCGTTTTTACCCATACTGATAATGCAGGTGAAAAAGCCtttaacg                                                                               >  2:1075526/1‑139 (MQ=255)
                                                                 gAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCg                                                                          >  1:531197/1‑143 (MQ=255)
                                                                  aaaaTGAATGATGCAATCCGTCGTTTGGGCCTCCCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCAATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCgg                                                                         <  1:685885/143‑1 (MQ=255)
                                                                     aTGAATAATGAAAACCGTCGTTTTTGCATACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTTTTTTCTTTTTACTTATTTAGTGCCATTTTTACCAAAACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTgg                                                                      <  1:434828/143‑1 (MQ=255)
                                                                       gAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAACATTTTATGGTTCGGTGGcc                                                                    >  2:683911/1‑142 (MQ=255)
                                                                              gaccaCCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTTGTGCCATTTTTACCCATTCTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTgg                                                             <  1:598434/139‑1 (MQ=255)
                                                                               taaaTCGTCGGTTTTGCCGACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCGGTTATGAAATTAGCGTTTTTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTggc                                                            <  1:658946/142‑1 (MQ=255)
                                                                               aaaaCCGTCGTTTTTGCCTACCACGATATGGGACGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGATATTAGCGCTATTTTTACCCAGTCTGATAAGCCCGGGGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTggc                                                            <  1:736111/143‑1 (MQ=255)
                                                                                  gtcgtcttTTTTGTCTACCACGATACGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAAGTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTggcggc                                                         <  1:56361/141‑1 (MQ=255)
                                                                                   ccGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCATTTATGGTACGGTGGCTCGTCTGgcggcg                                                        >  1:755348/1‑143 (MQ=255)
                                                                                     gtcgtTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTGGCGGCGGa                                                      >  2:6170/1‑143 (MQ=255)
                                                                                          tttgGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCTATTTTTACCCATACTTATAATCTTGGTGAAAAAGCAAATTATGGTTCGGTGGCTCGTCTGGCGGCGGAAagag                                                 <  1:582489/143‑1 (MQ=255)
                                                                                            ttGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCAGCGTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATAGTTCGGTGGCTCGTCTGGCGACGGAAagagac                                               >  2:826161/1‑141 (MQ=255)
                                                                                                  aCCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTGGCGGCGGAAAGAGGCATTCCg                                         <  2:521679/143‑1 (MQ=255)
                                                                                                      cGATATGGGATGCCTCGGTATTGAAGCCCGGCGAGCCGCCGGTGACGAAATTGGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAGAGCCTTTTATGGTTCGGTGGCTCGTCGGGCGGAGGAAAGAGGCATGCAGGttt                                     >  2:839349/1‑143 (MQ=255)
                                                                                                          atGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAACGACTTTTATGGTTCGGTGGTTCGTCTGGCGGCGGCCACAGGCATTCCGGTTTCTgc                                 >  2:736878/1‑143 (MQ=255)
                                                                                                           tGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTGGCGGCGGAAAGAGGCATTCCGGTTTATgcg                                >  1:1011106/1‑143 (MQ=255)
                                                                                                            gggATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTGGCGGCGGAAAGAGGCATTCCGGGTTATgagc                               >  1:342339/1‑143 (MQ=255)
                                                                                                              gATGCCTCGGTAGTTTAGCCTCGCTGTCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTGTTATGGTTCGGTGGCTCGTCTGGCGGCTGAAAGAGGCATTCCGGTTTATGCGCCg                             <  2:90800/143‑1 (MQ=255)
                                                                                                                   cTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCTGTGAAAAAGTCTTTTATGGTTGGGTGGCTCGTCTGGCGGCGGAAAGAGGCATTCCGGTTTATGCGCCGGATaa                        <  1:640213/143‑1 (MQ=255)
                                                                                                                     cGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTGGCGGCGGAAAGAGGCATTCCGGGTTATGCGGCGGATACCg                      >  2:614510/1‑143 (MQ=255)
                                                                                                                        tATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTGGCGGCGGAAAGAGGCATTGCGGGTTAGGCGCCGGAGCACGtgc                   >  2:442674/1‑141 (MQ=255)
                                                                                                                           tGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTGACGACGGAAAGAAACATTCCGGGTTATGCGCCGGATGAGGTTaagc                >  2:847730/1‑141 (MQ=255)
                                                                                                                            gAAGCCCTGCTGGCCGCCGGTTACGAAAATAGCGCCATTTTTACCCATACTGATAATCCCGGTGGAAAAGCCCTTTATGGTTCGGTGGCTCGTCTGGCGGCGCAAAGAGTCATTCCGGTTTATGCGCCGGATAACGATAaaca               >  1:848976/1‑143 (MQ=255)
                                                                                                                             aatcccTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTGGCGGCGGAAAGAGGCATTCCGGTTTATGCGCCGGATAACGTTAatcat              <  2:1011106/140‑1 (MQ=255)
                                                                                                                             aaccccTGCTGGCCGCCCGTGACTAAATTAGCGCCATTTTTACCCATACTGATATTCCCGGTGAAAAAGCCTTTTCTGGTTCGGTGGCTCGTCTGGCGGCGGAAAGAGGCATTCCGGTTTATGCGCCGGATAACGTTAatcat              <  2:1115487/140‑1 (MQ=255)
                                                                                                                             aaGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTGGCGGCGGAAAGAGGCATTCCGGTTTATGCGCCGGATAACGTTAatcat              >  1:441482/1‑143 (MQ=255)
                                                                                                                                 ccTGCTGGCCGCCTGTTACGAAATTAGCGCCTTTTTTACCCCTACTGATAACCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTGGCGGCGGAAAGAGGCATTCCGGTTTATGCGCCGGATAACGTTAATCATCCGc          <  2:1027778/143‑1 (MQ=255)
                                                                                                                                  ctggtgGCCGCCGGTTACGAAATTAGGGCCATTTTTTCCCAGACTGATAATCCCGATCACAAAGACGTCTATGGTTCGGTGGCTCGTCTGGCGGCGGAAAGAGGCATTCCGGTTTATGCGCCGGATAACGTTAATCATCCGCt         <  2:237359/143‑1 (MQ=255)
                                                                                                                                    gctggcagCCGGTTACTGAAGTAGCGCCATCTTTACTCATACTGATAATCCCGGTGCTATAGCCTTTTCTCAACCTGTCGCTTGTATTGCATCGGAAAGAGGCATTCCGGTCTATGCGCCGGATAACGTTAATCATCCGCtgt       <  2:1104329/143‑1 (MQ=255)
                                                                                                                                     ctggccgCCGGTGACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTTTTTTGTCTGGCGGCGGAAAGAGGCATTCCGGTTGAGGCGCCGGATAACGTTAATCATCCGCtgtg      <  1:12977/143‑1 (MQ=255)
                                                                                                                                         ctgCCAGTGACGAAATTTGCGCCATTCTTCCCCATAGTGCCAAGCCCGTTGAAAAAGCCTGTTATGGTTCGGGGGCTCGTGTGGCGGCGGAAAGAGGCATTCCGGTTTATGCGCCGGATAACGTTAATCATCCGCTgtgggtg  <  2:335018/143‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                             
TCTACACCGATGTCCGCGCCCGCCCCCGCTATTTTGTTCAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTGGCGGCGGAAAGAGGCATTCCGGTTTATGCGCCGGATAACGTTAATCATCCGCTGTGGGTG  >  NC_000913/2367963‑2368242

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