Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 991,557 G→A 8.7% G59D (GGT→GAT)  der → GTPase Der

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947991,5570GA 8.7% 126.0 / 2.8 46G59D (GGT→GAT) derGTPase Der
Reads supporting (aligned to +/- strand):  ref base G (21/21);  new base A (2/2);  total (23/23)
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.36e-01

CCCGAACGTCGGCAAATCCACCATGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTAT‑GTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAGATGATCGCCGAGCACCTGCGCAA  >  NC_002947/991411‑991704
                                                                                                                                                  |                                                                                                                                                    
cccGAACGTCGGCAAATCCACCATGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTac                                                                                                                                                   >  2:1302/1‑149 (MQ=255)
 ccGAACGTCGGCAAATCCACCATGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATc                                                                                                                                                  >  1:60996/1‑150 (MQ=255)
         cGGCAAATCCACCATGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTgt                                                                                                                                          <  1:395483/150‑2 (MQ=255)
            cAAATCCACCATGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGTCACCGGTGATACCACCGGTGTCGa                                                                                                                                       <  2:296619/150‑1 (MQ=255)
               aTCCACCATGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCAGCCTGATCGACACCGGTGGTATCACCGGTGACGAAGt                                                                                                                                    <  2:79629/150‑1 (MQ=255)
                tCCACCATGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTg                                                                                                                                   <  2:610543/150‑1 (MQ=255)
                tCCACCATGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTg                                                                                                                                   >  1:39159/1‑150 (MQ=255)
                 ccaccaTGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTgg                                                                                                                                  <  1:203161/150‑1 (MQ=255)
                      aTGGTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATg                                                                                                                             >  2:836351/1‑150 (MQ=255)
                                   tGACCAAGACCCGCGATGCCATCGTCGGTTACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATgg                                                                                                                >  1:871065/1‑150 (MQ=255)
                                   tGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATgg                                                                                                                >  2:145119/1‑150 (MQ=255)
                                      ccAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCg                                                                                                             <  2:39159/150‑1 (MQ=255)
                                        aaGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAg                                                                                                           <  1:514079/150‑1 (MQ=255)
                                          gACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCa                                                                                                         <  1:147668/150‑1 (MQ=255)
                                             ccGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGaaaa                                                                                                                     <  2:424853/135‑1 (MQ=255)
                                             ccGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGaaaa                                                                                                                     >  1:424853/1‑135 (MQ=255)
                                               gcgATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGc                                                                                                    >  1:799166/1‑150 (MQ=255)
                                                       aTCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAACATGGTGATGCCAGCAGGCAGCGACGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTCGGCATGGACGAAAAAATGGACGAGCAGTCGCTCATGGcc                                                                                            >  2:240117/1‑150 (MQ=255)
                                                        tcgtcgGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGGAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCa                                                                                           >  1:509498/1‑150 (MQ=255)
                                                                          gCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTc                                                                                                      <  2:49324/121‑1 (MQ=255)
                                                                          gCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTc                                                                                                      >  1:49324/1‑121 (MQ=255)
                                                                          gCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTAt‑g                                                                        >  1:199234/1‑150 (MQ=255)
                                                                              gACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTAT‑GTGCt                                                                    <  2:721645/150‑1 (MQ=255)
                                                                                 ccGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTAT‑GTGCTGtt                                                                 >  1:572474/1‑150 (MQ=255)
                                                                                 ccGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTAT‑GTGCTGtt                                                                 >  1:681409/1‑150 (MQ=255)
                                                                                        cGCAAGTATGGTGATGCCAGCTGGCAGAGTTGTTCCTTCACCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGNGCAGNCGCTCATGGCCATTGAAGAAGCCGACTAT‑GTGCTGTTCCTGGtt                                                          <  1:1302/150‑1 (MQ=255)
                                                                                                       gCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTAT‑GTGCTGTTCCTGGTTGATGCCCGTGCGGGt                                           <  1:162257/150‑1 (MQ=255)
                                                                                                        ccAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTAT‑GTGCTGTTCCTGGTTGATGCCCGTGCGGGTa                                          <  1:459886/150‑1 (MQ=255)
                                                                                                         cAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTgg                                                                                                                                  >  2:468711/1‑62 (MQ=255)
                                                                                                         cAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTgg                                                                                                                                  <  1:468711/62‑1 (MQ=255)
                                                                                                               gCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTAT‑GTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGc                                   <  2:871065/150‑1 (MQ=255)
                                                                                                                    gTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCAt                                                                                                                              <  1:309978/55‑1 (MQ=255)
                                                                                                                    gTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCAt                                                                                                                              >  2:309978/1‑55 (MQ=255)
                                                                                                                     tCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTAT‑GTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGa                             >  1:197798/1‑150 (MQ=255)
                                                                                                                      cGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTAT‑GTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGAc                            <  2:681409/150‑1 (MQ=255)
                                                                                                                        ttCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCCTGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTAC‑GTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCa                          <  2:572474/150‑1 (MQ=255)
                                                                                                                           cTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTAT‑GTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAgat                       >  2:979489/1‑150 (MQ=255)
                                                                                                                            ttCATCCTGATCGAAACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGGAGCCGACTAT‑GTGCTGTTCCTGGTTGATGCAAGTGCGGGTATGACTGCGGCAGAACAgatg                      >  2:994285/1‑150 (MQ=255)
                                                                                                                                tCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTAT‑GTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAGATGATCg                  <  2:197798/150‑1 (MQ=255)
                                                                                                                                   tGATCGACACCGGTGGTATCACCGGTGGAGCAGTGAGGCTGGACGAAAAAAGGTCCGAGCCGTCGCTCACTGACATCGAAGAAG‑AGAAAATCGTTCTGTTCCTGGTTGATGCCCGTGCGGGTATGAATGCGGCAGCCCGAATGATCGCCg               >  1:810352/1‑150 (MQ=255)
                                                                                                                                        gtcacCGGTGATACCACCGGTGTCGAAGTGGGCATGGACGAAAAAATGGCCGAg                                                                                                           <  2:447182/52‑1 (MQ=18)
                                                                                                                                        gtcacCGGTGATACCACCGGTGTCGAAGTGGGCATGGACGAAAAAATGGCCGAg                                                                                                           >  1:447182/3‑54 (MQ=18)
                                                                                                                                         tcacCGGTGATACCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTAT‑GTGCTGTTCCTGGTTGATGCCCGTGCGGGTCTGACTGCGGCAGACCAGATGATCGCCGAGCAcc         >  1:555151/2‑150 (MQ=255)
                                                                                                                                         acacCGGTGGTAGCACCGGTGACGAAGTGGTCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTAT‑GTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAGATGATCGCCGAGCAcc         <  1:240117/150‑1 (MQ=255)
                                                                                                                                           acCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTAT‑GTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAGATGATCGCCGAGCACCTg       <  2:810352/150‑1 (MQ=255)
                                                                                                                                             cGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTAT‑GTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAGATGATCGCCGCGAACCTgcg     >  2:139397/1‑150 (MQ=255)
                                                                                                                                                tggtATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTAT‑GTGCTGTTCCTGGTTGATGCCCCAGCGGGTAGGACTGCGGCATACCAGATGATCGCCGAGCCCCTgaacga  >  2:1030367/1‑145 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CCCGAACGTCGGCAAATCCACCATGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTAT‑GTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAGATGATCGCCGAGCACCTGCGCAA  >  NC_002947/991411‑991704

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