Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,541,456 A→G I75T (ATT→ACT)  cysW ← sulfate/thiosulfate ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,541,4560AG87.3% 147.9 / 12.0 55I75T (ATT→ACT) cysWsulfate/thiosulfate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (4/3);  new base G (25/23);  total (29/26)
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.79e-01

GTCCAGTAGCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTG  >  NC_000913/2541359‑2541550
                                                                                                 |                                                                                              
gTCCAGTAGCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTc                                                                                              <  2:25129/100‑1 (MQ=255)
   cAGTAGCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAgc                                                                                           >  1:1888658/1‑100 (MQ=255)
   cAGTAGCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGACAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAgc                                                                                           >  1:1244320/1‑100 (MQ=255)
    aGTAGCGTCATCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGCACACCAGGTTTACCGGTACGGCAGTCAgcg                                                                                          <  1:1279385/100‑1 (MQ=255)
         cGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGatc                                                                                     >  2:596855/1‑100 (MQ=255)
         cGTCAGCAGTAACTGGCGTCAAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGatc                                                                                     >  1:688197/1‑100 (MQ=255)
          gTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGatca                                                                                    <  2:1657463/100‑1 (MQ=255)
          gTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGatca                                                                                    >  1:880862/1‑100 (MQ=255)
          gTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGatca                                                                                    <  1:559976/100‑1 (MQ=255)
          gTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGatca                                                                                    <  2:1541550/100‑1 (MQ=255)
             agcagTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGAtcatca                                                                                 <  2:831490/100‑1 (MQ=255)
                  tAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTc                                                                            <  2:873731/100‑1 (MQ=255)
                       ggCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCa                                                                       >  1:1316207/1‑100 (MQ=255)
                         cGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGa                                                                     >  2:973951/1‑100 (MQ=255)
                             cAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGc                                                                 >  2:243835/1‑100 (MQ=255)
                              aGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCg                                                                >  2:1029128/1‑100 (MQ=255)
                               gggAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGt                                                               >  2:1409728/1‑100 (MQ=255)
                                ggAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTg                                                              >  1:468212/1‑100 (MQ=255)
                                          gcgcGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGGGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTcc                                                    <  1:1038290/100‑1 (MQ=255)
                                            gcgTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCgg                                                  >  2:550119/1‑100 (MQ=255)
                                             cgTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGa                                                 <  2:359377/100‑1 (MQ=255)
                                               tCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATc                                               >  2:1765116/1‑100 (MQ=255)
                                               tCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATc                                               <  2:326148/100‑1 (MQ=255)
                                               tCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATc                                               <  2:880863/100‑1 (MQ=255)
                                                 aCCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCgg                                             <  1:420576/100‑1 (MQ=255)
                                                 aCCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCgg                                             <  2:1375112/100‑1 (MQ=255)
                                                  ccagccagGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGc                                            >  1:233427/1‑100 (MQ=255)
                                                  ccagccagGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGc                                            >  2:1788438/1‑100 (MQ=255)
                                                   cagccagGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGcc                                           <  2:1300893/100‑1 (MQ=255)
                                                   cagccagGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGcc                                           <  2:1933319/100‑1 (MQ=255)
                                                   cagccagGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGcc                                           >  2:1372901/1‑100 (MQ=255)
                                                    agccagGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCa                                          >  2:1084036/1‑100 (MQ=255)
                                                     gccaggccagCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAg                                         <  2:39410/100‑1 (MQ=255)
                                                       caggccagCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGAt                                       >  1:256669/1‑100 (MQ=255)
                                                         ggccagCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTTAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTc                                     <  1:1362961/100‑1 (MQ=255)
                                                              tccgatctccGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTaa                                <  1:1157300/92‑1 (MQ=255)
                                                                 gAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAAc                             >  2:1154157/1‑100 (MQ=255)
                                                                  aaTGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAAcc                            >  1:531662/1‑100 (MQ=255)
                                                                      ccGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTaaagatcgga                        >  2:1157302/1‑93 (MQ=255)
                                                                        gAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATc                      <  1:596853/100‑1 (MQ=255)
                                                                           cacCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGc                   >  2:1063307/1‑100 (MQ=255)
                                                                             ccAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGccc                 >  2:1019247/1‑100 (MQ=255)
                                                                              cAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGcccc                >  1:1117160/1‑100 (MQ=255)
                                                                                ggTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCtt              <  1:973950/100‑1 (MQ=255)
                                                                                   ttACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCt           <  1:351807/100‑1 (MQ=255)
                                                                                    tACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTg          >  1:379355/1‑100 (MQ=255)
                                                                                    tACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTg          <  2:1240518/100‑1 (MQ=255)
                                                                                      ccGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGaa        >  2:487435/1‑100 (MQ=255)
                                                                                        ggTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATg      <  2:169803/100‑1 (MQ=255)
                                                                                        ggTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATg      >  2:828636/1‑100 (MQ=255)
                                                                                        ggTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGAGGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATg      >  1:1037455/1‑100 (MQ=255)
                                                                                           aCGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGcc    <  1:1765105/99‑1 (MQ=255)
                                                                                           aCGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCt   <  2:468214/100‑1 (MQ=255)
                                                                                            cGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTg  <  2:488551/100‑1 (MQ=255)
                                                                                            cGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTg  >  2:1348155/1‑100 (MQ=255)
                                                                                                 |                                                                                              
GTCCAGTAGCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTG  >  NC_000913/2541359‑2541550

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