Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,321,083 G→A A101V (GCG→GTG)  trpD ← fused glutamine amidotransferase (component II) of anthranilate synthase/anthranilate phosphoribosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,321,0830GA94.1% 122.5 / ‑6.1 34A101V (GCG→GTG) trpDfused glutamine amidotransferase (component II) of anthranilate synthase/anthranilate phosphoribosyl transferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (12/20);  minor base C (0/2);  total (12/22)
Fisher's exact test for biased strand distribution p-value = 5.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGGCGTTGATGGTTAAACCGGCCGGAATGTTACTGCCAACCAGCGAGTGATAACGCGCCACCGGCAGCGGGTTTGTTAATCCGGCAAACATCGCCTGACCGTCATGTTCAATGCTGGAGGCTTTACCGTGGAGAATTTCGCCCGCCTGACCGACATAGCCCCCGTAAGCTTCGACAATCGCCTGATGTCCGAGGCAAATGCCAATAATGGGCAGCTTGCCACGCAAGCGGGTGAGGAGTTCCGGCATACAACCGGCTTCGCTCGGCACACCGGGGCCAGG  >  NC_000913/1320939‑1321219
                                                                                                                                                |                                                                                                                                        
tGGGCGTTGATGGTTAAACCGGCCGGAATGTTACTGCCAACCAGCGAGTGATAACGCGCCACCGGCAGCGGGTTTGTTAATCCGGCAAACATCGCCTGACCGTCATGTTCAATGCTGGAGGCTTTACCGTGGAGAATTTCGCCCACCTGAc                                                                                                                                    >  2:253887/1‑151 (MQ=255)
        gATGGTTAAACCGGCCGGAATGTTACTGCCAACCAGCGAGTGATAACGCGCCACCGGCAGCGGGTTTGTTAATCCGGCAAACATCGCCTGACCGTCATGTTCAATGCTGGAGGCTTTACCGTGGAGAATTTCGCCCACCTGACCGACAt                                                                                                                              >  2:597431/1‑149 (MQ=255)
              tAAACCGGCCGGAATGTTACTGCCAACCAGCGAGTGATAACGCGCCACCGGCAGCGGGTTTGTTAATCCGGCAAACATCGCCTGACCGTCATGTTCAATGCTGGAGGCTTTACCGTGGAGAATTTCGCCCACCTGACCGACATAGccccc                                                                                                                       <  1:420704/150‑1 (MQ=255)
               aaaCCGGCCGGAATGTTACTGCCAACCAGCGAGTGATAACGCGCCACCGGCAGCGGGTTTGTTAATCCGGCAAACATCGCCTGACCGTCATGTTCAATGCTGGAGGCTTTACCGTGGAGAATTTCGCCCACCTGACCGACATAGCCCCCGt                                                                                                                     <  1:330550/151‑1 (MQ=255)
                   cggccggAATGTTACTGCCAACCAGCGAGTGATAACGCGCCACCGGCAGCGGGTTTGTTAATCCGGCAAACATCGCCTGACCGTCATGTTCAATGCTGGAGGCTTTACCGTGGAGAATTTCGCCCACCTGACCGACATAGCCCCCGTAAGc                                                                                                                 >  1:411858/1‑151 (MQ=255)
                     gccggAATGTTACTGCCAACCAGCGAGTGATAACGCGCCACCGGCAGCGGGTTTGTTAATCCGGCAAACATCGCCTGACCGTCATGTTCAATGCTGGAGGCTTTACCGTGGAGAATTTCGCCCACCTGACCGACATAGCCCCCGTAAGCtt                                                                                                               <  2:233408/151‑1 (MQ=255)
                      ccggAATGTTACTGCCAACCAGCGAGTGATAACGCGCCACCGGCAGCGGGTTTGTTAATCCGGCAAACATCGCCTGACCGTCATGTTCAATGCTGGAGGCTTTACCGTGGAGAATTTCGCCCACCTGACCGACATAGCCCCCGTAAGCTTc                                                                                                              <  1:297188/151‑1 (MQ=255)
                              ttACTGCCAACCAGCGAGTGATAACGCGCCACCGGCAGCGGGTTTGTTAATCCGGCAAACATCGCCTGACCGTCATGTTCAATGCTGGAGGCTTTACCGTGGAGAATTTCGCCCACCTGACCGACATAGCCCCCGTAAGCTTCGACAATCg                                                                                                      >  1:582552/1‑151 (MQ=255)
                                   gCCAACCAGCGAGTGATAACGCGCCACCGGCAGCGGGTTTGTTAATCCGGCAAACATCGCCTGACCGTCATGTTCAATGCTGGAGGCTTTACCGTGGAGAATTTCGCCCACCTGACCGACATAGCCCCCGTAAGCTTCGACAATCGCCTg                                                                                                  <  2:137135/150‑1 (MQ=255)
                                   gCCAACCAGCGAGTGATAACGCGCCACCGGCAGCGGGTTTGTTAATCCGGCAAACATCGCCTGACCGTCATGTTCAATGCTGGAGGCTTTACCGTGGAGAATTTCGCCCACCTGACCGACATAGCCCCCGTAAGCTTCGACAATCGCCTg                                                                                                  >  2:130227/1‑150 (MQ=255)
                                      aaCCAGCGAGTGATAACGCGCCACCGGCAGCGGGTTTGTTAATCCGGCAAACATCGCCTGACCGTCATGTTCAATGCTGGAGGCTTTACCGTGGAGAATTTCGCCCACCTGACCGACATAGCCCCTGTAAGCTTCGACAATCGCCTGATGt                                                                                              >  1:516443/1‑151 (MQ=255)
                                          aGCGAGTGATAACGCGCCACCGGCAGCGGGTTTGTTAATCCGGCAAACATCGCCTGACCGTCATGTTCAATGCTGGAGGCTTTACCGTGGAGAATTTCGCCCACCTGACCGACATAGCCCCCGTAAGCTTCGACAATCGCCTGATGTCCGa                                                                                          <  1:251290/151‑1 (MQ=255)
                                           gCGAGTGATAACGCGCCACCGGCAGCGGGTTTGTTAATCCGGCAAACATCGCCTGACCGTCATGTTCAATGCTGGAGGCTTTACCGTGGAGAATTTCGCCCACCTGACCGACATAGCCCCCGTAAGCTTCGACAATCGCCTGATGTCCg                                                                                           >  1:498820/1‑149 (MQ=255)
                                                 gATAACGCGCCACCGGCAGCGGGTTTGTTAATCCGGCAAACATCGCCTGACCGTCATGTTCAATGCTGGAGGCTTTACCGTGGAGAATTTCGCCCACCTGACCGACATAGCCCCCGTAAGCTTCGACAATCGCCTGATGTCCGAGGCaaa                                                                                    <  2:546861/150‑1 (MQ=255)
                                                    aaCGCGCCACCGGCAGCGGGTTTGTTAATCCGGCAAACATCGCCTGACCGTCATGTTCAATGCTGGAGGCTTTACCGTGGAGAATTTCGCCCACCTGACCGACATAGCCCCCGTGAGCTTCGACAATCGCCTGATGTCCGAGGCAAATGcc                                                                                <  1:148039/151‑1 (MQ=255)
                                                          ccACCGGCAGCGGGTTTGTTAATCCGGCAAACATCGCCTGACCGTCATGTTCAATGCTGGAGGCTTTACCGTGGAGAATTTCGCCCACCTGACCGACATAGCCCCTGTAAGCTTCGACAATCGCCTGATGTCCGAGGCAAATGCCaataat                                                                          <  2:516432/151‑1 (MQ=255)
                                                            aCCGGCAGCGGGTTTGTTAATCCGGCAAACATCGCCTGACCGTCATGTTCAATGCTGGAGGCTTTACCGTGGAGAATTTCGCCCACCTGACCGACATAGCCCCCGTAAGCTTCGACAATCGCCTGATGTCCGAGGCAAATGCCAATAATgg                                                                        >  2:196408/1‑151 (MQ=255)
                                                              cGGCAGCGGGTTTGTTAATCCGGGAAAGATCGCGTGGCCGGCATGGTCAATGGTGGAGGCTTTAGCGTGGAGAATTTCGCCCCCCTGACCGGCAAAGCGTACGGAAGCGTAGTCAAGCGCCTGATGTCGGAGGCAAATGGCAATAATgg                                                                        <  2:564569/149‑1 (MQ=255)
                                                                  aGCGGGTTTGTTAATCCGGCAAACATCGCCTGACCGTCATGTTCAATGCTGGAGGCTTTACCGTGGAGAATTTCGCCCACCTGACCGACATAGCCCCCGTAAGCTTCGACAATCGCCTGATGTCCGAg                                                                                         >  2:95690/1‑128 (MQ=255)
                                                                  aGCGGGTTTGTTAATCCGGCAAACATCGCCTGACCGTCATGTTCAATGCTGGAGGCTTTACCGTGGAGAATTTCGCCCACCTGACCGACATAGCCCCCGTAAGCTTCGACAATCGCCTGATGTCCGAg                                                                                         <  1:95690/128‑1 (MQ=255)
                                                                          tGTTAATCCGGCAAACATCGCCTGACCGTCATGTTCAATGCTGGAGGCTTTACCGTGGAGAATTTCGCCCACCTGACCGACATAGCCCCCGTAAGCTTCGACAATCGCCTGATGTCCGAGGCAAATGCCAATAATGGGCAGCTTGCCACg                                                           <  2:411851/150‑1 (MQ=255)
                                                                          tGTTAATCCGGCAAACATCGCCTGACCGTCATGTTCAATGCTGGAGGCTTTACCGTGGAGAATTTCGCCCACCTGACCGACATAGCCCCCGTAAGCTTCGACAATCGCCTGATGTCCGAGGCAAATGCCAATAATGGGCAGCTTGCCACg                                                           <  2:498810/150‑1 (MQ=255)
                                                                              aaTCCGGCAAACATCGCCTGACCGTCATGTTCAATGCTGGAGGCTTTACCGTGGAGAATTTCGCCCACCTGACCGACATAGCCCCCGTAAGCTTCGACAATCGCCTGATGTCCGAGGCAAATGCCAATAATGGGCAGCTTGCCACGCAAg                                                       <  2:625297/150‑1 (MQ=255)
                                                                              aaTCCGGCAAACATCGCCTGACCGTCATGTTCAATGCTGGAGGCTTTACCGTGGAGAATTTCGCCCACCTGACCGACATAGCCCCCGTAAGCTTCGACAATCGCCTGATGTCCGAGGCAAATGCCAATAATGGGCAGCTTGCCACGCAAGc                                                      >  1:28662/1‑151 (MQ=255)
                                                                                 ccGGCAAACATCGCCTGACCGTCATGTTCAATGCTGGAGGCTTTACCGTGGAGAATTTCGCCCACCTGACCGACATAGCCCCCGTAAGCTTCGACAATCGCCTGATGTCCGAGGCAAATGCCAATAATGGGCAGCTTGCCACGCAAGCgg                                                    <  1:384945/150‑1 (MQ=255)
                                                                                  cGGCAAACATCGCCTGACCGTCATGTTCAATGCTGGAGGCTTTACCGTGGAGAATTTCGCCCACCTGACCGACATAGCCCCCGTAAGCTTCGACAATCGCCTGATGTCCGAGGCAAATGCCAATAATGGGCAGCTTGCCACGCAAGCGGGt                                                  <  2:492966/151‑1 (MQ=255)
                                                                                     cAAACATCGCCTGACCGTCATGTTCAATGCTGGAGGCTTTACCGTGGAGAATTTCGCCCCCCTGACCGACATAGCCCCCGTAAGCTTCGCCAATCGCCTGATCTCCGAGGCAAATGCCAACCATGGGCAGCTTGCCACGCAAGCGGGTgag                                               <  2:557669/151‑1 (MQ=255)
                                                                                          aTCGCCTGACCGTCATGTTCAATGCTGGAGGCTTTACCGTGGAGAATTTCGCCCACCTGACCGACATAGCCCCCGTAAGCTTCGACAATCGCCTGATGTCCGAGGCAAATGCCAATAATGGGCAGCTTGCCACGCAAGCGGGTGAGGAGtt                                          <  2:110944/151‑1 (MQ=255)
                                                                                            cGCCTGACCGTCATGTTCAATGCTGGAGGCTTTACCGTGGAGAATTTCGCCCACCTGACCGACATAGCCCCCGTAAGCTTCGACAATCGCCTGATGTCCGAGGCAAATGCCAATAATGGGCAGCTTGCCACGCAAGCGGGTGAGGAGTTcc                                        <  2:229171/151‑1 (MQ=255)
                                                                                                 gACCGTCATGTTCAATGCTGGAGGCTTTACCGTGGAGAATTTCGCCCACCTGACCGACATAGCCCCCGTAAGCTTCGACAATCGCCTGATGTCCGAGGCAAATGCCAATAATGGGCAGCTTGCCACGCAAGCGGGTGAGGAGTTCCGGCAt                                   >  2:506631/1‑151 (MQ=255)
                                                                                                       cATGTTCAATGCTGGAGGCTTTACCGTGGAGAATTTCGCCCACCTGACCGACATAGCCCCCGTAAGCTTCGACAATCGCCTGATGTCCGAGGCAAATGCCAATAATGGGCAGCTTGCCACGCAAGCGGGTGAGGAGTTCCGGCATACAAcc                             <  2:90726/151‑1 (MQ=255)
                                                                                                                 gCTGGAGGCTTTACCGTGGAGAATTTCGCCCACCTGACCGACATAGCCCCCGTAAGCTTCGACAATCGCCTGATGTCCGAGGCAAATGCCAATAATGGGCAGCTTGCCACGCAAGCGGGTGAGGAGTTCCGGCATACAACCGGCTTCGCTc                   <  2:564780/151‑1 (MQ=255)
                                                                                                                               cGTGGAGAATTTCGCCCACCTGACCGACATAGCCCCCGTAAGCTTCGACAATCGCCTGATGTCCGAGGCAAATGCCAATAATGGGCAGCTTGCTACGCAAGCGGGTGAGGAGTTCCGGCATACAACCGGCTTCGCTCGGCACACCGGGGcc     <  2:253510/151‑1 (MQ=255)
                                                                                                                                   gagaATTTCGCCCACCTGACCGACATAGCCCCCGTAAGCTTCGACAATCGCCTGATGTCCGAGGCAAATGCCAATAATGGGCAGCTTGCCACGCAAGCGGGTGAGGAGTTCCGGCATACAACCGGCTTCGCTCGGCACACCGGGGCCAgg  >  1:93279/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                        
TGGGCGTTGATGGTTAAACCGGCCGGAATGTTACTGCCAACCAGCGAGTGATAACGCGCCACCGGCAGCGGGTTTGTTAATCCGGCAAACATCGCCTGACCGTCATGTTCAATGCTGGAGGCTTTACCGTGGAGAATTTCGCCCGCCTGACCGACATAGCCCCCGTAAGCTTCGACAATCGCCTGATGTCCGAGGCAAATGCCAATAATGGGCAGCTTGCCACGCAAGCGGGTGAGGAGTTCCGGCATACAACCGGCTTCGCTCGGCACACCGGGGCCAGG  >  NC_000913/1320939‑1321219

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