Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,223,984 (A)7→6 intergenic (+77/+295) ycgI → / ← minE protein YcgI/Z‑ring positioning protein MinE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,223,9780A.97.9% 189.5 / ‑2.1 48intergenic (+71/+301)ycgI/minEprotein YcgI/Z‑ring positioning protein MinE
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base . (19/28);  total (19/29)
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

GTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCAGGCATAAACAT  >  NC_000913/1223841‑1224109
                                                                                                                                         |                                                                                                                                    
gtATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCaaa                                                                                                                                    >  1:998666/1‑140 (MQ=255)
 tataGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCaaaa                                                                                                                                   <  2:998666/140‑1 (MQ=255)
    aGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCaaaaaag                                                                                                                                <  2:893647/140‑2 (MQ=255)
     gCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCaaaaaagc                                                                                                                               >  1:620231/1‑138 (MQ=255)
        aTGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAGCCCg                                                                                                                           >  1:59477/1‑140 (MQ=255)
          ggTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAGCCCGcc                                                                                                                         >  1:1046727/1‑140 (MQ=255)
          ggTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAGCCCGcc                                                                                                                         >  2:930463/1‑140 (MQ=255)
           gTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAGCCCGCCa                                                                                                                        >  2:836447/1‑140 (MQ=255)
           gTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAGCCCGCCa                                                                                                                        >  1:73179/1‑140 (MQ=255)
             gCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAGCCCGCCAGt                                                                                                                      >  2:283011/1‑140 (MQ=255)
              ccGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAGCCCGCCAGtt                                                                                                                     >  2:644136/1‑140 (MQ=255)
              ccGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAGCCCGCCAGtt                                                                                                                     >  2:1002186/1‑140 (MQ=255)
                ggTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAGCCCGCCAGTTAc                                                                                                                   >  2:878289/1‑140 (MQ=255)
                      gAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAGCCCGCCAGTTACGGCggg                                                                                                             <  1:878289/140‑1 (MQ=255)
                       aaTCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCGAAAAAGCCCGCCAGTTACGGCGGGa                                                                                                            <  1:149980/140‑1 (MQ=255)
                         tCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAGCCCGCCAGTTACGGCGGGaaa                                                                                                          <  2:833191/140‑1 (MQ=255)
                            ccGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTGGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAGCCCGCCAGTTACGGCGGGAAACCt                                                                                                       <  1:283011/140‑1 (MQ=255)
                            ccGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAGCCCGCCAGTTa                                                                                                                    <  1:39443/127‑1 (MQ=255)
                            ccGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAGCCCGCCAGTTa                                                                                                                    >  2:39443/1‑127 (MQ=255)
                            ccGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAGCCCGCCAGTTACGGCGGGAAACCt                                                                                                       <  1:1002186/140‑1 (MQ=255)
                            ccGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAGCCCGCCAGTTACGGCGGGAAACCt                                                                                                       <  2:73179/140‑1 (MQ=255)
                            ccGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAGCCCGCCAGTTACGGCGGGAAACCt                                                                                                       <  2:478333/140‑1 (MQ=255)
                              gTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCa                                                                                                     <  1:930463/140‑1 (MQ=255)
                                ggAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATc                                                                                                   <  1:830701/140‑1 (MQ=255)
                                ggAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATc                                                                                                   <  2:620231/140‑1 (MQ=255)
                                ggAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATc                                                                                                   <  2:59477/140‑1 (MQ=255)
                                      cgGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAGCCCGCCAGTTACGGCGGGAAAAcccatcccattg                                                                                             >  2:965628/1‑138 (MQ=39)
                                        gTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGga                                                                                           <  2:487455/140‑1 (MQ=255)
                                        gTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGga                                                                                           <  2:105389/140‑1 (MQ=255)
                                        gTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTAGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGga                                                                                           <  2:654355/140‑1 (MQ=255)
                                        agacagggTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTGCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGga                                                                                           <  1:627773/134‑1 (MQ=255)
                                                  tgAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAGCAATGaa                                                                                <  2:1046727/140‑1 (MQ=255)
                                                         gTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAGCCCGCCAGTTACGGCGGGAAAACCCATCCTATTGGGGAGGAAAGAATAAAGaa                                                                         >  2:356272/1‑140 (MQ=255)
                                                                 gACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAGCAATGAATAATGAAATTGCggg                                                                 <  2:675158/140‑1 (MQ=255)
                                                                      tCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAGCAATGAATAATGAAATTGCGGGGTTat                                                            <  1:59341/140‑1 (MQ=255)
                                                                      tCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAGCCCGCCAGTTACGGCGGGAAAACTCATCCTATTGGGGAGCAAAGAAAAAAGAAAATTCGGGGGTat                                                            >  2:771749/1‑140 (MQ=255)
                                                                        aaCGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAGCCCGCCAGTTACGGCGGGAAAACCCATCCTATGGGAGAGCAAAGAATAATGAAAATGCGGGGGTatca                                                          >  2:450148/1‑140 (MQ=255)
                                                                          cGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAGCAATGAATAATGAAATTGCGGGGTTatcatc                                                        <  1:722172/140‑1 (MQ=255)
                                                                           gATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAGCAATGAATAATGAAATTGCGGGGTTATCAtct                                                       <  2:322716/140‑1 (MQ=255)
                                                                            aTTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAGCAATGAATAATGAAATTGCGGGGTTATCAtctc                                                      <  2:177543/140‑1 (MQ=255)
                                                                            aTTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAGCAATGAATAATGAAATTGCGGGGTTATCAtctc                                                      <  1:274004/140‑1 (MQ=255)
                                                                                acaggcgcTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAGCAATGAATAATGAAATTGCGGGGTTATCATCTCCCAg                                                  <  1:267800/136‑1 (MQ=255)
                                                                                aaCTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAGCAATGAATAATGAAATTGCGGGGTTATCATCTCCCAg                                                  <  1:1095941/140‑1 (MQ=255)
                                                                                aaCTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAGCAATGAATAATGAAATTGCGGGGTTATCATCTCCCAg                                                  <  2:654639/140‑1 (MQ=255)
                                                                                 aCTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAGCAATGAATAATGAAATTGCGGGGTTATCACCTCCCAGt                                                 <  1:1114010/140‑1 (MQ=255)
                                                                                    gcgcTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAGCAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGtata                                              <  1:751325/140‑1 (MQ=255)
                                                                                               cAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAGCCCGCCAGTTACGGCGGGAAAACTCCTCCTATTGGGGAGGAAAGAAAAATGAAAATGCGGGGGTATTATCTCCCAGTATATTCATAcccgt                                   >  2:267800/1‑136 (MQ=255)
                                                                                                                   ttGAAAGGGAAAATTCTGGGGCAAAAAAGCCCGCCAGTTACGGCGGGGAAACTCCTCCTATGGGAGGGGAATGAATAATGAAATTGCGGGGGTATTATCTCCCCGTATATCCATACTAAAAAAAAAGGTATTTTTTCAAc               >  2:788599/1‑140 (MQ=255)
                                                                                                                   ttGAAAGGGAAAATTCTGGGGCAAAAAAGCCCGCCAGTTACGGCGGGAAACCCCCTCCCATGGGGGAGCAAAGAAAAAAGAAAATGCGGGGGTTTCATCTCCCCGCATATCCCTACTAACAATAAAGTTATTTTATCAAc               >  2:716755/1‑140 (MQ=255)
                                                                                                                           gAAAATTCTGGGGCAAAAAAGCCCGCCAGTTACGGCGGGAAAACTCATCCTTTTGGGGAGGAAAGAAAAATAAAAATTCGGGGGTATCATCTCCCCCTATATTCATACTAAAAAAAAAGATATTTATTTAAACAAGCATa       >  2:977967/1‑140 (MQ=255)
                                                                                                                           gAAAATTCTGGGGCAAAAAAGCCCGCCAGTTACGGCGGGAAAAAACCTTCTATTGGGGAGGAAAGGAAAAAAAAAATTGGGGGGTATGATCTCCCCATATATTCCAACTTAAAAAAAAATTTTTTTTTCAACCAGGCCTa       >  2:96163/1‑140 (MQ=255)
                                                                                                                                ttCTGGGGCAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAGCAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCAGGCATAAACAt  <  1:931023/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                    
GTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCAGGCATAAACAT  >  NC_000913/1223841‑1224109

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