Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,458,121 C→A 9.7% intergenic (+21/+19) PP_1276 → / ← algA hypothetical protein/mannose‑6‑phosphate isomerase/mannose‑1‑phosphate guanylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,458,1210CA 9.7% 118.1 / 2.8 41intergenic (+21/+19)PP_1276/algAhypothetical protein/mannose‑6‑phosphate isomerase/mannose‑1‑phosphate guanylyltransferase
Reads supporting (aligned to +/- strand):  ref base C (20/17);  new base A (2/2);  total (22/19)
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.76e-01

GAAATGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCA‑GCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCT‑CGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGACTGCACCTCGAT  >  NC_002947/1457975‑1458259
                                                                                                                                                   |                                                                                                                                           
gAAATGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCA‑Gcc                                                                                                                                            >  6:270980/1‑148 (MQ=255)
           cacaAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCA‑GCCCAGGGGCGAGt                                                                                                                                 >  5:129574/1‑148 (MQ=255)
              aaGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCA‑GCCCAGGGGc                                                                                                                                     <  6:243829/141‑1 (MQ=255)
              aaGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCA‑GCCCAGGGGc                                                                                                                                     >  5:243829/1‑141 (MQ=255)
               aGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCA‑GCCCAGGGGCGAGTCAGCa                                                                                                                            <  6:20016/149‑1 (MQ=255)
                  aTCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGGCT‑GACCAGGGGCGa                                                                                                                                   <  6:21094/139‑1 (MQ=255)
                  aTCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGGCT‑GACCAGGGGCGa                                                                                                                                   >  5:21094/1‑139 (MQ=255)
                     aGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTCCGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCA‑GCCCAGGGGCGAGTCAGCATTTACt                                                                                                                      <  3:279335/149‑1 (MQ=255)
                     aGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCA‑GCCCAGGGGCGAGTCAGCATTTACt                                                                                                                      <  5:270980/149‑1 (MQ=255)
                      gggCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCA‑GCCCAGGGGCGAGTct                                                                                                                               <  8:83800/139‑2 (MQ=255)
                      gggCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCA‑GCCCAGGGGCGAGTct                                                                                                                               >  7:83800/1‑138 (MQ=255)
                         cTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCA‑GCCCAGGGGCGAGTCAGCATTTACTgcgc                                                                                                                  <  3:274726/149‑1 (MQ=255)
                         cTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCA‑GCCCAGGGGCGAGTCAGCATTTACTgcgc                                                                                                                  >  2:31748/1‑149 (MQ=255)
                          tGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCA‑GCCCAGGGGCGAGTCAGCATTTACTgcgc                                                                                                                  >  6:9988/1‑148 (MQ=255)
                                cTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCA‑GCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTc                                                                                                           >  6:327566/1‑149 (MQ=255)
                                 ttCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCA‑GCCCAGGGGCGAGTCAGCAt                                                                                                                           <  6:361496/132‑1 (MQ=255)
                                 ttCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCA‑GCCCAGGGGCGAGTCAGCAt                                                                                                                           >  5:361496/1‑132 (MQ=255)
                                 ttCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGGCT‑GACCAGGGGCGAgc                                                                                                                                 <  8:232268/126‑2 (MQ=255)
                                 ttCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGGCT‑GACCAGGGGCGAgc                                                                                                                                 >  7:232268/1‑125 (MQ=255)
                                        ttGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCA‑GCCCAGGGGCGAGTCAGCATTTACTGCGCGa                                                                                                                <  8:206287/136‑1 (MQ=255)
                                        ttGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCA‑GCCCAGGGGCGAGTCAGCATTTACTGCGCGa                                                                                                                >  7:206287/1‑136 (MQ=255)
                                                gaggaTGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCGCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCA‑GCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAAcgccgc                                                                                           <  5:308138/149‑1 (MQ=255)
                                                    atgatTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCA‑GCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTc                                                                                       >  2:91662/1‑149 (MQ=255)
                                                                     ttCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCA‑GCCCAGGGGCGAGTCa                                                                                                                               <  7:331815/92‑1 (MQ=255)
                                                                     ttCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCA‑GCCCAGGGGCGAGTCa                                                                                                                               >  8:331815/1‑92 (MQ=255)
                                                                           gtggggGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCA‑GCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCAt                                                                <  8:131890/149‑1 (MQ=255)
                                                                                     gCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCA‑GCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCT‑c                                                     >  7:20085/1‑149 (MQ=255)
                                                                                      cccATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCA‑GCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCT‑Cg                                                    <  5:9988/149‑1 (MQ=255)
                                                                                         aTTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCA‑GCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCT‑CGAAg                                                 >  6:33251/1‑149 (MQ=255)
                                                                                                     cGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCT‑CGAAGCGCTCGATGt                                       >  1:262911/1‑148 (MQ=255)
                                                                                                           aGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCA‑GCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCT‑CGAAGCGCTCGATGTCATCCTCg                               >  3:151023/1‑149 (MQ=255)
                                                                                                              cccGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCA‑GCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATACT‑CGAAGCGCTCGATGTCATCCTCGccc                            <  8:20085/149‑1 (MQ=255)
                                                                                                                 gAGCGTTGAATAAGCAGGCTCGCCCCTGGTCA‑GCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCT‑CGAAGCGCTCGATGTCATCCTCGCCCAgg                         <  5:327566/149‑1 (MQ=255)
                                                                                                                        gAATAAGCAGGCTCGCCCCTGGTCA‑GCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCT‑CGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTcc                  >  6:95735/1‑149 (MQ=255)
                                                                                                                         aaTAAGCAGGCTCGCCCCTGGTCA‑GCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGCCGTGCGCCCATACACATCCT‑CGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTccc                 >  6:188694/1‑149 (MQ=255)
                                                                                                                         aaTAAGCAGGCTCGCCCCTGGTCA‑GCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCT‑CGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTccc                 >  7:275515/1‑149 (MQ=255)
                                                                                                                          aTAAGCAGGCTCGCCCCTGGTCA‑GCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCT‑CGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCg                <  5:33251/149‑1 (MQ=255)
                                                                                                                              gCAGGCTCGCCCCTGGTCA‑GCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCT‑CGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGCCt            >  3:15547/1‑149 (MQ=255)
                                                                                                                              gCAGGCTCGCCCCTGGTCA‑GCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCT‑CGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGACt            >  5:151903/1‑149 (MQ=255)
                                                                                                                                  gCTCGCCCCTGGTCA‑GCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCT‑CGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGACTGCAc        <  2:262911/149‑1 (MQ=255)
                                                                                                                                         ccTGGTCA‑GCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGTGCCCATACACATCCT‑CGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGACTGCACCTCgat  >  3:291251/1‑148 (MQ=255)
                                                                                                                                         ccTGGTCA‑GCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGACTGCACCTCgat  <  1:91662/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                           
GAAATGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCA‑GCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCT‑CGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGACTGCACCTCGAT  >  NC_002947/1457975‑1458259

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: ‑