Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,989,218 A→C 7.1% Y42D (TAT→GAT)  PP_2613 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,989,2180AC 7.1% 168.1 / 4.8 56Y42D (TAT→GAT) PP_2613hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (24/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.34e-01

CCCAAACTCATGACGGTCATCCAGGGCGCTGTTCTTCCATAACCATGCCTTTTGAGTATTTCGGCTGACAGCGTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAGTTGTGAGCGAAGTCTTCGCTGCGCGGAGGAGGTGATGACAAGGTGGTCATTGGACTTGCCCTCGGATACAGCGGCAATCGGGTAGGGGACAGGAGCCATC  >  NC_002947/2989078‑2989364
                                                                                                                                            |                                                                                                                                                  
cccAAACTCATGACGGTCATCCAGGGCGCTGTTCTTCCATAACCATGCCTTTTGAGTATTTCGGCTGACAGCGTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCTATCGGCCCatc                                                                                                                                            <  7:31683/149‑3 (MQ=255)
cccAAACTCATGACGGTCATCCAGGGCGCTGTTCTTCCATAACCATGCCTTTTGAGTATTTCGGCTGACAGCGTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGt                                                                                                                                            >  5:289418/1‑149 (MQ=255)
 ccAAACTCATGACGGTCATCCAGGGCGCTGTTCTTCCATAACCATGCCTTTTGAGTATTTCGGCTGACAGCGTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGtt                                                                                                                                           >  2:273964/1‑149 (MQ=255)
  cAAACTCATGACGGTCATCCAGGGCGCTGTTCTTCCATAACCATGCCTTTTGAGTATTTCGGCTGACAGCGTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTtc                                                                                                                                          >  4:193191/1‑149 (MQ=255)
  cAAACTCATGACGGTCATCCAGGGCGCTGTTCTTCCATAACCATGCCTTTTGAGTATTTCGGCTGACAGCGTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTtc                                                                                                                                          >  5:140799/1‑149 (MQ=255)
   aaaCTCATGACGGTCATCCAGGGCTCTGTTCTTCCATAACCATGCCTTTTGAGTATTTCGGCTGACAGCGTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTtct                                                                                                                                         <  1:235016/149‑1 (MQ=255)
   aaaCTCATGACGGTCATCCAGGGAGCTGTTCTTCCATAACCATGCCTTTTGAGTATTTCGGCTGAAAGAGTTTATAAGGTCTAGTAAAGATGGGGCATAAGAAGGAAACAGTTAACCTTCATGCGGGATGGGACGATAGGCACAGtta                                                                                                                                          >  4:6314/1‑147 (MQ=255)
   aaaCACATGACGGACATCCAGGGAGCTGTTCTTCCATAACCATGCCTTTTGAGTATTTAGGCTCACAGCGTTTATAAGGTAGAGTACAGAAGGGGCATAAGCACGAACCACTTCACCAAACTGCGGCAAGGACCGAAAGGCCCACttat                                                                                                                                         >  4:347542/1‑147 (MQ=255)
    aaCTCATGACGGTCATCCAGGGCGCTGTTCTTCCATAACCATGCCTTTTGAGTATTTCGGCTGACAGCGTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTtctc                                                                                                                                        <  8:250211/149‑1 (MQ=255)
      cTCATGACGGTCATCCAGGGCGCTGTTCTTCCATAACCATGCCTTTTGAGTATTTCGGCTGACAGCGTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGa                                                                                                                                      <  5:371085/149‑1 (MQ=255)
          tttcGGGCATCCAGGGCGGTGTTGTTGGATAACCATGCGTTTTCATTATTTCGGCTGACATCGTTTATAATGTCGAGTACGTATGGGTCATAAGCATGATCCAGTTCACGTTCTTGCGGGATGGGCCGATAGGCCCAGTTCTCGAttt                                                                                                                                   <  3:347542/145‑2 (MQ=255)
           gACGGTCATCCAGGGCGCTGTTCTTCCATAACCATGCCTTTTGAGTATTTCGGCTGACAGCGTTTATAAGGTCGAGTACCGATGGGGAATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGc                                                                                                                                  >  5:293209/1‑148 (MQ=255)
              ggTCATCCAGGGCGCTGTTCTTCCATAACCATGCCTTTTGAGTATTTCGGCTGACAGCGTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTg                                                                                                                              <  5:314079/149‑1 (MQ=255)
                tCATCCAGGGCGCTGTTCTTCCATAACCATGCCTTTTGAGTATTTCGGCTGACAGCGTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGt                                                                                                                            <  3:193191/149‑1 (MQ=255)
                 cATCCAGGGCGCTGTTCTTCCATAACCATGCCTTTTGAGTATTTCGGCTGACAGCGTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGt                                                                                                                            >  7:193962/1‑148 (MQ=255)
                    ccAGAGCGCTGTTCTTCCATAACCATGCCTTTTGAGTATTTCGGCTGACAGCGTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCAt                                                                                                                         >  8:178271/1‑148 (MQ=255)
                        ggCGCTGTTCTTCCATAACCATGCCTTTTGAGTATTTCGGCTGACAGCGTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGagca                                                                                                                    <  6:289418/149‑1 (MQ=255)
                         gcgcTGTTCTTCCATAACCATGCCTTTTGAGTATTTCGGCTGACAGCGTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGagcag                                                                                                                   <  4:281488/149‑1 (MQ=255)
                           gcTGTTCTTCCATAACCATGCCTTTTGAGTATTTCGGCTGACAGCGTTTATAAGGTCGAGTCCCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGAgcagca                                                                                                                 <  4:237455/149‑1 (MQ=255)
                           gcTGTTCTTCCATAACCATGCCTTTTGAGTATTTCGGCTGACAGCGTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGAgcagca                                                                                                                 <  8:44867/149‑1 (MQ=255)
                                   tcCATAACCATGCCTTTTGAGTATTTCGGCTGACAGCGTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCgg                                                                                                         >  7:412374/1‑149 (MQ=255)
                                    ccATAACCATGCCTTTTGAGTATTTCGGCTGACAGCGTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCgg                                                                                                         <  2:238207/148‑1 (MQ=255)
                                     cATAACCATGCCTTTTGAGTATTTCGGCTGACTGCTTTTATTAGGTCGTGTACCGATGGTGCATAAGCAGGAACCAGTTCACCTTCGTGCGTGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGag                                                                                                                      <  8:182325/134‑1 (MQ=255)
                                         aCCATGCCTTTTGAGTATTTCGGCTGACAGCGTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAGTtg                                                                                                   <  6:140799/149‑1 (MQ=255)
                                            ttGCCTTTTGAGTATTTCGGATGACAGCGTTTATAAGGTCGAGTACCGATTGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGTATGGGCCGATATGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAGTTGTGa                                                                                                <  3:6314/148‑1 (MQ=255)
                                                cTTTTGAGTATTTCGGCTGACAGCTTTTATTAGGTCTAGTACCGATGGGGGATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAGTTGTGAGCGa                                                                                            <  8:412374/149‑1 (MQ=255)
                                                  tttGAGTATTTCGGCTGACAGCGTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAGTTGTGAGCGAAg                                                                                          >  6:230806/1‑149 (MQ=255)
                                                                    cAGCGTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAGTTGTGAGCGAAGTCTTCGCTGCGCggagga                                                                        <  1:273964/149‑1 (MQ=255)
                                                                      gCGTTTATAAGGTCGAGTACCGATGGGGTATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAGTTGTGAGCGAAGTCTTCGCTGCGCGGaggagg                                                                      >  4:237039/1‑149 (MQ=255)
                                                                        gtctataCGGTCGAGTCCCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGTCATGAGCAGCAGAAAGCGGCAGTTGTGAGCGAAGTCTTCGCTGCGCGGAGGAGGtga                                                                   <  4:246192/146‑1 (MQ=255)
                                                                           tataAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAGTTGTGAGCGAAGTCTTCGCTGCGCGGAGGAGGtgatg                                                                 >  2:115190/1‑149 (MQ=255)
                                                                           tataAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAGTTGTGAGCGAAGTCTTCGCTGCGCGGAGGAGGtgat                                                                  <  3:380269/148‑1 (MQ=255)
                                                                           tataAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCACGAGCAGCAGAAAGCGGCAGTTGTGAGCGAAGTCTTCGCTGCGCGGAGGAGGtgat                                                                  <  8:193962/148‑1 (MQ=255)
                                                                           tataAGGTCGAGTACCGATGGGGAATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAGTTGTGAGCGAAGTCTTCGCTGCGCGGAGGAGGtgat                                                                  <  6:293209/148‑1 (MQ=255)
                                                                                ggTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAGTTGTGAGCGAAGTCTTCGCTGCGCGGAGGAGGTGATGACAAg                                                            >  1:277363/1‑149 (MQ=255)
                                                                                    gAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCa                                                                                                       >  6:240324/1‑102 (MQ=255)
                                                                                    gAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCa                                                                                                       <  5:240324/102‑1 (MQ=255)
                                                                                    gAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAATTGTGAGCGAAGTCTTCGCTGCGCggag                                                                          >  5:245387/1‑131 (MQ=255)
                                                                                    gAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAATTGTGAGCGAAGTCTTCGCTGCGCggag                                                                          <  6:245387/131‑1 (MQ=255)
                                                                                           gATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAGTTGTGAGCGAAGTCTTCGCTGCGCGGAGGAGGTGATGACAAGGTGGTCATTgg                                                 <  5:230806/149‑1 (MQ=255)
                                                                                                 gCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAGTTGTGAGCGAAGTCTTCGCTGCGCGGAGGAGGTGATGACAAGGTGGTCATTGGACtt                                             <  7:125509/147‑1 (MQ=255)
                                                                                                 gCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAGTTGTGAGCGAAGTCTTCGCTGCGCGGAGGAGGTGATGACAAGGTGGTCATTGGACtt                                             >  8:125509/1‑147 (MQ=255)
                                                                                                  cATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAGTTGTGAGCGAAGTCTTCGCTGCGCGGAGGAGGTGATGACAAGGTGGTCATTGGACTTGcc                                          >  7:220698/1‑149 (MQ=255)
                                                                                                   aTAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAGTTGTGAGCGAAGTCTTCGCTGCGCGGAGGAGGGGATGACAAGGTGGTCATTGGACTTGccc                                         >  3:284468/1‑149 (MQ=255)
                                                                                                       gCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAGTTGTGAGCGAAGTCTTCGCTGCGCGGAGGAGGTGATGACAAGGTGGTCATTGGACTTGCCCTCg                                      >  2:107213/1‑148 (MQ=255)
                                                                                                                    cACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAGTTGTGAGCGAAGTCTTCGCTGCGCGGAGGAGGTGATGACAAGGTGGTCATTGGACTTGCCCTCGGATACAGCGGCAAt                        >  2:229424/1‑149 (MQ=255)
                                                                                                                        ttCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAGTTGTGAGCGAAGTCTTCGCTGCGCGGAGGAGGTGATGACAAggtg                                                         >  7:176286/1‑112 (MQ=255)
                                                                                                                        ttCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAGTTGTGAGCGAAGTCTTCGCTGCGCGGAGGAGGTGATGACAAggtg                                                         <  8:176286/112‑1 (MQ=255)
                                                                                                                            tGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAGTTGTGAGCGAAGTCTTCGCTGCGCGGAGGAGGTGATGACAAGGTGGTCATTGGACTTGCCCTCGGATACAGCGGCAATCGGGTAgg                >  6:275081/1‑149 (MQ=255)
                                                                                                                              cGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAGTTGTGAGCGAAGTCTTCGCTGCGCGGAGGAGGTGATGACAAGGTGGTCATTGGACTTGCCCTCGGATACAGCGGCAATCGGGTAgg                <  2:277363/147‑1 (MQ=255)
                                                                                                                              agaactGGGCCTATCGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAGTTGTGAGCGAAGTCTTCGCTGCGCGGAGGAGGTGATGACAAGGTGGTCATTGGACTTGCCCTCGGATACAGCGGCAATCGGGTAggg               >  7:100103/6‑148 (MQ=255)
                                                                                                                                aactGGGCCTATCGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAGTTGTGAGCGAAGTCTTCGCTGCGCGGAGGAGGTGATGACAAGGTGGTCATTGGACTTGCCCTCGGATACAGCGGCAATCGGGTAgg                >  1:56608/4‑145 (MQ=255)
                                                                                                                                aactGGGCCTATCGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAGTTGTGAGCGAAGTCTTCGCTGCGCGGAGGAGGTGATGACAAGGTGGTCATTGGACTTGCCCTCGGATACAGCGGCAATCGGGTAgg                <  2:56608/142‑1 (MQ=255)
                                                                                                                                     ggCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAGTTGTGAGCGAAGTCTTCGCTGCGCGGAGGAGGTGATGACAAGGTGGTCATTGGACTTGCCCTCGGATACAGCGGCAATCGGGTAGGGGACAGGa        <  3:237039/148‑1 (MQ=255)
                                                                                                                                         gATAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAGTTGTGAGCGAAGTCTTCGCTGCGCGGAGGAGGTGATGACAAGGTGGTCATTGGACTTGCCCTCGGATACAGCGGCAATCGGGTAGGGGACAGGAGCCAt   >  7:261826/1‑149 (MQ=255)
                                                                                                                                          aTAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAGTTGTGAGCGAAGTCTTCGCTGCGCGGAGGAGGTGATGACAAGGTGGTCATTGGACTTGCCCTCGGATACAGCGGCAATCGGGTAGGGGACAGGAGCCATc  <  1:229424/149‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                  
CCCAAACTCATGACGGTCATCCAGGGCGCTGTTCTTCCATAACCATGCCTTTTGAGTATTTCGGCTGACAGCGTTTATAAGGTCGAGTACCGATGGGGCATAAGCAGGAACCAGTTCACCTTCCTGCGGGATGGGCCGATAGGCCCAGTTCTCGATTGCACTGGTCATGAGCAGCAGAAAGCGGCAGTTGTGAGCGAAGTCTTCGCTGCGCGGAGGAGGTGATGACAAGGTGGTCATTGGACTTGCCCTCGGATACAGCGGCAATCGGGTAGGGGACAGGAGCCATC  >  NC_002947/2989078‑2989364

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: