Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,740,857 T→C K98K (AAA→AAG)  aroF ← 3‑deoxy‑7‑phosphoheptulonate synthase, Tyr‑sensitive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,740,8570TC97.9% 158.9 / NA 47K98K (AAA→AAG) aroF3‑deoxy‑7‑phosphoheptulonate synthase, Tyr‑sensitive
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (19/27);  minor base A (1/0);  total (20/27)
Fisher's exact test for biased strand distribution p-value = 4.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCGTCGCCAGTGGCAGTCCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAGCCCGGCTTCTACATCAAAAGAGCCATCCATATGGGGATCGTTAATTAACCCTTTCCAGCCGACAGTGGTACGGGGTTTTTCAAAATAGACGCGCATTACCAGATAGAGGCTATCGCTGACCTCTGCGGCAAGGGCTTTAAATCGACGAGCATATTCCAGAGCAGTTTCCGGATCATGAATGGAACAAGGACCACATACTACCAGCAGACG  >  NC_000913/2740720‑2740991
                                                                                                                                         |                                                                                                                                      
tCCGTCGCCAGTGGCAGTCCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAGCCCGGCTTCTACATCAAAAGAGCCATCCATATGGGGATCGTTAATTAACCCTTTCCAGCCGACAGTGGTACGGGGctt                                                                                                                                      >  2:145039/1‑140 (MQ=255)
tCCGTCGCCAGTGGCAGTCCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAGCCCGGCTTCTACATCAAAAGAGCCATCCAGATGGGGATCGTTAATTTACCCTTTCAAGCCGACAGTGGTACGGGGctt                                                                                                                                      >  2:393016/1‑140 (MQ=255)
 ccGTCGCCAGTGGCAGTCCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAGCCCGGCTTCTACATCAAAAGAGCCATCCATATGGGGATCGTTAATTAACCCTTTCCAGCCGACAGTGGTACGGGGcttt                                                                                                                                     >  2:42451/1‑140 (MQ=255)
  cGTCGCCAGTGGCAGTCCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAGCCCGGCTTCTACATCAAAAGAGCCATCCATATGGGGATCGTTGATTAACCATTTGCATCCGACATTGGTACGGGGctttt                                                                                                                                    >  1:402148/1‑140 (MQ=255)
  cGTCGCCAGTGGCAGTCCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAGCCCGGCTTCTACATCAAAAGAGCCATCCATATGGGGATCGTGAATTAACCCTTTCCAGCCGGCAGGGGTACGGGGctttt                                                                                                                                    >  1:543019/1‑140 (MQ=255)
           tGGCAGTCCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAGCCCGGCTTCTACATCAAAAGAGCCATCCATATGGGGATCGTTAATTAACCCTTTCCAGCCGACAGTGGTACGGGGCTTTTCAAAATAGa                                                                                                                           >  1:516027/1‑140 (MQ=255)
               aGTCCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAGCCCGGCTTCTACATCAAAAGAGCCATCCATATGGGGATCGTTAATTAACCCTTTCCAGCCGACAGTGGTACGGGGCTTTTCAAAATAGAcgcg                                                                                                                       <  2:622392/140‑1 (MQ=255)
                 tCCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAGCCCGTCTTCTACATCAAAAGAGCCATCCATATGGGGATCGTTAATTAACCCTTTCCAGCCGACAGTGGTACGGGGCTTTTCAAAATAGACGCGCa                                                                                                                     <  2:854865/140‑1 (MQ=255)
                 gcccATATTCACCAGCTCAAACAGCAATTTACCCGCGATCTGCAGGCCGGCTTCTACATCAAAAGAGCCATCCATATGGGGATCGTTAATTAACCCTTTCCAGCCGACAGTGGTTCGGGGCTTTTCAAAATAGACGCGCa                                                                                                                     <  2:106323/139‑1 (MQ=37)
                      tatTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAGCCCGGCTTCTACATCAAAAGAGCCATCCATATGGGGATCGTTAATTAACCCTTTCCAGCCGACAGTGGTACGGGGCTTTTCAAAATAGACGCGCATTAcc                                                                                                                >  1:59151/1‑140 (MQ=255)
                           aCCAGCTCAAGCAGCAATTTACGCGCGATCTGCAGCCCGGCTTCTACATCAAAAGAGCCATCCATATGGGGATCGTTAATTAACCCTTTCCAGCCGACAGTGGTACGGGGATTTTCAAAATAGAGGCGCATTACCAGATa                                                                                                           >  1:247707/1‑140 (MQ=255)
                                cTCAAGCAGCAATTTTCGCGCGGTCGGCAACCCGGCTGCTACATCAAAAGAGCCATCCATATGGGGATCGTTAATTAACCCTTTCCAGCCGACAGTGGTACGGGGCTTTTCAAAATAGACGCGCATTACCAGATAGAGGc                                                                                                      <  1:393016/140‑1 (MQ=37)
                                  cAAGCAGCAATTTACGCGCGATCTGCAGCCCGGCTTCTACATCAAAAGAGCCATCCATATGGGGATCGTTAATTAACCCTTTCCAGCCGACAGTGGTACGGGGCTTTTCAGAATGGACGCGCATTACCAGATAGAGGCTa                                                                                                    >  1:225313/1‑140 (MQ=255)
                                          aaTTTACGCGCGATCTGCAGCCCGGCTTCTACATCAAAAGAGCCATCCATATGGGGATCGTTAATTAACCCTTTCCAGCCGACAGTGGTACGGGGCTTTTCAAAATAGACGCGCATTACCAGATAGAGGCTATCGCTGAc                                                                                            <  1:42451/140‑1 (MQ=255)
                                             ttACGCGCGATCTGCAGCCCGGCTTCTACATCAAAAGAGCCATCCATATGGGGATCGTTAATTAACCCTTTCCAGCCGACAGTGGTACGGGGCTTTTCAAAATAGACGCGCATTACCAGATAGAGGCTATCGCTGACctc                                                                                         <  2:516027/140‑1 (MQ=255)
                                                 gcgcgATCTGCAGCCCGGCTTCTACATCAAAAGAGCCATCCATATGGGGATCGTTAATTAACCCTTTCCAGCCGACAGTGGTACGGGGCTTTTCAAAATAGACGCGCATTACCAGATAGAGGCTATCGCTGACCTCTGCg                                                                                     >  1:465838/1‑140 (MQ=255)
                                                     gATCTGCAGCCCGGCTTCTACATCAAAAGAGCCATCCATATGGGGATCGTTAATTAACCCTTTCCAGCCGACAGTGGTACGGGGCTTTTCAAAATAGACGCGCATTACCAGATAGAGGCTATCGCTGACCTCTGCGGCaa                                                                                 <  1:667217/140‑1 (MQ=255)
                                                          gCAGCCCGGCTTCTACAGCAAACGAGCCATCCATACGGGCATCGTTAACTAACCCTTTCCAGCCGACAGTGGTACGGGGCTTTTCAAAATAGACGCGCATTACCAGATAGAGGCTATCGCTGACCTCTGCGGCAAGGGCt                                                                            <  2:248568/140‑1 (MQ=37)
                                                             gCCCGGCTTCTACATCAAAAGAGCCATCCATATGGGGATCGTTAATTAACCCTTTCCAGCCGACAGTGGTACGGGGCTTTTCAAAATAGACCCGCATTACGAGACAGAGGCTATCGCTGACCTCTGCGGCAAGGGCTTTa                                                                         >  2:337019/1‑140 (MQ=255)
                                                                 ggCTTCTACATCAAAAGAGCCATCCATATGGGGATCGTTAATTAACCCTTTCCAGCCGACAGTGGTACGGGGCTTTTCAAAATAGACGCGCATTACCAGATAGAGGCTATCGCTGACCTCTGCGGCAAGGGCTTTAAATc                                                                     >  2:537444/1‑140 (MQ=255)
                                                                   cTTCTACATCAAAAGAGCCATCCATATGGGGATCGTTAATTAACCCTTTCCAGCCGACAGTGGTACGGGGCTTTTCAAAATAGACGCGCATTACCAGATAGAGGCTATCGCTGACCTCTGCGGCAAGGGCTTTAAATcga                                                                   >  2:190878/1‑140 (MQ=255)
                                                                    ttCTACATCAAAAGATCCATCCCTACGGGGATCGTTAATTAACCCTTTCCAGCCGACAGTGGTACGGGGCTTTTCAAAATAGACGCGCATTACCAGATAGAGGCTATCGCTGACCTCTGCGGCAAGGGCTTTAAATcgac                                                                  <  1:365818/140‑1 (MQ=255)
                                                                    ttCTACATCAAAAGATCCATCCCTAAGGGAATCGTTATTTAACCCTTTCCAGCCGACAGTGGTACGGGGCTTTTCAAAATAGACGCGCATTACCAGATAGAGGCTATCGCTGACCTCTGTGGCAAGGGCTTTAAATcgac                                                                  <  1:224226/140‑1 (MQ=37)
                                                                    ttCTACATCAAAAGAGCCATCCATATGGGGATCGTTAATTAACCCTTTCCAGCCGACAGTGGTACGGGGCTTTTCAAAATAGACGCGCATTACCAGATAGAGGCTATCGCTGACCTCTGCGGCAAGGGCTTTAAATcgac                                                                  <  1:537444/140‑1 (MQ=255)
                                                                         cATCAAAAGAGCCATCCCTATGCGGATCGTTGATTTACCCTTTCCAGCAGACAGTGGTACGGGGCTTTTCAAAATAGACGCGCATTACCAGATAGAGGCTATCGCTGACCTCTGCGGCAAGGGCTTTAAATCGACGAGCa                                                             <  2:225313/140‑1 (MQ=37)
                                                                         cATCAAAAGAGCCATCCATATGGGGATCGTTAATTAACCCTTTCCAGCCGACAGTGGTACGGGGCTTTTCAAAATAGACGCGCATTACCAGATAGAGGCTATCGCTGACCTCTGCGGCAAGGGCTTTAAATCGACGAGCa                                                             <  2:465838/140‑1 (MQ=255)
                                                                            cAAAAGAGCCATCCATATGGGGATCGTTAATTAACCCTTTCCAGCCGACAGTGGTACGGGGCTTTTCAAAATAGACGCGCATTACCAGATAGAGGCTATCGCTGACCTCTGCGGCAAGGGCTTTAAATCGACGAGCatat                                                          >  2:448227/1‑140 (MQ=255)
                                                                                   gCCATCCATATGGGGATCGTTACTTAACCGTTTCCAGCCGACAGTGGTACGGGGCTTTTCAAAATAGACGCGCATTACCAGATAGAGGCTATCGCTGACCTCTGCGGCAAGGGCTTTAAATCGACGAGCATATTCCagag                                                   <  1:337019/140‑1 (MQ=255)
                                                                                       tccatATGGGGATCGTTAATTAACCCTTTCCAGCCGACAGTGGTACGGGGCTTTTCAAAATAGACGCGCATTACCAGATAGAGGCTATCGCTGACCTCTGCGGCAAGGGCTTTAAATCGACGAGCATATTCCAGAGCAGt                                               >  2:836914/1‑140 (MQ=255)
                                                                                                ggATCGTTAATTAACCCTTTCCAGCCGACGGTGGTACGGGGCTTTTCAAAATAGACGCGCATTACCAGATAGAGGCTATCGCTGACCTCTGCGGCAAGGGCTTTAAATCGACGAGCATATTCCAGAGCAGTTTCCGGATc                                      <  1:832736/140‑1 (MQ=255)
                                                                                                 gATCGTTAATTAACCCTTTCCAGCCGACAGTGGTACGGGGCTTTTCAAAATAGACGCGCATTACCAGATAGAGGCTATCGCTGACCTCTGCGGCAAGGGCTTTAAATCGACGAGCATATTCCAGAGCAGTTTCCGGATCa                                     >  1:123455/1‑140 (MQ=255)
                                                                                                    cGTTAATTAACCCTTTCCAGCCGACAGTGGTACGGGGCTTTTCAAAATAGACGCGCATTACCAGATAGAGGCTATCGCTGACCTCTGCGGCAAGGGCTTTAAATCGACGAGCATATTCCAGAGCAGTTTCCGGATCATGa                                  >  1:688753/1‑140 (MQ=255)
                                                                                                     gTTAATTAACCCTTTCCAGCCGACAGTGGTACGGGGCTTTTCAAAATAGACGCGCATTACCAGATAGAGGCTATCGCTGACCTCTGCGGCAAGGGCTTTAAATCGACGAGCATATTCCAGAGCAGTTTCCGGATCATGaa                                 <  2:688753/140‑1 (MQ=255)
                                                                                                     gTTAATTAACCCTTTCCAGCCGACAGTGGTACGGGGCTTTTCAAAATAGACGCGCATTACCAGATAGAGGCTATCGCTGACCTCTGCGGCAAGGGCTTTAAATCGACGAGCATATTCCAGAGCAGTTTCCGGATCATGaa                                 >  1:494056/1‑140 (MQ=255)
                                                                                                      ttaattaaCCCTTTCCAGCCGACAGTGGTACGGGGCTTTTCAAAATAGACGCGCATTACCAGATAGAGGCTATCGCTGACCTCTGCGGCAAGGGCTTTAAATCGACGAGCATATTCCAGAGCAGTTTCCGGATCATGAAt                                <  1:190878/140‑1 (MQ=255)
                                                                                                            aaCCCTTTCCAGCCGACAGTGGTACGGGGCTTTTCAAAATAGACGCGCATTACCAGATAGAGGCTATCGCTGACCTCTGCGGCAAGGGCTTTAAATCGACGAGCATATTCCAGAGCAGTTTCCGGATCATGAATGGAACa                          >  1:753963/1‑140 (MQ=255)
                                                                                                             aCCCTTTCCAGCCGACAGTGGTACGGGGCTTTTCAAAATAGACGCGCATTACCAGATAGAGGCTATCGCTGACCTCTGCGGCAAGGGCTTTAAATCGACGAGCATATTCCAGAGCAGTTTCCGGATCATGAATGGAACaa                         <  1:145039/140‑1 (MQ=255)
                                                                                                              cccTTTCCAGTCGACAGTGGTACGGGGCTTTTCAAAATAGACGCGCATTACCAGATAGAGGCTATCGCTGACCTCTGCGGCAAGGGCTTTAAATCGACGAGCATATTCCAGAGCAGTTTCCGGATCATGAATGGAACAAg                        >  1:852898/1‑140 (MQ=255)
                                                                                                              cccTTTCCAGCCGACAGCGGTACGGGGCTTTCCAACATACACGCGCATCACCAGATAGAGGCTATCGCTGACCTCTGCGGCAAGGGCTTTAAATCGACGAGCATATTCCAGAGCAGTTTCCGGATCATGAATGGAACAAg                        <  2:235531/140‑1 (MQ=255)
                                                                                                                cTTTCCAGCCGACAGTGGTACGGGGCTTTTCAAAATAGACGCGCATTACCAGATAGAGGCTATCGCTGACCTCTGCGGCAAGGGCTTTAAATCGACGAGCATATTCCAGAGCAGTTTCCGGATCATGAATGGAACAAGGa                      <  2:123455/140‑1 (MQ=255)
                                                                                                                cTTTCCAGCCGACAGTGGTACGGGGCTTTTCAAAATAGACGCGCATTACCAGATAGAGGCTATCGCTGACCTCTGCGACAAGGGCTTTAAATCGACGAGCATATCCCAAAGCAGTTTGCGGACCATGAATGGAACAAGGa                      >  2:318842/1‑140 (MQ=38)
                                                                                                                 tttCCAGCCGACAGTGGTACGGGGCTTTTCAAAATAGACGCGCATTACCAGATAGAGGCTATCGCTGACCTCTGCGGCAAGGGCTTTAAATCGACGAGCATATTCCAGAGCAGTTTCCGGATCATGAATGGAACAAGGAc                     >  2:184211/1‑140 (MQ=255)
                                                                                                                   tCCAGCCGACAGTGGTACGGGGCTTTTCAAGATAGACGCGCATTACCAGATAGAGGCTATCGCTGACCTCTGCGGCAAGGGCTTTAAATCGACGAGCATATTCCAGAGCAGTTTCCGGATCATGAATGGAACAAGGAc                     >  2:535734/1‑138 (MQ=255)
                                                                                                                   tCCAGCCGACAGTGGTACGGGGCTTTTCAAAATAGACGCGCATTACCAGATAGAGGCTATCGCTGACCTCTGCGGCAAGGGCTTTAAATCGACGAGCATATTCCAGAGCAGTTTCCGGATCATGAATGGAACAAGGAc                     >  2:162163/1‑138 (MQ=255)
                                                                                                                   tCCAGCCGACAGTGGTACGGGGCTTTTCAAAATAGACGCGCATTACCAGATAGAGGCTATCGCTGACCTCTGCGGCAAGGGCTTTAAATCGACGAGCATATTCCAGAGCAGTTTCCGGATCATGAATGGAACAAGGAc                     <  1:535734/138‑1 (MQ=255)
                                                                                                                   tCCAGCCGACAGTGGTACGGGGCTTTTCAAAATAGACGCGCATTACCAGATAGAGGCTATCGCTGACCTCTGCGGCAAGGGCTTTAAATCGACGAGCATATTCCAGAGCAGTTTCCGGATCATGAATGGAACAAGGAc                     <  1:162163/138‑1 (MQ=255)
                                                                                                                       gtcGACAGTGGTACGGGGCTTTTCAAAATAGACGCGCATTACCAGATAGAGGCTATCGCTGACCTCTGCGGCAAGGGCTTTAAATCGACGAGCATATTCCAGAGCAGTTTCCGGATCATGAATGGAACAAGGACCACAta               <  2:852898/138‑1 (MQ=255)
                                                                                                                       gCCGACAGTGGTACGGGGCTTTTCAAAATAGACGCGCATTACCAGATAGAGGCTATCGCTGACCTCTGCGGCAAGGGCTTTAAATCGACGAGCATATTCCAGAGCAGTTTCCGGATCATGAATGGAACAAGGACCACAta               <  2:494056/140‑1 (MQ=255)
                                                                                                                        ccGACAGTGGTACGGGGCTTTTCAAAATAGACGCGCATTACCAGATAGAGGCTATCGCTGACCTCTGCGGCAAGGGCTTTAAATCGACGAGCATATTCCAGAGCAGTTTCCGGATCATGAATGGAACAAGGACCACAtac              >  2:460493/1‑140 (MQ=255)
                                                                                                                              gTGGTACGGGGCTTTTCGAAATGGACGCGCAATACCAGATAGAGGCTAGCGCTGACCTCTGCGGCAAGGGCTTTAAATCGACGAGCATATTCCAGAGCAGTTTCCGGATCATGAATGGAACAAGGACCACATACTACcag        <  2:285309/140‑1 (MQ=255)
                                                                                                                              gTGGTACGGGGCTTTTCCAAATAGACGCGCATTACAAGATAGAGGCTAGCGCTGACCTCTGCGGCAAGGGCTTTAAATCGACGAGCATATTCCAGAGCAGTTTCCGGATCATGAATGGAACAAGGACCACATACTACcag        <  2:247707/140‑1 (MQ=255)
                                                                                                                                ggTACGGGGCTTTTCAAAATAGACGCGCATTACCAGATAGAGGCTATCGCTGACCTCTGCGGCAAGGGCTTTAAATCGACGAGCATATTCCAGAGCAGTTTCCGGATCATGAATGGAACAAGGACCACATACTACcagca      <  1:448227/140‑1 (MQ=255)
                                                                                                                                    cGGGGCTTTTCAAAATAGACGCGCATTACCAGATAGAGGCTATCGCTGACCTCTGCGGCAAGGGCTTTAAATCGACGAGCATATTCCAGAGCAGTTTCCGGATCATGAATGGAACAAGGACCACATACTACCAGCAGACg  <  2:59151/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                      
TCCGTCGCCAGTGGCAGTCCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAGCCCGGCTTCTACATCAAAAGAGCCATCCATATGGGGATCGTTAATTAACCCTTTCCAGCCGACAGTGGTACGGGGTTTTTCAAAATAGACGCGCATTACCAGATAGAGGCTATCGCTGACCTCTGCGGCAAGGGCTTTAAATCGACGAGCATATTCCAGAGCAGTTTCCGGATCATGAATGGAACAAGGACCACATACTACCAGCAGACG  >  NC_000913/2740720‑2740991

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