Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 503,187 C→T A18T (GCC→ACC)  ybaL ← putative transporter YbaL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913503,1870CT100.0% 140.2 / NA 42A18T (GCC→ACC) ybaLputative transporter YbaL
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (19/23);  total (19/23)

AGTTCCGGGGCAAGCTTGGTATCGGCAACAAAGCCCGGAGTGAATGGTCCTGCCAGCACACCCGCTAACAGATATCCCACCAGAGGAGAAATACGTAGTTTATTGGCCAGCATGCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCAAAATTACAGGTAATAGTTTATGACAATTTCA  >  NC_000913/503060‑503319
                                                                                                                               |                                                                                                                                    
aGTTCCGGGGCAAGCTTGGTATCGGCAACAAAGCCCGGAGTGAATGGTCCTGCCAGCACACCCGCTAACAGATATCCCACCAGAGGAGAAATACGTAGTTTATTGGCCAGCATGCCGAGGATAAAGGTGAGCACAAGgcc                                                                                                                          <  1:17674/140‑1 (MQ=255)
aGTTCCGGGGCAAGCTTGGTATCGGCAACAAAGCCCGGAGTGAATGGTCCTGCCAGCACACCCGCTAACAGATATCCCACCAGAGGAGAAATACGTAGTTTATTGGCCAGCATGCCGAGGATAAAGGTGAGCACAAGgcc                                                                                                                          <  2:650155/140‑1 (MQ=255)
  ttCCGGGGCAAGCTTGGTATCGGCAACAAAGCCCGGAGTGAATGGTCCTGCCAGCACACCCGCTAACAGATATCCCACCAGAGGAGAAATACGTAGTTTATTGGCCAGCATGCCGAGGATAAAGGTGAGCACAAGgccgc                                                                                                                        >  1:763728/1‑140 (MQ=255)
          cAAGCTTGGTATCGGCAACAAAGCCCGGAGTGAATGGTCCTGCCAGCACACCCGCTAACAGATATCCCACCAGAGGAGAAATACGTAGTTTATTGGCCAGCATGCCGAGGATAAAGGTGAGCACAAGGCCGCCAACAAtg                                                                                                                <  2:502376/140‑1 (MQ=255)
            aGCTTGGTATCGGCAACAAAGCCCGGAGTGAATGGTCCTGCCAGCACACCCGCTAACAGATATCCCACCAGAGGAGAAATACGTAGTTTATTGGCCAGCATGCCGAGGATAAAGGTGAGCACAAGGCCGCCAACAAtggt                                                                                                              >  2:696151/1‑140 (MQ=255)
            aGCTTGGTATCGGCAACAAAGCCCGGAGTGAATGGTCCTGCCAGCACACCCGCTAACAGATATCCCACCAGAGGAGAAATACGTAGTTTATTGGCCAGCATGCCGAGGATAAAGGTGAGCACAAGGCCGCCAACAAtggt                                                                                                              <  1:696151/140‑1 (MQ=255)
                       ggCAACAAAGCCCGGAGTGAATGGTCCTGCCAGCACACCCGCTAACAGATATCCCACCAGAGGAGAAATACGTAGTTTATTGGCCAGCATGCCGAGGATAAAGGTGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCg                                                                                                   >  2:281365/1‑140 (MQ=255)
                        gCAACAAAGCCCGGAGTGAATGGTCCTGCCAGCACACCCGCTAACAGATATCCCACCAGAGGAGAAATACGTAGTTTATTGGCCAGCATGCCGAGGATAAAGGTGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCgg                                                                                                  >  1:27379/1‑140 (MQ=255)
                         caacaaAGCCCGGAGTGAATGGTCCTGCCAGCACACCCGCTAACAGATATCCCACCAGAGGAGAAATACGTAGTTTATTGGCCAGCATGCCGAGGATAAAGGTGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCggg                                                                                                 >  1:13290/1‑140 (MQ=255)
                         caacaaAGCCCGGAGTGAATGGTCCTGCCAGCACACCCGCTAACAGATATCCCACCAGAGGAGAAATACGTAGTTTATTGGCCAGCATGCCGAGGATAAAGGTGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCggg                                                                                                 <  1:566923/140‑1 (MQ=255)
                          aacaaAGCCCGGAGTGAATGGTCCTGCCAGCACACCCGCTAACAGATATCCCACCAGAGGAGAAATACGTAGTTTATTGGCCAGCATGCCGAGGATAAAGGTGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCgggg                                                                                                <  1:501486/140‑1 (MQ=255)
                          aacaaAGCCCGGAGTGAATGGTCCTGCCAGCACACCCGCTAACAGATATCCCACCAGAGGAGAAATACGTAGTTTATTGGCCAGCATGCCGAGGATAAAGGTGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCgggg                                                                                                <  1:306512/140‑1 (MQ=255)
                              aaGCCCGGAGTGAATGGTCCTGCCAGCACACCCGCTAACAGATATCCCACCAGAGGAGAAATACGTAGTTTATTGGCCAGCATGCCGAGGATAAAGGTGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGc                                                                                            <  2:486241/140‑1 (MQ=255)
                                 cccGGAGTGAATGGTCCTGCCAGCACACCCGCTAACAGATATCCCACCAGAGGAGAAATACGTAGTTTATTGGCCAGCATGCCGAGGATAAAGGTGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTg                                                                                         >  2:445263/1‑140 (MQ=255)
                                              gTCCTGCCAGCACACCCGCTAACAGATATCCCACCAGAGGAGAAATACGTAGTTTATTGGCCAGCATGCCGAGGATAAAGGTGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATGCATTCCGtct                                                                            <  1:445263/140‑1 (MQ=255)
                                              gTCCTGCCAGCACACCCGCTAACAGATATCCCACCAGAGGAGAAATACGTAGTTTATTGGCCAGCATGCCGAGGATAAAGGTGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATGCATTCCGtct                                                                            <  1:281365/140‑1 (MQ=255)
                                                     cAGCACACCCGCTAACAGATATCCCACCAGAGGAGAAATACGTAGTTTATTGGCCAGCATGCCGAGGATAAAGGTGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTc                                                                     <  2:763728/140‑1 (MQ=255)
                                                        cacaCCCGCTAACAGATATCCCACCAGAGGAGAAATACGTAGTTTATTGGCCAGCATGCCGAGGATAAAGGTGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCtg                                                                  >  1:57923/1‑140 (MQ=255)
                                                                cTAACAGATATCCCACCAGAGGAGAAATACGTAGTTTATTGGCCAGCATGCCGAGGATAAAGGTGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTAt                                                          >  2:18634/1‑140 (MQ=255)
                                                                  aaCAGATATCCCACCAGAGGAGAAATACGTAGTTTATTGGCCAGCATGCCGAGGATAAAGGTGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTg                                                        <  2:27379/140‑1 (MQ=255)
                                                                       ataccccACCAGAGGAGAAATACGTAGTTTATTGGCCAGCATGCCGAGGATAAAGGTGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTc                                                                     >  1:393449/5‑122 (MQ=255)
                                                                       atacCCCACCAGAGGAGAAATACGTAGTTTATTGGCCAGCATGCCGAGGATAAAGGTGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTc                                                                     <  2:393449/122‑1 (MQ=255)
                                                                                   aggagAAATACGTAGTTTATTGGCCAGCATGCCGAGGATAAAGGTGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGa                                       >  1:107778/1‑140 (MQ=255)
                                                                                    ggagAAATACGTAGTTTATTGGCCAGCATGCCGAGGATAAAGGTGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGaa                                      <  2:107778/140‑1 (MQ=255)
                                                                                     gagaAATACGTAGTTTATTGGCCAGCATGCCGAGGATAAAGGTGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGaaa                                     >  1:614462/1‑140 (MQ=255)
                                                                                               tAGTTTATTGGCCAGCATGCCGAGGATAAAGGTGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCAAAAtt                            <  2:534652/139‑1 (MQ=255)
                                                                                               tAGTTTATTGGCCAGCATGCCGAGGATAAAGGTGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCAAAAtt                            >  1:534652/1‑139 (MQ=255)
                                                                                               tAGTTTATTGGCCAGCATGCCGAGGATAAAGGTGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCAAAATTa                           >  2:417024/1‑140 (MQ=255)
                                                                                                 gTTTATTGGCCAGCATGCCGAGGATAAAGGTGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCAAAATTACa                         <  2:8410/140‑1 (MQ=255)
                                                                                                     aTTGGCCAGCATGCCGAGGATAAAGGTGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCAAAATTACAGGTa                     >  1:537319/1‑140 (MQ=255)
                                                                                                      ttGGCCAGCATGCCGAGGATAAAGGTGAGCACAAGGCCGCCAACAATTGTGGTGATAAGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCaaaa                              >  1:165331/1‑130 (MQ=255)
                                                                                                      ttGGCCAGCATGCCGAGGATAAAGGTGAGCACAAGGCCGCCAACAATTGTGGTGATAAGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCaaaa                              <  2:165331/130‑1 (MQ=255)
                                                                                                      ttGGCCAGCATGCCGAGGATAAAGGTGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCAAAATTACAGGTaa                    <  1:498945/140‑1 (MQ=255)
                                                                                                      ttGGCCAGCATGCCGAGGATAAAGGTGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCAAAATTACAGGTaa                    <  1:417024/140‑1 (MQ=255)
                                                                                                      ttGGCCAGCATGCCGAGGATAAAGGTGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCAAAATTACAGGTaa                    <  1:244530/140‑1 (MQ=255)
                                                                                                      ttGGCCAGCATGCCGAGGATAAAGGTGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCAAAATTACAGGTaa                    <  2:715227/140‑1 (MQ=255)
                                                                                                          ccAGCATGCCGAGGATAAAGGTGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCAAAATTACAGGTAATAGt                >  2:327196/1‑140 (MQ=255)
                                                                                                              cATGCCGAGGATAAAGGTGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCAAAATTACAGGTaa                    <  1:611751/132‑1 (MQ=255)
                                                                                                              cATGCCGAGGATAAAGGTGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCAAAATTACAGGTaa                    >  2:611751/1‑132 (MQ=255)
                                                                                                                   cGAGGATAAAGGTGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCAAGATTACAGGTAATAGTTTTTGACaa       >  2:744561/1‑140 (MQ=255)
                                                                                                                       gATAAAGGTGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATGCATTCTGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCAAAATTACAGGTAATAGTTTATGACAATTTc   >  2:596910/1‑140 (MQ=255)
                                                                                                                        aTAAAGGTGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCAAAATTACAGGTAATAGTTTATGACAATTTCa  <  2:57923/140‑1 (MQ=255)
                                                                                                                               |                                                                                                                                    
AGTTCCGGGGCAAGCTTGGTATCGGCAACAAAGCCCGGAGTGAATGGTCCTGCCAGCACACCCGCTAACAGATATCCCACCAGAGGAGAAATACGTAGTTTATTGGCCAGCATGCCGAGGATAAAGGCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCAAAATTACAGGTAATAGTTTATGACAATTTCA  >  NC_000913/503060‑503319

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