Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 950,280 T→G 18.7% intergenic (+133/‑2) lrp → / → ftsK DNA‑binding transcriptional dual regulator, leucine‑binding/DNA‑binding membrane protein required for chromosome resolution and partitioning

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606950,2800TG18.7% 84.7 / 8.3 37intergenic (+133/‑2)lrp/ftsKDNA‑binding transcriptional dual regulator, leucine‑binding/DNA‑binding membrane protein required for chromosome resolution and partitioning
Reads supporting (aligned to +/- strand):  ref base T (15/14);  new base G (4/3);  total (20/17)
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

GACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACACTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGTCGCCTTCTGGAAGCGTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGA  >  REL606/950138‑950424
                                                                                                                                              |                                                                                                                                                
gACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGTCTTTCTTGAGc                                                                                                                                           <  1:301606/150‑1 (MQ=255)
           acGGAACAGGTGCAAAATCTGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATacac                                                                                                                                <  1:100500/150‑1 (MQ=255)
           acGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATacac                                                                                                                                <  2:599900/150‑1 (MQ=255)
              gAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACACTGa                                                                                                                             >  2:559452/1‑150 (MQ=255)
                                        gATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACACTGAAGACAAAGAAGTCACATTGACAAAGt                                                                                                   >  2:161031/1‑150 (MQ=255)
                                             cacTCCTGTTAATCCATACAGCAACAGTACTGGGGTAAACTGGTACTGTTGTCCGTTTTGGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACACTGAAGACAAAGAAGTCACATTGACAAAGTTAAGt                                                                                              <  2:670815/150‑1 (MQ=255)
                                               cTCCGGTTAATCCATACAGCAACAGTACTGGGGTCACCTGGTACTGTTGTCCGGTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCGTGCTTCAGCTAGGCATACACTGAAGACAAAGAAGTCGCATTGACAGAGTTCAGTac                                                                                            >  1:324513/1‑149 (MQ=255)
                                                 ccTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTCGTCCGTTGTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTCCTTGAGCCAGGAATACACTGAAGACAAAGAAGTCACATTGCCAAAGTTAAGTAGCg                                                                                          >  1:282804/1‑150 (MQ=255)
                                                 ccTGTTAATCCATACAGCAACACTACTGGGGTAACCTGGTACTGTTGTCCGTCTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGGGCCTTGCTTTAGCCAGGAATACACTTAAGACACAGAATTCACGTCGAGAAAGttctagtggc                                                                                          >  1:523047/1‑142 (MQ=255)
                                                  cTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACACTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCgg                                                                                         >  2:634381/1‑150 (MQ=255)
                                                  cTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCATTTCTTGAGCCAGGAATACACTGAAGACAAAGAAGTCACATTGACACAGTTAAGTAGCgg                                                                                         >  1:403334/1‑150 (MQ=255)
                                                                  cAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACACTGAAGACAAAGAAGTCACATTCACCAAGTTAAGTAGCGGCCGTCGCCTTCTGaaa                                                                         >  1:176582/1‑150 (MQ=255)
                                                                   aaCAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACACTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGTCGCCTTCTGGAAg                                                                        <  2:569730/150‑1 (MQ=255)
                                                                   aaCAGTACTGGAGTAACATGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATACACTGAAGAAAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGTCGCCTTCTGGAAg                                                                        <  1:132443/150‑1 (MQ=255)
                                                                     gtataagagacagtccTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGTAGAGCCTTTCTTGAGCCAGGAATACACTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGTCGCCTTCTGGAAGCg                                                                      <  2:592489/136‑1 (MQ=255)
                                                                     cAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACACTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGTCGCCTTCTGGAAGCg                                                                      <  2:206358/150‑1 (MQ=255)
                                                                       gTACTGGGGTCACCTGGTACTGTGGTCCGTTTTAGCATCGGGCCGGAAAAGCCTGTAACCGGGAGAGCCTTTCTTGCGCTAGGAATACACTGCAGACAAAGAAGTCACATTGACAATGGTAGGTCGCGGCCGTTGCTTTCTGGAAGCGtt                                                                    >  1:609001/1‑150 (MQ=255)
                                                                       gTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACACTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGTCGCCTTCTGTAAGCGtt                                                                    >  1:627784/1‑150 (MQ=255)
                                                                       gTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACACTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGTCGCCTTCTGGAAGCGtt                                                                    >  2:57761/1‑150 (MQ=255)
                                                                       gTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACACTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGTCGCCTTCTGGAAGCGtt                                                                    >  1:55682/1‑150 (MQ=255)
                                                                       gTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACACTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGTCGCCTTCTGGAAGCGtt                                                                    >  1:330037/1‑150 (MQ=255)
                                                                       gTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACACTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGTCGCCTTCTGGAAGCGtt                                                                    >  2:109662/1‑150 (MQ=255)
                                                                       gTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACACTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGTCGCCTTCTGGAAGCGtt                                                                    >  2:347741/1‑150 (MQ=255)
                                                                       gTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACATGGAGAGCCTTGCTTGAGCCAGGAATACACTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGTCGCCTTCTGGAAGCGtt                                                                    >  1:479912/1‑150 (MQ=255)
                                                                         aCTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATACACTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGTCGCCTTCTGGAAGCGTTGc                                                                  <  2:605651/150‑1 (MQ=255)
                                                                                 aaCCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACACTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGTCGCCTTCTGGAAGCGTTGCTGATCCtt                                                          <  2:330973/150‑1 (MQ=255)
                                                                                  aCCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACACTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGTCGCCTTCTGGAAGCGTTGCTGATCCTTa                                                         >  2:682128/1‑150 (MQ=255)
                                                                                  aCCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACACTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGTCGCCTTCTGGAAGCGCTgactcttatac                                                         >  1:590751/1‑140 (MQ=255)
                                                                                     tGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACACTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGTCGCCTTCTGGAAGCGTTGCTGATCCTTATTg                                                      <  1:148719/150‑1 (MQ=255)
                                                                                     tGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACACTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGTCGCCTTCTGGAAGCGTTGCTGATCCTTATTg                                                      <  2:335790/150‑1 (MQ=255)
                                                                                     tGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACACTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGTCGCCTTCTGGAAGCGTTGCTGATCCTTATTg                                                      <  1:57592/150‑1 (MQ=255)
                                                                                     tGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGGAACCTGGAGAGCCTTTCTTGAGCCAGGAATACACTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGTCGCCTTCTGGAAGCGTTGCTGATCCTTATTg                                                      <  1:632523/150‑1 (MQ=255)
                                                                                                             cGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACACTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGTCGCCTTCTGGAAGCGTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATgg                              <  1:680130/150‑1 (MQ=255)
                                                                                                             cGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACACTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGTCGCCTTCTGGAAGCGTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATgg                              <  1:160582/150‑1 (MQ=255)
                                                                                                                         gCCTGTAACCTGGAGAGCCTTGCTTGAGCCAGGAATACACTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGTCGCCTTCTGGAAGCGTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTaa                  >  1:371415/1‑150 (MQ=255)
                                                                                                                                  cTGGAGAGCCTTTCTTGAGCCAGGAATACACTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGTCGCCTTCTGGAAGCGTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAAcc         <  2:629633/150‑1 (MQ=255)
                                                                                                                                         gCCTTTCTTGAGCCAGGAATACACTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGTCGCCTTCTGGAAGCGTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGa  >  1:85511/1‑150 (MQ=255)
                                                                                                                                              |                                                                                                                                                
GACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACACTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGTCGCCTTCTGGAAGCGTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGA  >  REL606/950138‑950424

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