Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,104,698 C→T 10.7% E62E (GAG→GAA PP_1879 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,104,6980CT10.7% 147.5 / 6.9 57E62E (GAG→GAAPP_1879hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (25/25);  new base T (4/2);  total (30/27)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.25e-01

TAAGCTTGTCGACCGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCG  >  NC_002947/2104553‑2104839
                                                                                                                                                 |                                                                                                                                              
tAAGCTTGTCGACCGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTcc                                                                                                                                             <  2:124625/149‑1 (MQ=255)
  aGCTTGTCGACCGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCa                                                                                                                                            <  7:101419/148‑1 (MQ=255)
   gCTTGTCGACCGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAGTTCCAcc                                                                                                                                         >  1:353708/1‑149 (MQ=255)
   gCTTGTCGACCGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCAcc                                                                                                                                         >  8:10188/1‑149 (MQ=255)
        tCGACCGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAGTTCCACCCgtgg                                                                                                                                    <  2:353708/149‑4 (MQ=255)
          gACCGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAggg                                                                                                                                  <  7:10188/149‑1 (MQ=255)
            ccGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTcacc                                                                                                                                           >  5:294490/1‑136 (MQ=255)
            ccGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTcacc                                                                                                                                           <  6:294490/139‑4 (MQ=255)
            ccGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTACAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAggggg                                                                                                                                <  7:164284/149‑1 (MQ=255)
                    aaTATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTAAGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTTCCCACGGGTGGAAATCAACACGCAGGGGGGAATAGGt                                                                                                                        >  1:80283/1‑149 (MQ=255)
                         gAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAGTTACACCCgtg                                                                                                                                     >  3:145284/1‑129 (MQ=255)
                              gCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGaa                                                                                                               <  6:264506/148‑1 (MQ=255)
                                 ctTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTc                                                                                                            >  7:36982/1‑148 (MQ=255)
                                                    gCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTccc                                                                                          <  8:36982/147‑1 (MQ=255)
                                                      tGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCATCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCGGGAACTCGACCGTGTTGGGTTTCCCCCCt                                                                                      >  4:211937/1‑149 (MQ=255)
                                                       gAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCtt                                                                                     >  3:290186/1‑149 (MQ=255)
                                                        aaGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCtt                                                                                     >  3:224667/1‑148 (MQ=255)
                                                        aaGGGCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTc                                                                                    <  6:40775/149‑1 (MQ=255)
                                                             cttcttGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACAt                                                                               >  5:74109/1‑149 (MQ=255)
                                                               tcttGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTACCCCCCTGCAACATc                                                                              >  6:12536/1‑148 (MQ=255)
                                                                cttGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGt                                                                            >  4:293365/1‑149 (MQ=255)
                                                                    tCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTccac                                                                        >  7:325504/1‑149 (MQ=255)
                                                                    tCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTcca                                                                         <  1:173724/148‑1 (MQ=255)
                                                                                  aTTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACgg                                                          <  1:335020/149‑1 (MQ=255)
                                                                                  aTTCCAGTAGAGGAAAAGGTATTCCTCATCTAATTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGTGGGGAATAGCTTGCCCGGAACTCGACCGAGTTGGGTTTCCACCCTTCATCATAGTCCACAACGCTTGGTACgt                                                          >  1:218178/1‑148 (MQ=255)
                                                                                       aGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTggg                                                     <  4:224667/149‑1 (MQ=255)
                                                                                        gTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTgggg                                                    <  2:114486/149‑1 (MQ=255)
                                                                                         tAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTACACACTCAGGGGGGAATAGGTTGACCGGAACTCGACCGTGTTGGGTTTCAACCCTTCAACATAGTCAACCAGGCTTGTTACGGTTggggg                                                   >  5:312149/1‑149 (MQ=255)
                                                                                          agagGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGCTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCTCCATCGTCCACCCCGCTTGGTCCGGTTcgggg                                                   >  7:332346/1‑148 (MQ=255)
                                                                                           gagGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGa                                                 >  2:87003/1‑149 (MQ=255)
                                                                                              gAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGAc                                                >  1:54661/1‑147 (MQ=255)
                                                                                              gAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACaa                                              <  4:122150/149‑1 (MQ=255)
                                                                                               aaaaGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACaa                                              <  6:100103/148‑1 (MQ=255)
                                                                                                aaaGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGtt                                                                                                   <  7:324662/94‑1 (MQ=255)
                                                                                                aaaGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGtt                                                                                                   >  8:324662/1‑94 (MQ=255)
                                                                                                         ccTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTAAGGTTGGTGGGACAAGAGAGGACGAg                                   >  4:238571/1‑149 (MQ=255)
                                                                                                          cTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGTGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGt                                  <  6:312149/149‑1 (MQ=255)
                                                                                                          cTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGt                                  >  5:77202/1‑149 (MQ=255)
                                                                                                             aTCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCATCCCGCAGGGGGGAATAGGTTGCCCGGAACTCGAc                                                                                                         <  1:242741/76‑1 (MQ=255)
                                                                                                             aTCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAggggg                                                                                                                                <  4:15239/52‑1 (MQ=255)
                                                                                                             aTCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTc                                 <  1:87003/147‑1 (MQ=255)
                                                                                                               cTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGcc                                      >  6:165392/1‑140 (MQ=255)
                                                                                                               cTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGcc                                      <  5:165392/140‑1 (MQ=255)
                                                                                                                    tCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAg                         <  3:293365/148‑1 (MQ=255)
                                                                                                                       gggTGCGAGTTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCAACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGt                     >  5:359366/1‑149 (MQ=255)
                                                                                                                             gAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGt               >  5:180581/1‑149 (MQ=255)
                                                                                                                              aaTTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGtt              >  1:309729/1‑149 (MQ=255)
                                                                                                                                   gCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTTGAGTTCTGc          <  6:359366/148‑1 (MQ=255)
                                                                                                                                   gCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGc          <  6:77202/148‑1 (MQ=255)
                                                                                                                                    cctgcGGGTGGAGTTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCgg        >  7:128105/5‑149 (MQ=255)
                                                                                                                                    cctgcGGGTGGAGTTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCgg        >  4:343838/5‑149 (MQ=255)
                                                                                                                                    cccACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCg         <  4:290186/148‑1 (MQ=255)
                                                                                                                                    cccACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCg         <  2:309729/148‑1 (MQ=255)
                                                                                                                                     ccACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGt       >  1:138613/1‑149 (MQ=255)
                                                                                                                                       aCGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCGAGTTCAGGGGTTGGAGTTCTGCGGTTa     >  5:237470/1‑149 (MQ=255)
                                                                                                                                         gggTGGAGTTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTAcc   <  8:128105/149‑1 (MQ=255)
                                                                                                                                            tGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCg  >  7:260241/1‑147 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
TAAGCTTGTCGACCGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCG  >  NC_002947/2104553‑2104839

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: