Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006881 2,190,424 G→A H315H (CAC→CAT U712_10975 ← Tryptophan synthase beta chain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068812,190,4240GA97.4% 107.3 / NA 40H315H (CAC→CATU712_10975Tryptophan synthase beta chain
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base A (17/21);  minor base C (1/0);  total (18/22)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGCTTTTGCTAACGCATGGGCAGATTCGATTGCCGGCAAAATCCCCTCTTTTTCTGACAAGAGTTTTAATGCATCCACCGCTTCTTCATCTGTGATACTGTCATAAGTGACACGGCCGCTCTTATGCAAATATGCGTGCTCCGGACCGATTCCAGGATAGTCAAGTCCTGCTGAAATGGAGTAAGGCTCAATAATTTGCCCGAACTCATCTTGAATCAGATAAGTCAATGAACCGTGAATAACGCCTACGGTTCCTTTCGACATAGTA  >  CP006881/2190288‑2190556
                                                                                                                                        |                                                                                                                                    
aTGCTTTCGCTAACGCATGGGCAGATTCGATTGCCGGCAAAATCCCCTCTTTTTCTGACAAAAGCTTTAATGCATCCACCGCTTCTTCATCGGTTATACTGTCATAAGTGACACGGCCGCTTTTATTCAAATATGCATg                                                                                                                                    >  7:274911/1‑139 (MQ=255)
              gCATGGGCAGATTCGATTGCCGGCAAAATCCCCTCTTTTTCTGACAAGAGTTTTAATGCATCCACCGCTTCTTCATCTGTGATACTGTCATAAGTGACACGGCCGCTCTTATGCAAATATGCGTGCTCCGGACCGATTc                                                                                                                      <  7:122311/139‑1 (MQ=255)
                  gggCAGATTCGATTGCCGGCAAAATCCCCTCTTTTTCTGACAAAAGCTTTAATGCATCCACCGCTTCTTCATCTGTTATACTGTCATAAGTGACACGGCCGCTTTTATGCAAATATTCATTGTCCGGACCGATTCGAgg                                                                                                                  >  3:280094/1‑139 (MQ=255)
                  gggCAGATTCGATTGCCGGCAAAATCCCCTCTTTTTCTGACAAAAGCTTTAATGCATCCACCGCTTCTTCATCGGTTATACTGTCATAAGTGACACGGCCGCTTTTATGCAAATAATGAAGCGCCGGACCGATTCCAg                                                                                                                   >  5:146348/1‑138 (MQ=255)
                            gATTGCCGGCAAAATCCCCTCTTTTTCTGACAAAAGCTTTAATGCATCCACCGCTTCTTCATCGGTTATACTGTCATAAGTGACACGGCCGCTTTTATGCAAATATGCATGCTCCGGACCGATTCCAGGATAGTCGAga                                                                                                        <  5:167694/139‑2 (MQ=255)
                                ggcGGCAAAATCGCGTCTTTTTCTGACATAAGCTTTAATGAATCAACCGCTTCTTCATCGGTTATAGTGTAATAAGTGACACGGCCGCTTTTATGCAAATATGCATGCTCCGGACCGATTCCAGGATAGTCGAGAcc                                                                                                      <  1:205436/135‑1 (MQ=255)
                                 ccgcaaaaaGCCCCTCTTTTTCTGAAAAATGTTTTATAGCATCCCCCGGTTCTTCATCGGTTTTACTGTCATAAGTCACACGGCCGCTTTTATGCAAATATGCATGCTCCGGACCGATTCCAGGATAGTCGAGACcggc                                                                                                   <  5:52761/134‑1 (MQ=255)
                                 ccGGCAAAATCCCCTCTTTTTCTGACAAAAGCTTTAATGCATCCACCGCTTCTTCATCGGTTATACTGTCATAAGTGACACGGCCGCTTTTATGCAAATATGCATGCTCCGGACCGATTCCAGGATAGTCGAGAccgg                                                                                                    <  3:169507/138‑3 (MQ=255)
                                   ggCAAAATCCCCTCTTTTTCTGACAACAGCTTTAATCCATCCACCGCTTCTTCATCGGTAATACTGTCATCAGTGACACGGCCGCTTTTATGCAAATATGCATGCTCCGGACCGATTCCAGGATAGTCGAGACcggctg                                                                                                 <  8:57178/139‑1 (MQ=255)
                                   ggCAAAATCCCCTCTTTTTCTGACAAAAGCTTTAATGCATCCACCGCTTCTTCATCGGTTATACTGTCATAAGTGACACGGCCGCTTTTATGCAAATATGCATGCTCCGGACCGATTCCAGGATAGTCGAGACcggctg                                                                                                 <  2:306405/139‑1 (MQ=255)
                                      aaaaTCCCCTCTTTTTCTTACAAATCCTTTAATGCAACCCCCGCTTCTTCATGGTTTATACTGTCATACGTGCCACGGCCGCTTTTATGCAAATATGCATGCTCCGGACCGATTCCAGGATAGTCGAGACCGGCTGaaa                                                                                              <  5:269143/139‑1 (MQ=255)
                                      aaaaTCCCCTCTTTTTCTGACAAAAGCTTTAATGCATCCACCGCTTCTTCATCGGTTATACTGTCATAAGTGACACGGCCGCTTTTATGCAAATATGCATGCTCCGGACCGATTCCAGGATAGTCGAGACCGGCTGaaa                                                                                              <  2:12722/139‑1 (MQ=255)
                                             ccTCTTTTTCTGACAAAACCTTTAATGCATCCACCGCTTCTTCATCGGTTATACTGTCATAAGTGACACGGCCGCTTTTATGCAAATATGCATGCTCCGGACCGATTCCAGGATAGTCGAGACCGGCTGAAATAGAGTa                                                                                       <  3:210219/139‑1 (MQ=255)
                                             ccTCTTTATCTCACAAAAGCTTTATTGCCCCCACAGCTTATCCTTCGGTTAACCGGCCAAAAGTCACACGGCCGCTTTTATGCAAATATGCATGCTCCGGACCGATTCCAGGATAGTCGAGACCGGCTGAAATAGAGTa                                                                                       <  4:99989/139‑1 (MQ=255)
                                                         aCAAAAGCTTTAATGCATCCACCGCTTCTTCATCGGTTATACTGTCATAAGTGACACGGCCGCTTTTATGCAAATATGCATGCTCCGGACCGATTCCAGGATAGCCGAGACCGGCTGAAATAGAGTAGGGCTCaataat                                                                           >  3:218900/1‑139 (MQ=255)
                                                                          tCCACCGCTTCTTCATCGGTTATACTGTCATAAGTGACACGGCCGCTTTTATGCAAATATGGATGCTGCGGACCGATTGCAGGATTGTCGAGACCGGCTGAAATAGAGTAGGGCGCAATAATTTGCCCGAACTcctgcg                                                          >  4:412151/1‑134 (MQ=255)
                                                                          tCCACCGCTTCTTCATCGGTTATACTGTCATAAGTGACACGGCCGCTTTTATGCAAATATGCATGCTCCGGACCGATTCCAGGa                                                                                                                 >  1:100868/1‑84 (MQ=255)
                                                                          tCCACCGCTTCTTCATCGGTTATACTGTCATAAGTGACACGGCCGCTTTTATGCAAATATGCATGCTCCGGACCGATTCCAGGa                                                                                                                 <  2:100868/84‑1 (MQ=255)
                                                                          tCCACCGCTTCTTCATCGGTTATACTGTCATAAGTGACACGGCCGCTTTTATGCAAATATGCATGCTCCGGACCGATTCCAGGATAGTCGAGACCGGCTGAAATAGAGTAGGGCTCAATAATTTGCCCGAACTCATcct                                                          >  7:10417/1‑137 (MQ=255)
                                                                             aCCGCTTCTTCTTCGGTTAAACTGTATTAAGTGACACGGCCGCTTTTATGCAAATATGCATGCTCCGGACCGATTCCAGGATAGGCGAGACCGGCTGAAATAGAGTAGGGCTCAATAATTTGCCCGAACTCATCCTGaa                                                       <  4:199815/139‑1 (MQ=255)
                                                                                    cttcATCGGGAATACAGGCATAAGTGCCACGGCCGCTTTAATGCAAATATGCATGCTCCGGACCGATTCCAGCATAGGCGAGACCGGCTGAAATAGAGTAGGGCTCAATAATTTGCCCGAACTCATCCTGAATGAGATa                                                <  2:233334/139‑1 (MQ=255)
                                                                                                tACTGTCATAAGTGACACGGCCGCTTTTATGCAAATATGCATGCTCCGGACCGATTCCAGGATAGTCGAGACCGGCTGAAATAGAGTAGGGCTCAATAATTTGCCCGAACTCATCCTGAATGAGATAAGTCAATGAAc                                     <  8:10417/138‑1 (MQ=255)
                                                                                                tACTGTCATAAGTGACACGGCCGCTTTTATGCAAATATGCATGCGCCGGACCGATTCCAGGATAGTCGAGACCGGCTGAAATAGAGTAGGGCTCAATAATTTGCCCGAACTCATCCTGAATGAGATAAGTCAATGAAc                                     <  6:305967/138‑1 (MQ=255)
                                                                                                      cATAAGTGACACGGCCGCTTTTATGCAAATATGCCTGCTCCGGCCCCAATCCAGGATAGTAGAGCCCGGCGGAAATAGAGTAGTGCGCAAAAATTTGCCCGGACGCATCCGGCATGAGAGAAGTCAAAGAACCggtaa                               >  6:206862/1‑134 (MQ=255)
                                                                                                       aTACGTGCCACGGCCGCTTTTATGCAAATATGCATGCTCCGGCCCGATTCCAGGATAGTCGAGACCGGCTGAAATAGAGTAGGGCTCAATAATTTGCCCGAACTCATCCTGAATGAGATAAGTCAATGAACCGTGAAt                              <  5:206862/138‑1 (MQ=255)
                                                                                                       aTAAGTGACACGGCCGCTTTTATGCAAATATGCATGCTCCGGACCGATTCCAGGATAGTCGAGACCGGCTGAAATAGAGTAGGGCTCAATAATTTGCCCGAACTCATCCTGAATGAGATAAGTCAATGAACCGTGAAt                              >  8:82931/1‑138 (MQ=255)
                                                                                                       aTAAGTGACACGGCCGCTTTTAATCAAAAATTGATGCTCCGGACCGATACCACGCTAGTACCGACCGGCTGTAATATACTAGGGGTCAATAATATCAGCGCAAGCAACCAGAATGAGATAAGTCAATTAACCGAGaaa                              >  4:437794/1‑137 (MQ=255)
                                                                                                       aTAAGTGACACGGCCGCGTTTATGCAAATATGCATGCGACGGCCCGATTCCCGGATAGTCGAGACCGGCAGAAATAGAGTAGGGCGCTATATTTTGACCGTACTAAGCCTGCTTGAGAATAGTCAATGAACCGTGaca                              >  6:104618/1‑136 (MQ=255)
                                                                                                          aGTGACACGGCCGCTTTTATGCAAATATTCATGCTCCGGACCGATTCCAGGATATTCGAGACCGGCTGAAATAGAGTAGGGCTCAATAATTTTCCCGAACACATCCTGAATTAGATAAGTAAATTAAACCGGAATAAcc                          >  8:329988/1‑138 (MQ=255)
                                                                                                          aGTGACACGGCCGCTTTTATGCAAATATGCATGCTCCGGACCGATTCCAGGATAGTCGAGACCGGCTGAAATAGAGTAGGGCTCAATAATTTGCCCGAACTCATCCTGAATGAGATAAGTCAATGAACCGTGAATAAc                           <  7:165399/138‑1 (MQ=255)
                                                                                                               cacGGCCGCTTTTATGCAAATATGCATGCTCCGGACCGATTCCAGGATAGTCGAga                                                                                                        <  1:149878/56‑2 (MQ=255)
                                                                                                               cacGGCCGCTTTTATGCAAATATGCATGCTCCGGACCGATTCCAGGATAGTCGAga                                                                                                        >  2:149878/1‑55 (MQ=255)
                                                                                                               cacGGCCGCTTTTATGCAAATATGCATGCTCAGGTCCGATTGCAGGATAGTCGAGCCCGGCTGAAATAGAGTAGGGCTCAATAATTTGCCCGAACACATCCTGAATGACATAAGTAAATGAAACGTGAATAACCCCTAc                     >  8:152479/1‑139 (MQ=255)
                                                                                                                           tATGCAAATATGCATGCTCCGGACCGATTCCAGGATAGTCGAGACCGGCTGAAATAGAGTAGGGCTCAATAATTTGCCCGAACTCATCCTGAATGAGATAAGTCAATGAACCGTGAATAACCCCTACGGTTCCTTTCg          <  2:140144/138‑1 (MQ=255)
                                                                                                                                  atatGCATGCTCCGGACCGATTCCAGGATAGTCGAGCCCGGCTGAAATAGAGTAGGGCTCAATAATTTGCCCGAACTCATCCTGAATGAGATAAGTCAATGAACCGTGAATAACCCCTACGGTTCCTTTCGAAAtagtg  >  8:87431/1‑138 (MQ=255)
                                                                                                                                  atatGCATGCTCCGGACCGATTCCAGGATAGTCGAGAccgg                                                                                                    >  6:413888/1‑39 (MQ=255)
                                                                                                                                  atatGCATGCTCCGGACCGATTCCAGGATAGTCGAGAccgg                                                                                                    <  5:413888/41‑3 (MQ=255)
                                                                                                                                  atatGCATGCTCCGGACCGATTCCAGGATAGTCGAGACCGGCTGAAATGGAGTAGGGCTCAATAATTTGCCCGAACTCATCCTGAATGAGATAAGTCAATGAACCGTGAATAACCCCTACGGTTCCTTTCGAAAtagtg  >  3:281060/1‑138 (MQ=255)
                                                                                                                                  atatGCATGCTCCGGACCGATTCCAGGATAGTCGAGACCGGCTGAAATAGAGTAGGGCTCAATAATTTGCCCGAACTCATCCTGAATGAGATAAGTCAATGAACCg                                   <  2:127228/106‑1 (MQ=255)
                                                                                                                                  atatGCATGCTCCGGACCGATTCCAGGATAGTCGAGACCGGCTGAAATAGAGTAGGGCTCAATAATTTGCCCGAACTCATCCTGAATGAGATAAGTCAATGAACCg                                   >  1:127228/1‑106 (MQ=255)
                                                                                                                                        |                                                                                                                                    
ATGCTTTTGCTAACGCATGGGCAGATTCGATTGCCGGCAAAATCCCCTCTTTTTCTGACAAGAGTTTTAATGCATCCACCGCTTCTTCATCTGTGATACTGTCATAAGTGACACGGCCGCTCTTATGCAAATATGCGTGCTCCGGACCGATTCCAGGATAGTCAAGTCCTGCTGAAATGGAGTAAGGCTCAATAATTTGCCCGAACTCATCTTGAATCAGATAAGTCAATGAACCGTGAATAACGCCTACGGTTCCTTTCGACATAGTA  >  CP006881/2190288‑2190556

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