Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,164,022 C→T 7.7% R651H (CGT→CAT)  parC ← DNA topoisomerase IV, subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,164,0220CT 7.7% 153.9 / 3.1 52R651H (CGT→CAT) parCDNA topoisomerase IV, subunit A
Reads supporting (aligned to +/- strand):  ref base C (23/25);  new base T (2/2);  total (25/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATGGTCAGCGTGCTTTGCGGCGGCAGAACGTACAATTGCGCCAGACCATCTTCTCCACGCGCGGCTTCTGCCGATGGAATGTTGATAATCTTGTTGCCTTTGCCCTTCGACAGCTGCGGCAGATCACTTACCGGGAACATCAACATACGGCCTGCCTGAGTGATTGCCAGCAGCATATCGGAAGCATCTTCAATCACCACCGGCGGCATAACATGGGCATTTTCCGGTAAGGTGATCAAAGCCTTACCTGCACGGTTACGCGCCACCAGATCGTTAAAGGTGCAGAC  >  NC_000913/3163874‑3164161
                                                                                                                                                    |                                                                                                                                           
aaTGGTCAGCGTGCTTTGCGGCGGCAGAACGTACAATTGCGCCAGACCATCTTCTCCACGCGCGGCTTCTGCCGATGGAATGTTGATAATCTTGTTGCCTTTGCCCTTCGACAGCTGCGGCAGATCACTTACCGGGAACATCAACATAc                                                                                                                                             <  8:135705/149‑1 (MQ=255)
            gCTTTGCGGCGGCAGAACGTACAATTGCGCCAGACCATCTTCTCCACGGGCGGCTTCTGCCGATGGAATGTTGATAATCTTGTTGCCTTTGCCCTTCGACAGCTGCGGCAGATCACTTACCGGGAACATCAACATACGGCCTGCCTGAg                                                                                                                                 <  8:265266/149‑1 (MQ=255)
             cTTTGCGGCGGCAGAACGTACAATTGCGCCAGACCATCTTCTCCACGCGCGGCTTCTGCCGATGGAATGTTGATAATCTTGTTGCCTTTGCCCTTCGACAGCTGCGGCAGATCACTTACCGGGAACATCAACATACGGCCTGCCTGAg                                                                                                                                 <  3:236341/148‑1 (MQ=255)
             cTTGGCGGCGGCAGAACGTACAATTGCGCCAGACCATCTTCTCCACGCGCGGCTTCTGCCGATGGAATGTTGATAATCTTGTTGCCTTTGCCCTTCGACAGCTGCGGCAGATCACTTACCGGGAACATCAACATACGGCCTGCCTGAGt                                                                                                                                <  2:327427/149‑1 (MQ=255)
                      ggcAGAACGTACAATTGCGCCAGACCATCTTCTCCACGCGCGGCTTCTGCCGATGGAATGTTGATAATCTTGTTGCCTTTGCCCTTCGACAGCTGCGGCAGATCACTTACCGGGAACATCAACATACGGCCTGCCTGAGTGATTGCcag                                                                                                                       <  2:204962/149‑1 (MQ=255)
                             cGTACAATTGCGACAGACAATCTTCTCCACGCGCGGCTTCTGCCGATGGAATGTTGAAAATCTTGTTGCCTTTGCCCTTCGACAGCTTCGGCAGATCACTTACAGGGAACATCAACATACgg                                                                                                                                           >  5:88893/1‑122 (MQ=255)
                                   aTTGCGCCAGACGATCTTCTCCACGCGCGGCTTCTGCCGATGGAATGTTGATAATCTTGTTGCCTTTGCCCTTCGACAGCTGCGGCAGATCACTTACCGGGAACATCAACATACGGCCTGCCTGAGTGATTGCCAGCAGCATATCGGaa                                                                                                          <  1:122345/149‑1 (MQ=255)
                                    ttGCGCCAGACCATCTTCTCCACGCGCGGCTTCTGCCGATGGAATGTTGATAATCTTGTTGCCTTTGCCCTTCGACAGCTGCGGCAGATCACTTACCGGGAACATCAACATACGGCCTGCCTGAGTGATTGCCAGCAGCATATCGGAAg                                                                                                         <  7:328679/149‑1 (MQ=255)
                                         ccAGACCATCTTCTCCACGCGCGGCTTCTGCCGATGGAATGTTGATAATCTTGTTGCCTTTGCCCTTCGACAGCTGCGGCAGATCACTTACCGGGAACATCAACATACTGCCTGCCTGAGTGATTGCCAGCAGCATATCGGAAGCATCt                                                                                                    <  8:347553/149‑1 (MQ=255)
                                          cAGACCATCTTCTCCACGCGCGGCTTCTGCCGATGGAATGTTGATAATCTTGTTGCCTTTGCCCTTCGACAGCTGCGGCAGATCACTTACCGGGAACATCAACATACGGCCTGCCTGAGTGATTGCCAGCAGCATATCGGAAGCATCtt                                                                                                   <  1:110386/149‑1 (MQ=255)
                                               cATCTTCTCCACGCGCGGCTTCTGCCGATGGAATGTTGATAATCTTGTTGCCTTTGCCCTTCGACAGCTGCGGCAGATCACTTACCGGGAACATCAACATACGGCCTGCCTGAGTGATTGCCAGCAGCATATCGGAAGCATCTTCAATc                                                                                              >  5:44594/1‑149 (MQ=255)
                                                        cACGCGCGGCTTCTGCCGATGGAATGTTGATAATCTTGTTGCCTTTGCCCTTCGACAGCTGCGGCAGATCACTTACCGGGAACATCAACATATGGCCTGCCTGAGTGATTGCCAGCAGCATATCGGAAGCATCTTCAATCACCACcggc                                                                                     >  2:100974/1‑149 (MQ=255)
                                                        cACGCGCGGCTTCTGCCGATGGAATGTTGATAATCTTGTTGCCTTTGCCCTTCGACAGCTGCGGCAGATCACTTACCGGGAACATCAACATATGGCCTGCCTGAGTGATTGCCAGCAGCATATCGGAAGCATCTTCAATCACCACcggc                                                                                     >  2:106005/1‑149 (MQ=255)
                                                        cACGCGCGGCTTCTGCCGATGGAATGTTGATAATCTTGTTGCCTTTGCCCTTCGACAGCTGCGGCAGATCACTTACCGGGAACATCAACATACGGCCTGCCTGAGTGATTGCCAGCAGCATATCGGAAGCATCTTCAATCACCACcggc                                                                                     >  7:91953/1‑149 (MQ=255)
                                                        cACGCGCGGCTTCTGCCGATGGAATGTTGATAATCTTGTTGCCTTTGCCCTTCGACAGCTGCGGCAGATCACTTACCGGGAACATCAACATACGGCCTGCCTGAGTGATTGCCAGCAGCATATCGGAAGCATCTTCAATCACCACcggc                                                                                     >  3:145730/1‑149 (MQ=255)
                                                                             gAATGTTGATAATCTTGTTGCCTTTGCCCTTCGACAGCTGCGGCAGATCACTTACCGGGAACATCAACATACGGCCTGCCTGAGTGATTGCCAGCAGCATATCGGAAGCATCTTCAATCACCACCGGCGGCATAACATGGGCATTTTc                                                                 <  8:91953/148‑1 (MQ=255)
                                                                             gAATGTTGATAATCTTGTTGCCTTTGCCCTTCGACAGCTGCGGCAGATCACTTACCGGGAACATCAACATACGGCCTGCCTGAGTGATTGCCAGCAGCATATCGGAAGCATCTTCAATCACCACCGGCGGCATAACATGGGCATTTTc                                                                 <  6:44594/148‑1 (MQ=255)
                                                                               aTGTTGATAATCTTGTTGCCTTTGCCCTTCGACAGCTGCGGCAGATCACTTACCGGGAACATCAACATACGGCCTGCCTGAGTGATTGCCAGCAGCATATCGGAAGCATCTTCAATCACCACCGGCGGCATAACATGGGCATTTTCCgg                                                              <  4:145730/149‑1 (MQ=255)
                                                                                  ttGATAATCTTGTTGCCTTTGCCCTTCGACAGCTGCGGCAGATCACTTACCGGGAACATCAACATACGGCCTGCCTGAGTGATTGCCAGCAGCATATCGGAAg                                                                                                         <  3:86256/103‑1 (MQ=255)
                                                                                  ttGATAATCTTGTTGCCTTTGCCCTTCGACAGCTGCGGCAGATCACTTACCGGGAACATCAACATACGGCCTGCCTGAGTGATTGCCAGCAGCATATCGGAAg                                                                                                         >  4:86256/1‑103 (MQ=255)
                                                                                    gATAATCTTGTTGCCTTTGCCCTTCGACAGCTGCGGCAGATCACTTACCGGGAACATCAACATACGGCCTGCCTGAGTGATTGCCAGCAGCATATCGGAAGCATCTTCAATCACCACCGGCGGCATAACATGGGCATTTTCCGGTAAgg                                                         >  6:87899/1‑149 (MQ=255)
                                                                                      tAATCTTGTTGCCTTTGCCCTTCGACAGCTGCGGCAGATCACTTACCGGGAACATCAACATACGGCCTGCCTGAGTGATTGCCAGCAGCATATCGGaa                                                                                                          >  3:36150/1‑98 (MQ=255)
                                                                                      tAATCTTGTTGCCTTTGCCCTTCGACAGCTGCGGCAGATCACTTACCGGGAACATCAACATACGGCCTGCCTGAGTGATTGCCAGCAGCATATCGGaa                                                                                                          <  4:36150/98‑1 (MQ=255)
                                                                                         tCTTGTTGCCTTTGCCCTTCGACAGCTGCGGCAGATCACTTACCGGAAACATCAACATACGGCCTGCCTGAGTGATTGCCAGCAGCATATCGGAAGCATCTTCAATc                                                                                              <  6:116168/107‑1 (MQ=255)
                                                                                                gCCTTTGCCCTTCGACAGCTGCGGCAGATCACTTACCGGGAACATCAACATACGGCCTGCCTGAGTGATTGCCAGCAGCATATCGGAAGCATCTTCAATCACCACCGGCGGCATAACATGGGCATTTTCCGGTAAg                                                          >  3:303713/1‑136 (MQ=255)
                                                                                                gCCTTTGCACTTCGTCTGCTGCGGCAGATCACTTACCGGGAACATCAACATACTGCCTGCCTGATTGATTGCCAGCATCATATCGGTAGCATCTTCAATCACCAACGGCGGCATAAGATGGTCATTTTCCGGTAAg                                                          <  4:303713/136‑1 (MQ=255)
                                                                                                  cTTTGCCCTTCGACAGCTGTGGCAGATCACTTACCGGGAACATCAACATACGGCCTGCCTGAGTGATTGCCAGCAGCATATCGGAAGCATCTTCAATCACCACCGGCGGCATAACATGGGCATTTTc                                                                 >  2:56367/1‑127 (MQ=255)
                                                                                                  cTTTGCCCTTCGACAGCTGTGGCAGATCACTTACCGGGAACATCAACATACGGCCTGCCTGAGTGATTGCCAGCAGCATATCGGAAGCATCTTCAATCACCACCGGCGGCATAACATGGGCATTTTc                                                                 <  1:56367/127‑1 (MQ=255)
                                                                                                  cTTTGCCCTTCGACAGCTGCGGCAGATCACTTACCGGGAACATCAACATACGGCCTGCCTGAGTGATTGCCAGCAGCATATCGGAAGCATCTTCAATc                                                                                              <  2:210190/98‑1 (MQ=255)
                                                                                                  cTTTGCCCTTCGACAGCTGCGGCAGATCACTTACCGGGAACATCAACATACGGCCTGCCTGAGTGATTGCCAGCAGCATATCGGAAGCATCTTCAATc                                                                                              >  1:210190/1‑98 (MQ=255)
                                                                                                   tttGCCCTTCGACAGCTGCGGCAGATCACTTACCGGGAACATCAACATACGGCCTGCCTGAGTGATTGCCAGCAGCATATCGGAAGCATCTTCAATCACCACCGGCGGCATAACATGGGCATTTTCCGGTAAGGTGATCAAAGCCTTAc                                          >  4:263810/1‑149 (MQ=255)
                                                                                                      gCCCTTCGACAGCTGCGGCAGATCACTTACCGGGAACATCAACATACGGCCTGCCTGAGTGATTGCCAGCAGCATATCGGAAGCATCTTCAATCACCACCGGCGGCATAACATGGGCATTTTCCGGTAAGGTGATCAAAGCCTTACCt                                        <  3:263810/148‑1 (MQ=255)
                                                                                                        acTTCGACAGCTGCGGCAGATCACTTACCGGGAACATCAACATACGGCCTGCCTGAGTGATTGCCAGCAGCATATCGGAAGCATCTTCAATCACCACCGGCGGCATAACATGGGCATTTTCCGGTAAGGTGATCAAAGCCTTACCTGCa                                     >  7:82183/2‑149 (MQ=255)
                                                                                                                 gCTGCGGCAGATCACTTACCGGGAACATCAACATACGGCCTGCCTGAGTGATTGCCAGCAGCATATCGGAAGCATCTTCAATCACCACCGGCGGCATAACATGGGCATTTTCCGGTAAGGTGATCAAAGCCTTACCTGCACGGTTAcgc                            <  5:87899/149‑1 (MQ=255)
                                                                                                                 gCTGCGGCAGATCACTTACCGGGAACATCAACATACGGCCTGCCTGAGTGATTGCCAGCAGCATATCGGAAGCATCTTCAATCACCACCGGCGGCATAAAATGGGCATGTTCCGGTAAGGTGATAAAAGCCTTCCCTGCACGGTTAcgc                            >  6:63055/1‑149 (MQ=255)
                                                                                                                       gCAGATCACTTACCGGGAACATCAACATACGGCCTGCCTGAGTGATTGCCAGCAGCATATCGGAAGCATCTTCAATCACCACCGGCGGCATAACATGGGCATTTTCCGGTAAGGTGATCAAAGCCTTACCTGCACGGTTACGCGccac                       >  5:289020/1‑148 (MQ=255)
                                                                                                                          gATCACTTACCGGGAACATCAACATACGGCATGCCTGAGTGATTGCCACCAGCATATCGGAAGCATCTTCAAACACCACCGGCGGCATAACATGGGCAATTTACGGAAAGGTGATCAAAGCCTTAACTGCACGGTTACGCGACACAAg                    >  2:270645/1‑148 (MQ=255)
                                                                                                                               gttACCGGGTACATCAACATACGGCCTGCCTGATTGTTTGCCAGCAGCTTATCGGTTGCATCTTCATTCAGCAGCGGCGTCATTATATGGGCATTTTGCGGTAAGGTGATCAAAGCCTTACCTGCACGGTTAcgc                            <  4:145645/134‑1 (MQ=255)
                                                                                                                               cTTACCGGGAACATCAACATACGGCCTGCCTGAGTGATTGCCAGCAGCATATCGGAAGCATCTTCAATCACCACCGGCGGCATAACATGGGCATTTTCCGGTAAGGTGATCAAAGCCTTACCTGCACGGTTAcgc                            >  3:145645/1‑135 (MQ=255)
                                                                                                                                 tACCGGGAACATCAACATACGGCCTGCCTGAGTGATTGCCAGCAGCATATCGGAAGCATCTTCAATCACCACCGGCGGCATAACATGGGCATTTTCCGGTAAGGTGATCAAAGCCTTACCTGCACGGTTACGCGCCACCAGATCGTTa             >  5:297887/1‑148 (MQ=255)
                                                                                                                                    cGGGAACATCAACATACGGCCTGCCTGAGTGATTGCCAGCAGCATATCGGAAGCATCTTCAATCACCACCGGCGGCATAACATGTGCATTTTCCGGTAAGGTGATCAAAGCCTTACCTGCACGGTTACGCGCCACCAGATCGTTaa            >  1:340879/1‑146 (MQ=255)
                                                                                                                                    cGGGAACATCAACATACGGCCTGCCTGAGTGATTGCCAGCAGCATATCGGAAGCATCTTCAATCACCACCGGCGGCATAACATGTGCATTTTCCGGTAAGGTGATCAAAGCCTTACCTGCACGGTTACGCGCCACCAGATCGTTaa            <  2:340879/146‑1 (MQ=255)
                                                                                                                                    cGGGAACATCAACATACGGCCTGCCTGAGTGATTGCCAGCAGCATATCGGAAGCATCTTCAATCACCACCGGCGGCATAACATGGGCATTTTCCGGTAAGGTGATCAAAGCCTTACCTGCACGGTTACGCGCCACCAGATCGTTaaa           <  6:149938/147‑1 (MQ=255)
                                                                                                                                    cGGGAACATCAACATACGGCCTGCCTGAGTGATTGCCAGCAGCATATCGGAAGCATCTTCAATCACCACCGGCGGCATAACATGGGCATTTTCCGGTAAGGTGATCAAAGCCTTACCTGCACGGTTACGCGCCACCAGATCGTTaaa           >  5:149938/1‑147 (MQ=255)
                                                                                                                                      ggAACATCAACATACGGCCTGCCTGAGTGATTGCCAGCAGCATATCGGAAGCATCTTCAATCACCACcgg                                                                                      <  8:147168/70‑1 (MQ=255)
                                                                                                                                      ggAACATCAACATACGGCCTGCCTGAGTGATTGCCAGCAGCATATCGGAAGCAACTTCAATCACCACcgg                                                                                      >  7:147168/1‑70 (MQ=255)
                                                                                                                                         aCATCAACATATGGCCTGCCTGAGTGATTGCCAGCAGCATATCGGAAGCATCTTCAATCACCACCGGCGGCATAACATGGGCATTTTCCGGTAAGGTGATCAAAGCCTTACCTGCACGGTTACGCGCCACCAGATCGTTAAAGGTGCAg    <  1:100974/149‑1 (MQ=255)
                                                                                                                                         aCATCAACATATGGCATTCATGAGTGATTTCCAGCATCATATCGGAAGCATCTTCAATCAACAACTGCGGCATAACATGGGCATTTTCCTGTAAGGTGATCAAAGCCTTACCTGCACGGTTACGCGCCACCAGATCGTTAAAGGTGCAg    <  1:106005/149‑1 (MQ=255)
                                                                                                                                           aTCAACATACGGCCTGCCTGAGTGATTGCCAGCAGCATATCGGAAGCATCTTCAATCACCACCGGCGGCATAACATGGGCATTTTCCGGTAAGGTGATCAAAGCCTTACCTGCACGGTTACGCGCCACCAGATCGTTAAAGGTGc      >  7:285081/1‑145 (MQ=255)
                                                                                                                                           aTCAACATACGGCCTGCCTGAGTGATTGCCAGCAGCATATCGGAAGCATCTTCAATCACCACCGGCGGCATAACATGGGCATTTTCCGGTAAGGTGATCAAAGCCTTACCTGCACGGTTACGCGCCACCAGATCGTTAAAGGTGc      <  8:285081/145‑1 (MQ=255)
                                                                                                                                            tCAACATACGGCCTGCCTGAGTGATTGCCAGCAGCATATCGGAAGCATCTTCAATCACCACCGGCGGCATAACATGGGCATTTTCCGGTAAGGTGATCAAAGCCTTACCTGCACGGTTACGCGCCACCAGATCGTTAAAGGTGCAGAc  >  4:349534/1‑148 (MQ=255)
                                                                                                                                                   aCGGCCTGCCTGAGTGATTGCCAGCAGCATATCGGAAGCATCTTCAATCACCACCGGCGGCATAACATGGGCATTTTCCGGTAAGGTGATCAAAGCCTTACCTGCACGGTTACGCGCCACCAGATCGTTAAAg          <  2:220931/133‑1 (MQ=255)
                                                                                                                                                   aCGGCCTGCCTGAGTGATTGCCAGCAGCATATCGGAAGCATCTTCAATCACCACCGGCGGCATAACATGGGCATTTTCCGGTAAGGTGATCAAAGCCTTACCTGCACGGTTACGCGCCACCAGATCGTTAAAg          >  1:220931/1‑133 (MQ=255)
                                                                                                                                                    |                                                                                                                                           
AATGGTCAGCGTGCTTTGCGGCGGCAGAACGTACAATTGCGCCAGACCATCTTCTCCACGCGCGGCTTCTGCCGATGGAATGTTGATAATCTTGTTGCCTTTGCCCTTCGACAGCTGCGGCAGATCACTTACCGGGAACATCAACATACGGCCTGCCTGAGTGATTGCCAGCAGCATATCGGAAGCATCTTCAATCACCACCGGCGGCATAACATGGGCATTTTCCGGTAAGGTGATCAAAGCCTTACCTGCACGGTTACGCGCCACCAGATCGTTAAAGGTGCAGAC  >  NC_000913/3163874‑3164161

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑