Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 547,203 A→G 12.9% H18R (CAC→CGC)  fusA → elongation factor G 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947547,2030AG12.9% 134.5 / 11.1 54H18R (CAC→CGC) fusAelongation factor G 1
Reads supporting (aligned to +/- strand):  ref base A (19/28);  new base G (3/4);  total (22/32)
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.46e-01

AGAAGCGTGAAGACGTTCACCGTATGGCTGAAGCCAACAAAGCGTTCTCGCACTACCGCTTCTAATTCAAGCATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAATGGGCGAGGTGCATGACGGCGCCGCGACCACCGACTGGATGGTGCAGGAGCAGGAGCGGGGTATCACCATTACCTCCGCTG  >  NC_002947/547056‑547349
                                                                                                                                                   |                                                                                                                                                   
agaagCGTGAAGACGTTCACCGTATGGCTGAAGCCAACAAAGCGTTCTCGCACTATCGCTTCTATTTCAAGCATCAATCTTTTTGCGTGGGCTTTATGGCTCGTACTACAGGATTTTAGAGCTACAGTAACATTGTTATCTGTGCCCAc                                                                                                                                                    <  6:214127/149‑1 (MQ=255)
         aaGACGTTCACCGTATGGCTGAAGCCAACAAAGCGTTCTCGCACTACCGCTTCTAATTCAAGCATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGAcgcg                                                                                                                                           >  3:136704/1‑149 (MQ=255)
         aaGACGTTCACCGTATGGCTGAAGCCAACAAAGCGTTCTCGCACTACCGCTTCTAATTCAAGCATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGAcgc                                                                                                                                            <  3:119253/148‑1 (MQ=255)
               ttCAGCGTATGGCTGAAGCCTACAAAGCGTTCTCGCACTACCGCTTCTAATTCAAGCAGCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCGCGTCAACGTGGGCaca                                                                                                                                     <  3:235865/149‑3 (MQ=255)
               ttCACCGTATGGCTGAAGCCAACAAAGCGTTCTCGCACTACCGCTTCTAATTCAAGCATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCGCGTCAACGTGGGCac                                                                                                                                      >  6:211706/1‑147 (MQ=255)
               ttCACCGTATGGCTGAAGCCAACAAAGCGTTCTCGCACTACCGCTTCTAATTCAAGCATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCGCGTCAACGTGGGCac                                                                                                                                      <  5:211706/148‑2 (MQ=255)
                 cACCGTATGGCTGAAGCCAACAAAGCGTTCTCGCACTACCGCTTCTAATTCAAGCATCAATCTTTTTGCGAGGGCTTTACGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGa                                                                                                                                    <  6:60677/148‑1 (MQ=255)
                                gCCAACAAAGCGTTCTCGCACTACCGCTTCTAATTCAAGCATCAATCTTTTTGCGAGGTCTTTATTGCTAGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTTCCCACTTTGAAGCTGGCAAGACCACAACTACCGAGCg                                                                                                                    >  6:85173/1‑149 (MQ=255)
                                      aaaGCGTTCTCGCACTACCGCTTCTAATTCAAGCATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCGCGTCAACGTGGGCacag                                                                                                                                    <  6:12632/127‑4 (MQ=255)
                                       aaGCGTTCTCGCACTACCGCTTCTAATTCAAGCATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCGCGTCAACGTGGGCaca                                                                                                                                     <  4:96535/125‑3 (MQ=255)
                                       aaGCGTTCTCGCACTACCGCTTCTAATTCAAGCATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCGCGTCAACGTGGGCaca                                                                                                                                     >  3:96535/1‑123 (MQ=255)
                                         gCGTTCTCGCACTACCGCTTCTAATTCAAGCATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGtt                                                                                                           <  4:122132/149‑1 (MQ=255)
                                           gTTCTCGCACTACCGCTTCTAATTCAAGCATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGc                                                                                                                     >  8:84740/1‑137 (MQ=255)
                                           gTTCTCGCACTACCGCTTCTAATTCAAGCATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGc                                                                                                                     <  7:84740/137‑1 (MQ=255)
                                           gTTCTCGCACTACCGCTTCTAATTCAAGCATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCt                                                                                                         <  1:32015/149‑1 (MQ=255)
                                                  cACTACCGCTTCTAATTCAAGCATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAgg                                                                                                  >  6:46119/1‑149 (MQ=255)
                                                   aCTACCGCTTCTAATTCAAGCATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGt                                                                                                            <  7:10756/138‑1 (MQ=255)
                                                   aCTACCGCTTCTAATTCAAGCATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGt                                                                                                            >  8:10756/1‑138 (MQ=255)
                                                       ccGCTTCTAATTCAGGCATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGa                                                                                             <  1:130944/149‑1 (MQ=255)
                                                                 ttCAAGCATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGc                                                                                                                                        <  1:63866/96‑1 (MQ=255)
                                                                 ttCAAGCATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGc                                                                                                                                        >  2:63866/1‑96 (MQ=255)
                                                                 ttCAAGCATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAAt                                                                                   >  6:10516/1‑149 (MQ=255)
                                                                    aaGCATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAATgg                                                                                 <  8:47317/148‑1 (MQ=255)
                                                                     aGCATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAATGGGc                                                                               <  2:146785/149‑1 (MQ=255)
                                                                     aGCATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAATGGGc                                                                               >  8:72453/1‑149 (MQ=255)
                                                                       cATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAATGGGCGa                                                                             <  7:72453/149‑1 (MQ=255)
                                                                        aTCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAGCATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAATGGGCGAg                                                                            >  6:177296/1‑149 (MQ=255)
                                                                        aTCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGtt                                                                                                           >  8:61689/1‑118 (MQ=255)
                                                                        aTCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGtt                                                                                                           <  7:61689/118‑1 (MQ=255)
                                                                        aTCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAATGGGCGAg                                                                            >  1:171772/1‑149 (MQ=255)
                                                                                 tttGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAATGGGCGAGGTGCAt                                                                     <  4:136704/146‑1 (MQ=255)
                                                                                              tATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAATGGGAGAGGTGCATGACGGCGCCGCGacc                                                      <  5:46119/148‑1 (MQ=255)
                                                                                               aTGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAATGGGCGAGGTGCATGACGGCGCCGCGaccac                                                    >  3:61240/1‑149 (MQ=255)
                                                                                                tGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAATGGGCGAGGTGCATGACGGCGCCGCGaccac                                                    <  2:7819/148‑1 (MQ=255)
                                                                                                    tCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAATGGGCGAGGt                                                                          >  8:127579/1‑123 (MQ=255)
                                                                                                    tCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAATGGGCGAGGt                                                                          <  7:127579/123‑1 (MQ=255)
                                                                                                            cAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGac                                                                                                                                   >  2:81378/1‑58 (MQ=255)
                                                                                                            cAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGac                                                                                                                                   <  1:81378/58‑1 (MQ=255)
                                                                                                              gCAATTAACCGCTACCGTAGCATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAATGGGCGAGGTGCATGACGGCGCCGCGACCACCGACTGGATGGTgca                                     <  5:177296/149‑1 (MQ=255)
                                                                                                                aaTTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAATGGGCGAGGTGCATGACGGCGCCGCGACCACCGACTGGATGGTgcagg                                   <  5:10516/149‑1 (MQ=255)
                                                                                                                aaTTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAATGGGCGAGGTGCATGACGGCGCCGCGACCACCGACTGGATGGTgcag                                    >  1:4990/1‑148 (MQ=255)
                                                                                                                     aCCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAATGGGCGAGGTGCATGACGGCGCCGCGaccacc                                                   <  7:6820/128‑1 (MQ=255)
                                                                                                                     aCCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAATGGGCGAGGTGCATGACGGCGCCGCGaccacc                                                   >  8:6820/1‑128 (MQ=255)
                                                                                                                            ccGTAACATTGGTATCTGTGCCCAGGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAATGGGCGAGGTGCATGACGGCGCCGCGACCACCGACTGGATGGTGCAGGAGCAGGAGCggg                       <  5:85173/149‑1 (MQ=255)
                                                                                                                            ccGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAATGGGCGAGGTGCATGACGGCGCCGCGACCACCGACTGGATGGTGCAGGAGCAGGAGCggg                       <  2:171772/149‑1 (MQ=255)
                                                                                                                                 aCATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAATGGGCGAGGTGCCTTGACGGCGCCGCGACCACCGACTGGATGGTGCAGGAGCAGGAGCGGGGTAt                   <  1:202910/149‑1 (MQ=255)
                                                                                                                                   aTTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAATGGGCGAGGTGCATGACGGCGCCGCGACCACCGACTGGATGGTGCAGGAGCAGGAGCGGGGTATCAc                <  4:61240/149‑1 (MQ=255)
                                                                                                                                       gTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAATGGGCGAGGTGCATGACGGCGCcgcg                                                         >  5:115380/1‑104 (MQ=255)
                                                                                                                                       gTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAATGGGCGAGGTGCATGACGGCGCcgcg                                                         <  6:115380/104‑1 (MQ=255)
                                                                                                                                       gTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAATGGGCGAGGTGCATGACGGCGCCGCGACCACCGACTGGa                                            <  1:267647/117‑1 (MQ=255)
                                                                                                                                       gTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAATGGGCGAGGTGCATGACGGCGCCGCGACCACCGACTGGa                                            >  2:267647/1‑117 (MQ=255)
                                                                                                                                           tcttgCCCGCGTCAACGTGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAATGGGCGAGGTGCATGACGGCGCCGCGACCACCGACTGGATGGTGCAGGAGCAGGAGCGGGGTATCACCATTACCt        >  4:118971/4‑149 (MQ=25)
                                                                                                                                                 ccACGTTGACGCGGGCAAGACCACCACTGCCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAATGGGCGAGGTGCATGACGGCGCCGCGACCACCGGCTGGATGGTGCAGGAGCAGGAGCGGGGTATCACCATTACCTCCGCTg  >  6:107862/1‑149 (MQ=255)
                                                                                                                                                 ccACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAATGGGCGAGGTGCATGACGGCGCCGCGACCACCGACTGGATGGTgcagga                                  <  6:147457/117‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AGAAGCGTGAAGACGTTCACCGTATGGCTGAAGCCAACAAAGCGTTCTCGCACTACCGCTTCTAATTCAAGCATCAATCTTTTTGCGAGGGCTTTATGGCTCGTACTACAGCAATTAACCGCTACCGTAACATTGGTATCTGTGCCCACGTTGACGCGGGCAAGACCACCACTACCGAGCGGATCCTGTTCTACACAGGTCTGAGCCACAAAATGGGCGAGGTGCATGACGGCGCCGCGACCACCGACTGGATGGTGCAGGAGCAGGAGCGGGGTATCACCATTACCTCCGCTG  >  NC_002947/547056‑547349

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: