Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,175,118 A→G 5.1% R93R (CGT→CGC parE ← DNA topoisomerase IV, subunit B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,175,1180AG 5.1% 395.7 / 7.9 117R93R (CGT→CGCparEDNA topoisomerase IV, subunit B
Reads supporting (aligned to +/- strand):  ref base A (55/56);  new base G (3/3);  total (58/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.48e-01

Only 100 of 117 total aligned reads displayed.

TCAAAGGCGATGTTATAAACCTGACCATCGCGGCGCACGTTAACTTCTACGCGCTTCGACAGGGCGTTAACCACCGAAATCCCCACGCCATGCAGGCCGCCAGAGAACTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGTGAATATCCACCGGCATCCCGCGCCCATCGTCAATAACTTCTAACGACTGGTCAGCATGTAAAATAACGTCCACGCGTTTTGCGTGACCCGCCAGTGCTTCATCC  >  NC_000913/3174972‑3175265
                                                                                                                                                  |                                                                                                                                                   
tCAAAGGCGATGTTATAAACCTGACCATCGCGGCGCACGTTAACTTCTACGCGCTTCGACAGGGCGTTAACCACCGAAATCCCCACGCCATGCAGGCCGCCAGAGAACTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACg                                                                                                                                                   >  2:181179/1‑149 (MQ=255)
 cAAAGGCGATGTTATAAACCTGACCATCGCGGCGCACGTTAACTTCTACGCGCTTCGACAGGGCGTTAACCACCGAAATCCCCACGCCATGCAGGCCGCCAGAGAACTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACgg                                                                                                                                                  <  1:1042328/149‑1 (MQ=255)
  aaaGGCGATGTTATAAACCTGACCATCGCGGCGCACGTTAACTTCTACGCGCTTCGACAGGGCGTTAACCACCGAAATCCCCACGCCATGCAGGCCGCCAGAGAACTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGc                                                                                                                                                 >  1:986369/1‑149 (MQ=255)
   aaGGCGATGTTATAAACCTGACCATCGCGGCGCACGTTAACTTCTACGCGCTTCGACAGGGCGTTAACCACCGAAATCCCCACGCCATGCAGGCCGCCAGAGAACTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGCa                                                                                                                                                >  1:1112260/1‑149 (MQ=255)
    aGGCGATGTTATAAACCTGACCATCGCGGCGCACGTTAACTTCTACGCGCTTCGACAGGGCGTTAACCACCGAAATCCCCACGCCATGCAGGCCGCCAGAAAGCTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGCaa                                                                                                                                               >  1:427784/1‑149 (MQ=255)
     ggCGATGTTATAAACCTGACCATCGCGGCGCACGTTAACTTCTACGCGCTTCGACAGGGCGTTAACCACCGAAATCCCCACGCCATGCAGGCCGCCAGAGAACTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGCaaa                                                                                                                                              <  1:1065184/149‑1 (MQ=255)
     ggCGATGTTATAAACCTGACCATCGCGGCGCACGTTAACTTCTACGCGCTTCGACAGGGCGTTAACCACCGAAATCCCCACGCCATGCAGGCCGCCAGAAAACTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGCaaa                                                                                                                                              <  2:427784/149‑1 (MQ=255)
      gCGATGTTATAAACCTGACCATCGCGGCGCACGTTAACTTCTACGCGCTTCGACAGGGCGTTAACCACCGAAATCCCCACGCCATGCAGGCCGCCAGAGAACTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGAc                                                                                                                                                    >  3:799954/1‑142 (MQ=255)
        gATGTTATAAACCTGACCATCGCGGCGCACGTTAACTTCTACGCGCTTCGACAGGGCGTTAACCACCGAAATCCCCACGCCATGCAGGCCGCCAGAGAACTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGCAAAGaa                                                                                                                                           <  2:986369/149‑1 (MQ=255)
        gATGTTATAAACCTGACCATCGCGGCGCACGTTAACTTCTACGCGCTTCGACAGGGCGTTAACCACCGAAATCCCCACGCCATGCAGGCCGCCAGAGAACTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGCAAAGaa                                                                                                                                           >  3:603434/1‑149 (MQ=255)
        gATGTTATAAACCTGACCATCGCGGCGCACGTTAACTTCTACGCGCTTCGACAGGGCGTTAACCACCGAAATCCCCACGCCATGCAGGCCGCCAGAGAACTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGCAAAGaa                                                                                                                                           >  3:616990/1‑149 (MQ=255)
         aTGTTATAAACCTGACCATCGCGGCGCACGTTAACTTCTACGCGCTTCGACAGGGCGTTAACCACCGAAATCCCCACGCCATGCAGGCCGCCAGAGAACTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGAc                                                                                                                                                    <  1:1010373/139‑1 (MQ=255)
         aTGTTATAAACCTGACCATCGCGGCGCACGTTAACTTCTACGCGCTTCGACAGGGCGTTAACCACCGAAATCCCCACGCCATGCAGGCCGCCAGAGAACTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGAc                                                                                                                                                    >  2:1010373/1‑139 (MQ=255)
           gTTATAAACCTGACCATCGCGGCGCACGTTAACTTCTACGCGCTTCGACAGGGCGTTAACCACCGAAATCCCCACGCCATGCAGGCCGCCAGAGAACTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGCAAAGAATCa                                                                                                                                        >  4:17036/1‑149 (MQ=255)
              ttaAACCTGACCATCGCGGCGCACGTTAACTTCTACGCGCTTCGACAGGGCGTTAACCACCGAAATCCCCACGCCATGCAGGCCGCCAGAGAACTGGTAATTTTTGTTAGAGAATTTACCGCCTGTATGCAGACGGCAAAGAATCAGtt                                                                                                                                     <  3:1263360/148‑1 (MQ=255)
               taAACCTGACCATCGCGGCGCACGTTAACTTCTACGCGCTTCGACAGGGCGTTAACCACCGAAATCCCCACGCCATGCAGGCCGCCAGAGAACTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGCAAAGAATCAGTTc                                                                                                                                    <  3:1113160/149‑1 (MQ=255)
                aaaCCTGACCATCGCGGCGCACGTTAACTTCTACGCGCTTCGACAGGGCGTTAACCACCGAAATCCCCACGCCATGCAGGCCGCCAGAGAACTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGCAAAGAATCAGTTCa                                                                                                                                   >  1:358319/1‑149 (MQ=255)
                aaaCCTGACCATCGCGGCGCACGTTAACTTCTACGCGCTTCGACAGGGCGTTAACCACCGAAATCCCCACGCCATGCAGGCCGCCAGAGAACTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGCAAAGAATCAGTTCa                                                                                                                                   >  1:559479/1‑149 (MQ=255)
                 aaCCTGACCATCGCGGCGCACGTTAACTTCTACGCGCTTCGACAGGGCGTTAACCACCGAAATCCCCACGCCATGCAGGCCGCCAGAGAACTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGCAAAGAATCAGTTCaa                                                                                                                                  <  4:616990/149‑1 (MQ=255)
                    cTGACCATCGCGGCGCACGTTAACTTCTACGCGCTTCGACAGGGCGTTAACCACCGAAATCCCCACGCCATGCAGGCCGCCAGAGAACTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGCAAAGAATCAGTTCAAccg                                                                                                                               >  1:32991/1‑149 (MQ=255)
                        ccATCGCGGCGCACGTTAACTTCTACGCGCTTCGACAGGGCGTTAACCACCGAAATCCCCACGCCATGCAGGCCGCCAGAGAACTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGCAAAGAATCAGTTCAACCGCCgg                                                                                                                           <  2:32991/149‑1 (MQ=255)
                          aTCGCGGCGCACGTTAACTTCTACGCGCTTCGACAGGGCGTTAACCACCGAAATCCCCACGCCATGCAGGCCGCCAGAGAACTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGCAAAGAATCa                                                                                                                                        <  3:327039/134‑1 (MQ=255)
                          aTCGCGGCGCACGTTAACTTCTACGCGCTTCGACAGGGCGTTAACCACCGAAATCCCCACGCCATGCAGGCCGCCAGAGAACTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGCAAAGAATCa                                                                                                                                        >  4:327039/1‑134 (MQ=255)
                           tCGCGGCGCACGTTAACTTCTACGCGCTTCGACAGGGCGTTAACCACCGAAATCCCCACGCCATGCAGGCCGCCAGAGAACTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGCAAAGAATCAGTTCAACCGCCGGTac                                                                                                                        >  3:338000/1‑149 (MQ=255)
                                 cgcACGTTAACTTCTACGCGCTTCGACAGGGCGTTAACCACCGAAATCCCCACGCCATGCAGGCCGCCAGAGAACTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGCAAAGAATCAGTTCAACCGCCGGTACACCctc                                                                                                                  <  4:788946/149‑1 (MQ=255)
                                    aCGTTAACTTCTACGCGCTTCGACAGGGCGTTAACCACCGAAATCCCCACGCCATGCAGGCCGCCAGAGAACTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACgg                                                                                                                                                  <  2:1457120/114‑1 (MQ=255)
                                    aCGTTAACTTCTACGCGCTTCGACAGGGCGTTAACCACCGAAATCCCCACGCCATGCAGGCCGCCAGAGAACTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACgg                                                                                                                                                  >  1:1457120/1‑114 (MQ=255)
                                     cGTTAACTTCTACGCGCTTCGACAGGGCGTTAACCACCGAAATCCCCACGCCATGCAGGCCGCCAGAGAACTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTcc                                                                                                              <  3:1149505/149‑1 (MQ=255)
                                      gTTAACTTCTACGCGCTTCGACAGGGCGTTAACCACCGAAATCCCCACGCCATGCAGGCCGCCAGAGAACTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCg                                                                                                             >  3:953892/1‑149 (MQ=255)
                                       ttAACTTCTACGCGCTTCGACAGGGCGTTAACCACCGAAATCCCCACGCCATGCAGGCCGCCAGAGAACTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCgg                                                                                                            <  2:714581/149‑1 (MQ=255)
                                         aaCTTCTACGCGCTTCGACAGGGCGTTAACCACCGAAATCCCCACGCCATGCAGGCCGCCAGAGAACTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGt                                                                                                          <  2:559479/149‑1 (MQ=255)
                                         aaCTTCTACGCGCTTCGACAGGGCGTTAACCACCGAAATCCCCACGCCATGCAGGCCGCCAGAGAACTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGt                                                                                                          <  4:953892/149‑1 (MQ=255)
                                         aaCTTCTACGCGCTTCGACAGGGCGTTAACCACCGAAATCCCCACGCCATGCAGGCCGCCAGAGAACTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGt                                                                                                          <  2:129855/149‑1 (MQ=255)
                                          aCTTCTACGCGCTTCGACAGGGCGTTAACCACCGAAATCCCCACGCCATGCAGGCCGCCAGAGAACTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGTg                                                                                                         <  3:421030/149‑1 (MQ=255)
                                                    gcTTCGACAGGGCGTTAACCACCGAAATCCCCACGCCATGCAGGCCGCCAGAGAACTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGCAAAGAAt                                                                                                                                          <  4:1241688/106‑1 (MQ=255)
                                                         gACAGGGCGTTAACCACCGAAATCCCCACGCCATGCAGGCCGCCAGAGAACTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGTGAATATCCACCGGCAt                                                                                          >  3:362417/1‑149 (MQ=255)
                                                         gACAGGGCGTTAACCACCGAAATCCCCACGCCATGCAGGCCGCCAGAGAACTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGTGAATATCCACCGGCAt                                                                                          >  1:525289/1‑149 (MQ=255)
                                                             gggCGTTAACCACCGAAATCCCCACGCCATGCAGGCCGCCAGAGAACTGGTAATTTTTGTTAGAGAATTTACCGTCTGCATGCAGGCGGTaaa                                                                                                                                              <  4:230314/93‑1 (MQ=255)
                                                             gggCGTTAACCACCGAAATCCCCACGCCATGCAGGCCGCCAGAGAACTGGTAATTTTTGTTAGAGAATTTACCGTCTGCATGCAGGCGGTaaa                                                                                                                                              >  3:230314/1‑93 (MQ=255)
                                                                cGTTAACCACCGAAATCCCCACGCCATGCAGGCCGCCAGAGAACTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGTGAATATCCACCGGCATCCCgcgc                                                                                   <  2:525289/149‑1 (MQ=255)
                                                                cGTTAACCACCGAAATCCCCACGCCATGCAGGCCGCCAGAGAACTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGTGAATATCCACCGGCATCCCgcgc                                                                                   <  1:181179/149‑1 (MQ=255)
                                                                       caccGAAATCCCCACGCCATGCAGGCCGCCAGAGAACTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGTGAATATCCACCGGCATCCCGCGCCCATCGt                                                                            <  4:362417/149‑1 (MQ=255)
                                                                       caccGAAATCCCCACGCCATGCAGGCCGCCAGAGAACTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGTGAATATCCACCGGCATCCCGCGCCCATCGt                                                                            >  4:524025/1‑149 (MQ=255)
                                                                        accGAAATCCCCACGCCATGCAGGCCGCCAGAGAACTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGCAATGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGTGAat                                                                                                      <  4:851833/122‑1 (MQ=255)
                                                                        accGAAATCCCCACGCCATGCAGGCCGCCAGAGAACTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTc                                                                                                               >  4:878311/1‑113 (MQ=255)
                                                                        accGAAATCCCCACGCCATGCAGGCCGCCAGAGAACTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTc                                                                                                               <  3:878311/113‑1 (MQ=255)
                                                                        accGAAATCCCCACGCCATGCAGGCCGCCAGAGAACTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGTGAat                                                                                                      >  3:851833/1‑122 (MQ=255)
                                                                        accGAAATCCCCACGCCATGCAGGCCGCCAGAGAACTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGTGAATATCCACCGGCATCCCGCGCCCATCGTc                                                                           >  3:263676/1‑149 (MQ=255)
                                                                         ccGAAATCCCCACGCCATGCAGGCCGCCAGAGAACTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGTGAATATCCACCGGCATCCCGCGCCCATCGTCa                                                                          >  4:869686/1‑149 (MQ=255)
                                                                             aaTCCCCACGCCATGCAGGCCGCCAGAGAACTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGTGAATATCCACCGGCATCCCGCGCCCATCGTCAATaa                                                                      <  2:1112260/149‑1 (MQ=255)
                                                                             aaTCCCCACGCCATGCAGGCCGCCAGAGAACTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGTGAATATCCACCGGCATCCCGCGCCCATCGTCAATaa                                                                      <  3:524025/149‑1 (MQ=255)
                                                                              aTCCCCACGCCATGCAGGCCGCCAGAGAACTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGTGAATATCCACCGGCATCCCGCGCCCATCGTCAATAAc                                                                     <  2:1174843/149‑1 (MQ=255)
                                                                                 cccACGCCATGCAGGCCGCCAGAGAACTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGTGAATATCCACCg                                                                                              <  1:140823/121‑1 (MQ=255)
                                                                                 cccACGCCATGCAGGCCGCCAGAGAACTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGTGAATATCCACCg                                                                                              >  2:140823/1‑121 (MQ=255)
                                                                                  ccACGCCATGCAGGCCGCCAGAGAACTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGTGAATATCCACCGGCATCCCGCGCCCATCGTCAATAACTTCt                                                                 >  4:233610/1‑149 (MQ=255)
                                                                                        cATGCAGGCCGCCAGAGAACTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGTGAATATCCACCGGCATCCCGCGCCCATCGTCAATAACTTCTAACGAc                                                           <  3:17036/149‑1 (MQ=255)
                                                                                        cATGCAGGCCGCCAGAGAACTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGTGAATATCCACCGGCATCCCGCGCCCATCGTCAATAACTTCTAACGAc                                                           >  4:600789/1‑149 (MQ=255)
                                                                                               gccgccAGAGAACTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGTGAATATCCACCGGCATCCCGCGCCCATCGTCAATAACtt                                                                   <  1:1001620/134‑1 (MQ=255)
                                                                                               gccgccAGAGAACTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGTGAATATCCACCGGCATCCCGCGCCCATCGTCAATAACtt                                                                   >  2:1001620/1‑134 (MQ=255)
                                                                                                     agagAACTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGCAAAGa                                                                                                                                            <  1:1305600/55‑1 (MQ=255)
                                                                                                     agagAACTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGCAAAGa                                                                                                                                            >  2:1305600/1‑55 (MQ=255)
                                                                                                     agagAACTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGTGAATATCCACCGGCATCCCGCGCCCATCGTCAATAACTTCTAACGACTGGTCAGCATGTa                                              <  3:600789/149‑1 (MQ=255)
                                                                                                          aCTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGTGAATATCCACCg                                                                                              >  4:264413/1‑96 (MQ=255)
                                                                                                          aCTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGTGAATATCCACCg                                                                                              <  3:264413/96‑1 (MQ=255)
                                                                                                            tGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGTGAATATCCACCGGCATCCCGCGCCCATCGTCAATAACTTCTAACGACTGGTCAGCATGTAAAATAAc                                       <  3:1325373/149‑1 (MQ=255)
                                                                                                            tGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGTGAATATCCACCGGCATCCCGCGCCCATCGTCAATAACTTCTAACGACTGGTCAGCATGTAAAATAAc                                       <  3:869686/149‑1 (MQ=255)
                                                                                                              gTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGTGAATATCCACCGGCATCCCGCGCCCATCGTCAATAACTTc                                                                  <  1:1069229/120‑1 (MQ=255)
                                                                                                              gTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGTGAATATCCACCGGCATCCCGCGCCCATCGTCAATAACTTc                                                                  >  2:1069229/1‑120 (MQ=255)
                                                                                                               tAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGTGAATATCCACCGGCATCCCGCGCCCATCGTCAATAAc                                                                     >  3:1135441/1‑116 (MQ=255)
                                                                                                               tAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGTGAATATCCACCGGCATCCCGCGCCCATCGTCAATAAc                                                                     <  4:1135441/116‑1 (MQ=255)
                                                                                                                 aTTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGTGAATATCCACCGGCATCCCGCGCCCATCGTCAATAACTTCTAACGACTGGTCAGCATGTAAAATAACGTCCa                                  >  3:1048783/1‑149 (MQ=255)
                                                                                                                   ttttGTTAGAGAATTTACCGCCTGCATGCAGACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGTGAATATCCACCGGCATCCCGCGCCCATCGTCAATAACTTCTAACGACTGGTCAGCATGTAAAATAACGTCCAcg                                >  4:457336/1‑149 (MQ=255)
                                                                                                                        ttAGAGAATTTACCGCCTGCATGCAGACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGTGAATATCCACCGGCATCCCGCGCCCATCGTCAATAACTTCTAACGACTGGTCAGCATGTaaa                                            >  1:96808/1‑132 (MQ=255)
                                                                                                                          agagAATTTACCGCCTGCATGCAGACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGTGAATATCCACCGGCATCCCGCGCCCATCGTCAATAACTTCTAACGACTGGTCAGCa                                                  >  4:1393638/1‑124 (MQ=255)
                                                                                                                          agagAATTTACCGCCTGCATGCAGACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGTGAATATCCACCGGCATCCCGCGCCCATCGTCAATAACTTCTAACGACTGGTCAGCa                                                  <  3:1393638/124‑1 (MQ=255)
                                                                                                                          agagAATTTACCGCCTGCATGCAGACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGTGAATATCCACCGGCATCCCGCGCCCATCGTCAATAACTTCTAACGACTGGTCAGCATGTAAAATAACGTCCACGCGTTTTg                         >  4:1019282/1‑149 (MQ=255)
                                                                                                                           gagaATTTACCGCCTGCATGCAGACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGTGAATATCCACCGGCATCCCGCGCCCATCGTCAATAACTTCTaa                                                               >  1:1379075/1‑110 (MQ=255)
                                                                                                                           gagaATTTACCGCCTGCATGCAGACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGTGAATATCCACCGGCATCCCGCGCCCATCGTCAATAACTTCTAACGACTGGTc                                                      <  1:975509/119‑1 (MQ=255)
                                                                                                                           gagaATTTACCGCCTGCATGCAGACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGTGAATATCCACCGGCATCCCGCGCCCATCGTCAATAACTTCTAACGACTGGTc                                                      >  2:975509/1‑119 (MQ=255)
                                                                                                                           gagaATTTACCGCCTGCATGCAGACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGTGAATATCCACCGGCATCCCGCGCCCATCGTCAATAACTTCTAACGACTGGTCAGc                                                   >  4:484645/1‑122 (MQ=255)
                                                                                                                           gagaATTTACCGCCTGCATGCAGACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGTGAATATCCACCGGCATCCCGCGCCCATCGTCAATAACTTCTAACGACTGGTCAGc                                                   <  3:484645/122‑1 (MQ=255)
                                                                                                                           gagaATTTACCGCCTGCATGCAGACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGTGAATATCCACCGGCATCCCGCGCCCATCGTCAATAACTTCTAACGACTGGTCAGCATGTAAAATAACGTCCACGCGTTTTGc                        >  1:926928/1‑149 (MQ=255)
                                                                                                                            agaATTTACCGCCTGCATGCAGACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGTGAATATCCACCGGCATCCCGCGCCCATCGTCAATAACTTCTAACGACTGGTCAGCATGTAAAATAACGTCCACGCGTTTTGc                        <  2:169181/148‑1 (MQ=255)
                                                                                                                             gaATTTACCGCCTGCATGCAGACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGTGAATATCCACCGGCATCCCGCGCCCATCGTCAATAACTTCTAACGACTGGTCAGCATGTAAAATAACGTCCACGCGTTTTGCGt                      <  2:627435/149‑1 (MQ=255)
                                                                                                                               aTTTACCGCCTGCATGCAGACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGTGAATATCCACCGGCATCCCGCGCCCATCGTCAATAACTTCTAACGACTGGTCAGCATGTAAAATAACGTCCACGCGTTTTGCGTGa                    <  4:901992/149‑1 (MQ=255)
                                                                                                                                 ttgccGTCTGCATGCAGGCGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGTGAATATCCACCGGCATCCCGCGCCCATCGTCaa                                                                         >  1:1455711/4‑94 (MQ=255)
                                                                                                                                 ttgccGTCTGCATGCAGGCGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGTGAATATCCACCGGCATCCCGCGCCCATCGTCaa                                                                         <  2:1455711/91‑1 (MQ=255)
                                                                                                                                  tACCGCCTGCATGCAGACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGTGAATATCCACCGGCATCCCGCGCCCATCGTCAATAACTTCTAACGACTGGTCAGCATGTAAAATAACGTCCACGCGTTTTGCGTGAccc                 <  4:639605/149‑1 (MQ=255)
                                                                                                                                    ccGTCTGCATGCAGGCGGTAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGTGAATATCCACCGGCATCCCGCGCCCATCGTCAATAACTTCTAACGACTGGTCAGCATGTAAAATAACGTCCACGCGttt                           <  1:1287078/137‑1 (MQ=255)
                                                                                                                                    ccGTCTGCATGCAGGCGGTAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGTGAATATCCACCGGCATCCCGCGCCCATCGTCAATAACTTCTAACGACTGGTCAGCATGTAAAATAACGTCCACGCGttt                           >  2:1287078/1‑137 (MQ=255)
                                                                                                                                    ccGCCTGCATGCAGACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGTGAATATCCACCGGCATCCCGCGCCCATAGTCAATAACTTCTAACGACTGGTCAGCATGTAAAATAACGTCCACGCGTTTTGCGTGACCCGc               >  4:540548/1‑149 (MQ=255)
                                                                                                                                      gCCTGCATGCAGACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGTGAATATCCACCGGCATCCCGCGCCCATCATCAATAACTTCTAACGACTGGTCAGCATGTAAAATAACGTCCACGCGTTTTGCGTGACCCGCCa             <  3:540548/149‑1 (MQ=255)
                                                                                                                                        cTGCATGCAGACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGTGAATATCCACCGGCATCCCGCGCCCATCGTCAATaa                                                                      <  3:734141/90‑1 (MQ=255)
                                                                                                                                        cTGCATGCAGACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGTGAATATCCACCGGCATCCCGCGCCCATCGTCAATaa                                                                      >  4:734141/1‑90 (MQ=255)
                                                                                                                                         tgcatgcaGACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGTGAATATCCACCGGCATCCCGCGCCCATCGTCAATAACTTCTAACGACTGGTCAGCATGTAAAATAACGTCCACGCGTTTTGCGTGACCCGCCAGTg          >  4:1076739/1‑149 (MQ=255)
                                                                                                                                          gcatgcaGACGGCAAAGAATCAGTTCACCCGCCGGTACACCCTCTTCCGGGTGAATATCCACCGGCATCCCGCGCCCATCGTCaa                                                                         >  4:742174/1‑85 (MQ=255)
                                                                                                                                          gcatgcaGACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGTGAATATCCACCGGCATCCCGCGCCCATCGTCaa                                                                         <  3:742174/85‑1 (MQ=255)
                                                                                                                                               caGACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGTGAATATCCACCGGCATCCCGCGCCCATCGTCAATAACTTCTAACGACTGGTCAGCATGTAAAATAACGTCCACGCGTTTTGCGTGACCCGCCAGTGCTTCAt    <  3:999274/149‑1 (MQ=255)
                                                                                                                                                 gACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGTGAATATCCACCGGCATCCCGCGCCCATCGTCAATAACTTCTAACGACTGGTCAGCATGTAAAATAACGTCCACGCGTTTTGCGTGACCCGCCAGTGCTTCATcc  >  2:699608/1‑149 (MQ=255)
                                                                                                                                                 gACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGTGAATATCCACCGGCATCCCGCGCCCATCGTCAATAACTTCTAACGACTGGTCAGCATGTAAAATAACGTCCACGCGTTTTGCGTGACCCGCCAGTGCTTCATcc  >  4:668444/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TCAAAGGCGATGTTATAAACCTGACCATCGCGGCGCACGTTAACTTCTACGCGCTTCGACAGGGCGTTAACCACCGAAATCCCCACGCCATGCAGGCCGCCAGAGAACTGGTAATTTTTGTTAGAGAATTTACCGCCTGCATGCAGACGGCAAAGAATCAGTTCAACCGCCGGTACACCCTCTTCCGGGTGAATATCCACCGGCATCCCGCGCCCATCGTCAATAACTTCTAACGACTGGTCAGCATGTAAAATAACGTCCACGCGTTTTGCGTGACCCGCCAGTGCTTCATCC  >  NC_000913/3174972‑3175265

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