Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 932,884 T→C V97A (GTA→GCA)  lrp → leucine‑responsive global transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913932,8840TC97.9% 170.5 / ‑4.1 48V97A (GTA→GCA) lrpleucine‑responsive global transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (26/21);  total (26/22)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

GTGCGTCGGCTGGAAAGACAAGGGTTTATTCAGGGCTATACGGCGCTGCTTAACCCCCATTATCTGGATGCATCACTTCTGGTATTCGTTGAGATTACTCTGAATCGTGGCGCACCGGATGTGTTTGAACAATTCAATACCGCTGTACAAAAACTTGAAGAAATTCAGGAGTGTCATTTAGTATCCGGTGATTTCGACTACCTGTTGAAAACACGCGTGCCGGATATGTCAGCCTACCGTAAGTTGCTGGGGGAAACCCTGCTGCGTCTGCCTGGCGTCAATGACAC  >  NC_000913/932739‑933025
                                                                                                                                                 |                                                                                                                                             
gtgCGTCGGCTGGAAAGACAAGGGTTTATTCAGGGCTATACGGCGCTGCTTAACCCCCATTATCTGGATGCATCACTTCTGGTATTCGTTGAGATTACTCTGAATCGTGGCGCACCGGATGTGTTTGAACAATTCAATACCGCTGCACaaa                                                                                                                                          <  1:179449/151‑1 (MQ=255)
  gCGTCGGCTGGAAAGACAAGGGTTTATTCAGGGCTATACGGCGCTGCTTAACCCCCATTATCTGGATGCATCACTTCTGGTATTCGTTGAGATTACTCTGAATCGTGGCGCACCGGATGTGTTTGAACAATTCAATACCGCTGTACaaaa                                                                                                                                         <  2:88184/150‑1 (MQ=255)
  gCGTCGGCTGGAAAGACAAGGGTTTATTCAGGGCTATACGGCGCTGCTTAACCCCCATTATCTGGATGCATCACTTCTGGTATTCGTTGAGATTACTCTGAATCGTGGCGCACCGGATGTGTTTGAACAATTCAATACCGCTGCACaaaa                                                                                                                                         >  2:370956/1‑150 (MQ=255)
    gTCGGCTGGAAAGACAAGGGTTTATTCAGGGCTATACGGCGCTGCTTAACCCCCATTATCTGGATGCATCACTTCTGGTATTCGTTGAGATTACTCTGAATCGTGGCGCACCGGATGTGTTTGAACAATTCAATACCGCTGCACAAAAACt                                                                                                                                      >  1:414924/1‑151 (MQ=255)
    gTCGGCTGGAAAGACAAGGGTTTATTCAGGGCTATACGGCGCTGCTTAACCCCCATTATCTGGATGCATCACTTCTGGTATTCGTTGAGATTACTCTGAATCGTGGCGCACCGGATGTGTTTGAACAATTCAATACCGCTGCACAAAAACt                                                                                                                                      >  1:71037/1‑151 (MQ=255)
        gCTGGAAAGACAAGGGTTTATTCAGGGCTATACGGCGCTGCTTAACCCCCATTATCTGGATGCATCACTTCTGGTATTCGTTGAGATTACTCTGAATCGTGGCGCACCGGATGTGTTTGAACAATTCAATACCGCTGCACAAAAACTTga                                                                                                                                   >  1:626041/1‑150 (MQ=255)
           ggAAAGACAAGGGTTTATTCAGGGCTATACGGCGCTGCTTAACCCCCATTATCTGGATGCATCACTTCTGGTATTCGTTGAGATTACTCTGAATCGTGGCGCACCGGATGTGTTTGAACAATTCAATACCGCTGCACAAAAACTTgaagaa                                                                                                                               <  2:503475/151‑1 (MQ=255)
           ggAAAGACAAGGGTTTATTCAGGGCTATACGGCGCTGCTTAACCCCCATTATCTGGATGCATCACTTCTGGTATTCGTTGAGATTACTCTGAATCGTGGCGCACCGGATGTGTTTGAACAATTCAATACCGCTGCACAAAAACTTgaag                                                                                                                                 <  2:112355/149‑1 (MQ=255)
            gAAAGACAAGGGTTTATTCAGGGCTATACGGCGCTGCTTAACCCCCATTATCTGGATGCATCACTTCTGGTATTCGTTGAGATTACTCTGAATCGTGGCGCACCGGATGTGTTTGAACAATTCAATACCGCTGCACAAAAACTTGAAGaaa                                                                                                                              <  1:19521/151‑1 (MQ=255)
                 aCAAGGGTTTATTCAGGGCTATACGGCGCTGCTTAACCCCCATTATCTGGATGCATCACTTCTGGTATTCGTTGAGATTACTCTGAATCGTGGCGCACCGGATGTGTTTGAACAATTCAATACCGCTGCACAAAAACTTGAAGAAATTCAg                                                                                                                         <  2:196840/151‑1 (MQ=255)
                   aaGGGTTTATTCAGGGCTATACGGCGCTGCTTAACCCCCATTATCTGGATGCATCACTTCTGGTATTCGTTGAGATTACTCTGAATCGTGGCGCACCGGATGTGTTTGAACAATTCAATACCGCTGCACAAAAACTTGAAGAAATTCAGGa                                                                                                                       <  1:179649/151‑1 (MQ=255)
                    aGGGTTTATTCAGGGCTATACGGCGCTGCTTAACCCCCATTATCTGGATGCATCACTTCTGGTATTCGTTGAGATTACTCTGAATCGTGGCGCACCGGATGTGTTTGAACAATTCAATACCGCTGCACAAAAACTTGAAGAAATTCAgg                                                                                                                        >  2:559889/1‑149 (MQ=255)
                    aGGGTTTATTCAGGGCTATACGGCGCTGCTTAACCCCCATTATCTGGATGCATCACTTCTGGTATTCGTTGAGATTACTCTGAATCGTGGCGCACCGGATGTGTTTGAACAATTCAATACCGCTGCACAAAAACTTGAAGAAATTCAGGAg                                                                                                                      <  1:597667/151‑1 (MQ=255)
                          tATTCAGGGCTATACGGCGCTGCTTAACCCCCATTATCTGGATGCATCACTTCTGGTATTCGTTGAGATTACTCTGAATCGTGGCGCACCGGATGTGTTTGAACAATTCAATACCGCTGCACAAAAACTTGAAGAAATTCAGGAGTGTCa                                                                                                                 <  1:262320/150‑1 (MQ=255)
                          tATTCAGGGCTATACGGCGCTGCTTAACCCCCATTATCTGGATGCATCACTTCTGGTATTCGTTGAGATTACTCTGAATCGTGGCGCACCGGATGTGTTTGAACAATTCAATACCGCTGCACAAAAACTTGAAGAAATTCAGGAGTGTCAt                                                                                                                >  2:676100/1‑151 (MQ=255)
                              cAGGGCTATACGGCGCTGCTTAACCCCCATTATCTGGATGCATCACTTCTGGTATTCGTTGAGATTACTCTGAATCGTGGCGCACCGGATGTGTTTGAACAATTCAATACCGCTGCACAAAAACTTGAAGAAATTCAGGAGTGTCATTTAg                                                                                                            <  1:264939/151‑1 (MQ=255)
                                gggCTATACGGCGCTGCTTAACCCCCATTATCTGGATGCATCACTTCTGGTATTCGTTGAGATTACTCTGAATCGTGGCGCACCGGATGTGTTTGAACAATTCAATACCGCTGCACAAAAACTTGAAGAAATTCAGGAGTGTCATTTAGt                                                                                                           >  2:544555/1‑150 (MQ=255)
                                gggCTATACGGCGCTGCTTAACCCCCATTATCTGGATGCATCACTTCTGGTATTCGTTGAGATTACTCTGAATCGTGGCGCACCGGATGTGTTTGAACAATTCAATACCGCTGCACAAAAACTTGAAGAAATTCAGGAGTGTCATTTAGTa                                                                                                          <  2:516032/151‑1 (MQ=255)
                                 ggCTATACGGCGCTGCTTAACCCCCATTATCTGGATGCATCACTTCTGGTATTCGTTGAGATTACTCTGAATCGTGGCGCACCGGATGTGTTTGAACAATTCAATACCGCTGCACAAAAACTTGAAGAAATTCAGGAGTGTCATTTAGTAt                                                                                                         <  1:657109/151‑1 (MQ=255)
                                       aCGGCGCTGCTTAACCCCCATTATCTGGATGCATCACTTCTGGTATTCGTTGAGATTACTCTGAATCGTGGCGCACCGGATGTGTTTGAACAATTCAATACCGCTGCACAAAAACTTGAAGAAATTCAGGAGTGTCATTTAGTATCCGGTg                                                                                                   >  1:597666/1‑151 (MQ=255)
                                        cGGCGCTGCTTAACCCCCATTATCTGGATGCATCACTTCTGGTATTCGTTGAGATTACTCTGAATCGTGGCGCACCGGATGTGTTTGAACAATTCAATACCGCTGCACAAAAACTTGAAGAAATTCAGGAGTGTCATTTAGTATCCGGTg                                                                                                   >  2:188893/1‑150 (MQ=255)
                                         ggCGCTGCTTAACCCCCATTATCTGGATGCATCACTTCTGGTATTCGTTGAGATTACTCTGAATCGTGGCGCACCGGATGTGTTTGAACAATTCAATACCGCTGCACAAAAACTTGAAGAAATTCAGGAGTGTCATTTAGTATCCGGTGAt                                                                                                 <  2:503030/151‑1 (MQ=255)
                                                ctTAACCCCCATTATCTGGATGCATCACTTCTGGTATTCGTTGAGATTACTCTGAATCGTGGCGCACCGGATGTGTTTGAACAATTCAATACCGCTGCACAAAAACTTGAAGAAATTCAGGAGTGTCATTTAGTATCCGGTGATTTCGACt                                                                                          >  2:126672/1‑151 (MQ=255)
                                                            tATCTGGATGCATCACTTCTGGTATTCGTTGAGATTACTCTGAATCGTGGCGCACCGGATGTGTTTGAACAATTCAATACCGCTGCACAAAAACTTGAAGAAATTCAGGAGTGTCATTTAGTATCCGGTGATTTCGACTACCTGTTGaaa                                                                               >  2:674739/1‑150 (MQ=255)
                                                            tATCTGGATGCATCACTTCTGGTATTCGTTGAGATTACTCTGAATCGTGGCGCACCGGATGTGTTTGAACAATTCAATACCGCTGCACAAAAACTTGAAGAAATTCAGGAGTGTCATTTAGTATCCGGTGATTTCGACTACCTGTTGaaa                                                                               <  2:212131/150‑1 (MQ=255)
                                                              tCTGGATGCATCACTTCTGGTATTCGTTGAGATTACTCTGAATCGTGGCGCACCGGATGTGTTTGAACAATTCAATACCGCTGCACAAAAACTTGAAGAAATTCAGGAGTGTCATTTAGTATCCGGTGATTTCGACTACCTGTTGAAAac                                                                             >  2:365873/1‑150 (MQ=255)
                                                                 ggATGCATCACTTCTGGTATTCGTTGAGATTACTCTGAATCGTGGCGCACCGGATGTGTTTGAACAATTCAATACCGCTGCACAAAAACTTGAAGAAATTCAGGAGTGTCATTTAGTATCCGGTGATTTCGACTACCTGTTGAAAacac                                                                           <  1:666428/149‑1 (MQ=255)
                                                                 ggATGCATCACTTCTGGTATTCGTTGAGATTACTCTGAATCGTGGCGCACCGGATGTGTTTGAACAATTCAATACCGCTGCACAAAAACTTGAAGAAATTCAGGAGTGTCATTTAGTATCCGGTGATTTCGACTACCTGTTGAAAACAcg                                                                          <  2:38769/150‑1 (MQ=255)
                                                                       aTCACTTCTGGTATTCGTTGAGATTACTCTGAATCGTGGCGCACCGGATGTGTTTGAACAATTCAATACCGCTGCACAAAAACTTGAAGAAATTCAGGAGTGTCATTTAGTATCCGGTGATTTCGACTACCTGTTGAAAACACGCGTGCCg                                                                   >  1:434903/1‑151 (MQ=255)
                                                                         cACTTCTGGTATTCGTTGAGATTACTCTGAATCGTGGCGCACCGGATGTGTTTGAACAATTCAATACCGCTGCACAAAAACTTGAAGAAATTCAGGAGTGTCATTTAGTATCCGGTGATTTCGACTACCTGTTGAAAACACGCGTGCCgg                                                                  <  2:611269/150‑1 (MQ=255)
                                                                          aCTTCTGGTATTCGTTGAGATTACTCTGAATCGTGGCGCACCGGATGTGTTTGAACAATTCAATACCGCTGCACAAAAACTTGAAGAAATTCAGGAGTGTCATTTAGTATCCGGTGATTTCGACTACCTGTTGAAAACACGCGTGCCGGat                                                                <  2:583462/151‑1 (MQ=255)
                                                                           cTTCTGGTATTCGTTGAGATTACTCTGAATCGTGGCGCACCGGATGTGTTTGAACAATTCAATACCGCTGCACAAAAACTTGAAGAAATTCAGGAGTGTCATTTAGTATCCGGTGATTTCGACTACCTGTTGAAAACACGCGTGCCGGat                                                                >  2:43946/1‑150 (MQ=255)
                                                                             tCTGGTATTCGTTGAGATTACTCTGAATCGTGGCGCACCGGATGTGTTTGAACAATTCAATACCGCTGCACAAAAACTTGAAGAAATTCAGGAGTGTCATTTAGTATCCGGTGATTTCGACTACCTGTTGAAAACACGCGTGCCGGatat                                                              <  1:214849/150‑1 (MQ=255)
                                                                               ctgTATTCGTTGAGATTACTCTGAATCGTGGCGCACCGGATGTGTTTGAATATTTAAATACCGCTGCACAAAAACTTGAAGAAATTCAGGAGTGTCATTTAGTATCCGGTGATTTCGACTACCTGTTGAAAACACGTGAGCCGGata                                                               >  2:579713/3‑147 (MQ=255)
                                                                                        ttGAGATTACTCTGAATCGTGGCGCACCGGATGTGTTTGAACAATTCAATACCGCTGCACAAAAACTTGAAGAAATTCAGGAGTGTCATTTAGTATCCGGTGATTTCGACTACCTGTTGAAAACACGCGTGCCGGATATGTCAGCCTACCg                                                  >  1:672290/1‑151 (MQ=255)
                                                                                               tACTCTGAATCGTGGCGCACCGGATGTGTTTGAACAATTCAATACCGCTGCACAAAAACTTGAAGAAATTCAGGAGTGTCATTTAGTATCCGGTGATTTCGACTACCTGTTGAAAACACGCGTGCCGGATATGTCAGCCTACCGTAAGtt                                            <  2:428825/150‑1 (MQ=255)
                                                                                                aCTCTGAATCGTGGCGCACCGGATGTGTTTGAACAATTCAATACCGCTGCACAAAAACTTGAAGAAATTCAGGAGTGTCATTTAGTATCCGGTGATTTCGACTACCTGTTGAAAACACGCGTGCCGGATATGTCAGCCTACCGTAAGTTGc                                          >  2:137240/1‑151 (MQ=255)
                                                                                                  tctGAATCGTGGCGCACCGGATGTGTTTGAACAATTCAATACCGCTGCACAAAAACTTGAAGAAATTCAGGAGTGTCATTTAGTATCCGGTGATTTCGACTACCTGTTGAAAACACGCGTGCCGGATATGTCAGCCTACCGTAAGTTGCTg                                        >  2:564024/1‑151 (MQ=255)
                                                                                                   ctGAATCGTGGCGCACCGGATGTGTTTGAACAATTCAATACCGCTGCACAAAAACTTGAAGAAATTCAGGAGTGTCATTTAGTATCCGGTGATTTCGACTACCTGTTGAAAACACGCGTGCCGGATATGTCAGCCTACCGTAAGTTGCt                                         <  2:507857/149‑1 (MQ=255)
                                                                                                       aTCGTGGCGCACCGGATGTGTTTGAACAATTCAATACCGCTGCACAAAAACTTGAAGAAATTCAGGAGTGTCATTTAGTATCCGGTGATTTCGACTACCTGTTGAAAACACGCGTGCCGGATATGTCAGCCTACCGTAAGTTGCTggggg                                    >  2:607992/1‑150 (MQ=255)
                                                                                                        tCGTGGCGCACCGGATGTGTTTGAACAATTCAATACCGCTGCACAAAAACTTGAAGAAATTCAGGAGTGTCATTTAGTATCCGGTGATTTCGACTACCTGTTGAAAACACGCGTGCCGGATATGTCAGCCTACCGTAAGTTGCTGGGGGa                                   >  1:538349/1‑150 (MQ=255)
                                                                                                           tGGCGCACCGGATGTGTTTGAACAATTCAATACCGCTGCACAAAAACTTGAAGAAATTCAGGAGTGTCATTTAGTATCCGGTGATTTCGACTACCTGTTGAAAACACGCGTGCCGGATATGTCAGCCTACCGTAAGTTGCTGGGGGAAAcc                               >  2:431044/1‑151 (MQ=255)
                                                                                                              cgcACCGGATGTGTTTGAACAATTCAATACCGCTGCACAAAAACTTGAAGAAATTCAGGAGTGTCATTTAGTATCCGGTGATTTCGACTACCTGTTGAAAACACGCGTGCCGGATATGTCAGCCTACCGTAAGTTGCTGGGGGAAACCctg                            >  1:340115/1‑151 (MQ=255)
                                                                                                                    ggATGTGTTTGAACAATTCAATACCGCTGCACAAAAACTTGAAGAAATTCAGGAGTGTCATTTAGTATCCGGTGATTTCGACTACCTGTTGAAAACACGCGTGCCGGATATGTCAGCCTACCGTAAGTTGCTGGGGGAAACCCTGCTGCGt                      <  2:146480/151‑1 (MQ=255)
                                                                                                                     gATGTGTTTGAACAATTCAATACCGCTGCACAAAAACTTGAAGAAATTCAGGAGTGTCATTTAGTATCCGGTGATTTCGACTACCTGTTGAAAACACGCGTGCCGGATATGTCAGCCTACCGTAAGTTGCTGGGGGAAACCCTGCTGCGTc                     >  1:243943/1‑151 (MQ=255)
                                                                                                                      aTGTGTTTGAACAATTCAATACCGCTGCACAAAAACTTGAAGAAATTCAGGAGTGTCATTTAGTATCCGGTGATTTCGACTACCTGTTGAAAACACGCGTGCCGGATATGTCAGCCTACCGTAAGTTGCTGGGGGAAACCCTGCTGCGTCt                    >  1:47351/1‑151 (MQ=255)
                                                                                                                          gtTTGAACAATTCAATACCGCTGCACAAAAACTTGAAGAAATTCAGGAGTGTCATTTAGTATCCGGTGATTTCGACTACCTGTTGAAAACACGCGTGCCGGATATGTCAGCCTACCGTAAGTTGCTGGGGGAAACCCTGCTGCGTCTGCCt                >  1:504212/1‑151 (MQ=255)
                                                                                                                                        aTACCGCTGCACAAAAACTTGAAGAAATTCAGGAGTGTCATTTAGTATCCGGTGATTTCGACTACCTGTTGAAAACACGCGTGCCGGATATGTCAGCCTACCGTAAGTTGCTGGGGGAAACCCTGCTGCGTCTGCCTGGCGTCAATGacac  >  2:639179/1‑151 (MQ=255)
                                                                                                                                                 |                                                                                                                                             
GTGCGTCGGCTGGAAAGACAAGGGTTTATTCAGGGCTATACGGCGCTGCTTAACCCCCATTATCTGGATGCATCACTTCTGGTATTCGTTGAGATTACTCTGAATCGTGGCGCACCGGATGTGTTTGAACAATTCAATACCGCTGTACAAAAACTTGAAGAAATTCAGGAGTGTCATTTAGTATCCGGTGATTTCGACTACCTGTTGAAAACACGCGTGCCGGATATGTCAGCCTACCGTAAGTTGCTGGGGGAAACCCTGCTGCGTCTGCCTGGCGTCAATGACAC  >  NC_000913/932739‑933025

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑