Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,589,790 A→C 6.4% L444R (CTG→CGG) ‡ parE ← DNA topoisomerase IV subunit B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,589,7900AC 6.4% 664.7 / 26.3 202L444R (CTG→CGG) ‡parEDNA topoisomerase IV subunit B
Reads supporting (aligned to +/- strand):  ref base A (100/89);  new base C (4/9);  total (104/98)
Fisher's exact test for biased strand distribution p-value = 1.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.21e-01

Only 100 of 205 total aligned reads displayed.

GGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTGGCCTGCACAATCGGCCAGCTTGCCGGGCAATGCCGGGCCCTGGGTGATG  >  NC_002947/5589642‑5589939
                                                                                                                                                     |                                                                                                                                                     
ggCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCa                                                                                                                                                       <  1:1986686/150‑1 (MQ=255)
   gAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCgg                                                                                                                                                   <  2:3058022/150‑1 (MQ=255)
     ggATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCa                                                                                                                                                 <  1:1418791/150‑1 (MQ=255)
      gATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCGGATAGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCCTGCCGCTGCGCGGCaa                                                                                                                                                <  1:3533500/150‑2 (MQ=255)
       aTGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAgg                                                                                                                                               <  1:3422805/150‑1 (MQ=255)
          cAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATg                                                                                                                                            <  2:4295508/150‑1 (MQ=255)
          cAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCCTGCCGCTGCGCGGCAAGAtc                                                                                                                                            <  1:3679825/150‑2 (MQ=255)
            gATCTTGCCATAGCGCAGCTGCGCAATGTCGGTTGCACCCGGATCGACGCCAATGGCCTCGGCAATGTTGTGCTCTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGTCAGGATGGc                                                                                                                                          <  2:4632498/150‑1 (MQ=255)
              tCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCt                                                                                                                                        <  2:966677/150‑1 (MQ=255)
                ttGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTgg                                                                                                                                      <  1:3307967/150‑1 (MQ=255)
                 tGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGa                                                                                                                                     <  1:450675/150‑1 (MQ=255)
                  gCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGaa                                                                                                                                    <  1:4878765/150‑1 (MQ=255)
                  gCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGaa                                                                                                                                    <  1:4395893/150‑1 (MQ=255)
                   ccATAGCGCAGCTGCGCAAGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGAt                                                                                                                                             <  2:4955861/140‑1 (MQ=255)
                   ccATAGCGCAGCTGCGCAAGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAAt                                                                                                                                   >  1:2686804/1‑150 (MQ=255)
                    cATAGCGCAGCTGCGCAAGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAAtt                                                                                                                                  >  1:4109563/1‑150 (MQ=255)
                      tAGCGCAGCTGCGCAAGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTcc                                                                                                                                >  1:4057281/1‑150 (MQ=255)
                       aGCGCAGCTGCGCAAGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCt                                                                                                                               <  1:3616788/150‑1 (MQ=255)
                        gcgcAGCTGCGCAAGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAg                                                                                                                                                      <  1:705408/126‑1 (MQ=255)
                        gcgcAGCTGCGCAAGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCtt                                                                                                                              <  2:2285533/150‑1 (MQ=255)
                         cgcAGCTGCGCAAGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTg                                                                                                                             <  1:3757425/150‑1 (MQ=255)
                             gCTGCGCAAGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTccc                                                                                                                         <  2:4470045/150‑1 (MQ=255)
                                 cgcAAGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCgcgc                                                                                                                     <  2:3647718/150‑1 (MQ=255)
                                  gcAAGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGcc                                                                                                                    >  2:1655626/1‑150 (MQ=255)
                                    aaGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTg                                                                                                                  >  1:1117248/1‑150 (MQ=255)
                                        tCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTg                                                                                                              <  1:952948/150‑1 (MQ=255)
                                          ggTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGc                                                                                                            <  1:2358287/150‑1 (MQ=255)
                                          ggTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCCAGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGc                                                                                                            >  1:4913826/1‑150 (MQ=255)
                                            ttGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCTGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCg                                                                                                          >  2:797335/1‑150 (MQ=255)
                                            ttGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGGGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGACCTCCGTGTCCCGCGCCTGCTTGGCCg                                                                                                          >  2:1744728/1‑150 (MQ=255)
                                              gCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGaa                                                                                                        <  2:3013415/150‑1 (MQ=255)
                                              gCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGaa                                                                                                        <  2:2278386/150‑1 (MQ=255)
                                              gCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGGCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGaa                                                                                                        <  2:3195483/150‑1 (MQ=255)
                                                  ccGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAaccac                                                                                                    >  1:1533270/1‑150 (MQ=255)
                                                  ccGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAaccac                                                                                                    >  2:4867416/1‑150 (MQ=255)
                                                  ccGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAaccac                                                                                                    >  1:4113927/1‑150 (MQ=255)
                                                   cGGATCGACGCCAATGGCCACGGCAATGGTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAaccacc                                                                                                   <  2:2830476/150‑1 (MQ=255)
                                                    ggATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCg                                                                                                  <  2:649619/150‑1 (MQ=255)
                                                     gATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCgg                                                                                                 <  1:1829822/150‑1 (MQ=255)
                                                         gACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGa                                                                                             <  1:1219757/150‑1 (MQ=255)
                                                         gACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGa                                                                                             >  1:1327361/1‑150 (MQ=255)
                                                               aaTGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCAcc                                                                                       <  1:155957/150‑1 (MQ=255)
                                                                aTGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCAccc                                                                                      >  1:4779490/1‑150 (MQ=255)
                                                                aTGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCAccc                                                                                      <  2:535194/150‑1 (MQ=255)
                                                                    ccACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTcga                                                                                  <  1:4647585/150‑1 (MQ=255)
                                                                      aCGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTcgacg                                                                                >  1:2596738/1‑150 (MQ=255)
                                                                       cGGCAATGTTGTGCACGTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTcgacga                                                                               <  2:1149471/150‑1 (MQ=255)
                                                                       cGGCAATGTTGTGCACGTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTcgacga                                                                               <  2:752128/150‑1 (MQ=255)
                                                                        ggCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAg                                                                              <  1:582377/150‑1 (MQ=255)
                                                                        ggCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAg                                                                              <  2:847506/150‑1 (MQ=255)
                                                                            aTGTTGTGCACGTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGtac                                                                          <  2:4352280/150‑1 (MQ=255)
                                                                             tgttgtGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGaaca                                                                         >  2:2117737/1‑150 (MQ=255)
                                                                                 ggcgaagtCCTGGCCAGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTc                                                                     >  1:2012595/8‑150 (MQ=255)
                                                                                  tgCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCg                                                                    <  1:3829580/150‑1 (MQ=255)
                                                                                   gCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCgg                                                                   >  1:3925699/1‑150 (MQ=255)
                                                                                    cACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGc                                                                  >  2:1975790/1‑150 (MQ=255)
                                                                                     aCTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCa                                                                 <  2:4568566/150‑1 (MQ=255)
                                                                                       ttCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACg                                                               >  2:3314683/1‑150 (MQ=255)
                                                                                       ttCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACg                                                               >  1:3132376/1‑150 (MQ=255)
                                                                                             gctggcCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCgg                                                         <  2:2300898/150‑1 (MQ=255)
                                                                                                   cAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGGCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTg                                                   >  2:3165014/1‑150 (MQ=255)
                                                                                                   cAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTg                                                   >  1:1019712/1‑150 (MQ=255)
                                                                                                   cAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTg                                                   >  2:7204/1‑150 (MQ=255)
                                                                                                   cAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTg                                                   >  2:3396558/1‑150 (MQ=255)
                                                                                                   cAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTg                                                   >  1:2148800/1‑150 (MQ=255)
                                                                                                      gACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTGGcc                                                >  1:1681301/1‑150 (MQ=255)
                                                                                                         ttCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTGGCCTGc                                             >  2:2974882/1‑150 (MQ=255)
                                                                                                         ttCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTGGCCTGc                                             >  2:965389/1‑150 (MQ=255)
                                                                                                         ttCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTGGCCTGc                                             >  1:4313849/1‑150 (MQ=255)
                                                                                                         ttCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTGGCCTGc                                             >  2:4747581/1‑150 (MQ=255)
                                                                                                         ttCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTGGCCTGc                                             >  2:350995/1‑150 (MQ=255)
                                                                                                          tCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTGGCCTGca                                            >  1:2458691/1‑150 (MQ=255)
                                                                                                               taccGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTGGCCTGCACAAt                                        >  2:2356901/3‑150 (MQ=255)
                                                                                                                cccGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTGGCCTGCACAATc                                       <  2:4254432/150‑1 (MQ=255)
                                                                                                                 ccGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTGGCCTGCACAATCg                                      <  2:1591511/150‑1 (MQ=255)
                                                                                                                    tCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTGGCCTGCACAATCGGcc                                   >  1:3096571/1‑150 (MQ=255)
                                                                                                                    tCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGAGGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTGGCCTGCACAATCGGcc                                   <  2:2153186/150‑1 (MQ=255)
                                                                                                                    tCCACGTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTGGCCTGCACAATCGGcc                                   <  2:3894698/150‑1 (MQ=255)
                                                                                                                      cACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTGGCCTGCACAATCGGCCAg                                 >  2:4733636/1‑150 (MQ=255)
                                                                                                                       aCTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTGGCCTGCACAATCGGCCAGc                                >  2:673067/1‑150 (MQ=255)
                                                                                                                         gtCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTGGCCTGCACAATCGGCCAGCtt                              <  1:3991132/149‑1 (MQ=255)
                                                                                                                            ccAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTGGCCTGCACAATCGGCCAGCTTGcc                           >  1:2987743/1‑150 (MQ=255)
                                                                                                                                gtgtTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTGGCCTGCACAATCGGCCAGCTTGCCGGGc                       >  1:1998527/1‑150 (MQ=255)
                                                                                                                                gtgtTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTGGCCTGCACAATCGGCCAGCTTGCCGGGc                       <  2:4109599/150‑1 (MQ=255)
                                                                                                                                   ttGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTGGCCTGCACAATCGGCCAGCTTGCCGGGCAAt                    >  2:1554691/1‑150 (MQ=255)
                                                                                                                                     ggccaTCCTGCCGCTGCGCGGCAAGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTGGCCTGCACAATCGGCCAGCTTGCCGGGCAATGc                  >  2:4398665/5‑150 (MQ=255)
                                                                                                                                     gAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTGGCCTGCACAATCGGCCAGCTTGCCGGGCAATGc                  >  1:686428/1‑150 (MQ=255)
                                                                                                                                     gAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTGGCCTGCACAATCGGCCAGCTTGCCGGGCAATGc                  >  2:586287/1‑150 (MQ=255)
                                                                                                                                      aGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTGGCCTGCACAATCGGCCAGCTTGCCGGGCAATGcc                 >  2:3767637/1‑150 (MQ=255)
                                                                                                                                      aGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTGGCCTGCACAATCGGCCAGCTTGCCGGGCAATGcc                 >  1:1170006/1‑150 (MQ=255)
                                                                                                                                        gATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTGGCCTGCACAATCGGCCAGCTTGCCGGGCAATGCCgg               >  1:3613422/1‑150 (MQ=255)
                                                                                                                                        gATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTGGCCTGCACAATCGGCCAGCTTGCCGGGCAATGCCgg               >  1:2674119/1‑150 (MQ=255)
                                                                                                                                        gATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTGGCCTGCACAATCGGCCAGCTTGCCGGGCAATGCCgg               >  2:331750/1‑150 (MQ=255)
                                                                                                                                         aTCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTGGCCTGCACAATCGGCCAGCTTGCCGGGCAATGCCggg              >  1:979640/1‑150 (MQ=255)
                                                                                                                                         aTCCTGCCGCTGCGCGGCAAGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTGGCCTGCACAATCGGCCAGCTTGCCGGGCAATGCCggg              >  2:1012261/1‑150 (MQ=255)
                                                                                                                                                cgcgCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTGGCCTGCACAATCGGACAGCTTGCAGGGCAATGCCGGGCCCTggg       >  1:133450/1‑150 (MQ=255)
                                                                                                                                                 gcgcAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTGGCCTGCACAATCGGCCAGCTTGCCGGGCAATGCCGGGCCCTGGGt      >  1:232764/1‑150 (MQ=255)
                                                                                                                                                 gcgcAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTGGCCTGCACAATCGGCCAGCTTGCCGGGCAATGCCGGGCCCTGGGt      <  2:4057281/150‑1 (MQ=255)
                                                                                                                                                   gcAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTGGCCTGCACAATCGGCCAGCTTGCCGGGCAATGCCGGGCCCTGGGTGa    <  1:2633628/150‑1 (MQ=255)
                                                                                                                                                    cAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTGGCCTGCACAATCGGCCAGCTTGCCGGGCAATGCCGGGCCCTGGGTGAt   <  1:756974/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCGGATCGACGCCAATGGCCACGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCACTTCCCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCACGCATCGGGTCCTGGCCTGCACAATCGGCCAGCTTGCCGGGCAATGCCGGGCCCTGGGTGATG  >  NC_002947/5589642‑5589939

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