Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,458,119 A→T 9.3% intergenic (+19/+21) 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,1190AT 9.3% 117.2 / 3.3 43intergenic (+19/+21)PP_1276/algAhypothetical protein/mannose‑6‑phosphate isomerase/mannose‑1‑phosphate guanylyltransferase
Reads supporting (aligned to +/- strand):  ref base A (21/18);  new base T (2/2);  total (23/20)
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.58e-01

CTCGAAATGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCA‑GCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCT‑CGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGACTGCACCTCGAT  >  NC_002947/1457972‑1458259
                                                                                                                                                   |                                                                                                                                              
cTCGAAATGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCA‑g                                                                                                                                              <  7:94308/149‑1 (MQ=255)
   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                                                                                                                           >  5:361496/1‑132 (MQ=255)
                                    ttCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCA‑GCCCAGGGGCGAGTCAGCAt                                                                                                                           <  6:361496/132‑1 (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                                                                                                                               >  8:331815/1‑92 (MQ=255)
                                                                        ttCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCA‑GCCCAGGGGCGAGTCa                                                                                                                               <  7:331815/92‑1 (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)
                                                                                                                                                   |                                                                                                                                              
CTCGAAATGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCA‑GCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCT‑CGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGACTGCACCTCGAT  >  NC_002947/1457972‑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: ‑