Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,950,420 G→A *33* (TAG→TAA ilvL → ilvG operon leader peptide

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,950,4200GA98.2% 180.3 / ‑6.5 57*33* (TAG→TAAilvLilvG operon leader peptide
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (26/30);  minor base C (0/1);  total (26/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAAC  >  NC_000913/3950284‑3950554
                                                                                                                                        |                                                                                                                                      
ttAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTaaag                                                                                                                                      >  1:755721/1‑139 (MQ=255)
 tAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAaaga                                                                                                                                     >  1:831094/1‑139 (MQ=255)
    gCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCa                                                                                                                                  <  1:109583/139‑1 (MQ=255)
       ttCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGTGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTAAAAGATc                                                                                                                                   <  1:1094864/135‑1 (MQ=255)
       ttCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATc                                                                                                                                   >  2:1094864/1‑135 (MQ=255)
       ttCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGc                                                                                                                               >  1:1241315/1‑139 (MQ=255)
            tCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTaa                                                                                                                          <  2:1241315/139‑1 (MQ=255)
               aaCAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGa                                                                                                                       <  2:271751/139‑1 (MQ=255)
               aaCAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGa                                                                                                                       >  1:271751/1‑139 (MQ=255)
                         aaGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGAccc                                                                                                             <  1:612161/139‑1 (MQ=255)
                         aaGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGAccc                                                                                                             <  1:209508/139‑1 (MQ=255)
                              aaGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCa                                                                                                        <  2:831094/139‑1 (MQ=255)
                                   aaaaTGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGaa                                                                                                   <  2:755721/139‑1 (MQ=255)
                                      aTGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAgg                                                                                                >  2:597843/1‑139 (MQ=255)
                                                       aGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGttttttttt                                                                               >  1:71860/1‑138 (MQ=255)
                                                            ttAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCtt                                                                          <  2:1176644/139‑1 (MQ=255)
                                                                    gTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGCTTTTTTTGACCTTAAAAACAt                                                                  >  1:1196413/1‑139 (MQ=255)
                                                                     tCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATa                                                                 >  2:1456092/1‑139 (MQ=255)
                                                                                  ggtggtgATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGa                                                    <  1:1492264/139‑1 (MQ=255)
                                                                                    tggtgATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGa                                                    >  1:1555871/1‑137 (MQ=255)
                                                                                    tggtgATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGa                                                    <  2:1555871/137‑1 (MQ=255)
                                                                                        gATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGa                                                                                                    <  2:1158349/85‑1 (MQ=255)
                                                                                        gATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGa                                                                                                    >  1:1158349/1‑85 (MQ=255)
                                                                                              tatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACa                                        <  1:1456092/139‑1 (MQ=255)
                                                                                                     ccGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGa                                                                                                    <  2:1415330/72‑1 (MQ=255)
                                                                                                     ccGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGa                                                                                                    >  1:1415330/1‑72 (MQ=255)
                                                                                                      cGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACaaaa                                >  2:204027/1‑139 (MQ=255)
                                                                                                       gTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTGACCTTAAAAACATAACCga                                                            >  1:246427/1‑109 (MQ=255)
                                                                                                       gTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTGACCTTAAAAACATAACCga                                                            <  2:246427/109‑1 (MQ=255)
                                                                                                       gTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAAtt                              >  1:1387993/1‑139 (MQ=255)
                                                                                                        tGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTg                                                                               >  1:2180/1‑90 (MQ=255)
                                                                                                        tGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTg                                                                               <  2:2180/90‑1 (MQ=255)
                                                                                                              gCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGttt                        >  2:1392963/1‑139 (MQ=255)
                                                                                                                 gCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCa                    <  2:1387993/139‑1 (MQ=255)
                                                                                                                  cACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGAc                                                                             <  1:1351748/82‑1 (MQ=255)
                                                                                                                  cACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGAc                                                                             >  2:1351748/1‑82 (MQ=255)
                                                                                                                  cACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTaaaaa                                                                     >  2:583607/1‑90 (MQ=255)
                                                                                                                  cACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTaaaaa                                                                     <  1:583607/90‑1 (MQ=255)
                                                                                                                  cACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTaa                                                                        <  2:1431015/87‑1 (MQ=255)
                                                                                                                  cACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTaa                                                                        >  1:1431015/1‑87 (MQ=255)
                                                                                                                  cACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAAcc                                                              >  1:1208533/1‑97 (MQ=255)
                                                                                                                  cACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAAcc                                                              <  2:1208533/97‑1 (MQ=255)
                                                                                                                       ggACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGaa                                                                                                   >  2:1510895/1‑55 (MQ=255)
                                                                                                                       ggACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGaa                                                                                                   <  1:1510895/55‑1 (MQ=255)
                                                                                                                       ggACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCtt                                                                          >  2:46714/1‑80 (MQ=255)
                                                                                                                       ggACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCtt                                                                          <  1:46714/80‑1 (MQ=255)
                                                                                                                       ggACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCa                                                      <  2:857150/100‑1 (MQ=255)
                                                                                                                       ggACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCa                                                      >  1:857150/1‑100 (MQ=255)
                                                                                                                       ggACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGAtt               <  1:204027/139‑1 (MQ=255)
                                                                                                                          cGAGGAAAGGCTTACAGATCAAGCCTTAACGAACTAAGACCACCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGATAATGAATAACAGca                                     <  2:1015862/114‑1 (MQ=255)
                                                                                                                          cGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGca                                     >  1:1015862/1‑114 (MQ=255)
                                                                                                                            aGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGAc                                                                             <  2:1107028/72‑1 (MQ=255)
                                                                                                                            aGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGAc                                                                             >  1:1107028/1‑72 (MQ=255)
                                                                                                                              gAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACg        <  1:1335890/139‑1 (MQ=255)
                                                                                                                                 aGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACgggg     >  1:761951/1‑139 (MQ=255)
                                                                                                                                   gCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCa                                                      <  2:1246902/88‑1 (MQ=255)
                                                                                                                                   gCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCa                                                      >  1:1246902/1‑88 (MQ=255)
                                                                                                                                    cTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaac  <  2:1196413/139‑1 (MQ=255)
                                                                                                                                    cTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaac  <  2:707967/139‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                      
TTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAAC  >  NC_000913/3950284‑3950554

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