Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,828,478 T→C 11.8% I192M (ATA→ATG)  csd2 ← cysteine desulfurase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,828,4780TC11.8% 98.8 / 4.8 41I192M (ATA→ATG) csd2cysteine desulfurase
Reads supporting (aligned to +/- strand):  ref base T (19/17);  new base C (3/2);  total (22/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 = 9.78e-01

ATTAACTGTTCCACTTCTAACGCCATATTCTTGACCACCACCATGGACAGTTGGTTCAACATTTTGAATGTGATTTACAAGTAAGACGCCTTGTCCTTTTAAACCATTAAACTTGTGTCCACTTAAACTAATACTATCTATGTTATTGAGATCCATTGAAATTTTGCCGAATGCTTGAACCGCATCTACATGAAAATGTGCCTTAGGATAATTTTTTATAACTTTAGCCATTTGTGGAATAGGCTGTATTTGTCCAGTTACATTATTTACATACATACA  >  CP000730/1828340‑1828618
                                                                                                                                          |                                                                                                                                            
aTTAACTGTTCCACTTCTAACGCCATATTCTTGACCACCACCATGGACAGTTGGTTCAACATTTTGAATGTGATTTACAAGTAAGACGCCTTGTCCTTTTAAACCATTAAACTTGTGTCCACTTAAACTAATActatctat                                                                                                                                            <  2:129380/141‑1 (MQ=255)
      tGTTCCACTTCTAACGCCATATTCTTGACCACCACCATGGACAGGTGGTTCAACATTTTGAATGTGATTTACAAGTAAGACGCCTTGTCCTTTTAATCCATTAAACTTGTGTCCACTTAAACTAATACTATCTAGGTTAtt                                                                                                                                      <  2:223264/141‑1 (MQ=255)
           cACTTCTAACGCCATATTCTTGACCACCACAATGGACAGTTGGTTCAACATTTTGAATGTGATTTACAAGTAAGACGCCTTGTCCTTTTAAACCATTAAACTTGTGTCCACTTAAACTAATACTATCTATGTTATTGAGAt                                                                                                                                 >  2:425719/1‑141 (MQ=255)
               tCTAACGCCATATTCTTGACCACCACCATGGACAGTTGGTTCAACATTTTGAATGTGATTTACAAGTAAGACGCCTTGTCCTTTTAAACCATTAAACTTGTGTCCACTTAAACTAATACTATCTATGTTATTGAGATCCAt                                                                                                                             <  2:51838/141‑1 (MQ=255)
                 tAACGCCATATTCTTGACCACCACCATGGACAGTTGGTTCAACATTTTGAATGTGATTTACAAGTAAGACGCCTTGTCCTTTTAAACCATTAAACTTGTGTCCACTTAAACTAATACTATCTATGTTATTGAGATCCATTg                                                                                                                           >  2:258440/1‑141 (MQ=255)
                  aaCGCCATATTCTTGACCACCACCATGGACAGTTGGTTCAACATTTTGAATGTGATTTACAAGTAAGACGCCTTGTCCTTTTAAACCATTAAACTTGTGTCCACTTAAACTAATACTATCTATGTTATTGAGATCCATTgc                                                                                                                          <  2:241513/141‑2 (MQ=255)
                  aaCGCCATATTCTTGACCACCACCATGGACAGTTGGTTCAACATTTTGAATGTGATTTACAAGTAAGACGCCTTGTCCTTTTAAACCATTAAACTTGTGTCCACTTAAACTAATACTATCTATGTTATTGAGATCCATTGa                                                                                                                          >  1:372345/1‑141 (MQ=255)
                   aCGCCATATTCTTGACCACCACCATGGACAGTTGGTTCAACATTTTGAATGTGATTTACAAGTAAGACGCCTTGTCCTTTTAAACCATTAAACTTGTGTCCACTTAAACTAATACTATCTATGTTATTGAGATCCATTGaa                                                                                                                         >  2:405113/1‑141 (MQ=255)
                    cGCCATATTCTTGACCACCACCATGGACAGTTGGTTCAACATTTTGAATGTGATTTACAAGTAAGACGCCTTGTCCTTTTAAACCATTAAACTTGTGTCCACTTAAACTAATACTATCTATGTTATTGAGATCCATTGaaa                                                                                                                        >  1:172479/1‑141 (MQ=255)
                         tatTCTTGACCACCACCATGGACAGTTGGTTCAACATTTTGAATGTGATTTACAAGTAAGACGCCTTGTCCTTTTAAACCATTAAACTTGTGTCCACTTAAACTAATACTATCTATGTTATTGAGATCCATTGAAATTTTg                                                                                                                   <  1:367660/141‑1 (MQ=255)
                          atTCTTGACCACCACCATGGACAGTTGGTTCAACATTTTGAATGTGATTTACAAGTAAGACGCCTTGTCCTTTTAAACCATTAAACTTGTGTCCACTTAAACTAATACTATCTATGTTATTGAGATCCATTGAAATTTTGc                                                                                                                  <  1:232032/141‑1 (MQ=255)
                              ttGACCACCACCATAATCACAGGGTTTGACACTTAGAATGTGAGTTACTCATAAGACACCATGTCATTTTAAATCATTCCACTCGTATCCTCTTAAACTAAAACTATCCATGTTATTGAGATCCATTGAATGTTTGCCGag                                                                                                              >  2:163072/1‑140 (MQ=255)
                                       accaTGGACAGTTGGTTCAACATTTTGAATGTGATTTACAAGTAAGACGCCTTGTCCTTTTAAACCATTAAACTTGTGTCCACTTAAACTAATACTATCTATGTTATTGAGATCCATTGAAATTTTGCCGAATGCTTGAAc                                                                                                     <  2:342992/141‑1 (MQ=255)
                                               cAGTTGGTTCAACATTTTGAATGTGATTTACAAGTAAGACGCCTTGTCCTTTTAAACCATTAAACTTGTGTCCACTTAAACTAATACTATCTATGTTATTGAGATCCATTGAAATTTTGCCGACTGCTTGAACCGCATCTa                                                                                             >  2:298061/1‑141 (MQ=255)
                                                         aaCATTTTGAATGTGATTTACAAGTAAGACGCCTTGTCCTTTTAAACCATTAAACTTGTGTCCACTTAAACTAATACTATCTATGTTATTGAGATCCATTGAAATTTTGCCTAATGCTTGAACCGCATCTACATGAAAAtg                                                                                   >  2:217649/1‑141 (MQ=255)
                                                                 gAATGTGATTTACAAGTAAGACGCCTTGTCCTTTTAAACCATTAAACTTGTGTCCACTTAAACTAATACTATCTATGTTATTGAGATCCATTGAAATTTTGCCGAATGCTTGTACCGCATCTACATGAAAATGTGCCTTAg                                                                           <  2:172479/141‑1 (MQ=255)
                                                                 gAATGTGATTTACAAGTAAGACGCCTTGTCCTTTTAAACCATTAAACTTGTGTCCACTTAAACTAATACTATCTATGTTATTGAGATCCATTGAAATTTTGCCGAATGCTTGAACCGCATCTACATGAAAATGTGCCTTAg                                                                           >  1:6163/1‑141 (MQ=255)
                                                                   aTGTGATTTACAAGTAAGACGCCTTGTCCTTTTAAACCATTAAACTTGTGTCCACTTAAACTAATACTATCTATGTTATTGAGATCCATTGAAATTTTGCCGAATGCTTGAACCGCATCTACATGAAAATGTGCCTTAGGa                                                                         <  2:169590/141‑1 (MQ=255)
                                                                    tgtgATTTACAAGTAAGACGCCTTGTCCTTTTAAACCATTAAACTTGTGTCCACTTAAACTAATACTATCTATGTTATTGAGATCCATTGAAATTTTGCCGAATGCTTGAACCGCATCTACATGAAAATGTGCCTTAGGAt                                                                        >  1:434540/1‑141 (MQ=255)
                                                                    tgtgATTTACAAGTAAGACGCCTTGTCCTTTTAAACCATTAAACTTGTGTCCACTTAAACTAATACTATCTATGTTATTGAGATCCATTGAAATTTTGCCGAATGCTTGAACCGCATCTACATGAAAATGTGCCTTAGGAt                                                                        >  1:348382/1‑141 (MQ=255)
                                                                    tgtgATTTACAAGTAAGACGCCTTGTCCTTTTAAACCATTAAACTTGTGTCCACTTAAACTAATACTATCTATGTTATTGAGATCCATTGAAATTTTGCCGAATGCTTGAACCGCATCTACATGAAAATGTGCCTTAGGAt                                                                        <  2:466059/141‑1 (MQ=255)
                                                                             cAAGTAAGACGCCTTGTCCTTTTAAACCATTAAACTTGTGTCCACTTAAACTAATACTATCTATGTTATTGAGATCCATTGAAATTTTGCCGAATGCTTGAACCGCATCTACATGAAAATGTGCCTTAGGATAATTTTTta                                                               >  1:280492/1‑141 (MQ=255)
                                                                                       gCCTTGTCCTTTTAAACCATTAAACTTGTGTCCACTTAAACTAATACTATCTATGTTATTGAGATCCATTGAAATTTTGCCGAATGCTTGAACCGCATCTACATGAAAATGTGCCTTAGGATAATTTTTTATAACTTTAGc                                                     >  2:117979/1‑141 (MQ=255)
                                                                                           tggccTTTTAGACCATTAAACTTGTGTCCACTTAAACTAATACTAACTTGGGTATTGAGATCCATTGAGATTTTGCCGAATGCTTGAACCGCATCTACATGAAAATGTGCGTTAGGATATGTGTTTATATCGGTAGCCAtt                                                 >  2:47220/4‑141 (MQ=255)
                                                                                             tCCTTTTAAACCATTAAACTTGTGTCCACTTAAACTAATACTATCTATGTTATTGAGATCCATTGAAATTTTGCCGAATGCTTGAACCGCATCTACATGAAAATGTGCCTTAGGATAATTTTTTATAACTTTAGCCATTtg                                               <  1:308117/141‑1 (MQ=255)
                                                                                                  ttAAACCATTAAACTTGTGTCCACTTAAACTAATACTATCTATGTTATTGAGATCCATTGAAATTTTGCCGAATGCTTGAACCGCATCTACATGAAAATGTGCCTTAGGATAATTTTTTATAACTTTAGCCATTTGTGGaa                                          >  1:79897/1‑141 (MQ=255)
                                                                                                   tAAACCATTAAACTTGTGTCCACTTAAACTAATACTATCTATGTTATTGAGATCCATTGAAATTTTGCCGAATGCTTGCACCGCATCTACATGAAAATGTGCCTTAGGATAATTTTTTATATCTTTTGCTATTTGTGGAAt                                         >  2:96379/1‑141 (MQ=255)
                                                                                                   tAAACCATTAAACTTGTGTCCACTTAAACTAATACTATCTATGTTATTGAGATCCATTGAAATTTTGCCGAATGCTTGAACCGCATCTACATGAAAATGTGCCTTAGGATAATTTTTTATAACTTTAGCCATTTGTGGAAt                                         <  2:62009/141‑1 (MQ=255)
                                                                                                     aaCCATTAAACTTGTGTCCACTTAAACTAATACTATCTATGTTATTGAGATCCATTGAAATTTTGCCGAATGCTTGAACCGCATCTACATGAAAATGTGCCTTAGGATAATTTTTTATAACTTTAGCCATTTGTGGAATAg                                       <  1:312318/141‑1 (MQ=255)
                                                                                                         aTTAAACTTGTGTCCACTTAAACTAATACTATCTATGTTATTGAGATCCATTGAAATTTTGCCGAATGCTTGAACCGCATCTACATGAAAATGTGCCTTAGGATAATTTTTTATAACTTTAGCCATTTGTGGAATAGGCTg                                   >  1:504056/1‑141 (MQ=255)
                                                                                                         aTTAAACTTGTGTCCACTTAAACTAATACTATCCATGTTATTGAGATCCATTGAAATTTTGCCGAATGCTTGAACCGCATCTACATGaa                                                                                       <  1:195108/89‑1 (MQ=255)
                                                                                                         aTTAAACTTGTGTCCACTTAAACTAATACTATCCATGTTATTGAGATCCATTGAAATTTTGCCGAATGCTTGAACCGCATCTACATGaa                                                                                       >  2:195108/1‑89 (MQ=255)
                                                                                                         aTTAAACTTGTGTCCACTTAAACTAATACTATCCATGTTATTGAGATCCATTGAAATTTTGCCGAATGCTTGAACCGCATCTACATGaa                                                                                       >  2:195067/1‑89 (MQ=255)
                                                                                                         aTTAAACTTGTGTCCACTTAAACTAATACTATCCATGTTATTGAGATCCATTGAAATTTTGCCGAATGCTTGAACCGCATCTACATGaa                                                                                       <  1:195067/89‑1 (MQ=255)
                                                                                                          ttAAACTTGTGTCCACTTATACTAATACTATCTATGTTATTGAGATCCATTGAAATTTTGCCGAATGCTTGAACCGCATCTATTTGAAAATGTGCCTTAGGATAATTTTTTATAACTTTAGCCATTTGTGGAATAGGCTGt                                  <  1:285066/141‑1 (MQ=255)
                                                                                                          ttAAACTTGTGTCCACTTAAACTAATACTATCTATGTTATTGAGATCCATTGAAATTTTGCCGAATGCTTGAACCGCATCTACATGAAAATGTGCCTTAGGATAATTTTTTATAACTTTAGCCATTTGTGGAATAGGCTGt                                  <  2:500793/141‑1 (MQ=255)
                                                                                                                      ccACTTAAACTAATACTATCTATGTTATTGAGATCCATTGAAATTTTGCCGAATGCTTGAACCGCATCTACATGAAAATGTGCCTTAGGATAATTTTTTATAACTTTAGCCATTTGTGGAATAGGCTGTATTTGTCCAGtt                      >  1:136694/1‑141 (MQ=255)
                                                                                                                      ccACTTAAACTAATACTATCTATGTTATTGAGATCCATTGAAATTTTGCCGAATGCTTGAACCGCATCTACATGAAAATGTGCCTTAGGATAATTTTTTATAACTTTAGCCATTTGTGGAATAGGCTGTATTTGTCCAGtt                      >  1:137329/1‑141 (MQ=255)
                                                                                                                       cACTTAAACTAATACTATCTATGTTATTGAGATCCATTGAAATTTTGCCGAATGCTTGAACCGCATCTACATGAAAATGTGCCTTAGGATAATTTTTTATAACTTTAGCCATTTGTGGAATAGGCTGTATTTGTCCAGTTa                     >  1:160066/1‑141 (MQ=255)
                                                                                                                         cTTAAACTAATACTATTTATGTTATTGAGATCCATTGAAAGTTTGCCGAATGCTTGAACCGCATCTACATGAAAATGTGCCTTAGGATAATTTTTTATAACTTTAGCCATTTGTGGAATAGGCTGTATTTGTCCAGTTACa                   <  2:39734/141‑1 (MQ=255)
                                                                                                                           tAAACTAATACTATCTATGTTATTGAGATCCATTGAAATTTTGCCGAATGCTTGAACCGCATCTACATGAAAATGTGCCTTAGGATAATTTTTTATAACTTTAGCCATTTGTGGAATAGGCTGTATTTGTCCAGTTACatt                 <  1:298061/141‑1 (MQ=255)
                                                                                                                                  aTACTATCTATGTTATTGAGATCCATTGAAATTTTGCCGAATGCTTGAACCGCATCTACATGAAAATGTGCCTTAGGATAATTTTTTATAACTTTAGCCATTTGTGGAATAGGCTGTATTTGTCCAGTTACATTATTtaca          <  2:451996/141‑1 (MQ=255)
                                                                                                                                     ctatctatGTTATTGAGATCCATTGAAATTTTGCCGAATGCTTGAACCGCATCTACATGAAAATGTGCCTTAGGATAATTTTTTATAACTTTAGCCATTTGTGGAATAGGCTGTATTTGTCCAGTTACATTATTtacatac       <  1:117979/141‑1 (MQ=255)
                                                                                                                                         ctatGTTATTGAGATCCATTGAAATTTTGCCGAATGCTTGAACCGCATCTACATGAAAATGTGCCTTAGGATAATTTTTTATAACTTTAGCCATTTGTGGAATAGGCTGTATTTGTCCAGTTACATTATTtacatacatac   >  2:505595/1‑141 (MQ=255)
                                                                                                                                         ctatGTTATTGAGATCCATTGAAATTTTGCCGAATGCTTGAACCGCATCTACATGAAAATGTGCCTTAGGATAATTTTTTATAACTTTAGCCATTTGTGGAATAGGCTGTATTTGTCCAGTTACATTATTtacatacatac   >  1:51576/1‑141 (MQ=255)
                                                                                                                                          tatGTTATTGAGATCCATTGAAATTTTGCCGAATGCTTGAACCGCATCTACATGAAAATGTGCCTTAGGATAATTTTTTATAACTTTAGCCATTTGTGGAATAGGCTGTATTTGTCCAGTTACATTATTtacatacataca  >  1:297431/1‑141 (MQ=255)
                                                                                                                                          |                                                                                                                                            
ATTAACTGTTCCACTTCTAACGCCATATTCTTGACCACCACCATGGACAGTTGGTTCAACATTTTGAATGTGATTTACAAGTAAGACGCCTTGTCCTTTTAAACCATTAAACTTGTGTCCACTTAAACTAATACTATCTATGTTATTGAGATCCATTGAAATTTTGCCGAATGCTTGAACCGCATCTACATGAAAATGTGCCTTAGGATAATTTTTTATAACTTTAGCCATTTGTGGAATAGGCTGTATTTGTCCAGTTACATTATTTACATACATACA  >  CP000730/1828340‑1828618

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