Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,720,383 (C)5→6 intergenic (+40/+7) slyB → / ← slyA outer membrane lipoprotein SlyB/DNA‑binding transcriptional dual regulator SlyA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,720,3781.C95.7% 68.9 / ‑2.9 23intergenic (+35/+12)slyB/slyAouter membrane lipoprotein SlyB/DNA‑binding transcriptional dual regulator SlyA
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (10/12);  total (10/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.76e-01

AATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAACAGC  >  NC_000913/1720239‑1720512
                                                                                                                                            |                                                                                                                                      
aaTACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACg                                                                                                                                         <  1:540235/140‑1 (MQ=255)
 aTACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGc                                                                                                                                       <  2:680283/140‑1 (MQ=255)
        catGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGcccccctt                                                                                                                                >  2:295017/1‑137 (MQ=255)
           ggtggtACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGcccccct                                                                                                                                 >  1:111936/1‑134 (MQ=255)
           ggtggtACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGcccccct                                                                                                                                 <  2:111936/136‑3 (MQ=255)
                  cAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCCTTCATTTCAccc                                                                                                                       >  1:608439/1‑140 (MQ=255)
                    gAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCCTTCATTTCACCCtt                                                                                                                     <  1:685660/140‑1 (MQ=255)
                     aaaCAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGccc                                                                                                                                     <  2:706910/122‑1 (MQ=255)
                     aaaCAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGccc                                                                                                                                     >  1:706910/1‑122 (MQ=255)
                          aGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCCTTCATTTCACCCTTTg                                                                                                                   <  1:313038/136‑1 (MQ=255)
                          aGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCCTTCATTTCACCCTTTg                                                                                                                   >  2:313038/1‑136 (MQ=255)
                               acacTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACt                                                                                                         <  1:215901/140‑1 (MQ=255)
                                      tttCTCTCCGGGCCAACGTATCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGa                                                                                                   >  1:705212/1‑140 (MQ=255)
                                             ccGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATg                                                                                            <  1:413336/140‑1 (MQ=255)
                                             ccGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATg                                                                                            >  2:289852/1‑140 (MQ=255)
                                                           ggTACTGGCCAGCAATGGGAGTCAGGTGGCCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTTGCCGCACGCCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGa                                                                              <  2:421758/139‑1 (MQ=255)
                                                           cGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGa                                                                              <  2:608439/140‑1 (MQ=255)
                                                                  gCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGt                                                                       >  1:102290/1‑140 (MQ=255)
                                                                         aTGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGt                                                                <  1:294120/140‑1 (MQ=255)
                                                                                  aGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGt                                                       >  1:84713/1‑140 (MQ=255)
                                                                                                ccGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGAt                                         <  2:705212/140‑1 (MQ=255)
                                                                                                 cgcgCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATg                                        <  1:289852/140‑1 (MQ=255)
                                                                                                 cgcgCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATg                                        <  2:84713/140‑1 (MQ=255)
                                                                                                   cgcTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGcc                                      >  2:25976/1‑140 (MQ=255)
                                                                                                    gcTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATTATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCa                                     <  1:295017/140‑1 (MQ=255)
                                                                                                    gcTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCa                                     <  1:220062/140‑1 (MQ=255)
                                                                                                                                aCTGACCACACGCCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGtt                <  2:124676/133‑1 (MQ=255)
                                                                                                                                aCTGACCACACGCCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGtt                >  1:124676/1‑133 (MQ=255)
                                                                                                                                 cTGACCACACGCCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTtaat        <  1:25976/140‑1 (MQ=255)
                                                                                                                                   gACCACACGCCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTtaataa      >  1:451687/1‑140 (MQ=255)
                                                                                                                                       acacGCCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAACAGc  >  1:123922/1‑140 (MQ=255)
                                                                                                                                            |                                                                                                                                      
AATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAACAGC  >  NC_000913/1720239‑1720512

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: