Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,950,468 G→T intergenic (+48/‑39) ilvL → / → ilvX ilvG operon leader peptide/uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,950,4680GT92.5% 199.4 / 9.9 67intergenic (+48/‑39)ilvL/ilvXilvG operon leader peptide/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (1/4);  new base T (23/39);  total (24/43)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.51e-01

TGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAA  >  NC_000913/3950371‑3950564
                                                                                                 |                                                                                                
tgATTATGATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGttt                                                                                                >  2:1024781/1‑100 (MQ=255)
   ttattaTCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGtttttt                                                                                             >  2:1298992/1‑100 (MQ=255)
      ttatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTtttttttt                                                                                          >  2:404118/1‑100 (MQ=255)
      ttatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTtttttttt                                                                                          <  1:884376/100‑1 (MQ=255)
        atCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTtttttttttg                                                                                        <  2:1778198/100‑3 (MQ=255)
          cccACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGAcc                                                                                      <  2:709072/100‑1 (MQ=255)
          cccACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGAcc                                                                                      <  1:1732804/100‑1 (MQ=255)
          cccACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGAcc                                                                                      <  2:1810465/100‑1 (MQ=255)
           ccACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGAc                                                                                       <  2:1745947/98‑1 (MQ=255)
           ccACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCt                                                                                     >  2:1286238/1‑100 (MQ=255)
           ccACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCt                                                                                     >  2:1075235/1‑100 (MQ=255)
           ccACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTGACCtt                                                                                    <  1:1691319/100‑1 (MQ=255)
            cACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCtt                                                                                    <  1:373820/100‑1 (MQ=255)
             aCCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTa                                                                                   >  1:1339593/1‑100 (MQ=255)
             aCCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCGTAACGAACTAAGACCCCCGCGCCGAAAGGTCCGGGGTTTTTTTTTGACCTTa                                                                                   <  1:525313/100‑1 (MQ=255)
              ccGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTaa                                                                                  <  2:1028235/100‑1 (MQ=255)
               cGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTaaa                                                                                 <  1:1573337/100‑1 (MQ=255)
                gTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTaaaa                                                                                >  2:1838861/1‑100 (MQ=255)
                 tGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTaaaaa                                                                               <  2:1117427/100‑1 (MQ=255)
                 tGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTaaaaa                                                                               <  2:1528496/100‑1 (MQ=255)
                  gCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAAc                                                                              <  1:449417/100‑1 (MQ=255)
                   cGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTaaaaa                                                                               <  1:234309/98‑1 (MQ=255)
                   cGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACa                                                                             <  1:973286/100‑1 (MQ=255)
                    ggggCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTaaaacata                                                                            <  2:1243733/100‑5 (MQ=255)
                       ggtGCACTTGGACGAGGAAGGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAAc                                                                         <  2:21582/98‑1 (MQ=255)
                       gCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAAc                                                                         >  2:1687472/1‑100 (MQ=255)
                       gCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACATTAAAAACATAAc                                                                         <  1:582489/100‑1 (MQ=255)
                        cTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAAcc                                                                        <  1:1030263/100‑1 (MQ=255)
                          gCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCg                                                                       >  1:814080/1‑99 (MQ=255)
                                ggACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCa                                                                >  1:1067149/1‑100 (MQ=255)
                                   cGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGAc                                                             <  1:784537/100‑1 (MQ=255)
                                    gAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACa                                                            <  1:842574/100‑1 (MQ=255)
                                        aaaGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGa                                                        <  2:80720/100‑1 (MQ=255)
                                             cTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAAc                                                   >  2:953352/1‑100 (MQ=255)
                                                agagATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGc                                                <  2:432997/100‑1 (MQ=255)
                                                 gagaTCAAGCCTTAACGACCTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGca                                               <  1:1075234/100‑1 (MQ=255)
                                                 gagaTCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGca                                               <  1:1286231/100‑1 (MQ=255)
                                                 gagaTCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGca                                               >  2:457923/1‑100 (MQ=255)
                                                   gaTCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGcaca                                             >  1:1448213/1‑100 (MQ=255)
                                                      cAAGCCTTAACGAACTAAGACCCCCGCACCAAAGGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATACCAGCACaaaa                                          >  1:1278156/1‑100 (MQ=255)
                                                        aGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTc                                       <  2:729028/100‑1 (MQ=255)
                                                         gCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTc                                       >  1:1142509/1‑100 (MQ=255)
                                                          ccTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCt                                      >  1:820408/1‑100 (MQ=255)
                                                           cTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTg                                     <  1:193326/100‑1 (MQ=255)
                                                            ttAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGt                                    <  1:596364/100‑1 (MQ=255)
                                                              aaCGAACTAAGACCCCCGCACCGAAAGGGCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCATACAATGAATAACAGCACAAAATTCTtttt                                  >  1:673528/1‑100 (MQ=255)
                                                               aCGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGCGCAGACAATGAATAACAGCACAAAATTCTGTTtc                                 >  2:1834592/1‑100 (MQ=255)
                                                               aCGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTtc                                 <  2:688949/100‑1 (MQ=255)
                                                                cGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTTTTTTTTTAAACAAAACCGAGGAGCAGACAATGAATAACAGCACGAATTTCTGTTtc                                 >  1:98828/1‑99 (MQ=255)
                                                                  aaCTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCa                              <  1:1244831/100‑1 (MQ=255)
                                                                   aCTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCaa                             <  2:1567292/100‑1 (MQ=255)
                                                                     tAAGCCCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGa                           <  2:1339601/100‑1 (MQ=255)
                                                                     tAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTTTTTTTTAAAACATAAACGAAGAGCAGACAATGAATAAGAGCACAGAATTCCTTTTCTCaa                             >  1:813877/1‑98 (MQ=255)
                                                                      aaGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGAt                          >  1:802712/1‑100 (MQ=255)
                                                                       aGACCCCCGCACCGAAATGGCCGGGGTTTTTTTTTGACATTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCa                              <  2:961515/95‑1 (MQ=255)
                                                                        gACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTTACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTCTCTTAAGATTc                        <  1:1024780/100‑1 (MQ=255)
                                                                        gACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTc                        <  2:1517181/100‑1 (MQ=255)
                                                                         aCCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCa                       <  1:1687462/100‑1 (MQ=255)
                                                                          cccccGCAGCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTc                        <  1:1277695/98‑1 (MQ=255)
                                                                          cccccGCACCGAAAGGTCCGGGGTTGTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAg                      <  1:1149626/99‑1 (MQ=255)
                                                                            cacGCATCGAAAGGTCCGGGGTTTTTTTTTGACATTAAAAACATAACCGAGGAGCAGACAATGACTAACAGCACAAAATTCTGTTTCTCAAGATTCAGGa                    <  2:1334853/98‑1 (MQ=255)
                                                                             ccGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGa                    <  2:1383250/99‑1 (MQ=255)
                                                                                cACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACggg                <  1:1349853/100‑1 (MQ=255)
                                                                                  ccGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGa              >  2:763715/1‑100 (MQ=255)
                                                                                   cGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaa             >  2:262423/1‑100 (MQ=255)
                                                                                       aGGTCCGGGGTTTTTTTTTTTTTTTTTAAACATAACCGAGGAGCAGACAATGAATTACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaactaa         >  1:591779/1‑100 (MQ=255)
                                                                                        ggTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaactaac        >  2:1618932/1‑100 (MQ=255)
                                                                                         gTCCGGGGTTTTTTTTTTTCTTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaactaa         >  1:1190672/1‑98 (MQ=255)
                                                                                          tCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAactaacta      <  2:969574/100‑1 (MQ=255)
                                                                                              gggtTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGaa  >  2:605411/1‑100 (MQ=255)
                                                                                                 |                                                                                                
TGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAA  >  NC_000913/3950371‑3950564

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: