Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,951,536 Δ1 bp pseudogene (977/984 nt) ilvG → pseudogene, acetolactate synthase 2 large subunit, valine‑insensitive; acetolactate synthase II, large subunit, cryptic, interrupted

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,951,5340A.97.6% 177.9 / ‑1.6 43pseudogene (975/984 nt)ilvGpseudogene, acetolactate synthase 2 large subunit, valine‑insensitive; acetolactate synthase II, large subunit, cryptic, interrupted
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base . (17/24);  minor base G (0/1);  total (18/25)
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

CGGCAAACTGAACACCTTCGCGCCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTG  >  NC_000913/3951393‑3951663
                                                                                                                                             |                                                                                                                                 
cGGCAAACTGAACACCTTCGCGCCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACGAGCTGCGTCCGGCACATGTGGCAGTAGAAGGTGATTTAAACGCCCTGTTACCAGCATTCAAGCAGCCGTCaa                                                                                                                                  >  2:146048/1‑143 (MQ=255)
 ggCAAACTGAACACCTTCGCGCCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTTGTTTAGGTGCTCTGTCACCAGCGTCACAGCAGCCGTTaat                                                                                                                                 >  1:281322/1‑142 (MQ=255)
  gCAAACTGAACACCTTCGCGCCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTaatc                                                                                                                                >  1:548131/1‑141 (MQ=255)
       cTGAACACCTTCGCGCCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTT‑CATCAATGa                                                                                                                          <  1:599195/143‑1 (MQ=255)
       cTGAACACCTTCGCGACACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTT‑AATCAATGa                                                                                                                          <  1:500603/143‑1 (MQ=255)
            cacCTTCGCGCCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTT‑AATCAATGACTGgc                                                                                                                     <  1:569735/143‑1 (MQ=255)
             acgtgcgcgCCAGACGCCAGTGTTCTCCAGATGGAGATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTT‑AATCAATGACTGgca                                                                                                                    <  2:727383/138‑1 (MQ=255)
              ccTTCGCGCCACACGCCAGTATTATCCATATGGTTATCGACTCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTT‑GATCAATAACTGgcgg                                                                                                                   >  1:598576/1‑141 (MQ=255)
               cTTCGCGCCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTT‑AATCAATGACTGgcagc                                                                                                                  >  1:209941/1‑143 (MQ=255)
                 ccgcgCCACCCGCCAGTGTGATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTT‑AATCAATGACTGGCAGCaa                                                                                                                <  1:644498/142‑1 (MQ=255)
                       cacaCGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCAGCGTCAGGCACATGTGGCATCACAAGGGGATTTAAATGCTCTGTTACCAGCATTACACCAGCCGTT‑AAGCAATGACTGGCAGCAACACGgc                                                                                                          >  2:516106/1‑143 (MQ=255)
                          acGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTT‑AATCAATGACTGGCAGCAACACTgcgcg                                                                                                       >  1:687279/1‑143 (MQ=255)
                                 tgtTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTT‑AATCAATGACTGGCAGCAACACTGCGCGCAGCTGc                                                                                                <  1:659247/143‑1 (MQ=255)
                                 tgtTATCACTGTGGATATCGAGCCGGCAGAAATGAACAAGCGGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTT‑AATCAATGACTGGCAGCAACACTGCGCGCAGCTGc                                                                                                <  2:195244/143‑1 (MQ=255)
                                 tgtTAGCCATATGGAGATCGACCTGGCACAAATGAACAAGCGCCGTCAGGCACATGTGGCATTACAAGATGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTT‑AATCAATGACTGGCAGCAACACTGCGCGCAGCTGc                                                                                                <  2:611318/143‑1 (MQ=255)
                                     aTCCATAGGTATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTT‑AATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGtga                                                                                            <  1:469344/143‑1 (MQ=255)
                                         atatGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACCGCACCCGCT‑ACTCAAGGACTGGCAGCACCACTGCCCGCAGCTGCGTGATGaa                                                                                        >  2:576412/1‑143 (MQ=255)
                                         atatGGATATCGACCCGGCAGAAACGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTT‑AATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGaa                                                                                        <  1:663571/143‑1 (MQ=255)
                                           atGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTT‑AATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACa                                                                                      <  2:687279/143‑1 (MQ=255)
                                              gATAGCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTCGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTT‑AATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTc                                                                                   <  2:209344/143‑1 (MQ=255)
                                                        cGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACCAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTT‑AATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTAc                                                                         <  1:537823/143‑1 (MQ=255)
                                                          gaggcaaTGAACAAGCTGCGTCAGGCACATGTGGCATTACCAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTT‑AATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGa                                                                       <  2:238978/138‑1 (MQ=255)
                                                            aGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTT‑AATCAATGACTGGCAGCACCACTGCGCGCAGCTGCATGATGAACATTCCTGGCGTTACGAcc                                                                     >  2:143319/1‑143 (MQ=255)
                                                                  gagcAAGCTGCGTCAGGTACATGTGGCATTAGAAGGTGATTTAAATGCTCTGGTACCAACATGATAGC‑GTTGGGCAATCAATGACTGCCACCAACACTGCGCGCAGCTGCGTCACTAACCTCGGCAGCGTAACGACCATCCCg                                                               >  2:191213/4‑143 (MQ=255)
                                                                  gAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTT‑AATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCg                                                               >  1:741952/1‑143 (MQ=255)
                                                                   aaCAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTT‑AATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCgg                                                              <  1:464761/143‑1 (MQ=255)
                                                                      aaGCTGCGGCAGGGTCCTGTGCGATCACGAGGTGATTTAAATGCTCTGGTACCAGCAGTACAGCAGCCG‑TGAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGa                                                           <  2:519294/143‑1 (MQ=255)
                                                                             gTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTT‑AATCAATGACTGGCAGCAACACCGCGCGCAGCGGCGTGATGACCATTCCTGGCGTTACGACCATCCCGGTGACGCTATc                                                    >  2:159428/1‑143 (MQ=255)
                                                                                             aTTAAAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTT‑AATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgttgt                                    <  2:31449/143‑1 (MQ=255)
                                                                                              ttACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTT‑AATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgtt                                      <  2:431244/140‑1 (MQ=255)
                                                                                              ttACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTT‑AATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgtt                                      >  1:431244/1‑140 (MQ=255)
                                                                                                agacggTGATTTAAATGCTCTGTTACCAGCATGACAGCAGCCGTT‑AATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTaa                                 <  1:575672/139‑1 (MQ=255)
                                                                                                 cAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTT‑AATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTaaa                                >  1:70105/1‑143 (MQ=255)
                                                                                                    ggTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTT‑AATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAgaca                             >  1:333684/1‑143 (MQ=255)
                                                                                                      tGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTT‑AATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAaacaac                           >  2:59303/1‑143 (MQ=255)
                                                                                                         tttAAATGCTCTGTTACCAGCATTACAGCAGCCGTT‑AATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGAGGCTATCTACGCGCCGTTGTTGTTAAAACAACTGt                        <  1:562402/143‑1 (MQ=255)
                                                                                                         tgtAAATGCTCTGTGACCAGCACTACAGCAGCCGTT‑AATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGCTGAACATTCCTCGCGTTACGAACATCCCGGTGCCGCTATCTACGCGCCGTTGTTGTTAAaacaaccaa                        <  1:686113/141‑4 (MQ=255)
                                                                                                          ttAAATGCTCTGTTACCAGCATTACAGCAGCCGTT‑AATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGGTAAAACAACTGTc                       >  1:148318/1‑143 (MQ=255)
                                                                                                          ttAAATGCTCTGTTACCAGCATTACAGCAGCCGTT‑AATCAATGACTGGCAGCAACACTGCGCGAAGCTGCGTGATGAACATTCCTGGCGGTACGACTATCCCGGTGACGCTAGCTACGCGCCGTTGTTGGTAAAACAACTGTc                       >  1:147587/1‑143 (MQ=255)
                                                                                                              aTGCTCTGTTACCAGCATCACAGCAGCCGTT‑AATCAATGACTGGCAGCAACACTGCGCGCCGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGAt                   <  1:670637/143‑1 (MQ=255)
                                                                                                                gCTCTGTTACCAGCATTACAGCAGCCGTT‑AATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCg                 >  1:495558/1‑143 (MQ=255)
                                                                                                                  tctGTTACCAGCATTACAGCAGCCGTT‑AATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTATTCTTTTGATAACGCTCGGATCGtg               >  2:566430/1‑142 (MQ=255)
                                                                                                                    ggTTACCAGCATTACAGCAGCCGTT‑AATCAATGCCTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTaaa             <  1:435558/142‑1 (MQ=255)
                                                                                                                     gTTACCAGCATTAAAGCAGCCGTT‑AATCAATGACTGGCAGCAACACTGCGCGCCGCTTCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAAc            <  1:431821/143‑1 (MQ=255)
                                                                                                                        aCCAGCATTACAGCAGCCGTT‑AATCAATGACTGGCAGCAACACGGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAGACCTg         <  2:12302/143‑1 (MQ=255)
                                                                                                                               ttACAGCAGCCGTT‑AATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGGTGTTGACACCACTCTCGGCTCGTAACCATGCggctca  >  2:144466/1‑139 (MQ=255)
                                                                                                                                             |                                                                                                                                 
CGGCAAACTGAACACCTTCGCGCCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTG  >  NC_000913/3951393‑3951663

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