Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,668,555 T→C 7.3% A44A (GCA→GCG)  PP_3231 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,668,5550TC 7.3% 156.6 / 2.2 56A44A (GCA→GCG) PP_3231hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (23/28);  new base C (2/2);  total (26/30)
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.75e-01

GTCGGCCTCGACGGCGACGATGTTGAAGCCATGCTGTTCGATCAGGCGCCGGGTGATCGCCACTCGGGCGCGATAGAACTCGCTGCTGCCGTGGCTTGCCTCGCCCAGCAGGATGACGCGCGCGTCAGCGTAGCGGTCGAACAGCGCTGCGAATGCCGGCGTATCCAGGTTCGGCAAAGGTTCGGCGTACTGGTAAAGCACCGGGGCTATATCAGCAGTGTCGACCTGATACTGCTTGCCGTAGAGTTTTTCCATCAATTTCTTCGAAGGGCCGTTCATGGGGAGCCTCC  >  NC_002947/3668408‑3668697
                                                                                                                                                   |                                                                                                                                              
gTCGGCCTCGACGGCGACGATGTTGAAGCCATGCTGTTCGATCAGGCGCCGGGTGATCGCCACTCGGGCGCGATAGAACTCGCTGCTGCCGTGGCTTGCCTCGCCCAGCAGGATGACGCGCGCGTCAGCGTAGCGGTCGAACAGCGCTg                                                                                                                                               <  5:369843/149‑1 (MQ=255)
   ggCCTCGACGGCGACGATGTTGAAGCCATGCTGTTCGATCAGGCGCCGGGTGATCGCCACTCGGGCGCGATAGAACTCGCTGCTGCCGTGGCTTGCCTCGCCCAGCAGGATGACGCGCGCGTCAGCGTAGCGGTCGAACAGCGCTGCg                                                                                                                                             >  4:289458/1‑148 (MQ=255)
    gCCTCGACGGCGACGATGTTGAAGCCATGCTGTTCGATCAGGCGCCGGGTGATCGCCACTCGGGCGCGATAGAACTCGCTGCTGCCGTGGCTTGCCTCGCCCAGCAGGATGACGCGCGCGTCAGCGTAGCGGTCGAACAGCGCTGCGaa                                                                                                                                           <  4:156804/149‑1 (MQ=255)
           cGGCGACGATGTTGAAGCCATGCTGTTCGATCAGGCGCCGGGTGATCGCCACTCGGGCGCGATAGAACTCGCTGCTGCCGTGGCTTGCCTCGCCCAGCAGGATGACGCGCGCGTCAGCGTAGCGGTCGAACAGCGCTGCGAATGCCGGc                                                                                                                                    <  7:53434/149‑1 (MQ=255)
             gCGACGATGTTGAAGCCATGCTGTTCGATCAGGCGCCGGGTGATCGCCACTCGGGCGCGATAGAACTCGCTGCTGCCGTGGCTTGCCTCGCCCAGCAGGATGACGCGCGCGTCAGCGTAGCGGTCGAACAGCGCTGCGAATGCCGGCg                                                                                                                                   >  8:344189/1‑148 (MQ=255)
              cgacgaTGTTGAAGCCATGCTGTTCGATCAGGCGCCGGGTGATCGCCACTCGGGCGCGATAGAACTCGCTGCTGCCGTGGCTTGCCTCGCCCAGCAGGATGACGCGCGCGACAGCGTAGCGGTCGAACAGCGCAGCGAATGCCGGCGtt                                                                                                                                 >  5:410541/1‑148 (MQ=255)
               gacgaTGTTGAAGCCATGCTGTTCGATCAGGCGCCGGGTGATCGCCACTCGGGCGCGATAGAACTCGCTGCTGCCGTGGCTTGCCTCGCCCAGCAGGATGACGCGCGCGTCAGCGTAGCGGTCGAACAGCGCTGCGAATGCCGGCGTAt                                                                                                                                >  7:229179/1‑149 (MQ=255)
                 cgaTGTTGAAGCCATGCTGTTCGATCAGGCGCCGGGTGATCGCCACTCGGGCGCGATAGAACTCGCTGCTGCCGTGGCTTGCCTCGCCCAGCAGGATGACGCGCGCGTCAGCGTAGCGGTCGAACAGCGCTGCGAATGCCGGCGTATcc                                                                                                                              <  4:68136/149‑1 (MQ=255)
                      ttGAAGCCATGCTGTTCGATCAGGCGCCGGGTGATCGCCACTCGGGCGCGATAGAACTCGCTGCTGCCGTGGCTTGCCTCGCCCAGCAGGATGACGCGCGCGTCAGCGTAGCGGTCGAACAGCGCTGCGAATGCCGGCGTATCCAGGt                                                                                                                          <  5:268736/148‑1 (MQ=255)
                            ccATGCTGTTCGATCAGGCGCCGGGTGATCGCCACTCGGGCGCGATAGAACTCGCTGCTGCCGTGGCTTGCCTCGCCCAGCAGGATGACGCGCGCGTCAGCGTAGCGGTCGAACAGCGCTGCGAATGCCGGCGTATCCAGGTTCGGCa                                                                                                                    >  5:138158/1‑148 (MQ=255)
                              aTGCTGTTCGATCAGGCGCCGGGTGATCGCCACTCGGGCGCGATAGAACTCGCTGCTGCCGTGGCTTGCCTCGCCCAGCAGGATGACGCGCGCGTCAGCGTAGCGGTCGAACAGCGCCGCGAATGCCGGCGTATCCAGGTTCGGCAAAg                                                                                                                 >  6:360976/1‑149 (MQ=255)
                              aTGCTGTTCGATCAGGCGCCGGGTGATCGCCACTCGGGCGCGATAGAACTCGCTGCTGCCGTGGCTTGCCTCGCCCAGCAGGATGACGCGCGCGTCAGCGTAGCGGTCGAACAGCGCCGCGAATGCCGGCGTATCCAGGTTCGGCAAAg                                                                                                                 >  6:367315/1‑149 (MQ=255)
                                 cTGTTCGATCGGGCGCCGGGTGATCGCCACTCGGGCGCGATAGAACTCGCTGCTGCCGTGGCTTGCCTCGCCCAGCAGGATGACGCGCGCGTCAGCGTAGCGGTCGAACAGCGCTGCGAATGCCGGCGTATCCAGGTTCGGCAAAGGtt                                                                                                              <  2:379406/149‑1 (MQ=255)
                                 cTGTTCGATCAGGCGCCGGGTGATCGCCACTCGGGCGCGATAGAACTCGCTGCTGCCGTGGCTTGCCTCGCCCAGCAGGATGACGCGCGCGTCAGCGTAGCGGTCGAACAGCGCTGCGAATGCCGGCGTATCCAGGTTCGGCAAAg                                                                                                                 >  2:7846/1‑146 (MQ=255)
                                 cTGTTCGATCAGGCGCCGGGTGATCGCCACTCGGGCGCGATAGAACTCGCTGCTGCCGTGGCTTGCCTCGCCCAGCAGGATGACGCGCGCGTCAGCGTAGCGGTCGAACAGCGCTGCGAATGCCGGCGTATCCAGGTTCGGCAAAg                                                                                                                 <  1:7846/146‑1 (MQ=255)
                                    ttCGATCAGGCGCCGGGTGATCGCCACTCGGGCGCGATAGAACTCGCTGCTGCCGTGGCTTGCCTCGCCCAGCAGGATGACGCGCGCGTCAGCGTAGCGGTCGAACAGCGCTGCGAATGCCGGCGTATCCAGGTTCGGCAAAGGTTCg                                                                                                            <  8:271285/148‑1 (MQ=255)
                                      cGATCAGGCGCCGGGTGATCGCCACTCGGGCGCGATAGAACTCGCTGCTGCCGTGGCTTGCCTCGCCCAGCAGGATGACGCGCGCGTCAGCGTAGCGGTCGAACAGCGCTGCGAATGCCGGCGTATCCAGGTTCGGCAAAGGTTCGGCg                                                                                                         >  2:27092/1‑149 (MQ=255)
                                        aTCAGGCGCCGGGTGATCGCCACTCGGGCGCGATAGAACTCGCTGCTGCCGTGGCTTGCCTCGCCCAGCAGGATGACGCGCGCGTCAGCGTAGCGGTCGAACAGCGCTGCGAATGCCGGCGTATCCAGGTTCGGCAAAGGTTCGGCGTa                                                                                                       <  8:229179/149‑1 (MQ=255)
                                            ggCGCCGGGTGATCGCCACTCGGGCGCGATAGAACTCGCTGCTGCCGTGGCTTGCCTCGCCCAGCAGGATGACTCGCGCGTCAGCGTAGCGGTCGAACAGCGCTGCGAATGCCGGCGTATCCAGGTTCGGCAAAGGTTCTGCGTACTg                                                                                                    >  6:359523/1‑148 (MQ=255)
                                             gcgcCGGGTGATCGCCACTCGGGCGCGATAGAACTCGCTGCTGCCGTGGCTTGCCTCGCCCAGCAGGATGACGCGCGCGTCAGCGTAGCGGTCGAACAGCGCTGCGAATGCCGGCGTATCCAGGTTCGGCAAAGGTTCGGCGTACTgg                                                                                                   >  5:351599/1‑148 (MQ=255)
                                              cgcCGGGTGATCGCCACTCGGGCGCGATAGAACTCGCTGCTGCCGTGGCTTGCCTCGCCCAGCAGGATGACGCGCGCGTCAGCGTAGCGGTCGAACAGCGCCGCGAATGCCGGCGTATCCAGGTTCGGCAAAGGTTCGGCGTACTGGTa                                                                                                 <  5:367315/149‑1 (MQ=255)
                                              cgcCGGGTGATCGCCACTCGGGCGCGATAGAACTCGCTGCTGCCGTGGCTTGCCTCGCCCAGCAGGATGACGCGCGCGTCAGCGTAGCGGTCGAACAGCGCCGCGAATGCCGGCGTATCCAGGTTCGGCAAAGGTTCGGCGTACTGGTa                                                                                                 <  5:360976/149‑1 (MQ=255)
                                                  gggTGATCGCCACTCGGGCGCGATAGAACTCGCTGCTGCCGTGGCTTGCCTCGCCCAGCAGGATGACGCGCGCGTCAGCGTAGCGGTCGAACAGCGCTGCGAATGCCGGCGTATCCAGGTTCGGCAAAGGTTCGGCGTACTGGTaa                                                                                                >  6:241947/1‑146 (MQ=255)
                                                  gggTGATCGCCACTCGGGCGCGATAGAACTCGCTGCTGCCGTGGCTTGCCTCGCCCAGCAGGATGACGCGCGCGTCAGCGTAGCGGTCGAACAGCGCTGCGAATGCCGGCGTATCCAGGTTCGGCAAAGGTTCGGCGTACTGGTaa                                                                                                <  5:241947/146‑1 (MQ=255)
                                                       aTCGCCACTCGGGCGCGATAGAACTCGCTGCGGCCGTGGCTTGCCTCGCCCAGCAGGATGACGCGCGCGTCAGCGTAGCGGTCGAACAGCGCTGCGAATGCCGGCGTATCCAGGTTCGGCAAAGGTTCGGCGTACTGGTAAAGCACCgg                                                                                        <  6:351599/149‑1 (MQ=255)
                                                         cGCCACTCGGGCGCGATAGAACTCGCTGCTGCCGTGGCTTGCCTCGCCCAGCAGGATGACGCGCGCGTCAGCGTAGCGGTCGAACAGCGCTGCGAATGCCGGCGTATCCAGGTTCGGCAAAGGTTCGGCGTACTGGTAAAGCACCgggg                                                                                      >  8:128964/1‑149 (MQ=255)
                                                         cGCCACTCGGGCGCGATAGAACTCGCTGCTGCCGTGGCTTGCCTCGCCCAGCAGGATGACGCGCGCGTCAGCGTAGCGGTCGAACAGCGCTGCGAATGCCGGCGTATCCAGGTTCGGCAAAGGTTCGGCGTACTGGTAAAGCACCgggg                                                                                      >  3:330811/1‑149 (MQ=255)
                                                          gCCACTCGGGCGCGATAGAACTCGCTGCTGCCGTGGCTTGCCTCGCCCAGCAGGATGACGCGCGCGTCAGCGTAGCGGTCGAACAGCGCTGCGAATGCCGGCGTATCCAGGTTCGGCAAAGGTTCGGCGTACTGGTAAAGCACCGGGGc                                                                                     <  1:27092/149‑1 (MQ=255)
                                                           ccACTCGGGCGCGATAGAACTCGCTGCTGCCGTGGCTTGCCTCGCCCAGCAGGATGACGCGCGCGTCAGCGTAGCGGTCGAACAGCGCTGCGAATGCCGGCGTATCCAGGTTCGGCAAAGGTTCGGCGTACTGGTAAAGCACCGGGGc                                                                                     >  8:165992/1‑148 (MQ=255)
                                                                  ggCGCGATAGAACTCGCTGCTGCCGTGGCTTGCCTCGCCCAGCAGGATGACGCGCGCGTCAGCGTAGCGGTCGAACAGCGCTGCGAATGCCGGCGTATCCAGGTTCGGCAAAGGTTCGGCGTACTGGTAAAGCACCGGGGCTATATcag                                                                             <  7:128964/149‑1 (MQ=255)
                                                                       gATAGAACTCGCTGCTGCCGTGGCTTGCCTCGCCCAGCAGGATGACGCGCGCGTCAGCGTAGCGGTCGAACAGCGCTGCGAATGCCGGCGTATCCAGGTTCGGCAAAGGTTCGGCGTACTGGTAAAGCACCGGGGCTATATcagca                                                                           >  3:353393/1‑146 (MQ=255)
                                                                       gATAGAACTCGCTGCTGCCGTGGCTTGCCTCGCCCAGCAGGATGACGCGCGCGTCAGCGTAGCGGTCGAACAGCGCTGCGAATGCCGGCGTATCCAGGTTCGGCAAAGGTTCGGCGTACTGGTAAAGCACCGGGGCTATATcagca                                                                           <  4:353393/146‑1 (MQ=255)
                                                                        aTAGAACTCGCTGCTGCCGTGGCTTGCCTCGCCCAGCAGGATGACGCACGCGTCAGCGTAGCGGTCGAACAGCGCTGCGAATGCCGGCGTATCCAGGTTCGGCAAAGGTTCGGCGTACTGGTAAAGCACCGGGGCTATATCAGCAgtgt                                                                       >  4:2968/1‑149 (MQ=255)
                                                                         tAGAACTCGCTGCTGCCGTGGCTTGCCTCGCCCAGCAGGATGACGCGCGCGTCAGCGTAGCGGTCGAACAGCGCTGCGAATGCCGGCGTATCCAGGTTCGGCAAAGGTTCGGCGTACTGGTAAAGCACCGGGGCTATATCAGCAGTGTc                                                                      <  7:165992/149‑1 (MQ=255)
                                                                                    gctgcCGTGGCTTGCCTCGCCCAGCAGGATGACGCGCGCGTCAGCGTAGCGGTCGAACAGCGCTGCGAATGCCGGCGTATcc                                                                                                                              >  3:224927/1‑82 (MQ=255)
                                                                                    gctgcCGTGGCTTGCCTCGCCCAGCAGGATGACGCGCGCGTCAGCGTAGCGGTCGAACAGCGCTGCGAATGCCGGCGTATcc                                                                                                                              <  4:224927/82‑1 (MQ=255)
                                                                                     ctgcCGTGGCTTGCCTCGCCCAGCAGGATGACGCGCGCGTCAGCGTAGCGGTCGAACAGCGCTGCGAATGCCGGCGTATCCAGGTTCGGCAAAGGTTCGGCGTACTGGTAAAGCACCGGGGCTATATCAGCAGTGTCGACCTGATAc                                                            <  6:138158/147‑1 (MQ=255)
                                                                                        ccGTGGCTTGCCTCGCCCAGCAGGATGACGCGCGCGTCAGCGTAGCGGTCGAACAGCGCTGCGAATGCCGGCGTATCCAGGTTCGGCAAAGGTTc                                                                                                             >  5:365065/1‑95 (MQ=255)
                                                                                        ccGTGGCTTGCCTCGCCCAGCAGGATGACGCGCGCGTCAGCGTAGCGGTCGAACAGCGCTGCGAATGCCGGCGTATCCAGGTTCGGCAAAGGTTc                                                                                                             <  6:365065/95‑1 (MQ=255)
                                                                                        ccGTGGCTTGCCTCGCCCAGCAGGATGACGCGCGCGTCAGCGTAGCGGTCGAACAGCGCTGCGAATGCCGGCGTATCCAGGTTCGGCAAAGGTTCGGCGTACTGGTAAAGCACCGGGGCTATATCAGCAGTGTCGACCTGATACTGCtt                                                       <  7:344189/149‑1 (MQ=255)
                                                                                         cGTGGCTTGCCTCGCCCAGCAGGATGACGCGCGCGTCAGCGTAGCGGTCGAACAGCGCTGCGAATGCCGGCGTATCCAGGTTCGGCAAAGGTTCGGCGTACTGGTAAAGCACCGGGGCTATATcagca                                                                           <  2:235401/128‑1 (MQ=255)
                                                                                         cGTGGCTTGCCTCGCCCAGCAGGATGACGCGCGCGTCAGCGTAGCGGTCGAACAGCGCTGCGAATGCCGGCGTATCCAGGTTCGGCAAAGGTTCGGCGTACTGGTAAAGCACCGGGGCTATATcagca                                                                           >  1:235401/1‑128 (MQ=255)
                                                                                          gTGGCTTGCCTCGCCCAGCAGGATGACGCGCGCGTCAGCGTAGCGGTCGAACAGCGCTGCGAATGCCGGCGTATCCAGGTTCGGCAAAGGTTCGGCGTACTGGTaaa                                                                                               <  4:249110/107‑1 (MQ=255)
                                                                                          gTGGCTTGCCTCGCCCAGCAGGATGACGCGCGCGTCAGCGTAGCGGTCGAACAGCGCTGCGAATGCCGGCGTATCCAGGTTCGGCAAAGGTTCGGCGTACTGGTaaa                                                                                               >  3:249110/1‑107 (MQ=255)
                                                                                             gCTTGCCTCGCCCAGCAGGATGACGCGCGCGTCAGCGTAGCGGTCGAACAGCGCTGCGAATGCCGGCGTATCCAGGTTCGGCAAAGGTTCGGCGTACTGGTAAAGCACCGGGGCTATATCAGCAGTGTCGACCTGATACTGCTTGCCg                                                   >  5:176693/1‑148 (MQ=255)
                                                                                              cTTGCCTCGCCCAGCAGGATGACGCGCGCGTCAGCGTAGCGGTCGAACAGCGCTGCGAATGCCGGCGTATCCAGGTTCGGCAAAGGTTCGGCGTACTGGTAAAGCACCGGGGCTATATCAGCAGTGTCGACCTGATACTGCTTGCCGTa                                                 <  3:289458/149‑1 (MQ=255)
                                                                                               ttGCCTCGCCCAGCAGGATGACGCGCGCGTCAGCGTAGCGGTCGAACAGCGCTGCGAATGCCGGCGTATCCAGGTTCGGCAAAGGTTCGGCGTACTGGTAAAGCACCGGGGCTATATCAGCAGTGTCGACCTGATACTGCTTGCCg                                                   >  3:391276/1‑146 (MQ=255)
                                                                                               ttGCCTCCCCCAGCAGGATGACGCGCGCGTCAGCGTAGCGGTCGAACAGCGATGCGAAAGCCCGCGTATCCAGGTTCGGCAAAGGTTCGGCCTACTGGTAAAGCACCCGGGCTATATCAGCAGTGTCGACCTGATACTGCTTCCCg                                                   >  1:288885/1‑146 (MQ=255)
                                                                                                 gCCTCGCCCAGCAGGATGACGCACGCGTCAGCGTAGCGGTCGAACAGCGCTGCGAATGCCGGCGTATCCAGGTTCGGCAAAGGTTCGGCGTACTGGTAAAGCACCGGGGCTATATCAGCAGTGTCGACCTGATACTGCTTGCCGTaga                                               <  3:2968/148‑1 (MQ=255)
                                                                                                                   aCGCGCGCGTCAGCGTAGCGGTCGAACAGCGCTGCGAATGCCGGCGTATCCAGGTTCGGCAAAGGTTCGGCGTACTGGTAAAGCTCCGGGGCTATATCAGCAGTGTCGACCTGATACTGCTTGCCGTAGAGTTTTTCCATCAATttct                             <  5:359523/148‑1 (MQ=255)
                                                                                                                   aCGCGCGCGTCAGCGTAGCGGTCGAACAGCGCTGCGAATGCCGGCGTATCCAGGTTCGGCAAAGGTTCGGCGTACTGGTAAAGCACCGGGGCTATATCAGCAGTGTCGA‑CTGATACTGCTTGCCGTAGAGTTTTTCCATCAATttctt                            <  4:330811/148‑1 (MQ=255)
                                                                                                                       gcgcgTCAGCGTAGCGGTCGAACAGCGCTGCGAATGCCGGCGTATCCAGGTTCGGCAAAGGTTCGGCGTACTGGTAAAGCACCGGGGCTATATCAGCAGTGTCGACCTGATACTGCTTGCCGTAGAGTTTTTCCATCAATTTCTTCGaa                        <  4:391276/149‑1 (MQ=255)
                                                                                                                                     cGGTCGAACAGCGCTGCGAATGCCGGCGTATCCAGGTTCGGCAAAGGTTCGGCGTACTGGTAAAGCACCGGGGCTATATCAGCAGTGTCGACCTGATACTGCTTGCCGTAGAGTTTTTCCATCAATTTCTTCGAAGGGCCGTTCATggg          <  6:132638/149‑1 (MQ=255)
                                                                                                                                     cGGTCGAACAGCGCTGCGAATGCCGGCGTATCCAGGTTCGGC‑AAGGTTCGGCGTACTGGTAAAGCACCGGGGCTATATCAGCAGTGTCGACCTGATACTGCTTGCCGTAGAGTTTTTCCa                                      <  6:118014/120‑1 (MQ=255)
                                                                                                                                     cGGTCGAACAGCGCTGCGAATGCCGGCGTATCCAGGTTCGGC‑AAGGTTCGGCGTACTGGTAAAGCACCGGGGCTATATCAGCAGTGTCGACCTGATACTGCTTGCCGTAGAGTTTTTCCa                                      >  5:118014/1‑120 (MQ=255)
                                                                                                                                             cAGCGCTGCGAATGCCGGCGTATCCAGGTTCGGCAAAGGTTCGGCGTACTGGTAAAGCACCGGGGCTATATCAGCAGTGTCGACCTGATACTGCTTGCCGTAGAGTTTTTCCATGAATTTCTTCGAAGGGCCGTTCATGCGGAGCCtac  >  5:300597/1‑147 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
GTCGGCCTCGACGGCGACGATGTTGAAGCCATGCTGTTCGATCAGGCGCCGGGTGATCGCCACTCGGGCGCGATAGAACTCGCTGCTGCCGTGGCTTGCCTCGCCCAGCAGGATGACGCGCGCGTCAGCGTAGCGGTCGAACAGCGCTGCGAATGCCGGCGTATCCAGGTTCGGCAAAGGTTCGGCGTACTGGTAAAGCACCGGGGCTATATCAGCAGTGTCGACCTGATACTGCTTGCCGTAGAGTTTTTCCATCAATTTCTTCGAAGGGCCGTTCATGGGGAGCCTCC  >  NC_002947/3668408‑3668697

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