Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,676,076 A→T 8.0% L20* (TTG→TAG)  PP_2344 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,676,0760AT 8.0% 144.5 / 3.5 50L20* (TTG→TAG) PP_2344hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (22/24);  new base T (2/2);  total (24/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.54e-01

AGGGCACGCTGACGGGGTGCCGAGGGTAGCCACTCGCCGAACAGCGCCAGGTCCAGGCCC‑TGGATGCCCAGCGGCTTTGGCCCACCCGCCGCGCCCGGCAATACCCCGACATCGGCCGCCTGCAGGCCACGCTCGCGGATGTGTCGCAAGGCACGCTGGCCGGCCTTGAAGGTGAGGGCGGGGTATTTGATGTGGATGGCGCTCATGGCGGGTCTCTTGGGCATAGGAGCAGCCTTGTGCTGCGAGGGGCCATTACAGATAACACACCTGTAGGTGGCTGTATTGGCCCTTTCGCA  >  NC_002947/2675929‑2676224
                                                                                                                                                    |                                                                                                                                                    
aGGGCACGCTGACGGGGTGCCGAGGGTAGCCACTCGCCGAACAGCGCCAGGTCCAGGCCC‑TGGATGCCCAGCGGCTTTGGCCCACCCGCCGCGCCCGGCAATACCCCGACATCGGCCGCCTGCAGGCCACGCTCGCGGATGTGTCGCaa                                                                                                                                                     <  5:251532/149‑1 (MQ=255)
 gggCACGCTGACGGGGTGCCGAGGGTAGCCACTCGCCGAACAGCGCCAGGTCCAGGCCC‑TGGATGCCCAGCGGCTTTGGCCCACCCGCCGCGCCCGGCAATACCCCGACATCGGCCGCCTGCAGGCCACGCTCGCGGATGTGTCGCAAg                                                                                                                                                    <  1:51910/149‑1 (MQ=255)
     aCGCTGACGGGGTGCCGAGGGTAGCCACTCGCCGAACAGCGCCAGGTCCAGGCCC‑TGGATGCCCAGCGGCTTTGGCCCACCCGCCGCGCCCGGCAATACCCCGCCATCGGCCGCCTGCAGGCCACGCTCGCGGATGTGTCGCAAGGCAc                                                                                                                                                >  8:175953/1‑149 (MQ=255)
           aCGGGGTGCCGAGGGTAGCCACTCGCCGAACAGCGCCAGGTCCAGGCCC‑TGGATGCCCAGCGGCTTTGGCCCACCCGCCGCGCCCGGCAATACCCCGACATCGGCCGCCTGCAGGCCACGCTCGCGGATGTGTCGCAAGGCACGCTggc                                                                                                                                          <  5:94468/149‑1 (MQ=255)
           aCGGGGTGCCGAGGGTAGCCACTCGCCGAACAGCGCCAGGTCCAGGCCC‑TGGATGCCCAGCGGCTTTGGCCCACCCGCCGCGCCCGGCAATACCCCGACATCGGCCGCCTGCAGGCCACGCTCGCGGATGTGTCGCAAGGCACGCTg                                                                                                                                            <  5:43597/147‑1 (MQ=255)
               ggTGCCGAGGGTAGCCACTCGCCGAACAGCGCCAGGTCCAGGCCC‑TGGATGCCCAGCGGCTTTGGCCCACCCGCCGCGCCCGGCAATACCCCGACATCGGCCGCCTGCAGGCCACGCTCGCGGATGTGTCGCAAGGCACGCTggccggc                                                                                                                                      >  5:170507/1‑149 (MQ=255)
                    cGAGGGTAGCCACTCGCCGAACAGCGCCAGGTCCAGGCCC‑TGGATGCCCAGCGGCTTTGGCCCACCCGCCGCGCCCGGCAATACCCCGACATCGGCCGCCTGCAGGCCACGCTCGCGGATGTGTCGCAAGGCACGCTGGCCGGCCTTg                                                                                                                                  >  4:123191/1‑148 (MQ=255)
                      aGGGTAGCCACTCGCCGAACAGCGCCAGGTCCAGGCCC‑TGGATGCCCAGCGGCTTTGGCCCACCCGCCGCGCCCGGCAATACCCCGACATCGGCCGCCTGCAGGCCACGCTCGCGGATGTGTCGCTAGGCACGCTGGCCGGCCTTGAAg                                                                                                                               >  4:140394/1‑149 (MQ=255)
                      aGGGTAGCCACTCGCCGAACAGCGCCAGGTCCAGGCCC‑TGGATGCCCAGCGGCTTTGGCCCACCCGCCGCGCCCGGCAATACCCCGACATCGGCCGCCTGCAGGCCACGCTCGCGGATGTGTCGCTAGGCACGCTGGCCGGCCTTGAAg                                                                                                                               >  8:80095/1‑149 (MQ=255)
                              cACTCGCCGAACAGCGCCAGGTCCAGGCCCTTGGATGCCCAGCGGCTTTGGCCCACCCGCCGCGCCCGGCAATACCCCGACATCGGCCGCCTGCAGGCCACGCTCGCGGATGTGTCGCAAGGCACGCTGGCCGGCCTTGAAGGTGAggg                                                                                                                        <  8:261231/149‑1 (MQ=255)
                              cACTCGCCGAACAGCGCCAGGTCCAGGCCC‑TGGATGCCCAGCGGCTTTGGCCCACCCGCCGCGCCCGGCAATACCCCGACATCGGCCGCCTGCAGGCCACGCTCGCGGATGTGTCGCAAGGCACGCTGGCCGGCCTTGAAGGTGAGGGc                                                                                                                       >  1:339751/1‑149 (MQ=255)
                               aCTCGCCGAACAGCGCCAGGTCCAGGCCC‑TGGATGCCCAGCGGCTTTGGCCCACCCGCCGCGCCCGGCAATACCCCGACATCGGCCGCCTGCAGGCCACGCTCGCGGATGTGTCGCAAGGCACGCTGGCCGGCCTTGAAGGTGAGGGCg                                                                                                                      >  2:309250/1‑149 (MQ=255)
                                   gCCGAACAGCGCCAGGTCCAGGCCC‑TGGATGCCCAGCGGCTTTGGCCCACCCGCCGCGCCCGGCAATACCCCGACATCGGCCGCCTGCAGGCCACGCTCGCGGATGTGTCGCAAGGCACGCTGGCCGGCCTTGAAGGTGAGGGCggg                                                                                                                    <  2:169677/147‑1 (MQ=255)
                                   gCCGAACAGCGCCAGGTCCAGGCCC‑TGGATGCCCAGCGGCTTTGGCCCACCCGCCGCGCCCGGCAATACCCCGACATCGGCCGCCTGCAGGCCACGCTCGCGGATGTGTCGCAAGGCACGCTGGCCGGCCTTGAAGGTGAGGGCggg                                                                                                                    <  3:322729/147‑1 (MQ=255)
                                     cGAACAGCGCCAGGTCCAGGCCC‑TGGATGCCCAGCGGCTTTGGCCCACCCGCCGCGCCCGGCAATACCCCGACATCGGCCGCCTGCAGGCCACGCTCGCGGATGTGTCGCAAGGCACGCTGGCCGGCCTTGAAGGTGAGGGCGGGGTAt                                                                                                                <  3:53848/149‑1 (MQ=255)
                                       aaCAGCGCCAGGTCCAGGCCC‑TGGATGCCCAGCGGCTTTGGCCCACCCGCCGCGCCCGGCAATACCCCGACATCGGCCGCCTGCAGGCCACGCTCGCGGATGTGTCGCAAGGCACGCTGGCCGGCCTTGAAGGTGAGGGCGGGGTAttt                                                                                                              <  1:171765/149‑1 (MQ=255)
                                        aCAGCGCCAGGTCCAGGCCC‑TGGATGCCCAGCGGCTTTGGCCCACCCGCCGCGCCCGGCAATACCCCGACATCGGCCGCCTGCAGGCCACGCTCGCGGATGTGTCGCAAGGCACGCTGGCCGGCCTTGAAGGTGAGGGCGGGGTATTTg                                                                                                             <  1:309250/149‑1 (MQ=255)
                                           gcgcCAGGTCCAGGCCC‑TGGATGCCCAGCGGCTTTGGCCCACCCGCCGCGCCCGGCAATACCCCGACATCGGCCGCCTGCAGGCCACGCTCGCGGATGTGTCGCAAGGCACGCTGGCCGGCCTTGAAGGTGAGGGCggg                                                                                                                    >  2:257342/1‑139 (MQ=255)
                                           gcgcCAGGTCCAGGCCC‑TGGATGCCCAGCGGCTTTGGCCCACCCGCCGCGCCCGGCAATACCCCGACATCGGCCGCCTGCAGGCCACGCTCGCGGATGTGTCGCAAGGCACGCTGGCCGGCCTTGAAGGTGAGGGCggg                                                                                                                    <  1:257342/139‑1 (MQ=255)
                                            cgcCAGGTCCAGGCCC‑TGGATGCCCAGCGGCTTTGGCCCACCCGCCGCGCCCTGCAATACCCCGACATCGGCCGCCTGCAGGCCACGCTCGCGGATGTGTCGCAAGGCACGCTGGCCGGCCTTGAAGGTGAGGGCGGGGTATTTGAtgt                                                                                                         <  3:220053/149‑1 (MQ=255)
                                                          cc‑TGGATGCCCAGCGGCTTTGGCCCACCCGCCGCGCCCGGCAATACCCCGACATCGGCCGCCTGCAGGCCACGCTCGCGGATGTGTCGCAAGGCACGCTGGCCGGCCTTGAAGGTGAGGGCGGGGTATTTGATGTGGATGGCGCTCATg                                                                                           <  5:279764/149‑1 (MQ=255)
                                                          cc‑TGGATGCCCAGCGGCTTTGGCCCACCCGCCGCGCCCGGCAATACCCCGACATCGGCCGCCTGCAGGCCACGCTCGCGGATGTGTCGCAAGGCACGCTGGCCGGCCTTGAAGGTGAGGGCGGGGTATTTGATGTGGATGGCGCTCATg                                                                                           >  3:51595/1‑149 (MQ=255)
                                                          cc‑TGGATGCCCAGCGGCTTTGGCCCACCCGCCGCGCCCGGCAATACCCCGACATCGGCCGCCTGCAGGCCACGCTCGCGGATGTGTCGCAAGGCACGCTGGCCGGCCTTGAAGGTGAGGGCGGGGTATTTGATGTGGATGGCGCTCATg                                                                                           >  4:16681/1‑149 (MQ=255)
                                                               gATGCCCAGCGGCTTTGGCCCACCCGCCGCGCCCGGCAATACCCCGACATCGGCCGCCTGCAGGCCACGCTCGCGGATGTGTCGCAAGGCACGCTGGCCGGCCTTGAAGGTGAGGGCGGGGTATTTGATGTGGATGGCGCTCATGGc                                                                                         <  6:170507/147‑1 (MQ=255)
                                                                         ggCTTTGGCCCACCCGCCGCGCCCGGCAATACCCCGACATCGGCCGCCTGCAGGCCACGCGCGCGGATGTGTCGCAAGGCACGCTGGCCGGCCTTGAAGGTGAGGGCGGGGTATTTGATGTGGATGGCGCTCATGGCGGGTCTCTTggg                                                                             <  3:222898/149‑1 (MQ=255)
                                                                            tttGGCCCACCCGCCGCGCCCGGCAATACCCCGACATCGGCCGCCTGCAGGCCACGCTCGCGGATGTGTCGCAAGGCACGCTGGCCGGCCTTGAAGGTGAGGGCGGGGTATTTGATGTGGATGGCGCTCATGGCGGGTCTCTTggtgt                                                                           >  2:360917/1‑145 (MQ=255)
                                                                                 cccACCCGCCGCGCCCGGCAATACCCCGACATCGGCCGCCTGCAGGCCACGCGCGCGGATGTGTCGCAAGGCACGCTGGCCGGCCTTGAAGGTGAGGGCGGGGTATTTGATGTGGATGGCGCGCATGGCGGGTCTCTTGGGCATAGGag                                                                     <  7:175953/149‑1 (MQ=255)
                                                                                  ccACCCGCCGCGCCCGGCAATACCCCGACATCGGCCGCCTGCAGGCCACGCTCGCGGATGTGTCGCAAGGCACGCTGGCCGGCCTTGAAGGTg                                                                                                                            >  8:135254/1‑93 (MQ=255)
                                                                                  ccACCCGCCGCGCCCGGCAATACCCCGACATCGGCCGCCTGCAGGCCACGCTCGCGGATGTGTCGCAAGGCACGCTGGCCGGCCTTGAAGGTg                                                                                                                            <  7:135254/93‑1 (MQ=255)
                                                                                  ccACCCGCCGCGCCCGGCAATACCCCGACATCGGCCGCCTGCAGGCCACGCTCGCGGATGTGTCGCAAGGCACGCTGGCCGGCCTTGAAGGTGAGGGCGGGGTATTTGATGTGGATGGCGCTCATGGCGGGTCTCTTGGGCATAGGag                                                                     >  8:251015/1‑148 (MQ=255)
                                                                                        gccgcGCCCGGCAATACCCCGACATCGGCCGCCTGCAGGCCACGCTCGCGGATGTGTCGCAAGGCACGCTGGCCGGCCTTGAAGGTGAGGGCGGGGTATTTGATGTGGATGGCGCTCATGGCGGGTCTCTTGGGCATAGGAGCAGCCtt                                                              >  2:343832/1‑149 (MQ=255)
                                                                                         ccgcGCCCGGCAATACCCCGACATCGGCCGCCTGCAGGCCACGCTCGCGGATGTGTCGCTAGGCACGCTGGCCGGCCTTGAAGGTGAGGGCGGGGTATTTGATGTGGATGGCGCTCATGGCGGGTCTCTTGGGCATAGGAGCAGCCTtg                                                             <  7:80095/149‑1 (MQ=255)
                                                                                         ccgcGCCCGGCAATACCCCGACATCGGCCGCCTGCAGGCCACGCTCGCGGATGTGTCGCTAGGCACGCTGGCCGGCCTTGAAGGTGAGGGCGGGGTATTTGATGTGGATGGCGCTCATGGCGGGTCTCTTGGGCATAGGAGCAGCCTtg                                                             <  3:140394/149‑1 (MQ=255)
                                                                                                      tACCCCGACATCGGCCGCCTGCAGGCCACGCTCGCGGATGTGTCGCAAGGCACGCTGGCCGGCCTTGAAGGTGAGGGCGGGGTATTTGATGTGGATGGCGCTCATGGCGGGTCTCTTGGGCATAgg                                                                       >  7:45855/1‑126 (MQ=255)
                                                                                                      tACCCCGACATCGGCCGCCTGCAGGCCACGCTCGCGGATGTGTCGCAAGGCACGCTGGCCGGCCTTGAAGGTGAGGGCGGGGTATTTGATGTGGATGGCGCTCATGGCGGGTCTCTTGGGCATAgg                                                                       <  8:45855/126‑1 (MQ=255)
                                                                                                         cccGACATCGGCCGCCTGCAGGCCACGCTCGCGGATGTGTCGCAAGGCACGCTGGCCGGCCTTGAAGGTGAGGGCGGGGTATTTGATGTGGATGGCGCTCATGGCGGGTCTCTTGGGCATAGGAGCAGCCTTGTGCTGCGAgggg                                                 >  7:10866/1‑145 (MQ=255)
                                                                                                         cccGACATCGGCCGCCTGCAGGCCACGCTCGCGGATGTGTCGCAAGGCACGCTGGCCGGCCTTGAAGGTGAGGGCGGGGTATTTGATGTGGATGGCGCTCATGGCGGGTCTCTTGGGCATAGGAGCAGCCTTGTGCTGCGAgggg                                                 <  8:10866/145‑1 (MQ=255)
                                                                                                         cccGACATCGGCCGCCTGCAGGCCACGCTCGCGGATGTGTCGCAAGGCACGCTGGCCGGCCTTGAAGGTGAGGGCGGGGTATTTGATGTGGATGGCGCTCATGGCGGGTCTCTTGGGCATAGGAGCAGCCTTGTGCTGCGAGGGGGCAt                                             >  1:11802/1‑149 (MQ=255)
                                                                                                         cccGACATCGGCCGCCTGCAGGCCACGCTCGCGGATGTGTCGCAAGGCACGCTGGCCGGCCTTGAAGGTGAGGGCGGGGTATTTGATGTGGATGGCGCTCATGGCGGGTCTCTTGGGCATAGGAGCAGCCTTGTGCTGCGAGGGGCCAt                                             >  8:16038/1‑149 (MQ=255)
                                                                                                             aCATCGGCCGCCTGCAGGCCACGCTCGCGGATGTGTCGCAAGGCACGCTGGCCGGCCTTGAAGGTGAGGGCGGGGTATTTGATGTGGATGGCGCTCATGGCGGGTCTCTTGGGCATAGGAGCAGCCTTGTGCTGCGAGGGGCCATTACa                                         <  4:51595/149‑1 (MQ=255)
                                                                                                                tCGGCCGCCTGCAGGCCACGCTCGCGGATGTGTCGCAAGGCACGCTGGCCGGCCTTGAAGGTGAGGGCGGGGTATTTGATGTGGATGGCGCTCATGGCGGGTCTCtt                                                                                >  2:95610/1‑107 (MQ=255)
                                                                                                                tCGGCCGCCTGCAGGCCACGCTCGCGGATGTGTCGCAAGGCACGCTGGCCGGCCTTGAAGGTGAGGGCGGGGTATTTGATGTGGATGGCGCTCATGGCGGGTCTCtt                                                                                <  1:95610/107‑1 (MQ=255)
                                                                                                                      gccTGCAGGCCACGCTCGCGGATGTGTCGCAAGGCACGCTGGCCGGCCTTGAAGGTGAGGGCGGGGTATTTGATGTGGATGGCGCTCATGGCGGGTCTCTTGGGCATAGGAGCAGCCTTGTGCTGCGAGGGGCCATTACAGATAacaca                                >  7:187856/1‑149 (MQ=255)
                                                                                                                         tGCAGGCCACGCTCGCGGATGTGTCGCAAGGCACGCTGGCCGGCCTTGAAGGTGAGGGCGGGGTATTTGATGTGGATGGCGCTCATGGCGGGTCTCTTGGGCATAGGAGCAGCCTTGTGCTGCGAGGGGCCATTACAGATAACACACCt                             >  2:5725/1‑149 (MQ=255)
                                                                                                                               ccACGCTCGCGGATGTGTCGCAAGGCACGCTGGCCGGCCTTGAAGGTGAGGGCGGGGTATTTGATGTGGATGGCGCTCATGGCGGGTCTCTTGGGCATAGGAGCAGCCTTGTGCTGCGAGGGGCCATTACAGATAACACACCTGTAGGt                       <  2:339751/149‑1 (MQ=255)
                                                                                                                                cACGCTCGCGGATGTGTCGCAAGGCACGCTGGCCGGCCTTGAAGGTGAGGGCGGGGTATTTGATGTGGATGGCGCTCATGGCGGGTCTCTTGGGCATAGGAGCAGCCTTGTGCTGCGAgg                                                   >  4:325457/1‑120 (MQ=255)
                                                                                                                                cACGCTCGCGGATGTGTCGCAAGGCACGCTGGCCGGCCTTGAAGGTGAGGGCGGGGTATTTGATGTGGATGGCGCTCATGGCGGGTCTCTTGGGCATAGGAGCAGCCTTGTGCTGCGAgg                                                   <  3:325457/120‑1 (MQ=255)
                                                                                                                                       gcgGATGTGTCGCAAGGCACGCTGGCCGGCCTTGAAGGTGAGGGCGGGGTATTTGATGTGGATGGCGCTCATGGCGGGTCTCTTGGGCATAGGAGCAGCCTTGTGCTGCGAGGGGCCATTACAGATAACACACCTGTAGGTGGCTGt                 <  1:5725/147‑1 (MQ=255)
                                                                                                                                                  gCAAGGCACGCTGGCCGGCCTTGAAGGTGAGGGCGGGGTATTTGATGTGGATGGCGCTCATGGCGGGTCTCTTGGGCATAGGAGCAGCCTTGTGCTGCGAGGGGCCATTACAGATAACACACCTGTAGGTGGCTGTATTGGCCCTTTCg    <  7:16038/149‑1 (MQ=255)
                                                                                                                                                  gCAAGGCACGCTGGCCGGCCTTGAAGGTGAGGGCGGGGTATTTGATGTGGATGGCGCTCATGGCGGGTCTCGTGGGCATAGGAGCAGCCTTGTGCTGCGAGGGGCGATTACAGATAACACACCTGTAGGTGGCTGTATTTGCCCTTTCg    >  3:265095/1‑149 (MQ=255)
                                                                                                                                                   cAAGGCACGCTGGCCGGCCTTGAAGGTGAGGGCGGGGTATTTGATGTGGATGGCGCTCATGGCGGGTCTCTTGGGCATAGGAGCAGCCTTGTGCTGCGAGGGGCCATTACAGATAACACACCTGTAGGTGGCTGTATTGGCCCTTTCg    >  5:266820/1‑148 (MQ=255)
                                                                                                                                                    aaGGCACGCTGGCCGGCCTTGAAGGTGAGGGCGGGGTATTTGATGTGGATGGCGCTCATGGCGGGTCTCTTGGGCATAGGAGCAGCCTTGTGCTGCGAGGGGCCATTACAGATAACACACCTGTAGGTGGCTGTATTGGCCCTTTCgca  <  1:343832/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AGGGCACGCTGACGGGGTGCCGAGGGTAGCCACTCGCCGAACAGCGCCAGGTCCAGGCCC‑TGGATGCCCAGCGGCTTTGGCCCACCCGCCGCGCCCGGCAATACCCCGACATCGGCCGCCTGCAGGCCACGCTCGCGGATGTGTCGCAAGGCACGCTGGCCGGCCTTGAAGGTGAGGGCGGGGTATTTGATGTGGATGGCGCTCATGGCGGGTCTCTTGGGCATAGGAGCAGCCTTGTGCTGCGAGGGGCCATTACAGATAACACACCTGTAGGTGGCTGTATTGGCCCTTTCGCA  >  NC_002947/2675929‑2676224

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: ‑