Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,833,354 C→T 11.2% C395Y (TGC→TAC)  yicI ← alpha‑D‑xyloside xylohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,833,3540CT11.2% 75.7 / 3.3 36C395Y (TGC→TAC) yicIalpha‑D‑xyloside xylohydrolase
Reads supporting (aligned to +/- strand):  ref base C (14/18);  new base T (2/2);  total (16/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.70e-01

GTACGCATAATGGTTATGCATTTTCTGCGGATCGGAACCGTCAAACCACTGAACATCAGTTGGGATACGTTCGCCAAAGTCGGTCTTAAAGCAATCAACGCCCATCGCGACCAGACCTTTCAGTTTGTCGGCGTACCATTTGCAGGCATCCGGATTGGTAAAGTCATAAATCGCCAGACCTGGCTGCCATTTATCCCACTGCCATAGCGAACCGTCCGGGCGTTTGAGTAAATAGCCTTTCTCTTGTAACTCTTTAAAGACGGGGGATTTTTGACCGATATAGGG  >  NC_000913/3833212‑3833496
                                                                                                                                              |                                                                                                                                              
gTACGCATAATGGTTATGCATTTTCTGCGGATCGGAACCGTCAAACCCCTGAACATCAGTTGGGATACGTTCGCCAAAGTCGGTCTTAAAGCAAGCAACGCCCATCGCGACCAGACCTTTCAGTTTGTCGGCGTACCATTTGc                                                                                                                                                <  1:447997/143‑1 (MQ=255)
 ttcGCATAATGGTTATTTATTTTCTGCGGATCGGAACCGTCAAACCACTGAACATCAGTTGGGATACGTTCGCCAAAGTCGGTCTTAAAGCAATCAACGCCCATCGCGACCAGACCTTTCAGTTTGTCGGCGTACCATTTGCa                                                                                                                                               <  1:325693/141‑1 (MQ=255)
 tACGCATAATGGTTATGCATTTTCTGCGGATCGGAACCGTCAAACCACTGAACATCAGTTGGGATCCGCTCGCCGCAGTCGGTCTTACAGCACTCAACGCCCATCGCGACCAGCCCTTTCAGTTTGTCGGCGTACAATTTGCa                                                                                                                                               >  1:175734/1‑143 (MQ=255)
  aCGCATAATGGTTATGCATTTTCTGCGGATCGGAACCGTCAAACCACTGAACATCAGTTGGGATACGTTCGCCAAAGTCGGTCTTAAAGCAATCAACGCCCATCGCGACCAGACCTTTCAGTTTGTCGGCGTACCATTTGCAg                                                                                                                                              >  1:72437/1‑143 (MQ=255)
    gCATAATGGTTATGCATTTTCTGCGGATCGGAACCGTCAAACCACTGAACATCAGTTGGGATACGTTCGCCCAAGTCGGTCTTAAAGCAATCAACGCCCATCGCGACCAGCCCTTCCAGATCGTCGGCGTCACATTCGccagg                                                                                                                                            >  2:152895/1‑139 (MQ=255)
      aTAATGGTTATGCATTTTCTGCGGATCGGAACCGTCAAACCACTGAACATCAGTTGGGATACGTTCGCCAAAGTCGGTCTTAAAGCAATCAACGCCCATCGCGACCAGACCTTTCAGTTTGTCGGCGTACCATTTGCAGGCAt                                                                                                                                          <  1:360639/143‑1 (MQ=255)
             ttATGCATTTTCTGCGGATCGGAACCGTCAAACCACTGAACATCAGTTGGGATACGTTCGCCAAAGTCGGTCTTAAAGCAATCAACGCCCATCGCGACCAGACCTTTCAGTTTGTCGGCGTACCACTTGCAGGCATCCGGAtt                                                                                                                                   >  1:64850/1‑143 (MQ=255)
              gaggCATTTTCTGCGGATCGGAACCGTCAAACCACTGAACATCAGTTGGGATACGTTCGCCAAAGTCGGTCTTAAAGCAATCAACGCCCATCGCGACCAGACCTTTCAGTTTGTCGGCGTACCATTTGCAGGCATCCGGATTg                                                                                                                                  <  1:72350/140‑1 (MQ=255)
                 gCATTTTCTGCGGATCGGAACCGTCAAACCACTGAACATCAGGTGGGATACGTTCGCCAAAGTCGGTCTTAAAGCAATCAACGCCCATCGCGACCAGACCTTTCAGTTGGTCGGCGTACCATTTGCAGGCATCCGGATAGGTa                                                                                                                               >  2:60478/1‑143 (MQ=255)
                  cATTTTCTGCGGATCGGAACCGTCAAACCACTGAACATCAGTTGGGATACGTTCGCCAAAGTCGGTCTTAAAGCAATCAACGCCCATCGCGACCAGACATATCAGTTTGTAGGCGTACAATTTGCAGGCATCCGGATTGGTaa                                                                                                                              >  1:70408/1‑143 (MQ=255)
                   atatgcggCGGAACGGAACCGGCAAACCACGGAACATCAGTTGGGATACGTTCGCCAAAGGCGGTCTTAAAGCAATCAACGCCCATCGCGACCAGACCTTTCAGTTTGTCGGCGTACCATTTGCAGGCATCCGGATTGGTaaa                                                                                                                             <  1:146739/136‑1 (MQ=255)
                        gttcgtaTCGGTACCGCCACACCCTTGAATACTATTTGGCATACGTCCGCCAAAGTCGTTCTTAAAGCAATCAATGTCCATCTCGACCAGACCTTTCAGTTTGTCGGCGTACCATTTGCAGGCATCCGGATTGGTAAAGTCAt                                                                                                                        <  2:181014/137‑1 (MQ=255)
                            ggatcggatcCGTCAAACCACTGAACATCAGTTGGGATACGTTCGCCAAAGTCGGTCTTAAAGCAATCAACGCCCATCGCGACCAGACCTTTCAGTTTGTCGGCGTACCATTTGCAGGCATCCGGATTGGTAAAGTCATAAAt                                                                                                                    <  2:61677/143‑1 (MQ=255)
                               tCGGAACCGTCAAACCACTGAACATCAGTTGGGATACGTTCGCCAAAGTCGGTCTTAAAGCAATCAACGCCCATCGCGACCAGACCTTTCAGTTTGTCGGCGTACCATTTGCAGGCATCCGGATTGGTAAAGTCATAAATCGc                                                                                                                 <  2:210441/143‑1 (MQ=255)
                                    aCCGTCAAACCACTGAACATCAGTTGGGATACGTTCGCCAAAGTCGGTCTTAAAGCAATCAACGCCCATCGCGACCAGACCTTTCAGTTTGTCGGCGTACCATTTGCAGGCATCCGGATTGGTAAAGTCATAACCCGGCAGAc                                                                                                            <  2:439655/143‑1 (MQ=255)
                                              cACTGAACATCAGTTGGGATACATTCGCCAAAGACGGTCTTAAAGCAATCAACGCCCATCGCGACCAGACCTTTCAGTTTGTCGGCGAACCATTTGCAGGCATCCGGATTGGAAAAGTCATAAATCGCCAGACCTGGCTGccg                                                                                                  >  2:169922/1‑142 (MQ=255)
                                                   gaCATCAGTTGGGATACGTTCGCCAAAGGCGGTCTTAAAGCAATCAACGCCCATCGCGACCAGACCTTTCAGTTTGTCGGCGTACCATTTGCAGGCATCCGGATTGGTAAAGTCATAAATCGCCAGACCTGGCTGCCATTTAt                                                                                             <  2:467700/142‑1 (MQ=255)
                                                         aGTTGGGATACGTTCGCCAAAGTCGGTCTGAAAGCAATCAACGCCCATCGCGACCAGACCTTTCAGTTTGTCGGCGTACCATTTGCAGGCATCCGGATTGGTAAAGTCATAAATCACCAGACCTGGCTGCCATTTATCCCACt                                                                                       <  2:366522/143‑1 (MQ=255)
                                                          gTTGGGATACGTTCGCCAAAGTCGGTCTTAAAGCAATCAACGCCCATCGCGACCAGACCTTTCAGTTTGTCGGCGTACCATTTGCAGGCATCCGGATTGGTAAAGTCATAAATCGCAAGCGTGGGTTGCCATTTATCCCAcca                                                                                      >  2:183655/1‑141 (MQ=255)
                                                              ggctACGTTCGCCAAAGTCGGTCTTAAAGCAATCAACGCCCATCGCGACCAGACCTTTCAGTTTGTCGGCGTACCATTTGCAGGCATCCGGATTGGTAAAGTCATAAATCGCCAGACCTGGCTGCCATTTATCCCACTGCCAt                                                                                  >  2:190776/4‑143 (MQ=255)
                                                                  aCGGTCGCCAAAGTCGGTCTTAAAGCAATCACCGCCCCTCGCGACCATACCGTTCAGTTTGTCGGCGTACCATTTGCAGGCATCCGGATTGGTAAAGTCATAAATCGCCAGACCTGGCTGCCATTTATCCCACTGCCATAGCg                                                                              <  2:406375/143‑1 (MQ=255)
                                                                     ttCGCCAAAGTCGGTCTTAAAGCAATCAACGCCCATCGCGACCAGACCTTTCAGTTTGTCGGCGTACCATTTGCAGGCATCCGGATTGGTAAAGTCATAAATCGCCAGACCTGGCTGCCATTTATCCCACTGCCATAGCGAAc                                                                           >  1:449391/1‑143 (MQ=255)
                                                                           acactCGGTCTTAAAGAAATCAACGCCCATCGCGACCAAACCTTTCAGTTTGTAGACGTACCATTTGCAGGCATCCGGATTGGTAAAGTCATAAATCGCCAGACCAAGCTACCATTTATCCCACTGCCAAAGCGACCAGTCCg                                                                     >  2:81844/5‑143 (MQ=255)
                                                                                       aaaGCAATCAACGCCCATCGCGACCAGACCTTTCAGTTTGTCGGCGTACCATTTGTAGGCATCCGGATTGGTAAAGTCATAAATCGCCAGACCTGGCTGCCATTTATCCCACTGCCATAGCGAACCGTCCGGGCGTTTGAGTa                                                         >  1:178754/1‑143 (MQ=255)
                                                                                       aaaGCAATCAACGCCCATCGCGACCAGACCTTGCAGTTTGTCGGCGTACCATTTGTAGGCATCCGGATTGGTAAAGTCATAAATCGCCAGACCTGGCTGCCATTTATCCCACTGCCATAGCGAACCGTCCGGGCGTTTGAGTa                                                         >  1:424212/1‑143 (MQ=255)
                                                                                          gCAATCAACGCCCATCGCGACCAGACCTTTCAGTTTGTCGGCGTACCATTTGCAGGCATCCGGATTGGTAAAGTCATAAATCGCCAGACCTGGCTGCCATTTATCCCACTGCCATAGCGAACCGTCCGGGCGTTTGAGTaact                                                      >  1:284374/1‑141 (MQ=255)
                                                                                           cAATCAACGCCCATCGCGACCAGACCTTTCAGTTTGTCGGCGTACCATTTGCAGGCATCCGGATTGGTAAAGTCATAAATCGCCAGACCTGGCTGCCATTTCTGCGAATACCTAATCGAACCGTCCGGGCGTTTGAGTAAATa                                                     >  1:443915/1‑143 (MQ=255)
                                                                                                      cATCGCGACCAGACCTTTCAGTTTGTCGGCGTACCATTTGCAGGCATCCGGATTGGTAAAGTCATAAATCGCCAGACCTGGCTGCCATTTATCCCACTGCCATAGCGAACCGTCCGGGCGTTTGAGTAAATAGCCTTTCTCtt                                          <  1:156879/143‑1 (MQ=255)
                                                                                                            gACCAGACCTTTCAGTTTGTCGGCGTACCATTTGCAGGCATCCGGATTGGTAAAGTCATAAATCGCCAGACCTGGCTGCCATTTATCCCACTGCCATAGCGAACCGTCCGGGCGTTTGAGTAAATAGCCTTTCTCTTGTAAct                                    <  2:137479/143‑1 (MQ=255)
                                                                                                            caCCAGACCTTTCAGTTTGTCGGCGTACCATTTGCAGGCGTCCGGATTGGTACAGTCATACATCGCCAGACCTGGCTGCCATTTATCCCACTGCCATAGCGAACCGTCCGGGCGTTTGAGTAAATAGCCTTTCTCTTGTAAct                                    <  2:154596/142‑1 (MQ=255)
                                                                                                                 gACCTTTCAGTTTGTCGGCGTACCATTTGCAGGCATCCGGATTGGTAAAGTCATAAATCGCCAGACCTGGCTGCCATTTATCCCACTGCCATAGAGAACCGTACGGGCGTATGAGTAAATAGACTTTCTCTTGTAACACTTTa                               >  1:3673/1‑143 (MQ=255)
                                                                                                                      tgcAGTTTGTCGGCGTACCATTTGCAGGCATCCGGATTGGTAAAGTCATCAATCGCCAGACCTGGCTGCCATTTATCCCACTGCCATAGCGAACCGTCCGGGCGTTTGAGTAAATAGCCTTTCTCTTGTAACTCTTTAAAGAc                          <  1:29708/141‑1 (MQ=255)
                                                                                                                           tttGTCGGCGTACCATTTGCAGGCATCCGGATTGGTAAAGTCATAAATCGCCAGACCTGGCTGCCATTTATCCCACTGCCATAGCGAACCGTCCGGGCGTTTGAGTAAATAGCCTGTCTCTTGTAACTCTTTAAAGACggggg                     >  2:418424/1‑143 (MQ=255)
                                                                                                                             tGTCGGCGTACCATTTGCAGGCATCCGGATTGGTAAAGTCATAAATCGCCAGACCTGGCTGCCATTTATCCCACTGCCATAGCGAACCGTCCGGGCGTTTGAGTAAATAGCCTTTCTCTTGTAACTCTTTAAAGACGGGGGAt                   <  1:160394/143‑1 (MQ=255)
                                                                                                                                 ggCGTACCGTTTGCAGGCATCCGGCTTGTTACAGTCATCAAGCGCCAGACCTGGCTGCCATTTATCCCACTGCCATAGCGAACCGTCCGGGCGTTTGAGTAAATAGCCTTTCTCTTGTAACTCTTTAAAGACGGGGGAttttt               <  2:369283/143‑1 (MQ=255)
                                                                                                                                    gTACCATTTGCAGGCATCCGGATTGGTAAAGTCATAAATCGCCAGACCTGGCTGCCATTTATCCCACTGCCATAGCGAACCGTCCGGGCGTTTGAGTAAATAGCCTTTCTCTTGTAACTCTTTAAAGACGGGGGATTTTTGAc            <  2:449391/143‑1 (MQ=255)
                                                                                                                                           ttGTAGGCATCCGGATTGGTAAAGTCATAAATCGCCAGACCTGGCTGCCATTTATCCCACTGCCATAGCGAACCGTCCGGGCGTTTGAGTAAATAGCCTTTCTCTTGTAACTCTTTAAAGACGGGGGATTTTTGACCGAtata     <  2:178754/143‑1 (MQ=255)
                                                                                                                                           ttGTAGGCATCCGGATTGGTAAAGTCATAAATCGCCAGACCTGGCTGCCATTTATCCCACTGCCATAGCGAACCGTCCGGGCGTTTGAGTAAATAGCCTTTCTCTTGTAACTCTTTAAAGACGGGGGATTTTTGACCGAtata     <  2:424212/143‑1 (MQ=255)
                                                                                                                                              cAGGCGTCCGGCTTGGTAAAGTCGTAAAGCGCCAGACCTGGCTGCCATTTATCCCACTGCCAGAGCGAACTGTCGGGGCGTTTGAGTAAATAGCTTGTCTCTTGTAACGCTTTAAAGACGGGGGATTTTTGACCGATATAggg  <  2:41022/143‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                              
GTACGCATAATGGTTATGCATTTTCTGCGGATCGGAACCGTCAAACCACTGAACATCAGTTGGGATACGTTCGCCAAAGTCGGTCTTAAAGCAATCAACGCCCATCGCGACCAGACCTTTCAGTTTGTCGGCGTACCATTTGCAGGCATCCGGATTGGTAAAGTCATAAATCGCCAGACCTGGCTGCCATTTATCCCACTGCCATAGCGAACCGTCCGGGCGTTTGAGTAAATAGCCTTTCTCTTGTAACTCTTTAAAGACGGGGGATTTTTGACCGATATAGGG  >  NC_000913/3833212‑3833496

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