Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,769,786 A→G 7.1% intergenic (+33/+36) fpvA → / ← PP_4218 TonB‑dependent outer membrane ferripyoverdine receptor FpvA/lipase/esterase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,769,7860AG 7.1% 956.7 / 38.8 315intergenic (+33/+36)fpvA/PP_4218TonB‑dependent outer membrane ferripyoverdine receptor FpvA/lipase/esterase family protein
Reads supporting (aligned to +/- strand):  ref base A (146/140);  new base G (12/10);  total (160/155)
Fisher's exact test for biased strand distribution p-value = 8.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.85e-01

Only 100 of 325 total aligned reads displayed.

AACAACGTATTCGACAAGAAGTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTG  >  NC_002947/4769637‑4769935
                                                                                                                                                     |                                                                                                                                                     
 acaacGTATTCGACAAGAAGTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAg                                                                                                                                                      >  2:7936798/1‑150 (MQ=255)
 acaacGTATTCGACAAGAAGTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAg                                                                                                                                                      <  2:7472187/150‑1 (MQ=255)
 acaacGTATTCGACAAGAAGTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAg                                                                                                                                                      <  2:6256288/150‑1 (MQ=255)
   aacGTATTCGACAAGAAGTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGt                                                                                                                                                    >  2:2246662/1‑150 (MQ=255)
           cGACAAGAAGTACTACCAGGGCATAGCCAGCAACAACTCCCGGAGCCCGTACGACGTGTATGGTGACCCACGCAACATCACCATCACCGCCAAGTACAGCTTCTGATCGCCAGACGCTGAACGCAACAAACCGCACCCAGCAGCGGtatt                                                                                                                                            <  2:4378205/150‑1 (MQ=255)
            gACAAGAAGTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGttttt                                                                                                                                           >  1:2629911/1‑150 (MQ=255)
             aCAAGAAGTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGtttttt                                                                                                                                          <  1:7673256/150‑1 (MQ=255)
               aagaagTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCTGGGTGCGGttttttgc                                                                                                                                        >  2:6723698/1‑148 (MQ=255)
               aagaagTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCa                                                                                                                                        >  2:6910733/1‑150 (MQ=255)
               aagaagTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCa                                                                                                                                        >  2:234182/1‑150 (MQ=255)
                 gaagTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATg                                                                                                                                      <  1:1166075/150‑1 (MQ=255)
                        tacCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCt                                                                                                                               >  1:4731588/1‑150 (MQ=255)
                         acCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTg                                                                                                                              >  1:2140562/1‑150 (MQ=255)
                           cAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGaa                                                                                                                            >  2:5765286/1‑150 (MQ=255)
                            aGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGaaa                                                                                                                           <  2:2370206/150‑1 (MQ=255)
                              gtcATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCATCTTCTCCATCACCGACAAGTACTGATTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGTTTTTTTCATGGATGCCTTAAAgg                                                                                                                         <  2:5700165/148‑1 (MQ=255)
                              ggCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAgg                                                                                                                         <  2:8916926/150‑1 (MQ=255)
                                 aTTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTc                                                                                                                      <  1:4996814/150‑1 (MQ=255)
                                  ttGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTcc                                                                                                                     <  2:3536214/150‑1 (MQ=255)
                                   tGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCACACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTccc                                                                                                                    <  1:915752/150‑1 (MQ=255)
                                    gCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCt                                                                                                                   <  1:8607479/150‑1 (MQ=255)
                                     ccAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCtt                                                                                                                  >  1:729967/1‑150 (MQ=255)
                                      cAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTa                                                                                                                 <  1:4244706/150‑1 (MQ=255)
                                       aGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGATCGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTAc                                                                                                                <  2:2986409/150‑1 (MQ=255)
                                         caacaaCTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAg                                                                                                              >  2:7333024/1‑150 (MQ=255)
                                           acaacTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGAt                                                                                                            >  2:3143266/1‑150 (MQ=255)
                                             aacTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGAtgc                                                                                                          >  2:8803362/1‑150 (MQ=255)
                                             aacTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGAtgc                                                                                                          >  2:592206/1‑150 (MQ=255)
                                               cTCCTGGAGCCCGAACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTATTTTCATGGATGCCTGAAAGGCTCCCTTACAGAtgctg                                                                                                        <  2:4788303/150‑1 (MQ=255)
                                                tCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGAtgctgc                                                                                                       >  2:7153460/1‑150 (MQ=255)
                                                  cTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGTGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGAtgctgctg                                                                                                     >  1:3146651/1‑150 (MQ=255)
                                                  cTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGAtgctgctg                                                                                                     <  2:1811034/150‑1 (MQ=255)
                                                    ggAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCa                                                                                                   <  2:8240160/150‑1 (MQ=255)
                                                         ccGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGcc                                                                                              <  1:4656256/150‑1 (MQ=255)
                                                             acgacgTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGa                                                                                          <  2:4958009/150‑1 (MQ=255)
                                                                acgTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGacca                                                                                       <  1:3523149/150‑1 (MQ=255)
                                                                acgTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGacca                                                                                       <  1:1032917/150‑1 (MQ=255)
                                                                    gtATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAATCCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCAt                                                                                   >  1:1639720/1‑150 (MQ=255)
                                                                     tATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATc                                                                                  <  2:7521269/150‑1 (MQ=255)
                                                                          tGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAg                                                                             <  1:548329/150‑1 (MQ=255)
                                                                          tGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAg                                                                             <  1:2792076/150‑1 (MQ=255)
                                                                              ccACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCt                                                                         >  2:2176668/1‑150 (MQ=255)
                                                                                aCGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTc                                                                       >  2:4703478/1‑150 (MQ=255)
                                                                                aCGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTc                                                                       >  2:7997236/1‑150 (MQ=255)
                                                                                    aaCTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAACAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCg                                                                   <  2:4006581/150‑1 (MQ=255)
                                                                                    aaCTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAACAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACCGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCg                                                                   <  2:3214375/150‑1 (MQ=255)
                                                                                    aaCTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCg                                                                   <  2:7253718/150‑1 (MQ=255)
                                                                                    aaCTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCg                                                                   <  2:2140654/150‑1 (MQ=255)
                                                                                       ttCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGcc                                                                >  1:5577576/1‑150 (MQ=255)
                                                                                         cACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGcctc                                                              >  2:2486472/1‑150 (MQ=255)
                                                                                           ccATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGcctcct                                                            <  1:2665127/150‑1 (MQ=255)
                                                                                           ccATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGcctcct                                                            <  2:6667401/150‑1 (MQ=255)
                                                                                           ccATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGcctcct                                                            <  1:3413971/150‑1 (MQ=255)
                                                                                              tCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAACGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCAGCGGCCTCCTCga                                                         >  2:8721601/1‑150 (MQ=255)
                                                                                              tCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCga                                                         <  2:5288563/150‑1 (MQ=255)
                                                                                               cACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCgat                                                        >  1:9102826/1‑150 (MQ=255)
                                                                                                 ccGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTATTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCgatga                                                      <  2:5482985/150‑1 (MQ=255)
                                                                                                      aaGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGaa                                                 >  2:7374896/1‑150 (MQ=255)
                                                                                                      aaGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGaa                                                 <  1:247041/150‑1 (MQ=255)
                                                                                                      aaGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGaa                                                 <  2:3686604/150‑1 (MQ=255)
                                                                                                      aaGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGaa                                                 <  1:8980533/150‑1 (MQ=255)
                                                                                                       aGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAAc                                                <  2:1164550/150‑1 (MQ=255)
                                                                                                        gTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAAcc                                               <  1:327384/150‑1 (MQ=255)
                                                                                                        gTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAAcc                                               <  1:3748986/150‑1 (MQ=255)
                                                                                                          aCAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCa                                             <  2:3591829/150‑1 (MQ=255)
                                                                                                            aGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCCTCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATa                                           >  1:6720901/1‑150 (MQ=255)
                                                                                                            aGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATa                                           >  1:2551633/1‑150 (MQ=255)
                                                                                                             gCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATAc                                          <  2:2680118/150‑1 (MQ=255)
                                                                                                               ttCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGc                                        >  2:765529/1‑150 (MQ=255)
                                                                                                                  tGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGACTCCTCGATGATCGAACCCATACGCATg                                     >  1:4148578/1‑150 (MQ=255)
                                                                                                                    aTCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGaa                                   <  2:6787909/150‑1 (MQ=255)
                                                                                                                       gCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTc                                >  2:1689318/1‑150 (MQ=255)
                                                                                                                        ccTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGGCa                               >  1:546838/1‑150 (MQ=255)
                                                                                                                          tGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGc                                                                                <  1:3645921/99‑1 (MQ=255)
                                                                                                                          tGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTc                                                                       >  2:4226825/1‑108 (MQ=255)
                                                                                                                           gACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATgg                            <  2:1290352/150‑1 (MQ=255)
                                                                                                                                tGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCAc                       <  1:7225413/150‑1 (MQ=255)
                                                                                                                                  aaCGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTc                     >  1:1986350/1‑150 (MQ=255)
                                                                                                                                     gCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCg                  >  2:6007147/1‑150 (MQ=255)
                                                                                                                                     gCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCg                  <  1:1784682/150‑1 (MQ=255)
                                                                                                                                     gCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCg                  <  1:699283/150‑1 (MQ=255)
                                                                                                                                      cAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCgg                 <  1:728070/150‑1 (MQ=255)
                                                                                                                                       aaaaaaCCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCAt                              >  1:6555291/1‑136 (MQ=255)
                                                                                                                                       aaaaaaCCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCAt                              <  2:6555291/136‑1 (MQ=255)
                                                                                                                                       aaaaaaCCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCggg                >  2:2017542/1‑150 (MQ=255)
                                                                                                                                       aaaaaaCCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCggg                >  1:1424375/1‑150 (MQ=255)
                                                                                                                                       aaaaaaCCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCggg                >  1:4808690/1‑150 (MQ=255)
                                                                                                                                         aaaaCCGCACCTGGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGAt              >  1:8025293/1‑150 (MQ=255)
                                                                                                                                         aaaaCCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGag              >  1:2963017/1‑149 (MQ=255)
                                                                                                                                         aaaaCCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGAt              >  1:7580598/1‑150 (MQ=255)
                                                                                                                                           aaCCGCACCTGGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCa            >  1:839801/1‑150 (MQ=255)
                                                                                                                                           aaCCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCa            >  1:5721642/1‑150 (MQ=255)
                                                                                                                                             ccGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCAcc          >  1:4458614/1‑150 (MQ=255)
                                                                                                                                             ccGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCAcc          <  2:819352/150‑1 (MQ=255)
                                                                                                                                              cGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCg         >  2:8561679/1‑150 (MQ=255)
                                                                                                                                               gCACCCTGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGc        <  1:5847557/150‑1 (MQ=255)
                                                                                                                                                 aCCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCaa      <  1:234182/150‑1 (MQ=255)
                                                                                                                                                 aCCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCaa      >  1:6460678/1‑150 (MQ=255)
                                                                                                                                                    cAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCt   >  1:5186978/1‑150 (MQ=255)
                                                                                                                                                     aGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCTAGCTg  >  2:3830744/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AACAACGTATTCGACAAGAAGTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTG  >  NC_002947/4769637‑4769935

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: