Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 991,560 T→C 8.6% I60T (ATC→ACC)  der → GTPase Der

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947991,5600TC 8.6% 123.7 / 3.1 47I60T (ATC→ACC) derGTPase Der
Reads supporting (aligned to +/- strand):  ref base T (21/21);  new base C (2/2);  total (23/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.46e-01

CCGAACGTCGGCAAATCCACCATGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTAT‑GTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAGATGATCGCCGAGCACCTGCGCAAGCGG  >  NC_002947/991412‑991708
                                                                                                                                                    |                                                                                                                                                     
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                                                                                                                                       >  1:39159/1‑150 (MQ=255)
               tCCACCATGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTg                                                                                                                                       <  2:610543/150‑1 (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                                                                                                                         >  1:424853/1‑135 (MQ=255)
                                            ccGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGaaaa                                                                                                                         <  2:424853/135‑1 (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:681409/1‑150 (MQ=255)
                                                                                ccGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTAT‑GTGCTGtt                                                                     >  1:572474/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)
                                                                                                                                                   aTCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTAT‑GTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAGATGATCGCCGAGCACCTGCGCAAGCgg  <  1:1088040/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCGAACGTCGGCAAATCCACCATGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTAT‑GTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAGATGATCGCCGAGCACCTGCGCAAGCGG  >  NC_002947/991412‑991708

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