Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 216,233 T→C 12.4% I201V (ATT→GTT)  argD ← ornithine‑‑oxo‑acid transaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730216,2330TC12.4% 90.4 / 2.7 38I201V (ATT→GTT) argDornithine‑‑oxo‑acid transaminase
Reads supporting (aligned to +/- strand):  ref base T (17/16);  new base C (2/3);  total (19/19)
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

CAAACCATTTCCCAGTTCTACCAAGACCAACTTGAATTTCATCTGCAATCAATAATATTTGATGTTTATCACATAGTTGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTACTGTTGTAGTGCCTTGAAGTAGGGGGTGAAATCCTGCTTTATATGCGTCGTGGTTAGATAGTGATAATG  >  CP000730/216098‑216373
                                                                                                                                       |                                                                                                                                            
cAAACCATTTCCCAGTTCTACCAAGACCAACTTGAATTTCATCTGCAATCAATAATATTTGATGTTTATCACATAGTTGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTCAATTGGtt                                                                                                                                         >  2:149215/1‑141 (MQ=255)
  aaCCATTTCCCACTTCTACCAAGACTAACTTGATTTTCATCTGCAATCAATAATATTTGATGTTTATCACATAGTTGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTcc                                                                                                                                       <  2:167745/141‑1 (MQ=255)
       tttCCCAGTTCTACCAAGACCAACTTGAATTTCATCTGCAATCAATAATATTTGATGTTTATCACATAGTTGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAaat                                                                                                                                  >  2:304983/1‑141 (MQ=255)
         tCCCAGTTCTACCAAGACCAACTTTAAATTCATCTGCAATCAATAATATTGGATGTTTATCACATAGTTGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAaataa                                                                                                                                >  2:274574/1‑141 (MQ=255)
               ttCTACCAAGACCAACTTGAATTTCATCTGCAATCAATAGCATCTGATGTTTATCACATAGTTGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATtgctg                                                                                                                          <  1:349116/141‑1 (MQ=255)
                 cTACCAAGACCAACTTGAATTTCATCTGCAATCAATAATATTTGATGTTTATCACATAGTTGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATtgctgcg                                                                                                                        >  1:116347/1‑140 (MQ=255)
                 cTACCAAGACCAACTTGAATTTCATCTGCAATCAATAATATTTGATGTTTATCACATAGTTGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTgctgct                                                                                                                        <  1:88759/141‑1 (MQ=255)
                      aaGACCAACTTGAATTTCATCTGCAATCAATAATATTTGATGTTTATCACATAGTTGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAACTGGTTCCAAAATAATTGCTGCTGTAtt                                                                                                                   >  1:116434/1‑141 (MQ=255)
                      aaGACCAACTTGAATTTCATCTGCAATCAATAATATTTGATGTTTATCACATAGTTGACGCACAGCTTGAATATATCCTTTCGGGGGTATATTAACGCCACCTTCACCTTGACCGGGTGCCACAATACTTGCTGCTGGAtt                                                                                                                   >  1:106130/1‑141 (MQ=255)
                         aCCAACTTGAATTTCATCTGCAATCAATAATATTTGATGTTTATCACATAGTTGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCgg                                                                                                                <  1:304983/141‑1 (MQ=255)
                         aCCAACTTGAATTTCATCTGCAATCAATAATATTTGATGTTTATCACATAGTTGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCgg                                                                                                                <  2:269342/141‑1 (MQ=255)
                             aCTTGAATTTCATCTGCAATCAATAATATTTGATGTTTATCACATAGTTGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATATTTGCTGCTGTATTCGGTGaa                                                                                                            >  2:100441/1‑141 (MQ=255)
                                tGAATTTCATCTGCAATCAATAATATTTGATGTTTATCACATAGTTGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAAtg                                                                                                         >  1:43738/1‑140 (MQ=255)
                                    tttCATCTGCAATCAGTAATATTTGATGTTTATCACATAGTTGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCtt                                                                                                     >  1:276693/1‑141 (MQ=255)
                                        aTCTGCAATCAATAATATTTGATGTTTATCACATAGTTGACGCACAGCTTGCATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAGAATAATTGCTGCTGTATTCGGTGAAATAGCTTgtgt                                                                                                 >  1:311390/1‑141 (MQ=255)
                                        aTCTGCAATCAATAATATTTGATGTTTATCACATAGTTGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTCCTGTATTCGGTGAAATAGCCTgtgt                                                                                                 >  2:90363/1‑141 (MQ=255)
                                             caagcaacAATCTTTGATGTTTATCACATAGTTGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAAtt                                                                                            <  1:5492/141‑1 (MQ=255)
                                               atcaatAATATTTGATGTTTATCACATAGTTGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGGATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGTAATAGCTTGTGTTAATTGt                                                                                          <  1:242441/141‑1 (MQ=255)
                                                          ttGATGTTTATCACATAGTTGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAACTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTcc                                                                               <  2:116434/141‑1 (MQ=255)
                                                          tgtgtgtttcgcccaTAGTTCACGCACAGCGTGAATATATCCTTCCGCTGGTATATTAACCCCACCTTCCCCTTGAACTGGTTCCAAAATAATTGCTGCTGTATTCCGTGAAATAGCTTGTTTTAATTGTTCAATGTCTcc                                                                               <  2:106130/128‑1 (MQ=255)
                                                            gATGTTTATCACATAGTTGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCaa                                                                             >  2:276077/1‑141 (MQ=255)
                                                                    tCACATAGTTGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAATATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTGATTGTTCAATGTCTCCAAAATCTACt                                                                     <  2:223274/141‑1 (MQ=255)
                                                                        aTAGTTGACGCCCAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTACtgttg                                                                 <  2:284953/141‑1 (MQ=255)
                                                                            ttGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTTGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTACTGTTGTAGt                                                             <  1:4832/141‑1 (MQ=255)
                                                                                cGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTACTGTTGTAGTGCCt                                                         >  1:147782/1‑141 (MQ=255)
                                                                                 gCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTACTGTTGTAGTGCCtt                                                        >  1:209009/1‑141 (MQ=255)
                                                                                    caGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATTATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTAACTGTCTCCAAAATCTCCTGTTGTAATGCCTTGaa                                                     >  2:238911/1‑141 (MQ=255)
                                                                                    caGCTTGAATATATCCTTTCGGGTGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTACTGTTGTAGTGCCTTGaa                                                     <  2:10386/141‑1 (MQ=255)
                                                                                     aGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTACTGTTGTAGTGCCTTGAAg                                                    <  1:212261/141‑1 (MQ=255)
                                                                                                  ccTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTACTGTTGTAGTGCCTTGAAGTAGGGGGTGAAAt                                       <  2:325960/141‑1 (MQ=255)
                                                                                                  ccTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTACTGTTGTAGTGCCTTGAAGTAGGGGGTGAAAt                                       >  2:330486/1‑141 (MQ=255)
                                                                                                   cgtgcGGTGGTATATTTCCGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGATGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTTTCCAAAATCTACTGTTTTAGTGCCTTGAAGTAGTGGTTGAAATc                                      <  2:370797/137‑1 (MQ=255)
                                                                                                     ttCGGTGGTATATTAACGCCACCTTCACCTAGCACTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTACTGTTGTAGTGCCTTGAAGTAGGGGGTGAAATCCt                                    <  1:192742/141‑1 (MQ=255)
                                                                                                         gtggtATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTACTGTTGTAGTGCCTTGAAGTAGGGGGTGAAATCCTGCtt                                >  2:104848/1‑141 (MQ=255)
                                                                                                         gtggtATATTAACGCCACCATCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTACTGTTGTAGTGCCTTGAAGTAGGGGGAGAAATCCTGCtt                                >  2:208878/1‑141 (MQ=255)
                                                                                                                  tAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTACTGTTGTAGTGCCTTGAAGTAGGGGGTGAAATCCTGCTTTATATGcgt                       <  1:111913/141‑1 (MQ=255)
                                                                                                                            ttcaccttGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTACTGTTGTAGTGCCTTGAAGTAGGGGGTGAAATCCTGCTTTATATGCGTCGTGGTTAGa             >  1:135445/1‑141 (MQ=255)
                                                                                                                                      aTTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTACTGTTGTAGTGCCTTGAAGTAGGGGGTGAAATCCTGCTTTATATGCGTCGTGGTTAGATAGTGATAAt   <  1:19730/141‑1 (MQ=255)
                                                                                                                                       ttGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTACTGTTGTAGTGCCTTGAAGTAGGGGGTGAAATCCTGCTTTCTATGCGTCGTGGTTAGATAGTGATAATg  >  1:263091/1‑141 (MQ=255)
                                                                                                                                       |                                                                                                                                            
CAAACCATTTCCCAGTTCTACCAAGACCAACTTGAATTTCATCTGCAATCAATAATATTTGATGTTTATCACATAGTTGACGCACAGCTTGAATATATCCTTTCGGTGGTATATTAACGCCACCTTCACCTTGAATTGGTTCCAAAATAATTGCTGCTGTATTCGGTGAAATAGCTTGTGTTAATTGTTCAATGTCTCCAAAATCTACTGTTGTAGTGCCTTGAAGTAGGGGGTGAAATCCTGCTTTATATGCGTCGTGGTTAGATAGTGATAATG  >  CP000730/216098‑216373

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