Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,020,961 A→G intergenic (‑42/‑177) sulA ← / → sxy SOS cell division inhibitor/CRP‑S‑dependent promoter expression factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,020,9610AG95.2% 129.3 / ‑6.3 43intergenic (‑42/‑177)sulA/sxySOS cell division inhibitor/CRP‑S‑dependent promoter expression factor
Reads supporting (aligned to +/- strand):  ref base A (1/0);  major base G (24/16);  minor base T (2/0);  total (27/16)
Fisher's exact test for biased strand distribution p-value = 5.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACGGATCCGTTAAT  >  NC_000913/1020822‑1021097
                                                                                                                                           |                                                                                                                                        
atgggCCCGGCTGTAGTGGTTTCGGTAGAGAGACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTGCAt                                                                                                                                       <  2:497057/139‑1 (MQ=255)
aTAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGGCGAAGAACGATGTGCATAGCCTGGAGTGTCCATAATAAATACCCCCCTTCGGGGGTACTATATGGATGTACAGTTCAt                                                                                                                                       >  2:117567/1‑143 (MQ=255)
 tAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGTTGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATACTCAATCCAGCCCCTGTGAGTTACTGTATGGACgtacggtgcaga                                                                                                                                      >  1:380432/1‑141 (MQ=255)
   aGCCCGGCTGTTGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTGCATCCa                                                                                                                                    <  2:226382/143‑1 (MQ=255)
     cccGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTGCATCCAGt                                                                                                                                  >  1:395730/1‑143 (MQ=255)
        ggCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTGCATCCAGTGac                                                                                                                               >  2:376728/1‑143 (MQ=255)
          cTGGAGTGTTCTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGGAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTGCATCCAGTGacaa                                                                                                                             <  2:340919/143‑1 (MQ=255)
                 gtTTTCCGTAGAGACACGCGCGATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTGCATCCAGTGACAACAAAGAt                                                                                                                      <  2:491252/143‑1 (MQ=255)
                          agagACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTGCATCCAGTGACAACAAAGATCAACCCtct                                                                                                             >  2:98375/1‑141 (MQ=255)
                            agaCACGCGCAATTTTACTTGCTGTGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGGACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTGCATACAGTGAGAACAAAGATCACCCCTAtgt                                                                                                           >  2:496175/1‑141 (MQ=255)
                            agaCACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTGAATGCAGGGGCGACCAAGATTAACACTAttt                                                                                                           >  2:628537/1‑143 (MQ=255)
                                acGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTGCATCCAGTGACAACAAAGATCAACCCTATTTTCgg                                                                                                       >  1:582483/1‑143 (MQ=255)
                                        ttttACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTGCATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGcc                                                                                               >  1:660947/1‑143 (MQ=255)
                                              ttGCTGCGGATGAGAACGACGCAGAACGATGTGCATAGCCGGAAGTGTACATACTTATTAAAACCCCGGCGAAATGGCAGAGGGATGTACAGTTCATCCAGTGACAGCAAGGATCAAGCCTATTTTCGGAAAGCGCCTCGccc                                                                                         >  2:711189/1‑141 (MQ=255)
                                              ttGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTGCATCCAGTGACAACAAAGATCAACCCTATTTTCGGCAACAGCCCCGCaa                                                                                         >  2:714649/1‑143 (MQ=255)
                                              ttGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTGCATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCaa                                                                                         >  2:590479/1‑143 (MQ=255)
                                                  tgcGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTGCATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAAttt                                                                                     >  2:580775/1‑143 (MQ=255)
                                                   gcGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTTTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGTATGTACAGTGTATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAAtttt                                                                                    <  2:566341/143‑1 (MQ=255)
                                                                  gaagaaCGATGTGCATGTCCGGAAGTGTAAACAAACAAACCAGCCCCTGTGAGTTACTGTATTGATATACAGTACATCCACTGACAACAAAGATCAACCCTACTTTCGCCAAGACCCTCGCAAATCTTATCGTTGTTGACggg                                                                     >  2:522050/1‑143 (MQ=255)
                                                                  gaagaaCGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTGCATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACggg                                                                     >  1:387041/1‑143 (MQ=255)
                                                                        cGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTGCATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACa                                                               <  1:376728/143‑1 (MQ=255)
                                                                           tgtgCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTGCATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATaa                                                            <  1:272417/143‑1 (MQ=255)
                                                                                  aGCCTGAAGTCTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTGCATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATc                                                     <  1:209395/143‑1 (MQ=255)
                                                                                      tGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTGCATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGc                                                 >  1:341683/1‑143 (MQ=255)
                                                                                      tGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTGCATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGc                                                 >  1:68876/1‑143 (MQ=255)
                                                                                      tGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTGCATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGc                                                 <  1:394717/143‑1 (MQ=255)
                                                                                      tGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTGCATCCAGTGACAACAAAGATCAACCCGATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAACACATCAATTAATCTTGc                                                 >  1:367098/1‑143 (MQ=255)
                                                                                            ggaCAAAATGAACCCAAGCCCTGTGAGTTACTGTACGGACGTCCAGTGCATCCAGTGACAACAAAGATCAACCCACAGATCGGATAGAGCCGCGCAAATTTTGGCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTa                                           <  1:103302/141‑1 (MQ=255)
                                                                                              aCATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTGCATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAg                                         >  1:591612/1‑143 (MQ=255)
                                                                                               cATAATCTATCCAGCCCCTGTGAGTTACTGTATGTATGTACAGGGCATCCAGTGACAACAAAGATCAACCCTATTTTCGGAACGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGa                                        <  2:501892/143‑1 (MQ=255)
                                                                                               cATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTGCATCCAGTGACAACAAAGATCAACACTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGa                                        <  2:569196/143‑1 (MQ=255)
                                                                                                   atcaatcCAGCCCCTGTGAGTTACTGTATGGATGTACAGTGCATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTCGTGACCGGAAAACATAAATTAATCTTGCCCCTAAAGAGTaa                                    >  2:205479/1‑143 (MQ=255)
                                                                                                      aatcCAGCCCCTGTGAGTTACTGTATGGATGTACAGTGCATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGtt                                 <  1:259223/143‑1 (MQ=255)
                                                                                                                cTGTGAGTTACTGTATGGATGTACAGTGCATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGGTGCCTATTTTc                       >  1:519308/1‑143 (MQ=255)
                                                                                                                cTGCGAGTTCCTGGATGGATGTACAGTGCATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTc                       <  2:387041/143‑1 (MQ=255)
                                                                                                                 tgtgAGTTACTGTATGGATGTACAGTGCATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGCACCTCGCAAATAATGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGACTAAGTTGCCTATTTTCg                      >  2:440135/1‑143 (MQ=255)
                                                                                                                   tgAGTTACTGTATGGATGTACAGTGCATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTa                    <  1:580775/143‑1 (MQ=255)
                                                                                                                         aCTGTATGGATGTACAGTGCATCCAGTGACACCAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCCAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCATTAAGAATAAGTTGCATATTTTCGTAGTTAAc              >  2:267131/1‑143 (MQ=255)
                                                                                                                         aCTGTATGGATGTACAGTGCATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCATTAAGAATACGTTGCCTATTTTCATTGTTAAc              >  2:578956/1‑143 (MQ=255)
                                                                                                                              aTGGATGTACAGTGCATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACGGATc         >  2:199357/1‑143 (MQ=255)
                                                                                                                               tGGATGTACAGTGCATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAGGTTTTAGTTGCCTATTTTCGTAGTTAACGGATcc        >  1:492959/1‑143 (MQ=255)
                                                                                                                               tGGATGTACAGTGCATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACGGATcc        >  2:261734/1‑143 (MQ=255)
                                                                                                                                ggATGTACAGTGCATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACGGATCCg       >  1:634011/1‑143 (MQ=255)
                                                                                                                                   tGTACAGTGCATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAGCGGATCCGTTa    <  1:371133/143‑1 (MQ=255)
                                                                                                                                     tacagtgcaTCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTACGAATAAGTTGCCTATTTACGTAGTTAACGGGTCCGTTAAt  >  2:98373/1‑143 (MQ=255)
                                                                                                                                     tacagtgcaTCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACGGATCCGTTAAt  >  2:490243/1‑143 (MQ=255)
                                                                                                                                     aacagtgcaTCCAGCGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTTACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACGGATCCGTTAAt  <  1:213957/142‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                        
ATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACGGATCCGTTAAT  >  NC_000913/1020822‑1021097

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