Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,532,103 G→C 11.5% C280W (TGC→TGG cysM ← cysteine synthase B (O‑acetylserine sulfhydrolase B)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,532,1030GC11.5% 98.2 / 2.5 34C280W (TGC→TGGcysMcysteine synthase B (O‑acetylserine sulfhydrolase B)
Reads supporting (aligned to +/- strand):  ref base G (15/14);  new base C (2/2);  total (18/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCC  >  CP009273/2531849‑2532321
                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                          
cATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGccaccacca                                                                                                                                                                                                         >  1:92769/1‑274 (MQ=255)
                         aaaaaaacaggggtGTTTAATGCCGGTGGAGGAGTAAACGCCTTATCCGGCCCACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCAc                                                                                                                                                                                <  2:348576/262‑1 (MQ=38)
                                              gCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAACGATCGCCACGATCCCAGAAGATCGCCACCACCACCCCCTCAAGGTTAACTTTTGCCACCCCCAGTGccactggccctt                                                                                                                                                            >  1:348582/1‑263 (MQ=37)
                                                                                             tacAAAATCTTGCAGCAACCTTGTAGGCCTGATCAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGtt                                                                                                            <  2:92768/273‑1 (MQ=255)
                                                                                                                          tGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCGCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAg                                                                                                                                                 >  1:102491/1‑209 (MQ=255)
                                                                                                                              aaGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCGCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGa                                                                                                                                             <  2:102490/209‑1 (MQ=255)
                                                                                                                                           cAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGGTCGCCACGATCTCAGATTATCGCCACCACCACCGCGCCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCCCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCCCACCACCAGTTCGCGCACGGTGTTTTCCGCATCGCTCTGATTAACCTCCATCAcc                                                                >  2:326206/1‑273 (MQ=255)
                                                                                                                                           cAGGCAGTTTTGCGCTGGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATTCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATAGCCACGATCGCAgat                                                                                                                                                                                                                       >  1:404571/1‑122 (MQ=255)
                                                                                                                                               cAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCACCGCCCCCTGGCTAAAATGCTCTTCCCCAAGCACCCCGGTAGAAAGGTAGCGATCGCCACGATCCCATATGATc                                                                                                                                                                                                                   <  2:404564/122‑1 (MQ=255)
                                                                                                                                                  ttttGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCaccaccacc                                                                                                                                                                                                       >  2:720022/1‑131 (MQ=255)
                                                                                                                                                      gCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGt                                                                                                                                                                                                   <  1:720032/131‑1 (MQ=255)
                                                                                                                                                                          aTGCTATTAAACCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCGCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGc                               <  1:326212/275‑1 (MQ=255)
                                                                                                                                                                              tATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGccacc                                                                                                                                                                                                             >  2:21107/1‑97 (MQ=255)
                                                                                                                                                                                  aaTCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGccaccacca                                                                                                                                                                                                         <  1:21107/97‑1 (MQ=255)
                                                                                                                                                                                  aaTCCTTAAACCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGANCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATACTc                       >  1:601716/1‑275 (MQ=255)
                                                                                                                                                                                    tCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCGCCTCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGc                               <  1:446928/265‑1 (MQ=255)
                                                                                                                                                                                        tAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCTAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCATGGTTAGCTTTTGCCACCCGCAGTGCTCCCGCAACCGCGCCTCCGGAGCTGACTCCACAGACTATTCCGTCCCGCACCCCCAGTTCGCGCCTGGTGTTTTCCGCATTGCCCTGACGAATACCCAGCACCTAACCCACCAGGGAAGCGTAGAAACTCACCGCCAAATATTccnagc                  >  2:150553/1‑270 (MQ=255)
                                                                                                                                                                                        tAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACAACCCCCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGAc                                                                                                                                            >  2:326483/1‑152 (MQ=255)
                                                                                                                                                                                        tAAATCCCCGCCCCCTGGCGAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTACTTCCCGAACCGCCAGTTCGCGCCTGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGAAATTCCgacg                  >  1:585685/1‑271 (MQ=255)
                                                                                                                                                                                            tCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCca                                                                                                                                        <  1:326489/152‑1 (MQ=255)
                                                                                                                                                                                                       gggaacaaaCGCTATTCCCCACACCCCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCCCCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGcc  <  2:601707/269‑1 (MQ=255)
                                                                                                                                                                                                                   tctTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCca                                                                                                                                        >  2:464025/1‑129 (MQ=255)
                                                                                                                                                                                                                    ctTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACAGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCAGTAGGCCAgcg           >  1:401496/1‑253 (MQ=255)
                                                                                                                                                                                                                       ccccAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGa                                                                                                                                    <  1:464032/129‑1 (MQ=255)
                                                                                                                                                                                                                                 cccGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGc                                                                                         >  2:366603/1‑162 (MQ=255)
                                                                                                                                                                                                                                 cccGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGCTcgcg                                                                                                        >  1:574310/1‑147 (MQ=255)
                                                                                                                                                                                                                                     gTAGAAAGGTAGCGATCGCCACGTTCGCCGATGCCCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCCCCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATg                                                                                                    <  2:574301/147‑1 (MQ=255)
                                                                                                                                                                                                                                     gTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATcg                                                                                     <  1:366609/162‑1 (MQ=255)
                                                                                                                                                                                                                                                gCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTc                                                                                            >  1:238823/1‑144 (MQ=255)
                                                                                                                                                                                                                                                gCGATCGCCACGATCCCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAAc                                                                                                                                                              >  2:652110/1‑78 (MQ=255)
                                                                                                                                                                                                                                                    tCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCa                                                                                        <  2:238818/144‑1 (MQ=255)
                                                                                                                                                                                                                                                    tCGCCACGATCCCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACcgcg                                                                                                                                                          <  1:652119/78‑1 (MQ=255)
                                                                                                                                                                                                                                                        cACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCa                                                                                                                                                                                 >  2:418943/1‑51 (MQ=255)
                                                                                                                                                                                                                                                            aTCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCg                                                                                                                                                                             <  1:418950/51‑1 (MQ=255)
                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                          
CATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCCCAAACACCCCGGTAGAAAGGTAGCGATCGCCACGATCGCAGATGATCGCCACCACCACCGCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCC  >  CP009273/2531849‑2532321

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