Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,740,735 A→G L139P (CTG→CCG)  aroF ← 3‑deoxy‑D‑arabino‑heptulosonate‑7‑phosphate synthase, tyrosine‑repressible

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,740,7350AG97.7% 144.9 / NA 43L139P (CTG→CCG) aroF3‑deoxy‑D‑arabino‑heptulosonate‑7‑phosphate synthase, tyrosine‑repressible
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (22/20);  minor base T (0/1);  total (22/21)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAACCAACCGGCATGGAAAGCCCGGAGGCCATTTCACGGTGAGTTTGCGATTCCGTTGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCAGTCCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAGCCCGGCTTCTACATCAAAAGAGCCATCCATATGGGGATCGTTAATTAACCCTTTCCAGCCGACAGTGGTACGGGGTTTTTCAAAATAGACGCG  >  NC_000913/2740593‑2740874
                                                                                                                                              |                                                                                                                                               
aaaaCCAACCGGCATGGAAAGCCCGGAGGCCATTTCACGGTGAGTTTGCGATTCCGTTGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGTCCCat                                                                                                                                         <  2:441447/150‑1 (MQ=255)
 ataccaaccGCCATGGAAAGCACGGAAGCCATTTCACGGTGAGTTTGCGATTCCGTTGTACGAGCCCCAATTGCTGGCCAGCTAAACAGACCGCCCAGGTTTTGCGGGCTATTCGGATCTAACGGTTCCGTCACCAGTGGCGGTCCCat                                                                                                                                         <  2:550508/147‑1 (MQ=255)
 aaaCCAACCGGCATGGAAAGCCCGGAGGCCATTTCACGGTGAGTTTGCGATTCCGTTGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGTCCCatat                                                                                                                                       <  2:426911/151‑1 (MQ=255)
             aTGGAAAGCCCGGAGGCCATTTCACGGTGAGTTTGCGATTCCGTTGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGTCCCATATTCACCAGCTCaa                                                                                                                           <  2:830816/151‑1 (MQ=255)
             aTGGAAAGCCCGGAGGCCATTTCACGGTGAGTTTGCGATTCCGTTGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGTCCCACATTTACCCGcc                                                                                                                              >  1:640191/1‑147 (MQ=255)
                   aGCCCGGAGGCCATTTCACGGTGAGTTTGCGATTCCGTTGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGTCCCATATTCACCAGCTCAAgcagca                                                                                                                     <  1:15069/151‑1 (MQ=255)
                   aGCCCGGAGGCCATTTCACGGTGAGTTTGCGATTCCGTTGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGTCCCATATTCACCAGCTCAAgcagca                                                                                                                     <  1:78130/151‑1 (MQ=255)
                         gAGGCCATTTCACGGTGAGTTTGCGATTCCGTTGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGTCCCATATTCACCAGCTCAAGCAGCAATTTAc                                                                                                               <  2:4141/151‑1 (MQ=255)
                                   cACGGTGAGTTTGCGATTCCGTTGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGTCCCCTAATTACCCACCCCAGCAGCAAATTACGCGCGATc                                                                                                       >  1:13189/1‑149 (MQ=255)
                                          aGTTTGCGATTCCGTTGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGTCCCATATTTACCAGCTCAAACAGCAATTTACGCGCGAGCTGCAGccc                                                                                              >  2:290328/1‑150 (MQ=255)
                                                cGATTCCGTTGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGTCCCCTATTTACCAGCCCAAAAAGCAAATTTCGCGCGGTCTGCCGCCCGGc                                                                                           >  1:846125/1‑147 (MQ=255)
                                                     ccGTTGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGTCCCATATTTACCAGCCCAAGCAACAATTTTCGCGCGGTCTGCACCCCGGCtt                                                                                         >  2:237597/1‑144 (MQ=255)
                                                       gTTGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGTCCCCCTTTTAACCAACCCAACCGCAATTTTCTCGCGGTCTGCGGCCCCGCctct                                                                                       >  2:96814/1‑144 (MQ=255)
                                                       gTTGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGTCCCCTATTTACCAACCCAAACAACAATTTACGCGCGATCTTCAGCCCCGCTTCTACATc                                                                                  >  1:886360/1‑149 (MQ=255)
                                                         tGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGTCCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAGCCCGGCTTCTACATCaaa                                                                               <  1:479008/150‑1 (MQ=255)
                                                          gTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGTCCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAGCCCGGCTTCTACATCaaaa                                                                              <  1:654664/150‑1 (MQ=255)
                                                          gTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGTCCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAGCCCGGCTTCTACATCAAAag                                                                             <  2:808504/151‑1 (MQ=255)
                                                             cGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGTCCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAGCCCGGCTTCTACATCAAAAGAGc                                                                          <  1:394535/151‑1 (MQ=255)
                                                                   ccAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGTCCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAGCCCGGCTTCTACATCAAAAGAGccatcc                                                                     <  2:666519/150‑1 (MQ=255)
                                                                        tGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGTCCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAGCCCGGCTTCTACATCAAAAGAGCCATCCATATg                                                                <  1:607301/150‑1 (MQ=255)
                                                                          cTGACTAGCTCAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCTGTCCCATATTCGCCAGCTCAAGCAGCAATTTACGCGCGATCTGCAGCCCGGCTTCTACATCAAAAGAGCCATCCATATgg                                                               <  1:235775/149‑1 (MQ=255)
                                                                          cTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGTCCCATATTTCCCCGCTCCAGCAGCAACTTTCCCCCGATCTGCCGCCCCGCCTCTTCCTCACAAAAGCCCTCCATCTtggg                                                             >  2:179453/1‑151 (MQ=255)
                                                                            gACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGTCCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAGCCCGGCTTCTACATCAAAAGAGCCATCCATATGGGGa                                                            <  1:625709/150‑1 (MQ=255)
                                                                             aCCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGTCCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAGCCCGGCTTCTACATCAAAAGAGCCATCCATATGGGGAt                                                           <  2:109400/150‑1 (MQ=255)
                                                                              ccAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGTCCCCTATTCACCAGCACAAGCAGCAATTTACGCGCGGTCTGCACCCCGGCCTCTACATCAAAAAAACCCTCCATATgggggt                                                           >  1:664751/1‑147 (MQ=255)
                                                                                            cGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGTCCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAGCCCGGCTTCTACATCAAAAGAGCCATCCATATGGGGATCGTTAATTAACCCttt                                         <  2:698860/151‑1 (MQ=255)
                                                                                            cGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGTCCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAGCCCGGCTTCTACATCAAAAGAGCCATCCATATGGGGATCGTTAATTAACCCttt                                         <  2:595382/151‑1 (MQ=255)
                                                                                             gCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGTTCCATATTCAACAGCTCAAACAGCAATTTTCCCGCGGTCTCCACCCCGGCGACTAACTCATAACAACCACACATAATGGGGTCGGTTATTAACCCttt                                         >  2:180807/1‑150 (MQ=255)
                                                                                              cccAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGTCCCCTATTTACCAGCTCAAACAGCAATTTAAGCGCGATCTGCAGCCCCGCCTCTACATCAAAAAAGGCATCCCTATGGGGGTCGTTAATTAAACCctt                                         >  2:652395/1‑149 (MQ=255)
                                                                                              cccAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGTCCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAGCCCGGCTTCTACATCAAAAGAGCCATCCATATGGGGATCGTTAATTAACCCttt                                         <  1:652444/149‑1 (MQ=255)
                                                                                                cAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGTCCCATATTTCCCAGCTCAAGCAGCAATTTACGCGCGATCTTCAGCCCCGCTTCTACATCAAAAAAGCCATCCATATGGGGGTCGTTAATTAACCCTTTcc                                       >  1:592212/1‑149 (MQ=255)
                                                                                                  ggTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGTCCCCTATTTACCCGCTCAAGCAGCAATTTACCCGCGATCTGCAGCCCCGCCTCTTCATCAAAAAAGCCATCCATATTGGGGGCGGTAATTAACCCTTTTCAGcc                                   >  1:389241/1‑151 (MQ=255)
                                                                                                      ttGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGTCCCCTATTTACCAGCTCAAGCAGCAATTTTCGCGCGATCTGCAGCCCGGCTTCTACATCAAAAAAGACATCCATATTGGGGTCGTTAATTAAACCCTTTccgcccac                                 >  1:577035/1‑144 (MQ=255)
                                                                                                            gCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGTCCCCTATTTACCCGCCCAAGCAGGAATTTAAACGCGGTCCTCCGCCCGGCCTTTACATCAAAAAAGCCATCCAAATGGGGATCGTTAATTAACCCTTTTCAGCCGACAGTgg                           >  1:593443/1‑149 (MQ=255)
                                                                                                             cTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGTCCCATATTTACCAACCCCAGCAGCAATTTAAGCGCGGTCTGCAACCCGGCTTCTACATCAAAAAAGCCATCCATATGGGGGTCGGTAATTTAACCTTTTCAGCCGAAAGTgg                           >  1:56242/1‑148 (MQ=255)
                                                                                                              tATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGTTCCATATTTATCAGCTCAAACAACAATTTTTGCGCGCGCTGCGGCCCGGCCTCTTCATCAAAAAAGCCATCCCTATGGGGGTCGGTTAATTAAACCTTTCCGCCCACAGTGGGACg                       >  2:324268/1‑151 (MQ=255)
                                                                                                                         aaCGCTTCCGTCGCCAGTGGCGGTCCCCTATTTCCCCGCTCAAGCAGCAATTTTAGCGCGGTCTTCAGACCGGCTTCTACATCAAAAAAGCCCTCCCTATTGGGGTCGTTAAATTAACCTTTCCAGCCCACCGTGGGACGGGGgtttt               >  1:678662/1‑148 (MQ=255)
                                                                                                                          aCGCTTCCGTCGCCAGTGGCGGTCCCCTATTTACCAGCTCAAACAGCAATTTTCGCGCGCGCTGCCGCCCCGCCTCTTAATCCAAAAAGCCCTCCCCATGGGGGGCGGTTATTAAACCCTTCCCGCCGACAGTGGGAAGGGGgttttt              >  2:672279/1‑147 (MQ=255)
                                                                                                                          aCGCTTCCGTCGCCAGTGGCGGTCCCCTATTTACCAACCCCAACAGCAATTTATGCGCGGTCTGCAGCCCCGGCTTTTCATCAAAAAAACCATCCATATGGGGGTCGTTTATTTAAACTTTTCACCCGACAGGGGGACggggggttttaa            >  2:166906/1‑142 (MQ=255)
                                                                                                                            gCTTCCGTCGCCAGTGGCGGTCCCATATTCACCAGCTCAAGCAGCAAATTACGCGCGATCTGCAGCCCGGCTTCTACATCAAAAAAGCCATCCATATGGGGGTCGTTTATTAAACCTTTCCAGCCGACAGGGGGACGGGGGTTTTTaaaa          >  1:265887/1‑150 (MQ=255)
                                                                                                                             cTTCCGTCGCCAGTGGCGGTCCCCTTTTTCCCAGCCGAAACCGGCATTTATGCGCGGTCGGCCGCCCGGCTTTTTCATCAAAAAAGACCTCCCTTTTGGGGTTGTTAATTAACCCTTTCCAGCCGCCAGAGGGACGGGGgttttttaaa          >  2:345781/1‑144 (MQ=255)
                                                                                                                                    cGCCAGTGGCGGTCCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAGCCCGGCTTCTACATCAAAAGAGCCATCCATATGGGGATCGTTAATTAACCCTTTCCAGCCGACAGTGGTACGGGGTTTTTCAAAATAGAcgcg  <  1:288243/150‑1 (MQ=255)
                                                                                                                                        aGTCGCGGTCCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAGCCCGGCTTCTACATCAAAAGAGCCATCCACATGGGGATCGTTAATTAACCCTTTCCAGCCGACAGTGGTACGGAGTTTTTCAAAATAGAcg    <  2:832764/144‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                               
AAAACCAACCGGCATGGAAAGCCCGGAGGCCATTTCACGGTGAGTTTGCGATTCCGTTGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCAGTCCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAGCCCGGCTTCTACATCAAAAGAGCCATCCATATGGGGATCGTTAATTAACCCTTTCCAGCCGACAGTGGTACGGGGTTTTTCAAAATAGACGCG  >  NC_000913/2740593‑2740874

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: