Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,009,257 T→C G327G (GGT→GGC yedE → UPF0394 family sulphur transport domain‑containing inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,009,2570TC97.2% 117.8 / NA 37G327G (GGT→GGCyedEUPF0394 family sulphur transport domain‑containing inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (18/17);  minor base A (0/1);  total (18/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.94e-01

GTTGAACCCAAAATCATGTGGGCGGGACCAAACGCGGTAATTGGTGGTTTACTGTTTGGTTTTGGCATCGTGCTGGCTGGCGGCTGCGAAACCGGCTGGATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGTCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATTTCGCTCCGGCGCTGGCCACCGACTGGGACAAAATCAACCTGCTGAAAACCTTTGGTCCGATGGGTGGCCTGCTGGTGACATATTTGCTGT  >  NC_000913/2009114‑2009396
                                                                                                                                               |                                                                                                                                           
gTTGAACCCAAAATCATGTGGGCGGGACCAAACGCGGTAATTGGTGGTTTACTGTTTGGTTTTGGCATCGTGCTGGCTGGCGGCTGCGAAACCGGCTGGATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCa                                                                                                                                      <  1:103573/151‑1 (MQ=255)
   gAACCCAAAATCATGTGGGCGGGACCAAACGCGGTAATTGGTGGTTTACTGTTTGGTTTTGGCATCGTGCTGGCTGGCGGCTGCGAAACCGGCTGGATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAAt                                                                                                                                    <  1:90650/150‑1 (MQ=255)
     aCCCAAAATCATGTGGGCGGGACCAAACGCGGTAATTGGTGGTTTACTGTTTGGTTTTGGCATCGTGCTGGCTGGCGGCTGCGAAACCGGCTGGATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAAtgtg                                                                                                                                 >  2:171297/1‑151 (MQ=255)
       ccAAAATCATGTGGGCGGGACCAAACGCGGTAATTGGTGGTTTACTGTTTGGTTTTGGCATCGTGCTGGCTGGCGGCTGCGAAACCGGCTGGATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAATGTGAt                                                                                                                               >  1:158066/1‑151 (MQ=255)
                 gtgGGCGGGACCAAACGCGGTAATTGGTGGTTTACTGTTTGGTTTTGGCATCGTGCTGGCTGGCGGCTGCGAAACCGGCTGGATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAACg                                                                                                                     <  2:234225/151‑1 (MQ=255)
                  tgGGCGGGACCAAACGCGGTAATTGGTGGTTTACTGTTTGGTTTTGGCATCGTGCTGGCTGGCGGCTGCGAAACCGGCTGGATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAACg                                                                                                                     >  2:451531/1‑150 (MQ=255)
                   gggCGGGACCAAACGCGGTAATTGGTGGTTTACTGTTTGGTTTTGGCATCGTGCTGGCTGGCGGCTGCGAAACCGGCTGGATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAACGAt                                                                                                                   <  1:523892/151‑1 (MQ=255)
                        ggACCAAACGCGGTAATTGGTGGTTTACTGTTTGGTTTTGGCATCGTGCTGGCTGGCGGCTGCGAAACCGGCTGGATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTg                                                                                                               >  2:328990/1‑150 (MQ=255)
                        ggACCAAACGCGGCAATTGGTGGTTTACTGTTTGGTTTTGGCATCGTGCTGGCTGGCGGCTGCGAAACCGGCTGGATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTg                                                                                                               >  1:210830/1‑150 (MQ=255)
                         gACCAAACGCGGTAATTGGTGGTTTACTGTTTGGTTTTGGCATCGTGCTGGCTGGCGGCTGCGAAACCGGCTGGATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTGGc                                                                                                             >  1:423952/1‑151 (MQ=255)
                         gACCAAACGCGGTAATTGGTGGTTTACTGTTTGGTTTTGGCATCGTGCTGGCTGGCGGCTGCGAAACCGGCTGGATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTGGc                                                                                                             <  2:178507/151‑1 (MQ=255)
                          aCCAAACGCGGTAATTGGTGGTTTACTGTTTGGTTTTGGCATCGTGCTGGCTGGCGGCTGCGAAACCGGCTGGATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTGGCg                                                                                                            <  1:93080/151‑1 (MQ=255)
                                                  aCTGTTTGGTTTTGGCATCGTGCTGGCTGGCGGCTGCGAAACCGGCTGGATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGACTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATTTCGCTCCg                                                                                    <  1:517826/151‑1 (MQ=255)
                                                     gTTTGGTTTTGGCATCGTGCTGGCTGGCGGCTGCGAAACCGGCTGGATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATTTCGCTCCGgc                                                                                  >  2:311957/1‑150 (MQ=255)
                                                     gTTTGGTTTTGGCATCGTGCTGGCTGGCGGCTGCGAAACCGGCTGGATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATTTCGCTCCGgc                                                                                  <  1:311924/150‑1 (MQ=255)
                                                        tGGTTTTGGCATCGTGCTGGCTGGCGGCTGCGAAACCGGCTGGATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATTTCGCTCCGGCGCt                                                                               <  2:490495/150‑1 (MQ=255)
                                                          gTTTTGGCATCGTGCTGGCTGGCGGCTGCGAAACCGGCTGGATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATTTCGCTCCGGCGCTGGc                                                                            >  1:671713/1‑151 (MQ=255)
                                                          gTTTTGGCATCGTGCTGGCTGGCGGCTGCGAAACCGGCTGGATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATTTCGCTCCGGCGCTGGc                                                                            >  2:337498/1‑151 (MQ=255)
                                                                gCATCGTGCTGGCTGGCGGCTGCGAAACCGGCTGGATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATTTCGCTCCGGCGCTGGCCTcc                                                                        <  1:488996/149‑1 (MQ=255)
                                                                  aTCGTGCTGGCTGGCGGCTGCGAAACCGGCTGGATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATTTCGCTCCGGCGCTGGCCACCGACt                                                                    <  1:10046/151‑1 (MQ=255)
                                                                          ggctggcGGCTGCGAAACCGGCTGGATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATTTCGCTCCGGCGCTGGCCACCGACTGGGACaa                                                             >  1:318759/1‑150 (MQ=255)
                                                                           gctggcGGCTGCGAAACCGGCTGGATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAATGATTCTGGCGTATTACTGGGATGATTTCGCTCCGGCGCTGGCCACCGACTGGGACaaa                                                            >  2:613769/1‑150 (MQ=255)
                                                                                   cTGCGAAACCGGCTGGATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATTTCGCTCCGGCGCTGGCCACCGACTGGGACAAAATCAACctg                                                   >  2:342609/1‑151 (MQ=255)
                                                                                          aCCGGCTGGATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATTTCGCTCCGGCGCTGGCCACCGACTGGGACAAAATCAACCTGCTGaaaa                                            <  2:424003/151‑1 (MQ=255)
                                                                                           ccGGCTGGATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATTTCGCTCCGGCGCTGGCCACCGACTGGGACAAAATCAACCTGCTGAAAAc                                           <  2:708582/151‑1 (MQ=255)
                                                                                             ggCTGGATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATTTCGCTCCGGCGCTGGCCACCGACTGGGACAAAATCAACCTGCTGAAAACCt                                         >  2:692903/1‑151 (MQ=255)
                                                                                               cTGGATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATTTCGCTCCGGCGCTGGCCACCGACTGGGACAAAATCAACCTGCTGAAAACCttt                                       >  1:661205/1‑151 (MQ=255)
                                                                                                     gTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATTTCGCTCCGGCGCTGGCCACCGACTGGGACAAAATCAACCTGCTGAAAACCTTTGGTcc                                  <  1:266654/150‑1 (MQ=255)
                                                                                                       aCCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATTTCGCTCCGGCGCTGGCCACCGACTGGGACAAAATCAACCTGCTGAAAACCTTTGGTCCGAt                               >  1:93695/1‑151 (MQ=255)
                                                                                                         cgcgcgGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGGCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATTTCGCTCCGGCGCTGGCCACCGACTGGGACAAAATCATCCTGCTGAAAACCTTTGGTCCGATgg                             <  1:465674/151‑1 (MQ=255)
                                                                                                             cgGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATTTCGCTCCGGCGCTGGCCACCGACTGGGACAAAATCAACCTGCTGAAAACCTTTGGTCCGATGGGTgg                         >  2:735445/1‑151 (MQ=255)
                                                                                                                 aGAAGGCCAGGTGCACTACTGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATTTCGCTCCGGCGCTGGCCACCGACTGGGACAAAATCAACCTGCTGAAAACCTTTGGTCCGATGGGTGGCctg                     >  2:728226/1‑151 (MQ=255)
                                                                                                                       ccAGGTGCACTACTGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATTTCGCTCCGGCGCTGGCCACCGACTGGGACAAAATCAACCTGCTGAAAACCTTTGGTCCGATGGGTGGCCTGCTGGTg               <  2:10772/151‑1 (MQ=255)
                                                                                                                               actactGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATTTCGCTCCGGCGCTGGCCACCGACTGGGACAAAATCAACCTGCTGAAAACCTTTGGTCCGATGGGTGGCCTGCTGGTGACATAttt       <  2:523623/151‑1 (MQ=255)
                                                                                                                                  actGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATTTCGCTCCGGCGCTGGCCACCGACTGGGACAAAATCAACCTGCTGAAAACCTTTGGTCCGATGGGTGGCCTGCTGGTGACATATTTGCt    >  2:369204/1‑151 (MQ=255)
                                                                                                                                  actGGTGGGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATTTCGCTCCGGCGCTGGCCACCGACTGGGACAAAATCAACCTGCTGAAAACCTTTGGTCCGATGGGTGGCCTGCTGGTGACATATTTGCt    <  1:375502/151‑1 (MQ=255)
                                                                                                                                     ggtggGTCGGCCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATTTCGCTCCGGCGCTGGCCACCGACTGGGACAAAATCAACCTGCTGAAAACCTTTGGTCCGATGGGTGGCCTGCTGGTGACATATTTGCtgt  <  1:620610/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                           
GTTGAACCCAAAATCATGTGGGCGGGACCAAACGCGGTAATTGGTGGTTTACTGTTTGGTTTTGGCATCGTGCTGGCTGGCGGCTGCGAAACCGGCTGGATGTACCGCGCGGTAGAAGGCCAGGTGCACTACTGGTGGGTCGGTCTGGGCAATGTGATCGGCTCAACGATTCTGGCGTATTACTGGGATGATTTCGCTCCGGCGCTGGCCACCGACTGGGACAAAATCAACCTGCTGAAAACCTTTGGTCCGATGGGTGGCCTGCTGGTGACATATTTGCTGT  >  NC_000913/2009114‑2009396

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