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)
       cTGAACACCTTCGCGCCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTCATCAATGa                                                                                                                          <  1:599195/143‑1 (MQ=255)
       cTGAACACCTTCGCGACACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAATCAATGa                                                                                                                          <  1:500603/143‑1 (MQ=255)
            cacCTTCGCGCCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAATCAATGACTGgc                                                                                                                     <  1:569735/143‑1 (MQ=255)
             acgtgcgcgCCAGACGCCAGTGTTCTCCAGATGGAGATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAATCAATGACTGgca                                                                                                                    <  2:727383/138‑1 (MQ=255)
              ccTTCGCGCCACACGCCAGTATTATCCATATGGTTATCGACTCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTGATCAATAACTGgcgg                                                                                                                   >  1:598576/1‑141 (MQ=255)
               cTTCGCGCCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAATCAATGACTGgcagc                                                                                                                  >  1:209941/1‑143 (MQ=255)
                 ccgcgCCACCCGCCAGTGTGATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAATCAATGACTGGCAGCaa                                                                                                                <  1:644498/142‑1 (MQ=255)
                       cacaCGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCAGCGTCAGGCACATGTGGCATCACAAGGGGATTTAAATGCTCTGTTACCAGCATTACACCAGCCGTTAAGCAATGACTGGCAGCAACACGgc                                                                                                          >  2:516106/1‑143 (MQ=255)
                          acGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAATCAATGACTGGCAGCAACACTgcgcg                                                                                                       >  1:687279/1‑143 (MQ=255)
                                 tgtTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAATCAATGACTGGCAGCAACACTGCGCGCAGCTGc                                                                                                <  1:659247/143‑1 (MQ=255)
                                 tgtTATCACTGTGGATATCGAGCCGGCAGAAATGAACAAGCGGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAATCAATGACTGGCAGCAACACTGCGCGCAGCTGc                                                                                                <  2:195244/143‑1 (MQ=255)
                                 tgtTAGCCATATGGAGATCGACCTGGCACAAATGAACAAGCGCCGTCAGGCACATGTGGCATTACAAGATGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAATCAATGACTGGCAGCAACACTGCGCGCAGCTGc                                                                                                <  2:611318/143‑1 (MQ=255)
                                     aTCCATAGGTATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGtga                                                                                            <  1:469344/143‑1 (MQ=255)
                                         atatGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACCGCACCCGCTACTCAAGGACTGGCAGCACCACTGCCCGCAGCTGCGTGATGaa                                                                                        >  2:576412/1‑143 (MQ=255)
                                         atatGGATATCGACCCGGCAGAAACGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGaa                                                                                        <  1:663571/143‑1 (MQ=255)
                                           atGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACa                                                                                      <  2:687279/143‑1 (MQ=255)
                                              gATAGCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTCGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTc                                                                                   <  2:209344/143‑1 (MQ=255)
                                                        cGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACCAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTAc                                                                         <  1:537823/143‑1 (MQ=255)
                                                          gaggcaaTGAACAAGCTGCGTCAGGCACATGTGGCATTACCAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGa                                                                       <  2:238978/138‑1 (MQ=255)
                                                            aGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAATCAATGACTGGCAGCACCACTGCGCGCAGCTGCATGATGAACATTCCTGGCGTTACGAcc                                                                     >  2:143319/1‑143 (MQ=255)
                                                                  gagcAAGCTGCGTCAGGTACATGTGGCATTAGAAGGTGATTTAAATGCTCTGGTACCAACATGATAGCGTTGGGCAATCAATGACTGCCACCAACACTGCGCGCAGCTGCGTCACTAACCTCGGCAGCGTAACGACCATCCCg                                                               >  2:191213/4‑143 (MQ=255)
                                                                  gAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCg                                                               >  1:741952/1‑143 (MQ=255)
                                                                   aaCAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCgg                                                              <  1:464761/143‑1 (MQ=255)
                                                                      aaGCTGCGGCAGGGTCCTGTGCGATCACGAGGTGATTTAAATGCTCTGGTACCAGCAGTACAGCAGCCGTGAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGa                                                           <  2:519294/143‑1 (MQ=255)
                                                                             gTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAATCAATGACTGGCAGCAACACCGCGCGCAGCGGCGTGATGACCATTCCTGGCGTTACGACCATCCCGGTGACGCTATc                                                    >  2:159428/1‑143 (MQ=255)
                                                                                             aTTAAAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgttgt                                    <  2:31449/143‑1 (MQ=255)
                                                                                              ttACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgtt                                      <  2:431244/140‑1 (MQ=255)
                                                                                              ttACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgtt                                      >  1:431244/1‑140 (MQ=255)
                                                                                                agacggTGATTTAAATGCTCTGTTACCAGCATGACAGCAGCCGTTAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTaa                                 <  1:575672/139‑1 (MQ=255)
                                                                                                 cAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTaaa                                >  1:70105/1‑143 (MQ=255)
                                                                                                    ggTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAgaca                             >  1:333684/1‑143 (MQ=255)
                                                                                                      tGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAaacaac                           >  2:59303/1‑143 (MQ=255)
                                                                                                         tttAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGAGGCTATCTACGCGCCGTTGTTGTTAAAACAACTGt                        <  1:562402/143‑1 (MQ=255)
                                                                                                         tgtAAATGCTCTGTGACCAGCACTACAGCAGCCGTTAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGCTGAACATTCCTCGCGTTACGAACATCCCGGTGCCGCTATCTACGCGCCGTTGTTGTTAAaacaaccaa                        <  1:686113/141‑4 (MQ=255)
                                                                                                          ttAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGGTAAAACAACTGTc                       >  1:148318/1‑143 (MQ=255)
                                                                                                          ttAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAATCAATGACTGGCAGCAACACTGCGCGAAGCTGCGTGATGAACATTCCTGGCGGTACGACTATCCCGGTGACGCTAGCTACGCGCCGTTGTTGGTAAAACAACTGTc                       >  1:147587/1‑143 (MQ=255)
                                                                                                              aTGCTCTGTTACCAGCATCACAGCAGCCGTTAATCAATGACTGGCAGCAACACTGCGCGCCGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGAt                   <  1:670637/143‑1 (MQ=255)
                                                                                                                gCTCTGTTACCAGCATTACAGCAGCCGTTAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCg                 >  1:495558/1‑143 (MQ=255)
                                                                                                                  tctGTTACCAGCATTACAGCAGCCGTTAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTATTCTTTTGATAACGCTCGGATCGtg               >  2:566430/1‑142 (MQ=255)
                                                                                                                    ggTTACCAGCATTACAGCAGCCGTTAATCAATGCCTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTaaa             <  1:435558/142‑1 (MQ=255)
                                                                                                                     gTTACCAGCATTAAAGCAGCCGTTAATCAATGACTGGCAGCAACACTGCGCGCCGCTTCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAAc            <  1:431821/143‑1 (MQ=255)
                                                                                                                        aCCAGCATTACAGCAGCCGTTAATCAATGACTGGCAGCAACACGGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAGACCTg         <  2:12302/143‑1 (MQ=255)
                                                                                                                               ttACAGCAGCCGTTAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGGTGTTGACACCACTCTCGGCTCGTAACCATGCggctca  >  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: