Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006881 2,190,409 C→T K320K (AAG→AAA U712_10975 ← Tryptophan synthase beta chain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068812,190,4090CT95.1% 112.4 / ‑2.3 41K320K (AAG→AAAU712_10975Tryptophan synthase beta chain
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (18/21);  total (19/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

TGGCGAGTTTGAATGCTTTTGCTAACGCATGGGCAGATTCGATTGCCGGCAAAATCCCCTCTTTTTCTGACAAGAGTTTTAATGCATCCACCGCTTCTTCATCTGTGATACTGTCATAAGTGACACGGCCGCTCTTATGCAAATATGCGTGCTCCGGACCGATTCCAGGATAGTCAAGTCCTGCTGAAATGGAGTAAGGCTCAATAATTTGCCCGAACTCATCTTGAATCAGATAAGTCAATGAACCGTGAATAACGCCTAC  >  CP006881/2190276‑2190537
                                                                                                                                     |                                                                                                                                
tGGCGAGTTTGAATGCTTTCGCTAACGCATGGGCAGATTCGATTGCCGGCAAAATCCCCTCTTTTTCTGACAAAAGCTTTAATGCATCCACCGCTTCTTCATCGGTTATACTGTCATAAGTGACACGGCCGCTTTTATg                                                                                                                             >  5:305967/1‑139 (MQ=255)
 ggCGAGTTTGAATGCTTTCGCTAACGCATGGGCAGATTCGATTGCCGGCAAAATCCCCTCTTTTTCTGAAAAAAGCTTTAATGCATCCACCGCTTCTTCATCGGTTATACTGTCATAAGTGACACGGCCGCTTTTATg                                                                                                                             <  2:40423/138‑1 (MQ=255)
  gCGAGTTTGAATGCTTTCGCTACCGCATGGGCAGATTCGATTGCCGGCAAAATCCCCTCTTTTTCTGCCAAAAGCTTTAATGCATCCACCGCTTCTTCATCGGTTATACTGTCATAAGTGACACGGCCGCTTTTATGCa                                                                                                                           <  4:132443/139‑1 (MQ=255)
  gCGAGTTTGAATGCTTTCGCTAACGCATGGGCAGATTCGATTGCCGGCAAAATCCCCTCTTTTTCTGACAAAAGCTTTAATGCATCCACCGCTTCTTCATCGGTTATACTGTCATAAGTGACACGGCCGCTTTTATGc                                                                                                                            >  1:306405/1‑138 (MQ=255)
  gCGAGGTTGAATGCGTTCGCTAACGCTTGGGCAGATGCGATTGCCGGCAAAATCCCCTCTTTTTCTGACAAAAGCTTTAATCCATCCACCGCAGCTTCATCGGTTCTACTGTCATAAGTGACACGGCCGCTTTTATGCa                                                                                                                           <  3:226745/139‑1 (MQ=255)
      gTTTGAATGCTTTTGCTAACGCATGGGCAGATTCGATTGCCGGCAAAATCCCCTCTTTTTCTGACAAGAGTTTTAATGCATCCACCGCTTCTTCATCTGTGATACTGTCATAAGTGACACGGCCGCTCTTATGCAAata                                                                                                                       >  8:122311/1‑139 (MQ=255)
       tttGAATGCTTTCGCTAACGCATGGGCAGATTCGATTGCCGGCAAAATCCCCTCTTTTTCTGACAAAAGCTTTAATGCATCCACCGCTTCTTCATCGGTTATACTGTCATAAGTGACACGGCCGCTTTTTTGCAAaaat                                                                                                                      >  1:140144/1‑139 (MQ=255)
       tttGAATGCTTTCGCTAACGCATGGGCAGATTCGATTGCCGGCAAAATCCCCTCTTTTTCTGACAAAAGCTTTAATGCATCCACCGCTTCTTCATCGGTTATACTGTCATAAGTGACACGGCCGCTTTTATTGAAataa                                                                                                                      >  6:269143/1‑138 (MQ=255)
            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                                                                                     >  2:149878/1‑55 (MQ=255)
                                                                                                                           cacGGCCGCTTTTATGCAAATATGCATGCTCCGGACCGATTCCAGGATAGTCGAga                                                                                     <  1:149878/56‑2 (MQ=255)
                                                                                                                           cacGGCCGCTTTTATGCAAATATGCATGCTCAGGTCCGATTGCAGGATAGTCGAGCCCGGCTGAAATAGAGTAGGGCTCAATAATTTGCCCGAACACATCCTGAATGACATAAGTAAATGAAACGTGAATAACCCCTAc  >  8:152479/1‑139 (MQ=255)
                                                                                                                                     |                                                                                                                                
TGGCGAGTTTGAATGCTTTTGCTAACGCATGGGCAGATTCGATTGCCGGCAAAATCCCCTCTTTTTCTGACAAGAGTTTTAATGCATCCACCGCTTCTTCATCTGTGATACTGTCATAAGTGACACGGCCGCTCTTATGCAAATATGCGTGCTCCGGACCGATTCCAGGATAGTCAAGTCCTGCTGAAATGGAGTAAGGCTCAATAATTTGCCCGAACTCATCTTGAATCAGATAAGTCAATGAACCGTGAATAACGCCTAC  >  CP006881/2190276‑2190537

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: