Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,685,390 T→C 8.4% intergenic (+29/‑30) PP_5593 → / → PP_5594 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,685,3900TC 8.4% 168.1 / 4.7 59intergenic (+29/‑30)PP_5593/PP_5594hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (23/31);  new base C (2/3);  total (25/34)
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.97e-01

CCCGGGCAATTTCGCCAATGACCGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGAGCTTGGTCGGCGGCTGTTTCGACCGCGACAATGACCACCCTGGTAAAGATTCCGACCC  >  NC_002947/3685241‑3685517
                                                                                                                                                     |                                                                                                                               
cccGGGCAATTTCGCCAATGACCGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCt                                                                                                                                 >  2:1794331/1‑150 (MQ=255)
cccGGGCAATTTCGCCAATGACCGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCt                                                                                                                                 <  1:1836110/150‑1 (MQ=255)
 ccGGGCAATTTCGCCAATGACCGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTg                                                                                                                                <  2:1688412/150‑1 (MQ=255)
 ccGGGCAATTTCGCCAATGACCGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTg                                                                                                                                >  1:1457553/1‑150 (MQ=255)
 ccGGGCAATTTCGCCAATGACCGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTg                                                                                                                                <  1:659947/150‑1 (MQ=255)
  cGGGCAATTTCGCCAATGACCGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTgg                                                                                                                               <  1:1737561/150‑1 (MQ=255)
       aaTTTCGCCAATGACCGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCAGGTCGCCCCGGCACTg                                                                                                                          <  2:44124/150‑1 (MQ=255)
               cAATGACCGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATc                                                                                                                  >  2:725702/1‑150 (MQ=255)
               cAATGACCGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCAGGTCGCCCCGGCACTGcgc                                                                                                                       >  1:344882/1‑143 (MQ=255)
               cAATGACCGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCAGGTCGCCCCGGCACTGcgc                                                                                                                       <  2:344882/145‑3 (MQ=255)
                   gACCGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATcc                                                                                                                 <  2:444025/147‑1 (MQ=255)
                   gACCGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATcc                                                                                                                 >  1:444025/1‑147 (MQ=255)
                    aCCGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACTGAGGCGCAGTGCCGGGGCGACATGGCACTGCAATTATCCAAgg                                                                                                             <  2:275201/150‑1 (MQ=255)
                      cGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAgggg                                                                                                           >  1:1520962/1‑150 (MQ=255)
                      cGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAgggg                                                                                                           <  2:1925457/150‑1 (MQ=255)
                        tGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTgg                                                                                                                               <  1:1633268/128‑1 (MQ=255)
                        tGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTgg                                                                                                                               >  2:1633268/1‑128 (MQ=255)
                                   ccGACGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCAACAAGGACGGCCGTTCGTAACTGACCGAGTCTTAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATgcgc                                                                                              <  1:725702/150‑1 (MQ=255)
                                            gTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCaa                                                                                                                       >  2:886213/1‑116 (MQ=255)
                                            gTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCaa                                                                                                                       <  1:886213/116‑1 (MQ=255)
                                             tCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCAc                                                                                    <  1:1564557/150‑1 (MQ=255)
                                              cGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACt                                                                                   <  2:104050/150‑1 (MQ=255)
                                                   aGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGa                                                                                 >  1:247053/1‑147 (MQ=255)
                                                   aGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGa                                                                                 <  2:247053/147‑1 (MQ=255)
                                                    ggcggcCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGatcat                                                                             <  1:431262/150‑1 (MQ=255)
                                                     gcggcCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGatcatc                                                                            <  2:1457553/150‑1 (MQ=255)
                                                     gcggcCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGatcatc                                                                            <  1:357792/150‑1 (MQ=255)
                                                         ccAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGatcatc                                                                            >  2:254098/1‑146 (MQ=255)
                                                         ccAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGatcatc                                                                            <  1:254098/146‑1 (MQ=255)
                                                         ccAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCAGGTCGCCCCGGCACTGcg                                                                                                                        >  1:709126/1‑101 (MQ=38)
                                                         ccAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCAGGTCGCCCCGGCACTGcg                                                                                                                        <  2:709126/102‑2 (MQ=38)
                                                                      ggcAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTgcg                                                           >  2:1698902/1‑148 (MQ=255)
                                                                          aCATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGAGCTTg                                                       <  1:1794010/150‑1 (MQ=255)
                                                                             tGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCt                                                                                            >  1:1095197/1‑110 (MQ=255)
                                                                             tGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCt                                                                                            <  2:1095197/110‑1 (MQ=255)
                                                                                  cACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGAGCTTGGTcggcgg                                               >  2:1584896/1‑150 (MQ=255)
                                                                                        cgcgAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATGATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGAGCTTGGTCGGCGGCTGttt                                         >  1:1105209/1‑150 (MQ=255)
                                                                                          cgAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGAGCTTGGTcggcg                                                <  2:1465636/141‑1 (MQ=255)
                                                                                          cgAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGAGCTTGGTcggcg                                                >  1:1465636/1‑141 (MQ=255)
                                                                                              tCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGAGCTTGGTCGGCGGCTGTTTCGACcg                                   >  2:862278/1‑150 (MQ=255)
                                                                                                gggCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGa                                                            <  1:224630/123‑1 (MQ=255)
                                                                                                gggCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGa                                                            >  2:224630/1‑123 (MQ=255)
                                                                                                 ggCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGAGCTTGGTCGGCGGCTGTTTCGACCGCGa                                >  2:404611/1‑150 (MQ=255)
                                                                                                    cacaAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGatcatc                                                                            >  1:1345835/1‑103 (MQ=255)
                                                                                                    cacaAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGatcatc                                                                            <  2:1345835/103‑1 (MQ=255)
                                                                                                     acaAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTgg                                                                                                                               >  1:531596/1‑51 (MQ=255)
                                                                                                     acaAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTgg                                                                                                                               <  2:531596/51‑1 (MQ=255)
                                                                                                     acaAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCCTCCTGATGGGTTCGCTGAGCTTGGTCGGCGGCTGTTTCGACCGCGACAAt                            <  2:1520962/150‑1 (MQ=255)
                                                                                                     acaAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGAGCTTGGTCGGCGGCTGTTTCGACCGCGACAAt                            <  2:1426930/150‑1 (MQ=255)
                                                                                                       aaGGGCGGCCGTTCGTAACTGACCGAGGAGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGAGCTTGGTCGGCGGCTGTTTCGACCGCGACAATGa                          >  2:1677846/1‑150 (MQ=255)
                                                                                                        aGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGAGCTTGGTCGGCGGCTGTTTCGACCGCGACAATGac                         <  1:1794331/150‑1 (MQ=255)
                                                                                                                  ttCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGAGCTTGGTCGGCGGCTGTTTCg                                       >  2:1273805/1‑126 (MQ=255)
                                                                                                                  ttCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGAGCTTGGTCGGCGGCTGTTTCg                                       <  1:1273805/126‑1 (MQ=255)
                                                                                                                   tCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGAGCTTGGTCGGCGGCTGTTTCGACCGCGACAATGACCACCCTGGTaa              >  1:1931327/1‑150 (MQ=255)
                                                                                                                         cTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGAGCTTGGTCGGCGGCTGTTTCGACCGCGACAATGACCACCCTGGTAAAGATTc        <  1:712385/150‑1 (MQ=255)
                                                                                                                             ccGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGAGCTTGGTCGGCGGCTGTTTCGACCGCGACAATGACCACCCTGGTAAAGATTCCGAc    <  1:862278/150‑1 (MQ=255)
                                                                                                                               gAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGAGCTTGGTCGGCGGCTGTTTCGACCGCGACAATGACCACCCTGGTAAAGATTCCGAccc  <  2:1931327/150‑1 (MQ=255)
                                                                                                                                             gggCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGAg                                                           <  1:892836/79‑1 (MQ=255)
                                                                                                                                             gggCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGAg                                                           <  1:893037/79‑1 (MQ=255)
                                                                                                                                             gggCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGAg                                                           >  2:892836/1‑79 (MQ=255)
                                                                                                                                             gggCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGAg                                                           >  2:893037/1‑79 (MQ=255)
                                                                                                                                                     |                                                                                                                               
CCCGGGCAATTTCGCCAATGACCGTGAAAAAGCCTCCGAAGCCGGTCGCAAAGGCGGCCAGTCCTCGGGCGGCAACATGCCCCACGACCGCGAGTCGGGCCACAAGGGCGGCCGTTCGTAACTGACCGAGGCGCAGTGCCGGGGCGACCTGGCACTGCAATTATCCAAGGGGGTACTCCATGCGCCTGACTCCACTGATCATCCTGATGGGTTCGCTGAGCTTGGTCGGCGGCTGTTTCGACCGCGACAATGACCACCCTGGTAAAGATTCCGACCC  >  NC_002947/3685241‑3685517

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