Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,313,207 G→A 13.0% intergenic (‑29/+135) 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,2070GA13.0% 109.2 / 7.7 40intergenic (‑29/+135)rbfA/infB30s ribosome binding factor/translation initiation factor IF‑2
Reads supporting (aligned to +/- strand):  ref base G (14/20);  new base A (3/2);  total (17/23)
Fisher's exact test for biased strand distribution p-value = 6.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.43e-01

CCGGAAACGGTGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCTTAAGCAATGGTAC  >  NC_000913/3313060‑3313355
                                                                                                                                                   |                                                                                                                                                    
ccGGAAACGGTGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGt                                                                                                                                                     <  2:1134717/149‑1 (MQ=255)
 cGGAAACGGTGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGtt                                                                                                                                                    <  2:854085/149‑1 (MQ=255)
 cGGAAACGGTGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGtt                                                                                                                                                    <  2:705699/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                                              <  3:60761/132‑1 (MQ=255)
                                                                                                                        ataaATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTaaa                                              >  4:60761/1‑132 (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                                                               <  4:532256/90‑1 (MQ=255)
                                                                                                                                                 cTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCAt                                                               >  3:532256/1‑90 (MQ=255)
                                                                                                                                                   gTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCTTAAGCAATGGTAc  <  4:1011288/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CCGGAAACGGTGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCTTAAGCAATGGTAC  >  NC_000913/3313060‑3313355

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