Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 312,189 C→T G26G (GGC→GGT ykgL → uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913312,1890CT98.2% 167.3 / ‑5.3 54G26G (GGC→GGTykgLuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (28/25);  minor base G (1/0);  total (29/25)
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.97e-01

AGTACCAACAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCA  >  NC_000913/312064‑312315
                                                                                                                             |                                                                                                                              
aGTACCAACAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTcaac                                                                                                                   >  2:338283/1‑139 (MQ=255)
 gTACCAACAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTcaaca                                                                                                                  >  2:726623/1‑139 (MQ=255)
 gTACCAACAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGATGATGCCAGCATACAGGGGTGAAGATGGTcaaca                                                                                                                  >  1:711817/1‑139 (MQ=255)
     cAACAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTc                                                                                                              <  2:711817/139‑1 (MQ=255)
             ttACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTaa                                                                                                      >  1:306949/1‑139 (MQ=255)
                 cacaTTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCTGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATa                                                                                                  <  1:726623/139‑1 (MQ=255)
                      ttATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACAt                                                                                                         >  2:535821/1‑127 (MQ=255)
                      ttATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACAt                                                                                                         <  1:535821/127‑1 (MQ=255)
                             cTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTca                                                                                      >  1:532335/1‑138 (MQ=255)
                              ttAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCgg                                                                                     >  2:507880/1‑139 (MQ=255)
                                aaCAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGt                                                                                                               >  1:727019/1‑111 (MQ=255)
                                aaCAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGt                                                                                                               <  2:727019/111‑1 (MQ=255)
                                aaCAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACAtt                                                                                                        <  1:563536/118‑1 (MQ=255)
                                aaCAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACAtt                                                                                                        >  2:563536/1‑118 (MQ=255)
                                aaCAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGgag                                                                                   >  1:650219/1‑139 (MQ=255)
                                aaCAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGgag                                                                                   >  2:677590/1‑139 (MQ=255)
                                aaCAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGgag                                                                                   >  1:319719/1‑139 (MQ=255)
                                  cAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGATGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGgagag                                                                                 <  2:709257/139‑1 (MQ=255)
                                       tGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTc                                                                                       <  2:735272/128‑1 (MQ=255)
                                       tGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGgagag                                                                                 <  2:71596/134‑1 (MQ=255)
                                       tGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGgagag                                                                                 >  1:71596/1‑134 (MQ=255)
                                       tGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACg                                                                            >  1:1102000/1‑139 (MQ=255)
                                       tGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGGGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTc                                                                                       >  1:735272/1‑128 (MQ=255)
                                        gAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGaaa                                                                              <  2:9363/136‑1 (MQ=255)
                                        gAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGaaa                                                                              >  1:9363/1‑136 (MQ=255)
                                           tGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCa                                                                         >  1:509147/1‑138 (MQ=255)
                                           tGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCa                                                                         <  2:509147/138‑1 (MQ=255)
                                                       cTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTc                                                            >  1:781915/1‑139 (MQ=255)
                                                       cTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTc                                                            >  1:714622/1‑139 (MQ=255)
                                                               cAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGCTGCCAGCATACAGGGGTCAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATTAAc                                                    <  2:714622/139‑1 (MQ=255)
                                                                  aGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCATCAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAgg                                                 <  2:650219/139‑1 (MQ=255)
                                                                  aGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAgg                                                 <  2:319719/139‑1 (MQ=255)
                                                                     gaggagACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGc                                                               <  1:450969/122‑1 (MQ=255)
                                                                     gaggagACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGc                                                               >  2:450969/1‑122 (MQ=255)
                                                                              gAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATAAAGGGGTGAAACTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATcc                                     >  2:686418/1‑139 (MQ=255)
                                                                                  ggCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAg                                 <  1:405787/139‑1 (MQ=255)
                                                                                  ggCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAg                                 >  1:177877/1‑139 (MQ=255)
                                                                                   gCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCTACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAgg                                <  1:686418/139‑1 (MQ=255)
                                                                                        ttCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAAt                           <  2:1102000/139‑1 (MQ=255)
                                                                                             tCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGga                                                <  1:706308/113‑1 (MQ=255)
                                                                                             tCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGga                                                >  2:706308/1‑113 (MQ=255)
                                                                                                aaaaGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTg                   <  2:532335/139‑1 (MQ=255)
                                                                                                   agagGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGgtg                <  2:306949/139‑1 (MQ=255)
                                                                                                      gggTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGACTATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAAc                                                    <  1:678065/100‑1 (MQ=255)
                                                                                                      gggTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAAc                                                    >  2:678065/1‑100 (MQ=255)
                                                                                                      gggTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAAc                                                    >  2:879076/1‑100 (MQ=255)
                                                                                                      gggTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAAc                                                    <  1:879076/100‑1 (MQ=255)
                                                                                                       ggTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCa                                                                         <  1:802515/78‑1 (MQ=255)
                                                                                                       ggTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCa                                                                         >  2:802515/1‑78 (MQ=255)
                                                                                                       ggTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGc                                                               <  1:243895/88‑1 (MQ=255)
                                                                                                       ggTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGc                                                               >  2:243895/1‑88 (MQ=255)
                                                                                                       ggTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAAc                                                    >  2:1036927/1‑99 (MQ=255)
                                                                                                       ggTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAAc                                                    <  1:1036927/99‑1 (MQ=255)
                                                                                                                 aGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCa  >  1:573871/1‑139 (MQ=255)
                                                                                                                             |                                                                                                                              
AGTACCAACAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCA  >  NC_000913/312064‑312315

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