Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,245,412 (A)7→8 coding (253/441 nt) ycgY → uncharacterized protein YcgY

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,245,4051.A91.7% 151.0 / 6.0 49coding (246/441 nt)ycgYuncharacterized protein YcgY
Reads supporting (aligned to +/- strand):  ref base . (2/2);  new base A (19/25);  total (22/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGATTATTAATTTAACTCCAGAAAGCAAAAAAACGTTGATGAGTGAACACCGTAACAATATTCAGGTAGCAATTGATGCTTTATATAGCCAACGTCGGCGTAAATCGCCGGGTGATGAAAGTGAAACACTTACCCGTAGAAATGACGCTATCTTCGGCAA  >  NC_000913/1245267‑1245539
                                                                                                                                            |                                                                                                                                        
tCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGATTATTAATTTAACTCCAGAAAGCa                                                                                                                                         >  1:226306/1‑140 (MQ=255)
     cTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGATTATTAATTTAACTCCAGAAAGCaaaaaa                                                                                                                                   >  1:986249/1‑140 (MQ=255)
     cTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGATTATTAATTTAACTCCAGAAAGCaaaaaa                                                                                                                                   >  2:1364689/1‑140 (MQ=255)
      tGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGATTATTAATTTAACTCCAGAAAGCaaaaaaa                                                                                                                                  >  1:248041/1‑140 (MQ=255)
        gCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGATTATTAATTTAACTCCAGAAAGCaaa                                                                                                                                      >  1:257611/1‑134 (MQ=255)
        gCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGATTATTAATTTAACTCCAGAAAGCaaa                                                                                                                                      <  2:257611/134‑1 (MQ=255)
          aGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGATTATGAATCTAGCTCCAGTAAGCaaaaaaaacgt                                                                                                                              >  1:448239/1‑136 (MQ=255)
            aTTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGATTATTAATTTAACTCCAGAAAGCaaaaaaaaacgtt                                                                                                                            <  1:269048/140‑7 (MQ=255)
             ttAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGATTATTAATTTAACTCCAGAAAGCAAAAAAAACGTtga                                                                                                                           <  2:986249/140‑1 (MQ=255)
             ttAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGATTATTAATTTAACTCCAGAAAGCAAAAAAAACGTtga                                                                                                                           <  1:772746/140‑1 (MQ=255)
             ttAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGATTATTAATTTAACTCCAGAAAGCAAAAAAAACGTtga                                                                                                                           <  2:762800/140‑1 (MQ=255)
             ttAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGATTATTAATTTAACTCCAGAAAGCAAAAAAAACGTtga                                                                                                                           >  2:803220/1‑140 (MQ=255)
              tAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGATTATTAATTTAACTCCAGAAAGCAAAAAAAACGTtgat                                                                                                                          >  2:384559/1‑140 (MQ=255)
                   cAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGATTATTAATTTAACTCCAGAAAGCAAAAAAAACGTTGATGAGTg                                                                                                                     <  1:266104/140‑1 (MQ=255)
                     aaaTACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGGTTATTAATTTAACTCCAGAAAGCAAAAAAAACGTTGATGAGTGaa                                                                                                                   <  2:1292364/140‑1 (MQ=255)
                     aaaTACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGATTATTAATTTAACTCCAGAAAGCAAAAAAAACGTTGATGAGTGaa                                                                                                                   <  1:1083433/140‑1 (MQ=255)
                      aaTACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGATTATTAATGTAACTCCGGAAGGCGAAAAAAGCGTTGATGAGTGAac                                                                                                                  >  1:394082/1‑140 (MQ=255)
                      aaTACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGATTATTAATGTAACTCCGGAAAGCCAAAAAAACGTTGATGAGTGAac                                                                                                                  >  1:498654/1‑140 (MQ=255)
                        tACCCGTTCGGCAATGCCTCAGTTACGCGAACAGGCATGGTCTGAATTTTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGATTATTAATTTAACTCCAGAAAGCAAAAAAAACGTTGATGAGTGAacac                                                                                                                <  2:498654/140‑1 (MQ=255)
                          cccGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGATTATTAATTTAACTCCAGAAAGCAAAAAAAACGTTGATGAGTGAACACCg                                                                                                              <  2:248041/140‑1 (MQ=255)
                           ccGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGATTATTAATTTAACTCCAGAAAGCAAAAAAAACGTTGATGAGTGAACACCGt                                                                                                             >  1:932191/1‑140 (MQ=255)
                             gTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGATTATTAATTTAACTCCAGAAAGCAAAAAAAACGTTGATGAGTGAACACCGTaa                                                                                                           <  1:1364689/140‑1 (MQ=255)
                             gTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGATTATTAATTTAACTCCAGAAAGCAAAAAAAACGTTGAGGAGTGAACACCGTaa                                                                                                           >  1:734925/1‑140 (MQ=255)
                               tCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGATTATTAATTTAACTCCAGAAAGCAAGAAAAACGTTGATGGGTGAACACCGTAACa                                                                                                         >  1:423950/1‑140 (MQ=255)
                                  gCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGATTATTAATTTAACTCCAGAAAGCAAAAAAAACGTTGATGAGTGAACACCGTAACAata                                                                                                      <  2:932191/140‑1 (MQ=255)
                                           cAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGATTATTAATTTAACTCCAGAAAGCAAAAAAAACGTTGATGAGTGAACACCGTAACAATATTCAGGTAg                                                                                             <  1:384559/140‑1 (MQ=255)
                                                aCGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGATTATTAATTTAACTCCAGAAAGCAAAAAAAACGTTGATGAGTGAACACCGTAACAATATTCAGGTAGCAAtt                                                                                        >  1:1124229/1‑140 (MQ=255)
                                                 cgcgAAGATGAATGGTCTGAATCCTTTTTTTTGGTCCTAAATTCTCTGCAGTGTCTGGATTATGTGATTATTAATTTAACTCCAGAAAGCAAAAAAAACGTTGATGAGTGAACACCGTAACAATATTCAGGTAGCAATTg                                                                                       <  2:734925/140‑1 (MQ=255)
                                                 cgcgAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGATTATTAATTTAACTCCAGAAAGCAAAAAAAACGTTGATGAGTGAACACCGTAACAATATTCAGGTAGCAATTg                                                                                       <  2:1124229/140‑1 (MQ=255)
                                                 cgcgAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGATTATTAATTTAACTCCAGAAAGCAAAAAAAACGTTGATGAGTGAACACCGTAACAATATTCAGGTAGCAATTg                                                                                       <  1:1253370/140‑1 (MQ=255)
                                                 cgcgAAGAGGAATGGTCCGAATCCTTTTTTTAGTTGCTACATTCTCTGTAGTGTCTGGATTATGTGATTATTAATTTAACTCCAGAAAGCAAAAAAAACGTTGATGAGTGAACACCGTAACAATATTCAGGTAGCAATTg                                                                                       <  2:423950/140‑1 (MQ=255)
                                                 caggAAGAGGAATGGTCTGAATCCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGATTATTAATTTAACTCCAGAAAGCAAAAAAAACGTTGATGAGTGAACACCGTAACAATATTCAGGTAGCAATTg                                                                                       <  1:553444/137‑1 (MQ=255)
                                                   cgAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGATTATTAATTTAACTCCAGAAAGCAAAAAAAACGTTGATGAGTGAACACCGTAACAATATTCAGGTAGCAATTg                                                                                       <  1:1296549/138‑1 (MQ=255)
                                                   cgAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGATTATTAATTTAACTCCAGAAAGCAAAAAAAACGTTGATGAGTGAACACCGTAACAATATTCAGGTAGCAATTg                                                                                       >  2:1296549/1‑138 (MQ=255)
                                                     aaGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGATTATTAATTTAACTCCAGAAAGCAAAAAAAACGTTGATGAGTGAACACCGTAACAATATTCAGGTAGCAATTGATGc                                                                                   <  2:226306/140‑1 (MQ=255)
                                                      agagGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGATTATTAATTTAACTCCAGAAAGCAAAAAAAACGTTGATGAGTGAACACCGTAACAATATTCAGGTAGCAATTGATGCt                                                                                  <  1:803220/140‑1 (MQ=255)
                                                               gTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGATTATTAATTTAACTCCAGAAAGCAAAAAAAACGTTGATGAGTGAACAGCGTAACAGTATTCAGGCAGCAATTGATGCTTtatatctg                                                                         >  2:553444/1‑137 (MQ=255)
                                                                tCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGATTATTAATTTAACTCCAGAAAGCAAAAAAAACGTTGATGAGTGAACACCGTAACAATAGTCAAGTCGCAATTAATGCTTTATATAGcc                                                                        >  1:593944/1‑140 (MQ=255)
                                                                                  ttGCTAAATTCTCTGGAGTGTCTGGATTATGTGATTATTAATTTAACTCCAGAAAGCAAAAAAAACGTTGATGAGTGAACACCGTAACAATATTCAGGTAGCAATTGATGCTTTATATAGCCAACGTCGGCGTAAATCGc                                                      <  1:984266/140‑1 (MQ=255)
                                                                                     cTAAATTCTCTGGAGTGTCTGGATTATGTGATTATTAATTTAACTCCAGAAAGCAAAAAAAACGTTGATGAGTGAACACCGTAACAATATTCAGGTAGCAATTGATGCTTTATATAGCCAACGTCGGCGTAAATCGCCgg                                                   <  1:1033399/140‑1 (MQ=255)
                                                                                                  aGTGTCTGGATTATGTGATTATTAATTTAACTCCAGAAAGCAAAAAAAACGTTGATGAGTGAACACCGTAACAATATTCAGGTAGCAATTGATGCTTTATATAGCCAACGTCGGCGTAAATCGCCGGGTGATGAAAGTGa                                      <  1:831356/140‑1 (MQ=255)
                                                                                                    tgtCTGGATTATGTGATTATTAATTTAACTCCAGAAAGCAAAAAAAACGTTGATGAGTGAACACCGTAACAATATTCAGGTAGCAATTGATGCTTTATATAGCCAACGTCGGCGTAAATCGCCGGGTGATGAAAGTGaaa                                    >  1:234712/1‑140 (MQ=255)
                                                                                                    tgtCTGGATTATGTGATTATTAATTTAACTCCAGAAAGCAAAAAAAACGTTGATGAGTGAACACCGTAACAATATTCAGGTAGCAATTGATGCTTTATATAGCCAACGTCGGCGTAAATCGCCGAGTGATGAAAGTGaaa                                    >  1:231091/1‑140 (MQ=255)
                                                                                                    tgtCTGGATTATGTGATTATTAATTTAACTCCAGAAAGCAAAAAAAACGTTGATGAGTGAACACCGTAACAATATTCAGGTAGCAATTGATGCTTTATATAGCCAACGTCCGCGCAACTCGCCGGGTGATGAAAGTGaaa                                    >  1:556023/1‑140 (MQ=255)
                                                                                                      tCTGGATTATGTGATTATTAATTTAACTCCAGAAAGCAAAAAAAACGTTGATGAGTGAACACCGTAACAATATTCAGGTAGCAATTGATGCTTTATATAGCCAACGTCGGAGTAAATGGCCGGGGGATGAAAGTGAAaca                                  >  1:459456/1‑140 (MQ=255)
                                                                                                                gtgATTATTAATTTAACTCCAGAAAGCAAAAAAAACGTTGATGAGTGAACACCGTAACAATATTCAGGTAGCAATTGATGCTTTATATAGCCAACGTCGGCGTAAATCGCCGGGTGATGAAAGTGAAACACTTACCCGTa                        >  2:1376947/1‑140 (MQ=255)
                                                                                                                gtgATTATTAATTTAACTCCAGAAAGCAAAAAAAACGTTGATGAGTGAACACCGTAACAATATTCAGGTAGCAATTGATGCTTTATATAGCCAACGTCGGCGTAAATCGCCGGGTGATGAAAGTGAAACACTTACCCGTa                        >  1:1071517/1‑140 (MQ=255)
                                                                                                                gtgATTATTAATTTAACTCCAGAAAGCAAAAAAAACGTTGATGAGTGAACACCGTAACAATATTCAGGTAGCAATTGATGCTTTATATAGCCAACGTCGGCGTAAATCGCCGGGTGATGAAAGTGAAACACTTACCCGTa                        >  2:1541973/1‑140 (MQ=255)
                                                                                                                   attattAATTTAACTCCCGAAAACAAAAAAAGCGTCGATCAGTGATCATCGTAACAATATTCAGGTAGCAATTGATGCTTTATATAGCCAACGTCGGCGTAAATCGCCGGGTGATGAAAGTGAAACACTTACCCGTAGaa                     <  2:556023/140‑1 (MQ=255)
                                                                                                                    ttattaATTTAACTCCAGAAAGCAAAAAAAACGTTGATGAGTGAACACCGTAACAATATTCAGGTAGCAATTGATGCTTTATATAGCCAACGTCGGCGTAAATCGCCGGGTGATGAAAGTGAAACACTTACCCGTAGaaa                    >  1:1260441/1‑140 (MQ=255)
                                                                                                                    ttattaATTTAACTCCAGAAAGCAAAAAAAACGTTGATGAGTGAACACCGTAACAATATTCAGGTAGCAATTGATGCTTTATATAGCCAACGTCGGCGTAAATCGCCGGGTGATGAAAGTGAAACACTTACCCGTAGaaa                    >  1:1007912/1‑140 (MQ=255)
                                                                                                                      attaATTTAACTCCAGAAAGCAAAAAAAACGTTGATGAGTGAACACCGTAACAATATTCAGGTAGCAATTGATGCTTTATATAGCCAACGTCGGCGTAAATCGCCGGGTGATGAAAGTGAAACACTTACCCGTAGAAATg                  <  1:1541973/140‑1 (MQ=255)
                                                                                                                            ttAACTCCAGAAAGCAACAAAATCGTTGATGCGTGCACACCGTAACAACATTCAGGTAGCAATTGATGCTTTATATAGCCAACGTCGGCGTAAATCGCCGGGTGATGAAAGTGAAACACTTACCCGTAGAAATGACGCTa            <  2:459456/140‑1 (MQ=255)
                                                                                                                                cTCCAGAAAGCAAAAAAAACGTTGATGAGTGAACACCGTAACAATATTCAGGTAGCAATTGATGCTTTATATAGCCAACGTCGGCGTAAATCGCCGGGTGATGAAAGTGAAACACTTACCCGTAGAAATGACGCTa            >  1:1429550/1‑136 (MQ=255)
                                                                                                                                cTCCAGAAAGCAAAAAAAACGTTGATGAGTGAACACCGTAACAATATTCAGGTAGCAATTGATGCTTTATATAGCCAACGTCGGCGTAAATCGCCGGGTGATGAAAGTGAAACACTTACCCGTAGAAATGACGCTa            <  2:1429550/136‑1 (MQ=255)
                                                                                                                                 tCCAGAAAGCAAAAAAAACGTTGATGAGTGAACACCGTAACAATATTCAGGTAGCAATTGATGCTTTATATAGCCAACGTCGGCGTAAATCGCCGGGTGATGAAAGTGAACCACTTAGCCGTAGCAATGACGCTATCTTc       >  1:691480/1‑140 (MQ=255)
                                                                                                                                      aaaGCAAAAAAAACGTTGATGACTGAACACCATAACAATATTCAGGTAACAATTGATGCTTTATATAGCCAACGTCGGCGTAAATCGCCGGGTGATGAAAGTGAAACACTTACCCGTAGAAATGACGCTATCTTCGGCaa  <  2:691480/140‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                        
TCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGATTATTAATTTAACTCCAGAAAGCAAAAAAACGTTGATGAGTGAACACCGTAACAATATTCAGGTAGCAATTGATGCTTTATATAGCCAACGTCGGCGTAAATCGCCGGGTGATGAAAGTGAAACACTTACCCGTAGAAATGACGCTATCTTCGGCAA  >  NC_000913/1245267‑1245539

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