Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,834,428 (G)6→7 coding (110/2319 nt) yicI ← alpha‑D‑xyloside xylohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,834,4221.G95.0% 131.0 / ‑0.3 40coding (116/2319 nt)yicIalpha‑D‑xyloside xylohydrolase
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base G (23/15);  total (25/15)
Fisher's exact test for biased strand distribution p-value = 5.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

AGCGGATAATGAGGACCGTTATTCAGCGCCCCCTGAAAATGCTCAATCCGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTT  >  NC_000913/3834287‑3834552
                                                                                                                                        |                                                                                                                                   
aGCGGATAATGAGGACCGTTATTCAGCGCCCCCTGAAAATGCTCAATCCGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACgggg                                                                                                                                 >  2:428895/1‑140 (MQ=255)
   ggATAATGAGGACCGTTATTCAGCGCCCCCTGAAAATGCTCAATCCGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACggggggg                                                                                                                             <  1:480526/140‑2 (MQ=255)
   ggATAATGAGGACCGTTATTCAGCGCCCCCTGAAAATGCTCAATCCGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACggggggg                                                                                                                             <  1:465846/140‑2 (MQ=255)
     aTAATGAGGACCGTTATTCAGCGCCCCCTGAAAATGCTCAATCCGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACgggggggca                                                                                                                           >  1:578628/1‑137 (MQ=255)
        aTGAGGACCGTTATTCAGCGCCCCCTGAAAATGCTCAATCCGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGgcagca                                                                                                                         <  2:615427/140‑1 (MQ=255)
               ccGTTATTCAGCGCCCCCTGAAAATGCTCAATCCGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGTTACGTTCCCGCTCATCACGGGGGGGCAGCATAGacca                                                                                                                  >  1:265999/1‑140 (MQ=255)
               ccGTTATTCAGCGCCCCCTGAAAATGCTCAATCCGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGGCAGCATAGacca                                                                                                                  >  1:649531/1‑140 (MQ=255)
                cGTTATTCAGCGCCCCCTGAAAATGCTCAATCCGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGGCAGCATAGaccac                                                                                                                 >  2:779543/1‑140 (MQ=255)
                  ttATTCAGCGCCCCCTGCAAATCCTCAATCCCCACCCCGCCAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGGCAGCATAGAccacca                                                                                                               <  2:238065/140‑1 (MQ=255)
                  ttATTCAGCGCCCCCTGAAAATGCTCAATCCGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGGCAGCATAGAccacca                                                                                                               >  1:540737/1‑140 (MQ=255)
                  ttATTCAGCGCCCCCTGAAAATGCTCAATCCGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACgggg                                                                                                                                 <  2:624318/122‑1 (MQ=255)
                  ttATTCAGCGCCCCCTGAAAATGCTCAATCCGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACgggg                                                                                                                                 >  1:624318/1‑122 (MQ=255)
                          cgcCCCCTGAAAATGCTCAATCCGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGGACATCACGGGGGGGAAGCATAGACCACCATTTCATTa                                                                                                       >  1:238065/1‑140 (MQ=255)
                          cgcCCCCTGAAAATGCTCAATCCGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGTGGTAGCATAGACCACCATTTCATTa                                                                                                       >  1:231176/1‑140 (MQ=255)
                          cgcCCCCTGAAAATGCTCAATCCGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGCGGCAGCATAGACCACCATTTCATTa                                                                                                       >  2:423061/1‑140 (MQ=255)
                            cccccTGAAAATGCTCAATCCGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGGCAGCATAGACCACCATTTCATTATc                                                                                                     >  2:715676/1‑140 (MQ=255)
                              cccTGAAAATGCTCAATCCGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGGCAGCATAGACCACCATTTCATTATCct                                                                                                   <  2:776511/140‑1 (MQ=255)
                                       tGCTCAATCCGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGGCAGCATCGACAACGATTTCAGTATCCTGCTGTTCaa                                                                                          >  1:268818/1‑140 (MQ=255)
                                       tGCTCAATCCGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCaa                                                                                          >  1:531267/1‑140 (MQ=255)
                                            aaTCCGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCg                                                                                     <  1:428895/140‑1 (MQ=255)
                                            aaTCCGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCg                                                                                     <  2:649531/140‑1 (MQ=255)
                                            aaTCCGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCg                                                                                     <  1:779543/140‑1 (MQ=255)
                                              tCCGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGaa                                                                                   >  1:822238/1‑140 (MQ=255)
                                                            aTACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGt                                                                     <  1:715676/140‑1 (MQ=255)
                                                                cTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAAt                                                                 >  1:534551/1‑140 (MQ=255)
                                                                 ttCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATc                                                                >  2:688208/1‑140 (MQ=255)
                                                                    cTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCaaa                                                             >  1:107164/1‑140 (MQ=255)
                                                                           gagaAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGc                                                      >  1:38490/1‑140 (MQ=255)
                                                                            agaAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGcc                                                     >  2:65526/1‑140 (MQ=255)
                                                                            agaAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGcc                                                     >  1:42578/1‑140 (MQ=255)
                                                                                      cAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATc                                           <  2:531267/140‑1 (MQ=255)
                                                                                        aCGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACACACATCACGGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCaa                                         >  2:128350/1‑140 (MQ=255)
                                                                                                   ggCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCGGGTTGAATCAACCAGTTTCCAt                              >  2:271722/1‑140 (MQ=255)
                                                                                                   ggCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCAt                              <  1:65526/140‑1 (MQ=255)
                                                                                                   ggCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGGGAAGCGACCTGAGGCCGAGTTAAATCAACCAGTTTCCAt                              >  2:404518/1‑140 (MQ=255)
                                                                                                     cGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAAGCAAATTGAGGCAAGGTTGAATAAACCAGTTTCCATCg                            >  1:323951/1‑140 (MQ=255)
                                                                                                        aTCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTa                         <  2:540737/140‑1 (MQ=255)
                                                                                                              cTGCCAGGTACGTTCACGCACATCACGGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTc                   >  2:143150/1‑140 (MQ=255)
                                                                                                               tGCCAGGTACGTTCACGCACATCACGGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGTGGCGAGGGTGAATCAAGCAGTATCCATCGCTAATTTTCa                  >  2:233264/1‑140 (MQ=255)
                                                                                                                 ccAGGTACGTTCACGCACATCACGGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCAtt                >  1:553792/1‑140 (MQ=255)
                                                                                                                     gTACGTTCACGCACATCTCGGGGGGGCAGCAAAGACCACCATTTCATTGTCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCg            <  2:268818/140‑1 (MQ=255)
                                                                                                                     gTACGTTCACGCACATCACGGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCg            <  1:423061/140‑1 (MQ=255)
                                                                                                                     gTACGTTCACGCACATCACGGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCg            <  2:265999/140‑1 (MQ=255)
                                                                                                                     gTACGTTCACGCACATCACGGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCg            <  1:688208/140‑1 (MQ=255)
                                                                                                                     gTACGTTCACACACATCACGGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCg            <  1:128350/140‑1 (MQ=255)
                                                                                                                          ttCACGCACATCACGGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCt       <  2:107164/140‑1 (MQ=255)
                                                                                                                               gCACATCACGGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGtt  >  2:827608/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                                   
AGCGGATAATGAGGACCGTTATTCAGCGCCCCCTGAAAATGCTCAATCCGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAATCAACCAGTTTCCATCGCTAATTTTCATTGCCGTTCCTTAGTT  >  NC_000913/3834287‑3834552

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: