Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,951,535 Δ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,5330A.97.7% 186.4 / ‑1.6 45pseudogene (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/26);  minor base G (0/1);  total (18/27)
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/3951392‑3951662
                                                                                                                                             |                                                                                                                                 
cGGCAAACTGAACACCTTCGCGCCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACGAGCTGCGTCCGGCACATGTGGCAGTAGAAGGTGATTTAAACGCCCTGTTACCAGCATTCAAGCAGCCGTCaa                                                                                                                                  >  1:146048/1‑143 (MQ=255)
 ggCAAACTGAACACCTTCGCGCCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTTGTTTAGGTGCTCTGTCACCAGCGTCACAGCAGCCGTTaat                                                                                                                                 >  2:281322/1‑142 (MQ=255)
  gCAAACTGAACACCTTCGCGCCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTaatc                                                                                                                                >  2:548131/1‑141 (MQ=255)
       cTGAACACCTTCGCGCCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTT‑CATCAATGa                                                                                                                          <  2:599195/143‑1 (MQ=255)
       cTGAACACCTTCGCGACACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTT‑AATCAATGa                                                                                                                          <  2:500603/143‑1 (MQ=255)
            cacCTTCGCGCCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTT‑AATCAATGACTGgc                                                                                                                     <  2:569735/143‑1 (MQ=255)
             acgtgcgcgCCAGACGCCAGTGTTCTCCAGATGGAGATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTT‑AATCAATGACTGgca                                                                                                                    <  1:727383/138‑1 (MQ=255)
              ccTTCGCGCCACACGCCAGTATTATCCATATGGTTATCGACTCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTT‑GATCAATAACTGgcgg                                                                                                                   >  2:598576/1‑141 (MQ=255)
               cTTCGCGCCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTT‑AATCAATGACTGgcagc                                                                                                                  >  2:209941/1‑143 (MQ=255)
                 ccgcgCCACCCGCCAGTGTGATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTT‑AATCAATGACTGGCAGCaa                                                                                                                <  2:644498/142‑1 (MQ=255)
                       cacaCGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCAGCGTCAGGCACATGTGGCATCACAAGGGGATTTAAATGCTCTGTTACCAGCATTACACCAGCCGTT‑AAGCAATGACTGGCAGCAACACGgc                                                                                                          >  1:516106/1‑143 (MQ=255)
                          acGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTT‑AATCAATGACTGGCAGCAACACTgcgcg                                                                                                       >  2:687279/1‑143 (MQ=255)
                                 tgtTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTT‑AATCAATGACTGGCAGCAACACTGCGCGCAGCTGc                                                                                                <  2:659247/143‑1 (MQ=255)
                                 tgtTATCACTGTGGATATCGAGCCGGCAGAAATGAACAAGCGGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTT‑AATCAATGACTGGCAGCAACACTGCGCGCAGCTGc                                                                                                <  1:195244/143‑1 (MQ=255)
                                 tgtTAGCCATATGGAGATCGACCTGGCACAAATGAACAAGCGCCGTCAGGCACATGTGGCATTACAAGATGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTT‑AATCAATGACTGGCAGCAACACTGCGCGCAGCTGc                                                                                                <  1:611318/143‑1 (MQ=255)
                                 gttgactctgttGGTTACCGACCCGGCACCAATGACCAAGCTGCTTCAGGCACATGTGGCATTACAAGGTGATTTACATTCTCTGTTACTATCATTACAGCAGCCGTT‑AATCAATGACTGGCAGCAACACTGCGCGCAGCTGc                                                                                                <  1:611617/132‑1 (MQ=255)
                                     aTCCATAGGTATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTT‑AATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGtga                                                                                            <  2:469344/143‑1 (MQ=255)
                                         atatGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACCGCACCCGCT‑ACTCAAGGACTGGCAGCACCACTGCCCGCAGCTGCGTGATGaa                                                                                        >  1:576412/1‑143 (MQ=255)
                                         atatGGATATCGACCCGGCAGAAACGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTT‑AATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGaa                                                                                        <  2:663571/143‑1 (MQ=255)
                                           atGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTT‑AATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACa                                                                                      <  1:687279/143‑1 (MQ=255)
                                              gATAGCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTCGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTT‑AATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTc                                                                                   <  1:209344/143‑1 (MQ=255)
                                                        cGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACCAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTT‑AATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTAc                                                                         <  2:537823/143‑1 (MQ=255)
                                                          gaggcaaTGAACAAGCTGCGTCAGGCACATGTGGCATTACCAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTT‑AATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGa                                                                       <  1:238978/138‑1 (MQ=255)
                                                            aGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTT‑AATCAATGACTGGCAGCACCACTGCGCGCAGCTGCATGATGAACATTCCTGGCGTTACGAcc                                                                     >  1:143319/1‑143 (MQ=255)
                                                                  gagcAAGCTGCGTCAGGTACATGTGGCATTAGAAGGTGATTTAAATGCTCTGGTACCAACATGATAGC‑GTTGGGCAATCAATGACTGCCACCAACACTGCGCGCAGCTGCGTCACTAACCTCGGCAGCGTAACGACCATCCCg                                                               >  1:191213/4‑143 (MQ=255)
                                                                  gAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTT‑AATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCg                                                               >  2:741952/1‑143 (MQ=255)
                                                                   aaCAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTT‑AATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCgg                                                              <  2:464761/143‑1 (MQ=255)
                                                                      aaGCTGCGGCAGGGTCCTGTGCGATCACGAGGTGATTTAAATGCTCTGGTACCAGCAGTACAGCAGCCG‑TGAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGa                                                           <  1:519294/143‑1 (MQ=255)
                                                                             gTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTT‑AATCAATGACTGGCAGCAACACCGCGCGCAGCGGCGTGATGACCATTCCTGGCGTTACGACCATCCCGGTGACGCTATc                                                    >  1:159428/1‑143 (MQ=255)
                                                                                             aTTAAAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTT‑AATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgttgt                                    <  1:31449/143‑1 (MQ=255)
                                                                                              ttACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTT‑AATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgtt                                      <  1:431244/140‑1 (MQ=255)
                                                                                              ttACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTT‑AATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgtt                                      >  2:431244/1‑140 (MQ=255)
                                                                                                agacggTGATTTAAATGCTCTGTTACCAGCATGACAGCAGCCGTT‑AATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTaa                                 <  2:575672/139‑1 (MQ=255)
                                                                                                 cacagttgtttATATGCGCTGGTACCAGCATTACAGCAGCCGTT‑AATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTaaa                                <  1:180361/135‑1 (MQ=255)
                                                                                                 cAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTT‑AATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTaaa                                >  2:70105/1‑143 (MQ=255)
                                                                                                    ggTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTT‑AATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAgaca                             >  2:333684/1‑143 (MQ=255)
                                                                                                      tGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTT‑AATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAaacaac                           >  1:59303/1‑143 (MQ=255)
                                                                                                         tttAAATGCTCTGTTACCAGCATTACAGCAGCCGTT‑AATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGAGGCTATCTACGCGCCGTTGTTGTTAAAACAACTGt                        <  2:562402/143‑1 (MQ=255)
                                                                                                         tgtAAATGCTCTGTGACCAGCACTACAGCAGCCGTT‑AATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGCTGAACATTCCTCGCGTTACGAACATCCCGGTGCCGCTATCTACGCGCCGTTGTTGTTAAaacaaccaa                        <  2:686113/141‑4 (MQ=255)
                                                                                                          ttAAATGCTCTGTTACCAGCATTACAGCAGCCGTT‑AATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGGTAAAACAACTGTc                       >  2:148318/1‑143 (MQ=255)
                                                                                                          ttAAATGCTCTGTTACCAGCATTACAGCAGCCGTT‑AATCAATGACTGGCAGCAACACTGCGCGAAGCTGCGTGATGAACATTCCTGGCGGTACGACTATCCCGGTGACGCTAGCTACGCGCCGTTGTTGGTAAAACAACTGTc                       >  2:147587/1‑143 (MQ=255)
                                                                                                              aTGCTCTGTTACCAGCATCACAGCAGCCGTT‑AATCAATGACTGGCAGCAACACTGCGCGCCGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGAt                   <  2:670637/143‑1 (MQ=255)
                                                                                                                gCTCTGTTACCAGCATTACAGCAGCCGTT‑AATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCg                 >  2:495558/1‑143 (MQ=255)
                                                                                                                  tctGTTACCAGCATTACAGCAGCCGTT‑AATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTATTCTTTTGATAACGCTCGGATCGtg               >  1:566430/1‑142 (MQ=255)
                                                                                                                    ggTTACCAGCATTACAGCAGCCGTT‑AATCAATGCCTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTaaa             <  2:435558/142‑1 (MQ=255)
                                                                                                                     gTTACCAGCATTAAAGCAGCCGTT‑AATCAATGACTGGCAGCAACACTGCGCGCCGCTTCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAAc            <  2:431821/143‑1 (MQ=255)
                                                                                                                        aCCAGCATTACAGCAGCCGTT‑AATCAATGACTGGCAGCAACACGGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAGACCTg         <  1:12302/143‑1 (MQ=255)
                                                                                                                               ttACAGCAGCCGTT‑AATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGGTGTTGACACCACTCTCGGCTCGTAACCATGCggctca  >  1:144466/1‑139 (MQ=255)
                                                                                                                                             |                                                                                                                                 
CGGCAAACTGAACACCTTCGCGCCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTG  >  NC_000913/3951392‑3951662

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