Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,767,819 A→G 6.6% 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,767,8190AG 6.6% 222.2 / 5.0 75intergenic (+33/+36)fpvA/PP_4218TonB‑dependent outer membrane ferripyoverdine receptor FpvA/lipase/esterase family protein
Reads supporting (aligned to +/- strand):  ref base A (34/36);  new base G (2/3);  total (36/39)
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.77e-01

CGTATTCGACAAGAAGTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCG  >  NC_002947/4767675‑4767961
                                                                                                                                                |                                                                                                                                              
cGTATTCGACAAGAAGTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTg                                                                                                                                            >  7:221387/1‑149 (MQ=255)
    ttCGACAAGAAGTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCgg                                                                                                                                         <  2:378102/148‑1 (MQ=255)
        aCAAGAAGTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCTGGGTGCGGttttt                                                                                                                                    <  1:113733/149‑1 (MQ=255)
          aagaagTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAANTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCTGGGTGCGGttttttg                                                                                                                                  <  1:125029/149‑2 (MQ=255)
            gaagTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCAt                                                                                                                                <  5:324018/149‑1 (MQ=255)
            gaagTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCAt                                                                                                                                <  8:221387/149‑1 (MQ=255)
               gTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGa                                                                                                                             >  3:184801/1‑149 (MQ=255)
                tactacCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGAt                                                                                                                            >  2:115166/1‑149 (MQ=255)
                     ccAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGcc                                                                                                                         <  6:338924/147‑1 (MQ=255)
                     ccAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCt                                                                                                                        >  5:378808/1‑148 (MQ=255)
                      cAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGa                                                                                                                      <  1:115166/149‑1 (MQ=255)
                       aGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGa                                                                                                                      >  3:85316/1‑148 (MQ=255)
                         ggCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAg                                                                                                                   >  2:20276/1‑149 (MQ=255)
                         ggCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAg                                                                                                                   >  7:267429/1‑149 (MQ=255)
                            cttGCCAGCAACAACTCCTGGAGCCCGTACGACTTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCt                                                                                                                <  3:244727/148‑1 (MQ=255)
                             ttGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTc                                                                                                               <  4:85316/149‑1 (MQ=255)
                                    caacaaCTTCTGGACCCCGTACGACGTTTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTGTGATAGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACa                                                                                                        <  5:278067/149‑1 (MQ=255)
                                       gaacTCCTGGAGACGGTAGGACGTGTATGTTTACCCTCGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCATGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGAt                                                                                                     <  6:286856/148‑1 (MQ=255)
                                         acTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGAtgc                                                                                                   <  6:71083/149‑1 (MQ=255)
                                           tCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGAtgctg                                                                                                 >  5:130821/1‑149 (MQ=255)
                                             cTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGAtgctgct                                                                                               <  6:378808/149‑1 (MQ=255)
                                                 aGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAg                                                                                           <  1:8236/149‑1 (MQ=255)
                                                  gCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGt                                                                                          <  1:20276/149‑1 (MQ=255)
                                                  gCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGt                                                                                          <  1:228307/149‑1 (MQ=255)
                                                    ccGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGAtgctgctg                                                                                              >  2:73050/1‑143 (MQ=255)
                                                    ccGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGAtgctgctg                                                                                              <  1:73050/143‑1 (MQ=255)
                                                      gTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCg                                                                                                                                          >  5:311184/1‑97 (MQ=255)
                                                      gTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCg                                                                                                                                          <  6:311184/97‑1 (MQ=255)
                                                      gTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCt                                                                                      <  6:255228/149‑1 (MQ=255)
                                                       tACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCATGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTc                                                                                     <  3:192454/149‑1 (MQ=255)
                                                       tACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTc                                                                                     <  4:184801/149‑1 (MQ=255)
                                                          gacgTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGac                                                                                  <  7:355919/149‑1 (MQ=255)
                                                                   ggTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCTGGGTGCGGtttt                                                                                                                                     >  2:162562/1‑89 (MQ=255)
                                                                   ggTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCTGGGTGCGGtttt                                                                                                                                     <  1:162562/89‑1 (MQ=255)
                                                                     tGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTTCAGTGCCTCGACCACCATCGCCa                                                                       <  6:130821/149‑1 (MQ=255)
                                                                           aCGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGAtg                                                                                                    <  1:207753/114‑1 (MQ=255)
                                                                           aCGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGAtg                                                                                                    >  2:207753/1‑114 (MQ=255)
                                                                            cGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGtttttt                                                                                                                                   <  5:28683/82‑1 (MQ=255)
                                                                            cGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGtttttt                                                                                                                                   >  6:28683/1‑82 (MQ=255)
                                                                               aaCTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGAtgc                                                                                                   <  4:77277/111‑1 (MQ=255)
                                                                               aaCTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGAtgc                                                                                                   >  3:77277/1‑111 (MQ=255)
                                                                               aaCTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGAGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATc                                                             >  8:171614/1‑149 (MQ=255)
                                                                                  ttCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACa                                                                                                        <  6:258005/103‑1 (MQ=255)
                                                                                  ttCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACa                                                                                                        >  5:258005/1‑103 (MQ=255)
                                                                                  ttCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGc                                                          >  2:194486/1‑149 (MQ=255)
                                                                                     aCCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCTAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGcctc                                                       <  1:194486/149‑1 (MQ=255)
                                                                                     aCCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGcctc                                                       >  3:373515/1‑149 (MQ=255)
                                                                                         tCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCg                                                   <  4:373515/149‑1 (MQ=255)
                                                                                           aCCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCga                                                  <  2:264735/148‑1 (MQ=255)
                                                                                           aCCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCga                                                  <  8:267429/148‑1 (MQ=255)
                                                                                             cGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGc                                                                         <  2:332578/123‑1 (MQ=255)
                                                                                             cGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGc                                                                         >  1:332578/1‑123 (MQ=255)
                                                                                               ccAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGGGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATc                                             <  5:47150/149‑1 (MQ=255)
                                                                                                          ttCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAAccc                                       >  3:306025/1‑144 (MQ=255)
                                                                                                          ttCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAAccc                                       <  4:306025/144‑1 (MQ=255)
                                                                                                          ttCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATAc                                   >  7:4389/1‑148 (MQ=255)
                                                                                                              gATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCAt                               >  1:34909/1‑148 (MQ=255)
                                                                                                                   ccTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGt                          >  5:257962/1‑148 (MQ=255)
                                                                                                                    cTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCa                        >  4:86828/1‑149 (MQ=255)
                                                                                                                       aCGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATgg                     <  3:86828/149‑1 (MQ=255)
                                                                                                                        cGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATggg                    <  7:171614/149‑1 (MQ=255)
                                                                                                                              aCGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTc              >  7:38231/1‑149 (MQ=255)
                                                                                                                               cGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCAt                                     >  3:60573/1‑125 (MQ=255)
                                                                                                                               cGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCAt                                     <  4:60573/125‑1 (MQ=255)
                                                                                                                               cGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTcc             >  6:2291/1‑149 (MQ=255)
                                                                                                                                gCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTccc            >  4:81122/1‑149 (MQ=255)
                                                                                                                                 gaaaaaaCCGCACCTGGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCGGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCg           >  5:316211/2‑149 (MQ=255)
                                                                                                                                  aaaaaaCCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCAt                       <  2:39355/136‑1 (MQ=255)
                                                                                                                                  aaaaaaCCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCAt                       >  1:39355/1‑136 (MQ=255)
                                                                                                                                    aaaaCCGCACCCAGGTGCGGTTTTTTCATGGATGCGTGAAAGGCTCCCTTTCAGATGCTGCTGCAGTGCCTCGACCACCAACGCCAGGCCTTCATGGGCCTCCTAGATGATCGATCCCATACGCATGAAGTCATGGGTCACTCCCGGGa        >  7:42075/1‑149 (MQ=255)
                                                                                                                                    aaaaCCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGa        >  8:300850/1‑149 (MQ=255)
                                                                                                                                         cGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCAcc   >  2:212413/1‑149 (MQ=255)
                                                                                                                                          gCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCCCCg  >  3:303009/1‑149 (MQ=255)
                                                                                                                                          gCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCg  >  3:193546/1‑149 (MQ=255)
                                                                                                                                          gCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCg  >  5:56707/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                              
CGTATTCGACAAGAAGTACTACCAGGGCATTGCCAGCAACAACTCCTGGAGCCCGTACGACGTGTATGGTGACCCACGCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTCATGGATGCCTGAAAGGCTCCCTTACAGATGCTGCTGCAGTGCCTCGACCACCATCGCCAGGCCTTCATCGGCCTCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCG  >  NC_002947/4767675‑4767961

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