Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,340,081 A→G 20.7% intergenic (+66/‑265) ribAB‑II → / → PP_3814 3,4‑dihydroxy‑2‑butanone‑4‑phosphate synthase/polyamine ABC transportersubstrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,340,0810AG20.7% 97.9 / 15.0 49intergenic (+66/‑265)ribAB‑II/PP_38143,4‑dihydroxy‑2‑butanone‑4‑phosphate synthase/polyamine ABC transportersubstrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (14/24);  new base G (7/3);  total (22/27)
Fisher's exact test for biased strand distribution p-value = 8.05e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑TTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGCAGGC  >  NC_002947/4339938‑4340228
                                                                                                                                                |                                                                                                                                                   
cTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑TTTACCCGCGAATGCGAGGGTGAcacca                                                                                                                                                <  2:349266/149‑1 (MQ=255)
   aaaTACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGCGGCCCTGTAGGAGCGGG‑TTTACCCGCGAATGCGAGGGTGACaccacc                                                                                                                                              <  3:316537/148‑1 (MQ=255)
     aTACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑TTTTCCCGCGCATGGGAGGGTGCCCCCCCCggc                                                                                                                                           >  8:166345/1‑147 (MQ=255)
     aTACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑TTTACCCGCGAATGCGAGGGTGACACCACCGAc                                                                                                                                           <  3:356526/149‑1 (MQ=255)
             ggCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑TTTACCCGCGAATGCGAGTGTGACACCACCGACGCATTcgc                                                                                                                                   <  7:233900/149‑1 (MQ=255)
             ggCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑TTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTcgc                                                                                                                                   <  6:174514/149‑1 (MQ=255)
             ggCTGGCGGCCTATGAGCTGGAGGTGCTGCAGAGCAACCCGTTCGACCCGGGTATGGATGGCTCACAACACACAACCAATTTATTTAGCGGGGCCCTGTAGGAGCGGG‑TTTACCCGCGAATGCGATGGTGACACCACCGACGCATTcgc                                                                                                                                   <  7:232037/149‑1 (MQ=38)
              gCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGAATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑TTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTcgc                                                                                                                                   >  2:219851/1‑148 (MQ=255)
               cTGGCGGGCTATGAGCTGGAGGTGGTGCATAGCATCCCGTTCGTCCCGGGGTTGGATGGCTGAGAACACACAAGGAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑TTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCgg                                                                                                                                 <  8:283595/149‑1 (MQ=255)
               cTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGAGCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑TTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCgg                                                                                                                                 <  7:330574/149‑1 (MQ=25)
               cTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑TTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCgg                                                                                                                                 <  7:98182/149‑1 (MQ=255)
                      gCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATCGCTCAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑TTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAAcc                                                                                                                          <  5:24870/149‑1 (MQ=255)
                        tATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑TTTACCCGCGAATGCGAGGGTGACACCACCGCCGCATTCGCCGGTAAACCCg                                                                                                                        >  1:45492/1‑149 (MQ=255)
                          tGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑TTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCt                                                                                                                      <  3:91952/149‑1 (MQ=255)
                                ggAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑TTTACCCGCGAATGCGAGGGTGACACCCCCGGCGCCTTCGCGGGTAAACCCGCTCCTACa                                                                                                                >  5:177471/1‑149 (MQ=255)
                                 gAGGTGGGCCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTCAGAACAAACAAGCAATTTATTTCGCGGGGCCCTGTAGGAGCGGG‑TTTACCCGCTAATGCGAGGGTTACACCACCGACGCATTCGCGGTTAAACCCGCTCCTACa                                                                                                                <  2:55793/148‑1 (MQ=37)
                                                   ccGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑TTTAACCGCGAATGCGAGGGTGAAACCCCCGACGCATTAGCGGGTAAAACCCCAGCTACAGGGCCCTTGGTGCTTTTGa                                                                                             >  5:332120/1‑149 (MQ=38)
                                                     gTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑TTTAACCGCGACTGCGCGGGTGACACCACCGACGCATTAGCGGGGAAACCCGCGCCGACAGGGCCCGTGGTGGGTTTGa                                                                                             >  6:73418/1‑147 (MQ=37)
                                                            ccGGGGATGGATGGCTTAGAACACACAAGAAATTGTTTTAGCGGGGCCCTGTAGGAGCGGG‑TTTACCCGCTAATGCGAGGGTGACACCACCGACGCATTCGCGGTTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGa                                                                                    <  2:45492/149‑1 (MQ=38)
                                                                   tggatggCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑TTTAACCGCGAATGCGAGGGTGGCACCACCGACGCAATAGCGGGTAAACCCGCTCCTACAGGGCCCGTGGTGGTTTTGACACGGCCGAGTGAttg                                                                             >  5:107971/1‑147 (MQ=37)
                                                                     gatggCTGAGAACACACACGCAATTTATTTCGCGGGGCCCTGTAGGAGCGGG‑TTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGcc                                                                           <  5:73418/149‑1 (MQ=255)
                                                                          cTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTTACCCCCAAAT‑CGAGGGGGGCCCCACCCCCGCCTTAGCGGGGAAACCCCCGCCGACAAGGGCCGTGGGGGTTTTTAGAGGGGGGGGTGGTGGGCggg                                                                        >  8:182590/1‑146 (MQ=255)
                                                                          cTGACAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGCGCGGG‑TTTACCCGCTAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGt                                                                        <  7:182590/147‑1 (MQ=255)
                                                                            gataaCACACATGCATTTTATTTCGGGGGGCCGTGTTGGCGGGGG‑TTTACCGGCTATTGCTTGGGTTACACCACCGACGCTTTCGCGGTTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGc                                                                     <  6:332120/145‑1 (MQ=255)
                                                                                  acacaAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑TTTACCCGTGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAAcc                                                                                                                          >  8:319508/1‑89 (MQ=255)
                                                                                  acacaAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑TTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAAcc                                                                                                                          <  7:319508/89‑1 (MQ=255)
                                                                                    acaAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑TTTAACCGCGCATGGGAGGGTGGCACCCCCGCCGCATTAGCGGGGAAACCCCCGGCGACAGGGCCCGTGGTGGTTTTGACGGGGCGGAGTGATGGCCGGTGGCCCCCTaac                                                             >  4:111120/1‑146 (MQ=255)
                                                                                      aaGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑TTTAACCGCGAATGCGAGGGTGACACCCCCGACGCATTAGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACctct                                                          >  5:173908/1‑149 (MQ=255)
                                                                                                  accaagggCCCTGTAGGAGCGGG‑TTTAACCGCGAATGGGTCGGTGGTGTGACCCACGCCTTAGCGGGTGAACCCGCTCCGACAGGGCCCGTGGTGGTTTTGACAGGGCGGGGTGATGGGCGGTCGCCCCCAACCTCTTGCCAAAATAtt                                              >  2:262928/6‑149 (MQ=9)
                                                                                                   aGCGGGGCCCTGTAGGCGGGTG‑TTTACCCGCTAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCtt                                              <  7:283107/148‑1 (MQ=255)
                                                                                                    ggggggCCCTGTTGGTGCGGG‑TTTACCCGGTAATGGGAGGGTGACACCACCGACGAATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTg                                             <  6:107971/146‑1 (MQ=38)
                                                                                                     cGGGGCCGTTTAGGAGGGGT‑TTTACCCGCTAATGGGTGGGTTACCCCACCGACCAATTCGGGTGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGa                                           <  1:120607/149‑1 (MQ=255)
                                                                                                          ccgtGTAGGCGCGGG‑TTTACCCGCTAAGGCGTGGGTGACACCACCGACGCATTCGCGGTTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATAcc                                      <  3:123612/146‑1 (MQ=38)
                                                                                                                aGGAGCGGG‑TTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTAGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGc                                                                            >  8:61259/1‑105 (MQ=255)
                                                                                                                 ggAGCGGG‑TTAACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGTTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACGGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCataata                                <  3:111120/148‑1 (MQ=255)
                                                                                                                               ccGTGAATGCGTGGGTGACACCACCGACGCTTTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGa                  <  6:177471/149‑1 (MQ=255)
                                                                                                                               ccGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTg                                             >  6:374933/1‑122 (MQ=21)
                                                                                                                               ccGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTg                                             <  5:374933/122‑1 (MQ=21)
                                                                                                                               ccGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGa                  >  5:147229/1‑149 (MQ=255)
                                                                                                                                cgcgAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCCCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAAAACCATAATATGATAGACGGTAGAc                 >  1:128626/1‑149 (MQ=21)
                                                                                                                                cgcgAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGAAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAc                 <  1:262928/149‑1 (MQ=25)
                                                                                                                                 gcgAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAcc                >  5:228545/1‑149 (MQ=25)
                                                                                                                                 gcgAATGCGAGGGTGACACCACCGATGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAcc                >  1:104921/1‑149 (MQ=255)
                                                                                                                                 gcgAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCGTTTGCCAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATAAGATATCCCGTAGccc                >  6:11838/1‑149 (MQ=255)
                                                                                                                                  cgAATGCGAGGGTGACACCACCGATGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCt               <  2:104921/149‑1 (MQ=255)
                                                                                                                                       gCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATAc                                       >  2:95595/1‑120 (MQ=255)
                                                                                                                                       gCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATAc                                       <  1:95595/120‑1 (MQ=255)
                                                                                                                                       gCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGt          >  5:5961/1‑149 (MQ=255)
                                                                                                                                               gACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGCAGGc  >  4:57870/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
CTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑TTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGCAGGC  >  NC_002947/4339938‑4340228

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