Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,313,208 T→A 15.8% intergenic (‑30/+134) rbfA ← / ← infB 30s ribosome binding factor/translation initiation factor IF‑2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,313,2080TA15.8% 105.9 / 13.4 38intergenic (‑30/+134)rbfA/infB30s ribosome binding factor/translation initiation factor IF‑2
Reads supporting (aligned to +/- strand):  ref base T (14/18);  new base A (3/3);  total (17/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

CCGGAAACGGTGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCTTAAGCAATGGTACG  >  NC_000913/3313060‑3313356
                                                                                                                                                    |                                                                                                                                                    
ccGGAAACGGTGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGt                                                                                                                                                      <  2:1134717/149‑1 (MQ=255)
 cGGAAACGGTGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGtt                                                                                                                                                     <  2:705699/149‑1 (MQ=255)
 cGGAAACGGTGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGtt                                                                                                                                                     <  2:854085/149‑1 (MQ=255)
     aaCGGTGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGcc                                                                                                                                                 <  2:1079640/149‑1 (MQ=255)
          tggtCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCttt                                                                                                                                            <  1:1145041/149‑1 (MQ=255)
          tggtCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGACACCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGAACAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTCACAGCCCCttt                                                                                                                                            <  4:105018/149‑1 (MQ=255)
               atcatGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTaaa                                                                                                                                       >  2:176094/1‑149 (MQ=255)
                tcatGCCCAGGCGAGGATCTTTAATTTCACGCTGTAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTaaaa                                                                                                                                      <  1:939669/149‑1 (MQ=255)
                  atGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGCCCAAATTCTTTAGCCATAATAGATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAaatt                                                                                                                                    <  3:886818/149‑1 (MQ=255)
                     cccAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAaattaat                                                                                                                                 <  4:33851/149‑1 (MQ=255)
                          gCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAg                                                                                                                            <  1:395911/149‑1 (MQ=255)
                           cGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAgg                                                                                                                           <  3:50159/149‑1 (MQ=255)
                                   tttAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAggg                                                                                                                   >  3:761170/1‑149 (MQ=255)
                                      aaTTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTg                                                                                                                <  4:761170/149‑1 (MQ=255)
                                           cACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCa                                                                                                                             <  1:436859/131‑1 (MQ=255)
                                           cACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCa                                                                                                                             >  2:436859/1‑131 (MQ=255)
                                                               aTCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACgcgc                                                                                       >  3:790034/1‑149 (MQ=255)
                                                                           aTTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATc                                                                           <  1:176094/149‑1 (MQ=255)
                                                                            tttCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCa                                                                                                                             <  3:933846/98‑1 (MQ=255)
                                                                            tttCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCa                                                                                                                             >  4:933846/1‑98 (MQ=255)
                                                                                             tGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGAt                                                         >  4:46790/1‑149 (MQ=255)
                                                                                                    gACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAaattaat                                                                                                                                 <  1:46748/70‑1 (MQ=255)
                                                                                                    gACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAaattaat                                                                                                                                 >  2:46748/1‑70 (MQ=255)
                                                                                                        aaaTTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAAc                                              >  3:565837/1‑149 (MQ=255)
                                                                                                              tttCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTAACAGCCCCtttt                                                                                                                                           >  2:793606/1‑50 (MQ=25)
                                                                                                              tttCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTAACAGCCCCtttt                                                                                                                                           <  1:793606/50‑1 (MQ=21)
                                                                                                                        ataaATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTaaa                                               >  4:60761/1‑132 (MQ=255)
                                                                                                                        ataaATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTaaa                                               <  3:60761/132‑1 (MQ=255)
                                                                                                                             ttCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAAccc                         >  3:1011288/1‑149 (MQ=255)
                                                                                                                                    gACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAAc                  >  1:514594/1‑149 (MQ=255)
                                                                                                                                     tttaaaaaGGGGCTAACAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCgg                                                           >  2:88955/4‑107 (MQ=38)
                                                                                                                                     tttaaaaaGGGGCTAACAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCgg                                                           <  1:88955/104‑1 (MQ=39)
                                                                                                                                        aaaaaGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCTTa              <  4:565837/149‑1 (MQ=255)
                                                                                                                                          aaaGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCTTAAg            >  1:320688/1‑149 (MQ=255)
                                                                                                                                          aaaGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGGTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCTTAAg            >  2:985620/1‑149 (MQ=255)
                                                                                                                                          aaaGGGGCTAACAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGc                                            >  2:483835/1‑117 (MQ=255)
                                                                                                                                               ggCTGTTAGCCCCTTTTTAAAAGTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCCTAAGCAATg       <  1:985620/149‑1 (MQ=255)
                                                                                                                                                gCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCTTAAGCAATgg      <  3:46790/149‑1 (MQ=255)
                                                                                                                                                 cTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCAt                                                                >  3:532256/1‑90 (MQ=255)
                                                                                                                                                 cTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCAt                                                                <  4:532256/90‑1 (MQ=255)
                                                                                                                                                   gTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCTTAAGCAATGGTAc   <  4:1011288/149‑1 (MQ=255)
                                                                                                                                                    ttAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCTTaa             >  1:504822/1‑138 (MQ=255)
                                                                                                                                                    ttAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCTTaa             <  2:504822/138‑1 (MQ=255)
                                                                                                                                                    ttAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCTTAAGCAATGGTACg  >  4:1036730/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCGGAAACGGTGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCTTAAGCAATGGTACG  >  NC_000913/3313060‑3313356

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