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