Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,350,507 T→A 12.2% K11N (AAA→AAT) ‡ PP_1176 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,350,5070TA12.2% 109.1 / 8.1 41K11N (AAA→AAT) ‡PP_1176hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (16/20);  new base A (2/3);  total (18/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.44e-01

GAGAGCCGTTTCTGACGAGAACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATGGGAACGGGTTGGTAGTC  >  NC_002947/1350360‑1350654
                                                                                                                                                   |                                                                                                                                                   
gagaGCTGTTTCGGACGAGAACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCAttt                                                                                                                                                   <  2:488786/150‑1 (MQ=255)
     ccGTTTCTGACGAGAACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAg                                                                                                                                              >  1:578753/1‑150 (MQ=255)
     ccGTTTCTGACGAGAACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAg                                                                                                                                              >  1:555970/1‑150 (MQ=255)
      cGTTTCTGACGAGAACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTgagatc                                                                                                                                             >  1:188236/1‑145 (MQ=255)
         ttCTGACGAGAACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCAAATGGCAGCCtt                                                                                                                                          <  2:473626/150‑1 (MQ=255)
           cTGACGAGAACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGACCCACATGTATCTATTGCTAATAATGTTCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCttgt                                                                                                                                        <  1:527952/150‑1 (MQ=255)
           cTGACGAGAACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGACTGCCATTTGGCAGCCttgt                                                                                                                                        <  2:525439/150‑1 (MQ=255)
            tGACGAGAACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCttgtt                                                                                                                                       <  1:298639/150‑1 (MQ=255)
             gACGAGAACGAGATGCCGCTATTGGAGTGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCACCCTGGTCCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCAAATGGCAGCCttttta                                                                                                                                      <  1:692624/150‑2 (MQ=255)
              aCGAGAACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGa                                                                                                                                     >  1:584747/1‑150 (MQ=255)
                      gggATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGc                                                                                                                             >  1:370087/1‑150 (MQ=255)
                       ggATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGcc                                                                                                                            <  2:774224/150‑1 (MQ=255)
                               cTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGtcaat                                                                                                                    <  2:415567/150‑1 (MQ=255)
                               cTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGtcaat                                                                                                                    <  1:358845/150‑1 (MQ=255)
                                        cGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCa                                                                                                                                               <  1:124685/114‑1 (MQ=255)
                                        cGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCa                                                                                                                                               >  2:124685/1‑114 (MQ=255)
                                             ttAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCaa                                                                                                      >  2:424711/1‑150 (MQ=255)
                                             ttAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCaa                                                                                                      >  2:605387/1‑150 (MQ=255)
                                                  cgcCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGTCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGccc                                                                                                 <  1:387782/150‑1 (MQ=255)
                                                  cgcCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTg                                                                                                                                                  <  1:296936/101‑1 (MQ=255)
                                                  cgcCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTg                                                                                                                                                  >  2:296936/1‑101 (MQ=255)
                                                  cgcCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGccc                                                                                                 <  2:578753/150‑1 (MQ=255)
                                                          aaCGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTc                                                                                         <  2:370087/150‑1 (MQ=255)
                                                                           cAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCACTCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACAcg                                                                        <  2:584747/150‑1 (MQ=255)
                                                                            aaaaaCAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCt                                                                                        <  1:673179/133‑1 (MQ=255)
                                                                            aaaaaCAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCt                                                                                        >  2:673179/1‑133 (MQ=255)
                                                                            aaaaaCAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCAAATGGCAGc                                                                                                                                             >  2:395048/1‑80 (MQ=255)
                                                                            aaaaaCAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCAAATGGCAGc                                                                                                                                             <  1:395048/80‑1 (MQ=255)
                                                                             aaaaCAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTAACACAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACAcgcg                                                                      >  1:93552/1‑150 (MQ=255)
                                                                                 cAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGtctc                                                                  <  1:424711/150‑1 (MQ=255)
                                                                                    gCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGATGCCGGATAGCCCGGTCAATCAATCCCAGTTCAACGCCACACCCGTCTGATACTCGATCACACGCGTCTCaaa                                                               >  1:232726/1‑150 (MQ=255)
                                                                                           cGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGttc                                                        <  1:345304/150‑1 (MQ=255)
                                                                                                     aGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCcttct                                              <  1:439308/150‑1 (MQ=255)
                                                                                                              ccTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCa                                     >  1:361550/1‑150 (MQ=255)
                                                                                                                       cttctgcCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCACTCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAACTTCTTCTCCTTCTTCAAGTCCATGATCTCGc                            <  2:525219/150‑1 (MQ=255)
                                                                                                                              cAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCa                     >  2:754585/1‑150 (MQ=255)
                                                                                                                               aGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCAt                    <  1:605387/150‑1 (MQ=255)
                                                                                                                                ggCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATg                   <  1:754585/150‑1 (MQ=255)
                                                                                                                                      aacaaGGCTGCCAAATGGCAGCCTTTTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATGGGAACg             >  2:403410/4‑150 (MQ=255)
                                                                                                                                                ccATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATGGGAACGGGTTGGtggg   <  2:555970/150‑4 (MQ=255)
                                                                                                                                                ccATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATGGGAACGGGTTGGTAGt   >  2:733623/1‑150 (MQ=255)
                                                                                                                                                 cATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATGGGAACGGGTTGGTAGTc  >  1:38543/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GAGAGCCGTTTCTGACGAGAACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATGGGAACGGGTTGGTAGTC  >  NC_002947/1350360‑1350654

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