Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,533,049 C→T A214V (GCG→GTG)  cysK → cysteine synthase A, O‑acetylserine sulfhydrolase A subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,533,0490CT91.4% 115.4 / ‑3.1 35A214V (GCG→GTG) cysKcysteine synthase A, O‑acetylserine sulfhydrolase A subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (7/25);  minor base G (3/0);  total (10/25)
Fisher's exact test for biased strand distribution p-value = 1.83e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTACGCTGACTGGCGTCAGCCGCTACATTAAAGGCACCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCAGGCGCTGGCAGGTGAAGAGATTAAACCTGGCCCGCATAAAATTCAGGGTATTGGCGCTGGTTTTATCCCGGCTAACCTCGATCTCAAGCTGGTCGATAAAGTCATTGGCATCACCAATGAAGAAGCGATTTCTACCGCGCGTCGTCT  >  NC_000913/2532907‑2533193
                                                                                                                                               |                                                                                                                                                
ggTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTACGCTGACTGGCGTCAGCGGCTACATTAAAGGCACCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCAGGTGCTGGCa                                                                                                                                           <  2:621284/150‑1 (MQ=255)
ggTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTACGCTGACTGGCGTCAGCCGCTACATTAAAGGCACCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCAGGTGCTGGCAg                                                                                                                                          <  2:279924/151‑1 (MQ=255)
   cAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTACGCTGACTGGCGTCAGCCGCTACATTAAAGGCACCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCAGGTGCTGGCAGGt                                                                                                                                        <  1:350555/150‑1 (MQ=255)
     ggTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTAAGCTGACTTGCGGCAGACGCCAAATTTAAAGGACACAAAGCAAGACCCAAATTATTTCCGTCTCCGCTTGGGCCAACCATTTTCCCCGTATTGGCCACGGGCTGGCAgggg                                                                                                                                       >  2:92931/1‑147 (MQ=255)
            gTATTTATTGCTGGCGTTGGGACTGGCGGTACGCCGGCTGGCGGCAGCCGCTACATTAAAAGCACCAAAAGCAAAAACGATCTTTTTTCTTTCGCCGTTTAGCCAACCCATTCTCCAGTTATTTCCCCGGGGCTGGCAGGGGAAGAGAtt                                                                                                                               >  2:143515/1‑150 (MQ=255)
                 tattGCTGGCGTTGGGACTGGCGGTACGCTGACTGGCGTCAGCCGCTACATTAAAGGCACCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCAGGTGCTGGCAGGTGAAGAGATTAAAcc                                                                                                                          <  2:53/150‑1 (MQ=255)
                  attGCTGGCGTTGGGACTGGCGGTACGCTGACTGGCGTCAGCCGCTACATTAAAGGCACCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCAGGTGCTGGCAGGTGAAGAGATTAAACCt                                                                                                                         <  2:573977/150‑1 (MQ=255)
                     gCTGGCGTTGGGACTGGCGGTACGCTGACTGGCGTCAGCCGCTACATTAAAGGCACCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCAGGTGCTGGCAGGTGAAGAGATTAAACCTGGcc                                                                                                                     <  2:23002/151‑1 (MQ=255)
                      cTGGCGTTGGGACTGGCGGTACGCTGACTCGGGGCTGCCGCTTCATTTAAAGGACCAAAAGGAAGAAAGGGCTTATTTTTTTCGCCGGTGGGCGCAACGATTCTTCAGGTATTGCCCACGGGCTGGGAGGGGGAAAGAGTTAAACTGGccc                                                                                                                    >  1:383222/1‑151 (MQ=255)
                              gggACTGGCGGTACGCTGACTGGCGTCAGCCGCTACATTAAAGGCACCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCAGGTGCTGGCAGGTGAAGAGATTAAACCTGGCCCGCATaaaa                                                                                                            <  2:404283/151‑1 (MQ=255)
                              gggACTGGCGGTACGCTGACTGGCGTCAGCCGCTACATTAAAGGCACCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCAGGTGCTGGCAGGTGAAGAGATTAAACCTGGCCCGCATaaa                                                                                                             <  1:534104/150‑1 (MQ=255)
                                    ggCGGTACGCTGACTGGCGTCAGCCGCTACATTAAAGGCACCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCAGGTGCTGGCAGGTGAAGAGATTAAACCTGGCCCGCATAAAATTCa                                                                                                        <  1:638138/149‑1 (MQ=255)
                                    ggCGGTACGCTGACTGGCGTCAGCCGCTACATTAAAGGCACCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCAGGTGCTGGCAGGTGAAGAGATTAAACCTGGCCCGCATAAAATTCAgg                                                                                                      >  1:521968/1‑151 (MQ=255)
                                          aCGCTGACTGGCGTCAGCCGCTACATTAAAGGCACCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCAGGTGCTGGCAGGTGAAGAGATTAAACCTGGCCCGCATAAAATTCAGGGTAtt                                                                                                 <  2:544501/150‑1 (MQ=255)
                                                   ggCGTCAGCCGCTACATTAAAGGCACCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCAGGTGCTGGCAGGTGAAGAGATTAAACCTGGCCCGCATAAAATTCAGGGTATTGGCGCTGGt                                                                                        <  1:240115/150‑1 (MQ=255)
                                                                    tAAAGGCACCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCAGGTGCTGGCAGGTGAAGAGATTAAACCTGGCCCGCATAAAATTCAGGGTATTGGCGCTGGTTTTATCCCGGCTAACCt                                                                       <  2:60828/150‑1 (MQ=255)
                                                                     aaaGGCACCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCACAGGTGCTGGCAGGTGAAGAGATTAAACCTGGCCCGCATAAAATTCAGGGTATTGGCGCTGGTTTTATCCCGGCTAACCTCg                                                                     <  1:526797/151‑1 (MQ=255)
                                                                           aCCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCAGGTGCTGGCAGGTGAAGAGATTAAACCTGGCCCGCATAAAATTCAGGGTATTGGCGCTGGTTTTATCCCGGCTAACCTCGAtctc                                                                >  1:228586/1‑150 (MQ=255)
                                                                            ccAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCAGGTGCTGGCAGGTGAAGAGATTAAACCTGGCCCGCATAAAATTCAGGGTATTGGCGCTGGTTTTATCCCGGCTAACCTCGAtctc                                                                <  1:237479/149‑1 (MQ=255)
                                                                             cAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCAGGTGCTGGCAGGTGAAGAGATTAAACCTGGCCCGCATAAAATTCAGGGTATTGGCGCTGGTTTTATCCCGGCTAACCTCGATCTCAAg                                                             <  2:333778/151‑1 (MQ=255)
                                                                                aGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCAGGTGCTGGCAGGTGAAGAGATTAAACCTGGCCCGCATAAAATTCAGGGTATTGGCGCTGGTTTTATCCCGGCTAACCTCGATCTCAAGCTg                                                          >  1:499274/1‑151 (MQ=255)
                                                                                     aGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCAGGTGCTGGCAGGTGAAGAGATTAAACCTGGCCCGCATAAAATTCAGGGTATTGGCGCTGGTTTTATCCCGGCTAACCTCGATCTCAAGCTGGTCGa                                                     <  1:358753/151‑1 (MQ=255)
                                                                                          gATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCCCCCAGGTGCTGGCAGGTGAAGAGATTAAACCTGGCCCGCATAAAATTCAGGGTATTGGCGCTGGTTTTATCCCGGCTAACCTCGATCTCAAGCTGGTCGATaa                                                  >  1:576136/1‑149 (MQ=255)
                                                                                               tATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCAGGTGCTGGCAGGTGAAGAGATTAAACCTGGCCCGCATAAAATTCAGGGTATTGGCGCTGGTTTTATCCCGGCTAACCTCGATCTCAAGCTGGTCGATaaa                                                 <  1:143511/145‑1 (MQ=255)
                                                                                                aTCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCAGGTGCTGGCAGGTGAAGAGATTAAACCTGGCCCGCATAAAATTCAGGGTATTGGCGCTGGTTTTATCCCGGCTAACCTCGATCTCAAGCTGGTCGATAAAGTCAtt                                           <  2:41286/150‑1 (MQ=255)
                                                                                                       tCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCAGGTGCTGGCAGGTGAAGAGATTAAACCTGGCCCGCATAAAATTCAGGGTATTGGCGCTGGTTTTATCCCGGCTAACCTCGATCTCAAGCTGGTCGATAAAGTCATTGGCATCAc                                   >  2:526327/1‑151 (MQ=255)
                                                                                                        cGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCAGGTGCTGGCAGGTGAAGAGATTAAACCTGGCCCGCATAAAATTCAGGGTATTGGCGCTGGTTTTATCCCGGCTAACCTCGATCTCAAGCTGGTCGATAAAGTCATTGGCATCAc                                   <  2:304125/150‑1 (MQ=255)
                                                                                                            gTTGAGCCAACCGATTCTCCAGTTATCGCCCAGGTGCTGGCAGGTGAAGAGATTAAACCTGGCCCGCATAAAATTCAGGGTATTGGCGCTGGTTTTATCCCGGCTAACCTCGATCTCAAGCTGGTCGATAAAGTCATTGGCATCACCAATg                              <  1:130619/151‑1 (MQ=255)
                                                                                                              ttGAGCGTAACGATTCTCCAGTTATCGCCCATGTGGTTGCAGGTGAAGAGATTAAACCTGGACCGCATAAAATTCAGGGTAGTGGCGCTGCTTTTATCCCGGCTAACCTCGATCTCAAGCTGGTCGATAAAGTCATTGGCATCTCCaa                                <  1:649775/148‑1 (MQ=255)
                                                                                                               tGAGCCAACCGATTCTCCAGTTATCGCCCAGGTGCTGGCAGGTGAAGAGATTAAACCTGGCCCGCATAAAATTCAGGGTATTGGCGCTGGTTTTATCCCGGCTAACCTCGATCTCAAGCTGGTCGATAAAGTCATTGGCATCACCAATga                             >  1:21269/1‑150 (MQ=255)
                                                                                                                    cAACCGATTCTCCAGTTATCGCCCAGGTGCTGGCAGGTGAAGAGATTAAACCTGGCCCGCATAAAATTCAGGGTATTGGCGCTGGTTTTATCCCGGCTAACCTCGATCTCAAGCTGGTCGATAAAGTCATTGGCATCACCAATGAAGAAGc                       <  2:271315/151‑1 (MQ=255)
                                                                                                                       ccGATTCTCCAGTTATCGCCCAGGTGCTGGCAGGTGAAGAGATTAAACCTGGCCCGCATAAAATTCAGGGTATTGGCGCTGGTTTTATCCCGGCTAACCTCGATCTCAAGCTGGTCGATAAAGTCATTGGCATCACCAATGAAGAAGCGAt                    <  2:29213/151‑1 (MQ=255)
                                                                                                                        cGATTCTCCAGTTATCGCCCAGGTGCTGGCAGGTGAAGAGATTAAACCTGGCCCGCATAAAATTCAGGGTATTGGCGCTGGTTTTATCCCGGCTAACCTCGATCTCAAGCTGGTCGATAAAGTCATTGGCATCACCAATGAAGAAGCGAtt                   >  2:268590/1‑151 (MQ=255)
                                                                                                                                 aGTTATCGCCCAGGTGCTGGCAGGTGAAGAGATTAAACCTGGCCCGCATAAAATTCAGGGTATTGGCGCTGGTTTTATCCCGGCTAACCTCGATCTCAAGCTGGTCGATAAAGTCATTGGCATCACCAATGAAGAAGCGATTTCTACcgcg          <  2:58900/151‑1 (MQ=255)
                                                                                                                                         cccAGGTGCTGGCAGGTGAAGAGATTAATCCTGGCCCGCATAAAATTCAGGGTGTTGGCGCTGGTTTTATCCCGGCTAACCTCGATCTCAAGCTGGTCGATAAAGTCATTGGCATCACCAATGAAGAAGCTATTTCTACCGCGCGTCGTCt  <  1:143031/151‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                
GGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTACGCTGACTGGCGTCAGCCGCTACATTAAAGGCACCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCAGGCGCTGGCAGGTGAAGAGATTAAACCTGGCCCGCATAAAATTCAGGGTATTGGCGCTGGTTTTATCCCGGCTAACCTCGATCTCAAGCTGGTCGATAAAGTCATTGGCATCACCAATGAAGAAGCGATTTCTACCGCGCGTCGTCT  >  NC_000913/2532907‑2533193

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