Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,322,453 G→A 8.0% S165L (TCA→TTA)  trpE ← component I of anthranilate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,322,4530GA 8.0% 149.0 / 3.3 50S165L (TCA→TTA) trpEcomponent I of anthranilate synthase
Reads supporting (aligned to +/- strand):  ref base G (24/22);  new base A (2/2);  total (26/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.06e-01

TTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTCTTCCTTCGGTACATTCAACAGATTCTGCAATAAACGGAAAGCGTCAAAAACCGAAAGGGAGCATAAGCG  >  NC_000913/1322309‑1322598
                                                                                                                                                |                                                                                                                                                 
ttcAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAg                                                                                                                                               >  1:75243/1‑149 (MQ=255)
     gCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCg                                                                                                                                          >  5:318537/1‑149 (MQ=255)
     gCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCg                                                                                                                                          <  2:75243/149‑1 (MQ=255)
      cGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCgg                                                                                                                                         <  6:318537/149‑1 (MQ=255)
       gAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGTCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGt                                                                                                                                        <  3:11595/149‑1 (MQ=255)
       gAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCgg                                                                                                                                         >  3:206918/1‑148 (MQ=255)
        aGCAGTGAGACGGTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGt                                                                                                                                        >  3:10758/1‑148 (MQ=255)
                       tttttCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATccc                                                                                                                         <  8:351495/148‑1 (MQ=255)
                       tttttCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCg                                                                                                                        <  8:207685/149‑1 (MQ=255)
                             ttcttcATTCGGATCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTc                                                                                                                                >  5:85682/1‑135 (MQ=255)
                             ttcttcATTCGGATCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTc                                                                                                                                <  6:85682/135‑1 (MQ=255)
                             ttcttcATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCcac                                                                                                                    >  2:6350/1‑147 (MQ=255)
                                ttcATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGt                                                                                                               >  4:270215/1‑149 (MQ=255)
                                ttcATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGt                                                                                                               >  8:184812/1‑149 (MQ=255)
                                    ttCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTg                                                                                                                                            >  7:366933/1‑116 (MQ=255)
                                    ttCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTg                                                                                                                                            <  8:366933/116‑1 (MQ=255)
                                    ttCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCAt                                                                                                            >  1:188104/1‑148 (MQ=255)
                                          gCAAACCGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGaa                                                                                                     >  4:310808/1‑149 (MQ=255)
                                             aaCAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAg                                                                                                  <  4:206918/149‑1 (MQ=255)
                                                                      tttttttCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGt                                                                         >  6:361499/1‑149 (MQ=255)
                                                                       gaaacgcTGATGGTGATTGACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGtt                                                                        >  7:379239/7‑149 (MQ=255)
                                                                         ttttCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCTGCCGAAGAACATGGCTTCTCGttc                                                                       <  4:258421/148‑1 (MQ=255)
                                                                         ttttCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGttct                                                                      >  5:17338/1‑149 (MQ=255)
                                                                          tttCTGATTGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGttct                                                                      <  1:241355/148‑1 (MQ=255)
                                                                          tttCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGttct                                                                      <  7:184812/148‑1 (MQ=255)
                                                                              tGATGGTCAATCACCATCAGCTTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTcttcctt                                                                 <  6:146366/149‑1 (MQ=255)
                                                                                aTGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCt                                                                               <  4:76557/133‑1 (MQ=255)
                                                                                aTGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCt                                                                               >  3:76557/1‑133 (MQ=255)
                                                                                aTGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTcttcctt                                                                 <  5:289521/147‑1 (MQ=255)
                                                                                           ccATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTTCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTCTTCCTTCGGTACATTCAAc                                                    <  6:17338/149‑1 (MQ=255)
                                                                                           ccATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTCTTCCTTCGGTACATTc                                                       >  4:292806/1‑146 (MQ=255)
                                                                                           ccATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTCTTCCTTCGGTACATTc                                                       <  3:292806/146‑1 (MQ=255)
                                                                                            cATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTCTTCCTTCGGTACATTCAACa                                                   <  5:361499/149‑1 (MQ=255)
                                                                                            cATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCAGCTGACAGTTGCGGTAAATCTTCAAATCCAGCCACAAGGGCATAAGAGAACAGGCCGCCGACGAACATGGCTTCTCGTTCTTCCTTCGGTACATTCaa                                                     >  1:301787/1‑147 (MQ=255)
                                                                                                    tttCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTCTTCCTTCGGTACATTCAACAGATTCTg                                            >  6:177835/1‑148 (MQ=255)
                                                                                                    tttCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTCTTCCTTCGGTACATTCAACAGATTCTg                                            <  5:177835/148‑1 (MQ=255)
                                                                                                     ttCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTCTTACTTCGGTACATTCAACAGATTCTGc                                           >  4:351304/1‑148 (MQ=255)
                                                                                                      tCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTCTTCCTTCGGTACATTCAACAGATTCTGCa                                          >  7:166973/1‑148 (MQ=255)
                                                                                                         gCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTCTTCCTTCGGTACATTCAACAGATTCTGCAATa                                       >  7:392171/1‑148 (MQ=255)
                                                                                                             aaTAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTCTTCCTTCGGTACATTCAACAGATTCTGCAATAAACg                                   >  8:285998/2‑148 (MQ=255)
                                                                                                                    aCAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTCTTCCTTCGGTACATTCAACAGATTCTGCAATAAACGGAAAGCGt                           <  1:6350/149‑1 (MQ=255)
                                                                                                                               ggCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTCTTCCTGCGGTACATTCACCAGATTCTGCAATAAACGGAAAGCGTCAAAAACCGa                 >  5:320087/1‑148 (MQ=255)
                                                                                                                                gCAGTTATTTTCCGCTGCCAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGCGAACAGGCCGCCGAAGAACATGGCTTCTCGTTCTTCCTTCGGTACATTCAACAGATTCTGCAATAAACGGAAAGCGTCAAAAACCGaaa               <  4:10758/149‑1 (MQ=255)
                                                                                                                                    ttATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCCCAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTCTTCCTTCGGTACATTCAACAGATTCTGCAATAAACGGAAAGCGTCAAAAACCGAAAgg             >  2:210033/1‑147 (MQ=255)
                                                                                                                                    ttATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTCTTCCTTCGGTACATTCAACAGATTCTGCAATAAACGGAAAGCGTCAAAAACCGAAAGGGa           <  7:285998/149‑1 (MQ=255)
                                                                                                                                      aTTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTCTTCCTTCGGTACATTCAACAGATTCTGCAATAAACGGAAAGCGTCAAAAACCGAAAGGGAg          <  3:351304/148‑1 (MQ=255)
                                                                                                                                      aTTTTCCGCTAACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTCTTCCTTCGGTACTTTCAACAGATTCTGCAATAAACGGAAAGCGTCTAAAACCGATAGGGAg          >  2:367077/1‑148 (MQ=255)
                                                                                                                                      aTTTTCCGCTAACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTCTTCCTTCGGTACATTCAACAGATTCTGCAATAAACGGAAAGCGTCAAAAACCGAAAGGGAg          >  2:361740/1‑148 (MQ=255)
                                                                                                                                         ttCCGCTAACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTCTTCCTTCGGTACATTCAACAGATTCTGCAATAAACGGAAAGCGTCAAAAACCGAAAGGGAGCAt       <  1:367077/148‑1 (MQ=255)
                                                                                                                                         ttCCGCTAACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTCTTCCTTCGGTACATTCAACAGATTCTGCAATAAACGGAAAGCGTCAAAAACCGAAAGGGAGCAt       <  1:361740/148‑1 (MQ=255)
                                                                                                                                             gccgACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTCTTCCTTCGGTACATTCAACAGATTCTGCAATAAACGGAAAGCGTCAAAAACCGAAAGGGAGCATAAGCg  >  3:166620/4‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
TTCAGGCGAGCAGTGAGACGTTGTTTTTCTTCTTCATTCGGAGCAAACAGGCTGGCCTGAATACGGGTGCTTTTTTTCTGATGGTCAATCACCATCAGCGTTTCAGCGAGATAAAAACAGAAATCAGGGCAGTTATTTTCCGCTGACAGTTGCGGTAAATCTTCAAATCCCGCCACAAGGTCATAAGAGAACAGGCCGCCGAAGAACATGGCTTCTCGTTCTTCCTTCGGTACATTCAACAGATTCTGCAATAAACGGAAAGCGTCAAAAACCGAAAGGGAGCATAAGCG  >  NC_000913/1322309‑1322598

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑