Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,852,559 T→A 7.0% K76* (AAA→TAA) ‡ ilvB ← acetolactate synthase 2 large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,852,5590TA 7.0% 188.4 / 5.8 57K76* (AAA→TAA) ‡ilvBacetolactate synthase 2 large subunit
Reads supporting (aligned to +/- strand):  ref base T (23/30);  new base A (2/2);  total (25/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

GGAAGGCGTCGGTGCCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGATAGAACCGCCCGGAATGCCTGTCACAATCTTAATG  >  NC_000913/3852418‑3852704
                                                                                                                                             |                                                                                                                                                 
ggAAGGCGTCGGTGCCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCg                                                                                                                                            >  1:434803/1‑149 (MQ=255)
ggAAGGCGTCGGTGCCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGACGGTAAACCggcg                                                                                                                                            >  2:375797/1‑146 (MQ=255)
ggAAGGCGTCGGTGCCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGACGGTAAACCGGCg                                                                                                                                            >  2:519134/1‑149 (MQ=255)
  aaGGCGTCGGTGCCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGACGGTAAACCGGCGGt                                                                                                                                          <  1:519134/149‑1 (MQ=255)
  aaGGCGTCGGTGCCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGACGGTAAACCGGCGGt                                                                                                                                          <  1:375797/149‑1 (MQ=255)
      cGTCGGTGCCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTgcgc                                                                                                                                      <  1:140167/149‑1 (MQ=255)
              ccGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTcc                                                                                                                              <  1:33578/149‑1 (MQ=255)
              ccGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTcc                                                                                                                              >  1:59793/1‑149 (MQ=255)
                gATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTTACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCt                                                                                                                            >  4:371076/1‑149 (MQ=255)
                gATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCt                                                                                                                            >  1:80061/1‑149 (MQ=255)
                gATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCt                                                                                                                            >  4:379053/1‑149 (MQ=255)
                       gAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGAt                                                                                                                     >  3:395462/1‑149 (MQ=255)
                        aGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGATTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATa                                                                                                                    >  4:327300/1‑149 (MQ=255)
                                  cTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCgcgc                                                                                                          <  3:251546/149‑1 (MQ=255)
                                  cTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCgcgc                                                                                                          <  2:255647/149‑1 (MQ=255)
                                        aGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGtt                                                                                                    >  2:89550/1‑149 (MQ=255)
                                         gTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTc                                                                                                   <  3:17927/149‑1 (MQ=255)
                                           gATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCAt                                                                                                 <  2:434803/149‑1 (MQ=255)
                                            aTGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATg                                                                                                >  1:16640/1‑149 (MQ=255)
                                             tGCAAATCAGCGGGATTGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGa                                                                                               >  2:368533/1‑149 (MQ=255)
                                               cAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACg                                                                                             <  2:126260/149‑1 (MQ=255)
                                                          gATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATCCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAatat                                                                                  <  4:353306/149‑1 (MQ=255)
                                                           aTGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATg                                                                                 <  2:469293/149‑1 (MQ=255)
                                                            tGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATgg                                                                                <  1:89550/149‑1 (MQ=255)
                                                                gTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGTGGa                                                                            <  2:288532/149‑1 (MQ=255)
                                                                   cAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGAttt                                                                         >  4:124477/1‑149 (MQ=255)
                                                                          gcATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCt                                                                  >  3:67601/1‑149 (MQ=255)
                                                                            aTCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCttt                                                                >  2:384523/1‑149 (MQ=255)
                                                                             tCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTg                                                               >  1:202044/1‑149 (MQ=255)
                                                                             tCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTg                                                               <  3:483688/149‑1 (MQ=255)
                                                                              cGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTgg                                                              <  2:59793/149‑1 (MQ=255)
                                                                                gCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCt                                                            <  4:523998/149‑1 (MQ=255)
                                                                                   aTGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTaa                                                         <  2:80061/149‑1 (MQ=255)
                                                                                    tGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAg                                                        <  2:241530/149‑1 (MQ=255)
                                                                                    tGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAg                                                        >  3:247358/1‑149 (MQ=255)
                                                                                         gTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATc                                                   <  4:402433/149‑1 (MQ=255)
                                                                                            aCCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGACTTGCGTGCTTTGGCTTAAGGCATCGTa                                                <  3:379053/149‑1 (MQ=255)
                                                                                            aCCAGGTTAGTCGCACCCGGTCCGCCACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTa                                                <  3:371076/149‑1 (MQ=255)
                                                                                              cAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAg                                                                                                                 >  4:377862/1‑82 (MQ=255)
                                                                                              cAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAg                                                                                                                 <  3:377862/82‑1 (MQ=255)
                                                                                                     gTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAAcaggtag                                       <  4:67601/149‑1 (MQ=255)
                                                                                                       cGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGa                                     <  4:209765/149‑1 (MQ=255)
                                                                                                         cACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGATa                                   <  3:327300/149‑1 (MQ=255)
                                                                                                          aCCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGATAg                                  >  3:73836/1‑149 (MQ=255)
                                                                                                                tCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGATAGAACCGc                            <  3:141070/149‑1 (MQ=255)
                                                                                                                    cTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGATAGAACCGCCCgg                        >  3:244078/1‑149 (MQ=255)
                                                                                                                       cAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGATAGAACCGCCCGGAAt                     <  4:60771/149‑1 (MQ=255)
                                                                                                                         ggCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGATAGAACCGCCCGGAATGc                   >  3:114229/1‑149 (MQ=255)
                                                                                                                         ggCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGATAGAACCGCCCGGAATGc                   >  1:442184/1‑149 (MQ=255)
                                                                                                                           ccATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGATAGAACCGCCCGGAATGCCt                 <  2:104724/149‑1 (MQ=255)
                                                                                                                            cATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGATAGAACCGCCCGGAATGCCTg                >  4:349945/1‑149 (MQ=255)
                                                                                                                              tACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGATAGAACCGCCCGGAATGCCTGTc              <  3:349945/149‑1 (MQ=255)
                                                                                                                                  gACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGATAGAACCGCCCGGAATGCCTGTCACaa          >  3:186397/1‑149 (MQ=255)
                                                                                                                                   aCCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGATAGAACCGCCCGGAATGCCTGTCACAAt         <  4:395462/149‑1 (MQ=255)
                                                                                                                                      gccgGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGATAGAACCGCCCGGAATGCCTGTCACAATCtt      >  1:125341/1‑149 (MQ=255)
                                                                                                                                         ggTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGATAGAACCGCCCGGAATGCCTGTCACAATCTTAAt   <  4:247358/149‑1 (MQ=255)
                                                                                                                                          gTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGATAGAACCGCCCGGAATGCCTGTCACAATCTTAATg  <  3:107339/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                 
GGAAGGCGTCGGTGCCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGATAGAACCGCCCGGAATGCCTGTCACAATCTTAATG  >  NC_000913/3852418‑3852704

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