Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,728,698 A→G 20.9% intergenic (+50/+73) cspA‑I → / ← PP_1523 major cold shock protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,728,6980AG20.9% 39.5 / 7.1 24intergenic (+50/+73)cspA‑I/PP_1523major cold shock protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (10/9);  new base G (2/3);  total (12/12)
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.21e-01

GCATCGTACCCTGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCC  >  NC_002947/1728558‑1728844
                                                                                                                                            |                                                                                                                                                   
gCATCGTACCCTGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTg                                                                                                                                            >  2:268057/1‑150 (MQ=255)
gCATCGTACCCTGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTg                                                                                                                                            <  2:342616/150‑1 (MQ=255)
  aTCGTACCCTGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGcc                                                                                                                                          >  1:42169/1‑150 (MQ=255)
      tACCCTGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGTG‑‑GCAGCCACCGTGGTGGCTGcttcgcgg                                                                                                                                    <  2:379467/150‑8 (MQ=38)
        cccTGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCGAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGAGAAGCCGCCACCGTGGTGGCTGcttctcgg                                                                                                                                    <  2:30271/150‑8 (MQ=255)
         ccTGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCa                                                                                                                                  >  2:378875/1‑150 (MQ=255)
         ccTGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAACTCCCATGGCTTGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCa                                                                                                                                  <  1:378875/150‑1 (MQ=255)
                              ggAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTaaa                                                                                                             <  1:77913/150‑1 (MQ=255)
                                     gCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTaaa                                                                                                             >  1:462192/1‑143 (MQ=255)
                                     gCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTaaa                                                                                                             <  2:462192/143‑1 (MQ=255)
                                           gTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAACAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACaa                                                                                                >  2:478423/1‑149 (MQ=255)
                                                  gCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGaa                                                                                         >  1:232830/1‑149 (MQ=255)
                                                       aaaGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAGGCAGCCACCGTGGTGGCTGcttcgc                                                                                                                                      >  2:260271/1‑96 (MQ=38)
                                                       aaaGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAGGCAGCCACCGTGGTGGCTGcttcgc                                                                                                                                      <  1:260271/101‑6 (MQ=38)
                                                                    ccGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGCTTcgtc                                                                       <  2:192084/150‑4 (MQ=255)
                                                                                      tctGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGAttt                                                     >  1:35875/1‑150 (MQ=255)
                                                                                         gAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCct                                                  <  2:218055/150‑1 (MQ=255)
                                                                                          aGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCctc                                                 <  1:250078/150‑1 (MQ=255)
                                                                                                   tttCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGACACCACGGTGGATCCAACGGCGCAAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTc                                        >  1:60626/1‑150 (MQ=255)
                                                                                                      cAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCAACACGGTGGCTGCCAGGGGCATAGCCCGTATTCGCGGGTAGGCCGATTCGTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGgc                                     <  1:478423/150‑1 (MQ=255)
                                                                                                            gTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGAt                               >  1:395006/1‑150 (MQ=255)
                                                                                                                gAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCgg                           <  2:42169/150‑1 (MQ=255)
                                                                                                                      gCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCAAGAGGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAg                     >  1:374397/1‑150 (MQ=255)
                                                                                                                                   gcagcCACCGTGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCAtctt        >  2:356452/1‑150 (MQ=255)
                                                                                                                                       ccaccaCGGTGGCTGCCACGGGACTACCTTGTACCGGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCttcttc    <  2:374397/150‑1 (MQ=255)
                                                                                                                                         accacGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTccc  >  2:293222/1‑150 (MQ=255)
                                                                                                                                            |                                                                                                                                                   
GCATCGTACCCTGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGCCACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCC  >  NC_002947/1728558‑1728844

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