Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,619,266 A→G 7.8% F16S (TTT→TCT)  PP_2292 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,619,2660AG 7.8% 141.7 / 2.3 51F16S (TTT→TCT) PP_2292hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (22/25);  new base G (2/2);  total (24/27)
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.48e-01

TGAGGGGGTGTACACCTCTGCAAGAATTTCGCCGATTCGTTCCTCACTTGCTCCTACTGTGCCGCTCAGTTCTATCGCGATTTTGTCTAGGGTTTTCATGTTGGAGAATGGCGGAATAGCCAATCCTGCAGTTTCAAAGATCCCTTCAAAAGCTTTCAAGTCGTTCAATTCGTTACGCTCCATTTAGTGAATGCATTTCCGTAGGTAGCCTAGTCTAGCTGCTAATAAAATTTTTGCTCGTAAACAACAGGTGCCCCCAGGCGCGTCCTGGGCGGCCTTG  >  NC_002947/2619130‑2619409
                                                                                                                                        |                                                                                                                                               
tgAGGGGGTGTACACCTCTGCAAGAATTTCGCCGATTCGTTCCTCACTTGCTCCTACTGTGCCGCTCAGTTCTATCGCGATTTTGTCTAGGGTTTTCATGTTGGAGAATGGCGGAATAGCCAATCCTGCAGTTTCAAAGATCCCTTCa                                                                                                                                      <  6:273296/148‑1 (MQ=255)
  aGGGGGTGTACACCTCTGCAAGAATTTCGCCGATTCGTTCCTCACTTGCTCCTACTGTGCCGCTCAGTTCTATCGCGATTTTGTCTAGGGTTTTCATGTTGGAGAATGGCGGAATAGCCAATCCTGCAGTTTCAAAGATCCCTTCaa                                                                                                                                     >  2:13410/1‑147 (MQ=255)
    ggggTGTACACCTCTGCAAGAATTTTGCCGATTCGTTCCTCACTTGCTCCTACTGTGCCGCTCAGTTCTATCGCGATTTTGTCTAGGGTTTTCATGTTGGAGAATGGCGGAATAGCCAATCCTGCAGTTTCAAAGATCCCTTCAAAAGc                                                                                                                                 >  8:5316/1‑149 (MQ=255)
      ggTGTACACCTCTGCAAGAATTTCGCCGATTCGTTCCTCACTTGCTCCTACTGTGCCGCTCAGTTCTATCGCGATTTTGTCTAGGGTTTTCATGTTGGAGAATGGCGGAATAGCCAATCCTGCAGTTTCAAAGATCCCTTCAAAAGCtt                                                                                                                               <  7:200004/149‑1 (MQ=255)
            cacCTCTGCAAGAATTTCGCCGATTCGTTCCTCACTTGCTCCTACTGTGCCGCTCAGTTCTATCGCGATTTTGTCTAGGGTTTTCATGTTGGAGAATGGCGGAATAGCCAATCCTGCAGTTTCAAAGATCCCTTCAAAAGCTTTCAAGt                                                                                                                         <  6:306842/149‑1 (MQ=255)
                 ctGCAAGAATTTCGCCGATTCGTTCCTCACTTGCTCCTACTGTGCCGCTCAGTTCTATCGCGATTTTGTCTAGGGTTTTCATGTTGGAGAATGGCGGAATAGCCAATCCTGCAGTTTCAAAGATCCCTTCAAAAGCTTTCAAGTCGTTc                                                                                                                    <  5:341501/149‑1 (MQ=255)
                    cAAGAATTTCGCCGATTCGTTCCTCACTTGCTCCTACTGTGCCGCTCAGTTCTATCGCGATTTTGTCTAGGGTTTTCATGTTGGAGAATGGCGGAATAGCCAATCCTGCAGTTTCAAAGATCCCTTCAAAAGCTTTCAAGTCGTTc                                                                                                                    <  5:336123/146‑1 (MQ=255)
                       gAATTTCGCCGATTCGTTCCTCACTTGCTCCTACTGTGCCGCTCAGTTCTATCGCGATTTTGTCTAGGGTTTTCATGTTGGAGAATGGCGGAATAGCCAATCCTGCAGTTTCAAAGATCCCTTCAAAAGCTTTCAAGTCGTTCAATTCg                                                                                                              >  7:254410/1‑149 (MQ=255)
                                      gTTCCTCACTTGCTCCTACTGTGCCGCTCAGTTCTATCGCGATTTTGACTAGGGTTTTCATGTTGGAGAATGGCGGAATAGGCAATCCTGCAGTTTCAAAGATCCCTTCAAAAGCTTTCAAGGCGGTCAATTCGTTACGGTCCATGTAg                                                                                               >  8:167251/1‑149 (MQ=255)
                                       ttCCTCACTTGCTCCTACTGTGCCGCTCAGTTCTATCGCGATTTTGTCTAGGGTTTTCATGTTGGAGAATGGCGGAATAGCCAATCCTGCAGTTTCAAAGATCCCTTCAAAAGCTTTCAAGTCGTTCAATTCGTTACGCTCCATTTAGt                                                                                              <  8:280429/149‑1 (MQ=255)
                                            cACTTGCTCCTACTGTGCCGCTCAGTTCTATCGCGATTTTGTCTAGGGTTTTCATGTTGGAGAATGGCGGAATAGCCAATCCTGCAGTTTCAGAGATCCCTTCAAAAGCTTTCAAGTCGTTCAATTCGTTACGCTCCATTTAGTGAATg                                                                                         >  5:359472/1‑149 (MQ=255)
                                            cACTTGCTCCTACTGTGCCGCTCAGTTCTATCGCGATTTTGTCTAGGGTTTTCATGTTGGAGAATGGCGGAATAGCCAATCCTGCAGTTTCAGAGATCCCTTCAAAAGCTTTCAAGTCGTTCAATTCGTTACGCTCCATTTAGTGAATg                                                                                         >  5:364400/1‑149 (MQ=255)
                                             aCTTGCTCCTACTGTGCCGCTCAGTTCTATCGCGATTTTGTCTAGGGTTTTCATGTTGGAGAATGGCGGAATAGCCAATCCTGCAGTTTCAAAGATCCCtt                                                                                                                                        >  6:274275/1‑101 (MQ=255)
                                             aCTTGCTCCTACTGTGCCGCTCAGTTCTATCGCGATTTTGTCTAGGGTTTTCATGTTGGAGAATGGCGGAATAGCCAATCCTGCAGTTTCAAAGATCCCtt                                                                                                                                        <  5:274275/101‑1 (MQ=255)
                                             aCTTGCTCCTACTGTGCCGCTCAGTTCTATCGCGATTTTGTCTAGGGTTTTCATGTTGGAGAATGGCGGAATAGCCAATCCTGCAGTTTCAAAGATCCCTTCAAAAGCTTTCAAGTCGt                                                                                                                      >  7:61066/1‑119 (MQ=255)
                                             aCTTGCTCCTACTGTGCCGCTCAGTTCTATCGCGATTTTGTCTAGGGTTTTCATGTTGGAGAATGGCGGAATAGCCAATCCTGCAGTTTCAAAGATCCCTTCAAAAGCTTTCAAGTCGt                                                                                                                      <  8:61066/119‑1 (MQ=255)
                                             aCTTGCTCCTACTGTGCCGCTCAGTTCTATCGCGATTTTGTCTAGGGTTTTCATGTTGGAGAATGGCGGAATAGCCAATCCTGCAGTTTCAAAGATCCCTTCAAAAGCTTTCAAGTCGTTCAATTCGTTACGCTCCATTTAGTGAAt                                                                                          <  3:151096/147‑1 (MQ=255)
                                                        cTGTGCCGCTCAGTTCTATCGCGATTTTGTCTAGGGTTTTCATGTTGGAGAATGGCGGAATAGCCAATCCTGCAGTTTCAAAGATCCCTTCAAAAGCTTTCAAGTCGt                                                                                                                      <  6:88943/108‑1 (MQ=255)
                                                        cTGTGCCGCTCAGTTCTATCGCGATTTTGTCTAGGGTTTTCATGTTGGAGAATGGCGGAATAGCCAATCCTGCAGTTTCAAAGATCCCTTCAAAAGCTTTCAAGTCGt                                                                                                                      >  5:88943/1‑108 (MQ=255)
                                                             ccGCTCAGTTCTATCGCGATTTTGTCTAGGGTTTTCATGTTGGAGAATGGCGGAATAGCCAATCCTGCAGTTTCAAAGATCCCTTCAAAAGCTTTCAAGTCGTTCAATTCGTTACGCTCCATTTAGTGAATGCATTTCCGTAGGTAGcc                                                                        >  5:49137/1‑149 (MQ=255)
                                                               gCTCAGTTCTATCTCGATTTTGTCTAGGGTTTTCATGATGGAGAATGGCGGAATAGCCAATCCTGCAGTTTCAAAGATCCCTTCAAAAGCTTTCAAGTCGTTCAATTCGTTACGCTCCATTTAGTGAATGCATTTCCGTAGGTAGCCTa                                                                      <  4:67801/149‑1 (MQ=255)
                                                                cTCAGTTCTATCGCGATTTTGTCTAGGGTTTTCATGTTGGAGAATGGCGGAATAGCCAATCCTGCAGTTTCAAAGATCCCTTCAAAAGCTTTCAAGTCGTTCAATTCGTTACGCTCCATTTAGTGAATGCATTTCCGTAGGTAGCCTAg                                                                     >  4:278873/1‑149 (MQ=255)
                                                                   aGTTCTATCGCGATTTTGTCTAGGGTTTTCATGTTGGAGAATGGCGGAATAGCCAATCCTGCAGTTTCAAAGATCCCTTCAAAAGCTTTCAAGTCGTTCAATTCGTTACGCTCCATTTAGTGAATGCATTTCCGTAGGTAGCCTAGt                                                                    <  6:49137/147‑1 (MQ=255)
                                                                     ttCTATCGCGATTTTGTCTAGGGTTTTCATGTTGGAGAATGGCGGAATAGCCAATCCTGCAGTTTCAAAGATCCCTTCAAAAGCTTTCAAGTCGTTCAATTCGTTACGCTCCATTTAGTGAATGCATTTCCGTAGGTAGCCTAGTCTAg                                                                >  2:224777/1‑149 (MQ=255)
                                                                        tATCGCGATTTTGTCTAGGGTTTTCATGTTGGAGAATGGCGGAATAGCCAATCCTGCAGTTTCAAAGATCCCTTCAAAAGCTTTCAAGTCGTTCAATTCGTTACGCTCCATTTAGTGAATGCATTTCCGTAGGTAGCCTAGTCTAgctg                                                             >  8:132096/1‑149 (MQ=255)
                                                                           cgcgATTTTGTCTAGGGTTTTCATGTTGGAGAATGGCGGAATAGCCAATCCTGCAGTTTCAAAGATCCCTTCAAAAGCTTTCAAGTCGTTCAATTCGTTACGCTCCATTTAGTGAATGCATTTCCGTAGGTAGCCTAGTCTAgctgct                                                           >  3:272944/1‑148 (MQ=255)
                                                                            gcgATTTTGTCTAGGGTTTTCATGTTGGAGAATGGCGGAATAGCCAATCCTGCAGTTTCAAAGATCCCTTCAAAAGCTTTCAAGTCGTTCAATTCGTTACGCTCCATTTAGTGAATGCATTTCCGTAGGTAGCCTAGTCTAgctgc                                                            <  8:254410/146‑1 (MQ=255)
                                                                            gcgATTTTGTCTAGGGTTTTCATGTTGGAGAATGGCGGAATAGCCAATCCTGCAGTTTCAAAGATCCCTTCAAAAGCTTTCAAGTCGTTCAATTCGTTACGCTCCATTTAGTGAATGCATTTCCGTAGGTAGCCTAGTCTAGCTGCta                                                          <  7:5316/148‑1 (MQ=255)
                                                                            gcgATCTTGTCTAGGGTTTTCATGTTGGAGAATGGCGGAATACCCAATCCTGCAGTTTCAAAGATCCCTTCAAAAGCTTTCAAGTCGTTCAATTCGTTACGCTCCATTTAGTGAATGCATTTCCGTAGGTAGCCTAGTCTAGCTGCta                                                          <  4:42099/148‑1 (MQ=255)
                                                                             ccaTTTTCTCTAGGGTATTCATGTTGGAGAATGGCGGAAAAGCCAATCCTCCACTTTCAAAGATCCCTTCAAAAGCTTTCAAGTCGTACAATTCGATACGCTCCATTTAGTCAATGCATTTCCGTAGGTAGCATAGTCTAGATGCtaa                                                         >  8:277418/3‑148 (MQ=255)
                                                                                         gggTTTTCATGTTGGAGAATGGCGGAATAGGCAATCCTGCAGTTTCAAAGATCCCTTCAAAAGCTTTCAAGTCGTTCAATTCGTTACGCTCCATTTAGTGAATGCATTTCCGTAGGTAGCCTAGTCTAGCTGCTAATAAAATTTTTGc                                             <  7:167251/148‑1 (MQ=255)
                                                                                         gggTTTTCATGTTGGAGAATGGCGGAATAGCCAATCCTGCAGTTTCAAAGATCCCTTCAAAAGCTTTCAAGTCGTTCAATTCGTTACGCTCCATTTAGTGAATGCATTTCCGTAGGTAGCCTAGTCTAGCTGCTAATAAAAttttt                                               <  1:13410/146‑1 (MQ=255)
                                                                                              ttCATGTTGGAGAATGGCGGAATAGCCAATCCTGCAGTTTCAAAGATCCCTTCAAAAGCTTTCAAGTCGTTCAATTCGTTACGCTCCATTTAGTGAATGCATTTCCGTAGGTAGCCTAGTCTAGCTGCTAATAAAATTTTTGCTCGTa                                        <  4:272944/148‑1 (MQ=255)
                                                                                                   gTTGGAGAATGGCGGAATAGCCAATCCTGCAGTTTCAAAGATCCCTTCAAAAGCTTTCAAGTCGTTCAATTCGTTACGCTCCATTTAGTGAATGCATTTCCGTAGGTAGCCTAGTCTAGCTGCTAATAAAATTTTTGCTCGTAaacaa                                   >  6:107123/1‑148 (MQ=255)
                                                                                                    ttGGAGAATGGCGGAATAGCCAATCCTGCAGTTTCAAAGATCCCTTCAAAAGCTTTCAAGTCGTTCAATTCGTTACGCTCCATTTAGTGAATGCATTTCCGTAGGTAGCCTAGTCTAgctgc                                                            <  8:162729/122‑1 (MQ=255)
                                                                                                    ttGGAGAATGGCGGAATAGCCAATCCTGCAGTTTCAAAGATCCCTTCAAAAGCTTTCAAGTCGTTCAATTCGTTACGCTCCATTTAGTGAATGCATTTCCGTAGGTAGCCTAGTCTAgctgc                                                            >  7:162729/1‑122 (MQ=255)
                                                                                                    ttGGAGAATGGCGGAATAGCCAATCCTGCAGTTTCAAAGATCCCTTCAAAAGCTTTCAAGTCGTTCAATTCGTTACGCTCCATTTAGTGAATGCATTTCCGTAGGTAGCCTAGTCTAGCTGCTAATAAAATTTTTGCTCGTAaacaac                                  <  3:278873/148‑1 (MQ=255)
                                                                                                    ttGGAGAATGGCGGAATAGCCAATCCTGCAGTTTCAAAGATCCCTTCAAAAGCTTTCAAGTCGTTCAATTCGTTACGCTCCATTTAGTGAATGCATTTCCGTAGGTAGCCTAGTCTAGCTGCTAATAAAATTTTTGCTCGTAaacaa                                   >  6:97146/1‑147 (MQ=255)
                                                                                                      ggAGAATGGCGGAATAGCCAATCCTGCAGTTTCAAAGATCCCTTCAAAAGCTTTCAAGTCGTTCAATTCGTTACGCTCCATTTAGTGAATGCATTTCCGTAGGTAGCCTAGTCTAGCTGCTAATAAAATTTTTGCTCGTAAACAACAg                                <  1:224777/148‑1 (MQ=255)
                                                                                                       gagaATGGCGGAATAGCCAATCCTGCAGTTTCAGAGATCCCTTCAAAAGCTTTCAAGTCGTTCAATTCGTTACGCTCCATTTAGTGAATGCATTTCCGTAGGTAGCCTAGTCTAGCTGCTAATAAAATTTTTGCTCGTAAACAACAGGt                              <  6:364400/149‑1 (MQ=255)
                                                                                                       gagaATGGCGGAATAGCCAATCCTGCAGTTTCAGAGATCCCTTCAAAAGCTTTCAAGTCGTTCAATTCGTTACGCTCCATTTAGTGAATGCATTTCCGTAGGTAGCCTAGTCTAGCTGCTAATAAAATTTTTGCTCGTAAACAACAGGt                              <  6:359472/149‑1 (MQ=255)
                                                                                                                   aTAGCCAATCCTGCAGTTTCAAAGATCCCTTCAAAAGCTTTCAAGTCGTTCAATTCGTTACGCTCCATTTAGTGAATGCATTTCCGTAGGTAGCCTAGTCTAGCTGCTAATAAAATTTTTGCTCGTAAACAACAGGTGCCCCCAGgcgc                  <  5:97146/149‑1 (MQ=255)
                                                                                                                    tAGCCAATCCTGCAGTTTCAAAGATCCCTTCAAAAGCTTTCAAGTCGTTCAATTCGTTACGCTCCATTTAGTGAATGCATTTCCGTAGGTAGCCTAGTCTAGCTGCTAATAAAATTTTTGCTCGTAAACAACAGGTGCCCCCAGGcgcg                 <  5:107123/149‑1 (MQ=255)
                                                                                                                       ccAATCCTGCAGTTTCAATGATCCCTTCAAAAGCTTTCAAGTCGTTCAATTCGTTACGCTCCATTTAGTGAATGCATTTCCGTAGGTAGCCTAGTCTAGCTGCTAATAAAATTTTTGCTCGTAAACAACAGGTGCCCCCAGGCGCGTcc              >  8:104441/1‑149 (MQ=255)
                                                                                                                       ccAATCCTGCAGTTTCAAAGATCCCTTCAAAAGCTTTCAAGTCGTTCAATTCGTTACGCTCCATTTAGTGAATGCATTTCCGTAGGTAGCCTAGTCTAGCTGCTAATAAAATTTTTGCTCGTAAACAACAGGTGCCCCCAGGCGCGTcc              >  3:240837/1‑149 (MQ=255)
                                                                                                                         aaTCCTGCAGTTTCAAAGATCCCTTCAAAAGCTTTCAAGTCGTTCAATTCGTTACGCTCCATTTAGTGAATGCATTTCCGTAGGTAGCCTAGTCTAgctgc                                                            >  4:260058/1‑101 (MQ=255)
                                                                                                                         aaTCCTGCAGTTTCAAAGATCCCTTCAAAAGCTTTCAAGTCGTTCAATTCGTTACGCTCCATTTAGTGAATGCATTTCCGTAGGTAGCCTAGTCTAGCTGCTAATAAAATTTTTGCTCGTAAACAACAGGTGCCCCCAGGCGCGt                >  2:302948/1‑145 (MQ=255)
                                                                                                                         aaTCCTGCAGTTTCAAAGATCCCTTCAAAAGCTTTCAAGTCGTTCAATTCGTTACGCTCCATTTAGTGAATGCATTTCCGTAGGTAGCCTAGTCTAGCTGCTAATAAAATTTTTGCTCGTAAACAACAGGTGCCCCCAGGCGCGt                <  1:302948/145‑1 (MQ=255)
                                                                                                                          aTCCTGCAGTTTCAAAGATCCCTTCAAAAGCTTTCAAGTCGTTCAATTCGTTACGCTCCATTTAGTGAATGCATTTCCGTAGGTAGCCTAGTCTAGCTGCTAATAAAATTTTTGCTCGTAAACAACAGGTGCCCCCAGGCGCGTCCTgg           <  7:132096/149‑1 (MQ=255)
                                                                                                                           tCCTGCAGTTTCAAAGATCCCTTCAAAAGCTTTCAAGTCGTTCAATTCGTTACGCTCCATTTAGTGAATGCATTTCCGTAGGTAGCCTAGTCTAGCTGCTAATAAAATTTTTGCTCGTAAACAACAGGTGCCCCCAGGCGCGTCCTggg          >  6:266400/1‑149 (MQ=255)
                                                                                                                                   tttCAAAGATCCCTTCAAAAGCTTTCAAGTCGTTCAATTCGTTACGCTCCATTTAGTGAATGCATTTCCGTAGGTAGCCTAGTCTAGCTGCTAATAAAATTTTTGCTCGTAAACAACAGGTGCCCCCAGGCGCGTCCTGGGCGGCCTTg  >  7:143216/1‑149 (MQ=255)
                                                                                                                                        |                                                                                                                                               
TGAGGGGGTGTACACCTCTGCAAGAATTTCGCCGATTCGTTCCTCACTTGCTCCTACTGTGCCGCTCAGTTCTATCGCGATTTTGTCTAGGGTTTTCATGTTGGAGAATGGCGGAATAGCCAATCCTGCAGTTTCAAAGATCCCTTCAAAAGCTTTCAAGTCGTTCAATTCGTTACGCTCCATTTAGTGAATGCATTTCCGTAGGTAGCCTAGTCTAGCTGCTAATAAAATTTTTGCTCGTAAACAACAGGTGCCCCCAGGCGCGTCCTGGGCGGCCTTG  >  NC_002947/2619130‑2619409

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑