Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 614,850 A→C 5.1% T232P (ACC→CCC)  entF → enterobactin synthase multienzyme complex component, ATP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913614,8500AC 5.1% 274.1 / 6.1 79T232P (ACC→CCC) entFenterobactin synthase multienzyme complex component, ATP‑dependent
Reads supporting (aligned to +/- strand):  ref base A (41/34);  new base C (2/2);  total (43/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.67e-01

GCAATACCGCGAAAGCGAAGCCTGGCAGCGTGATGCGGCATTCTGGGCAGAACAGCGTCGTCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATGGACTACGCCGCCGGATTTATCTTTATGCGTCGACTGG  >  NC_000913/614702‑614994
                                                                                                                                                    |                                                                                                                                                
gCAATACCGCGAAAGCGAAGCCTGGCAGCGTGATGCGGCATTCTGGGCAGAACAGCGTCGTCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTa                                                                                                                                                  <  2:826066/149‑1 (MQ=255)
 cAATACCGCGAAAGCGAAGCCTGGCAGCGTGATGCGGCATTCTGGGCAGAACAGCGTCGTCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATCTAc                                                                                                                                                 >  1:122971/1‑149 (MQ=255)
      ccGCGAAAGCGAAGCCTGGCAGCGTGATGCGGCATTCTGGGCAGAACAGCGTCGTCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACg                                                                                                                                            >  1:525387/1‑149 (MQ=255)
       cgcgAAAGCGAAGCCTGGCAGCGTGATGCGGCATTCTGGGCAGAACAGCGTCGTCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACgg                                                                                                                                           <  2:414827/149‑1 (MQ=255)
        gcgAAAGCGAAGCCTGGCAGCGTGATGCGGCATTCTGGGCAGAACAGCGTCGTCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACggg                                                                                                                                          <  3:583216/149‑1 (MQ=255)
          gAAAGCGAAGCCTGGCAGCGTGATGCGGCATTCTGGGCAGAACAGCGTCGTCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGa                                                                                                                                        <  2:241817/149‑1 (MQ=255)
           aaaGCGAAGCCTGGCAGCGTGATGCGGCATTCTGGGCAGAACAGCGTCGTCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGaa                                                                                                                                       <  3:250430/149‑1 (MQ=255)
           aaaGCGAAGCCTGGCAGCGTGATGCGGCATTCTGGGCAGAACAGCGTCGTCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGaa                                                                                                                                       <  1:899266/149‑1 (MQ=255)
              gCGAAGCCTGGCAGCGTGATGCGGCATTCTGGGCAGAACAGCGTCGTCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTc                                                                                                                                    >  2:306231/1‑149 (MQ=255)
                      tGGCAGCGTGATGCGGCATTCTGGGCAGAACAGCGTCGTCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTCCCCGTCGGTAAATTccag                                                                                                                                 <  4:359510/144‑3 (MQ=255)
                      tGGCAGCGTGATGCGGCATTCTGGGCAGAACAGCGTCGTCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTCCCCGTCGGTAAATTccag                                                                                                                                 >  3:359510/1‑142 (MQ=255)
                        gCAGCGTGATGCGGCATTCTGGGCAGAACAGCGTCGTCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTgg                                                                                                                          >  4:139481/1‑149 (MQ=255)
                          aGCGTGATGCGGCATTCTGGGCAGAACAGTGTCGTCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCt                                                                                                                        <  3:78412/149‑1 (MQ=255)
                               gATGCGGCATTCTGGGCAGAACAGCGTCGTCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCa                                                                                                                   >  2:884327/1‑149 (MQ=255)
                                  gCGGCATTCTGGGCAGAACAGCGTCGTCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACt                                                                                                                <  2:525387/149‑1 (MQ=255)
                                   cGGCATTCTGGGCAGAACAGCGTCGTCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACgg                                                                                                                                           <  4:713756/121‑1 (MQ=255)
                                   cGGCATTCTGGGCAGAACAGCGTCGTCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACgg                                                                                                                                           >  3:713756/1‑121 (MQ=255)
                                    ggCATTCTGGGCAGAACAGCGTCGTCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACttt                                                                                                              >  3:746571/1‑149 (MQ=255)
                                       aTTCTGGGCAGAACAGCGTCGTCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAg                                                                                                           <  1:97726/149‑1 (MQ=255)
                                        ttCTGGGCAGAACAGCGTCGTCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAgg                                                                                                          >  2:50992/1‑149 (MQ=255)
                                         tCTGGGCAGAACAGCGTCGTCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGgt                                                                                                         <  3:353148/149‑1 (MQ=255)
                                            gggCAGAACAGCGTCGTCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGtgtg                                                                                                      <  4:495432/149‑1 (MQ=255)
                                              gCAGAACAGCGTCGTCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTCCCCGTCGGTAAATTcca                                                                                                                                  >  2:433076/1‑118 (MQ=255)
                                              gCAGAACAGCGTCGTCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTCCCCGTCGGTAAATTcca                                                                                                                                  <  1:433076/119‑2 (MQ=255)
                                                aGAACAGCGTCGTCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGc                                                                                                  >  4:97992/1‑149 (MQ=255)
                                                aGAACAGCGTCGTCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGc                                                                                                  <  2:339043/149‑1 (MQ=255)
                                                       cgtcgtCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGa                                                                                           >  1:761559/1‑149 (MQ=255)
                                                        gtcgtcAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGAt                                                                                          >  4:399786/1‑149 (MQ=255)
                                                           gtcAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTa                                                                                       >  4:11780/1‑149 (MQ=255)
                                                             cAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGc                                                                                     <  3:97992/149‑1 (MQ=255)
                                                             cAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGc                                                                                     >  4:68008/1‑149 (MQ=255)
                                                              aaCTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGcc                                                                                    >  3:703924/1‑149 (MQ=255)
                                                              aaCTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGcc                                                                                    >  2:692046/1‑149 (MQ=255)
                                                                 tGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCtt                                                                                 >  3:311761/1‑149 (MQ=255)
                                                                 tGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCtt                                                                                 <  2:289187/149‑1 (MQ=255)
                                                                 tGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCtt                                                                                 >  4:429211/1‑149 (MQ=255)
                                                                 tGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCtt                                                                                 <  4:363666/149‑1 (MQ=255)
                                                                  gccgccgccCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTg                                                                                <  2:536309/149‑1 (MQ=255)
                                                                        gcccgcGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTgg                                                                          <  3:68008/149‑1 (MQ=255)
                                                                        gcccgcGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTgg                                                                          <  3:545366/149‑1 (MQ=255)
                                                                               tCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCtt                                                                   >  3:359199/1‑149 (MQ=255)
                                                                                   tttCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTgg                                                               >  3:345224/1‑149 (MQ=255)
                                                                                       tcCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTgg                                                           <  4:435379/149‑1 (MQ=255)
                                                                                       tcCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTgg                                                           <  4:432426/149‑1 (MQ=255)
                                                                                       tcCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTgg                                                           >  3:32522/1‑149 (MQ=255)
                                                                                         cGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTgggg                                                         <  2:890105/149‑1 (MQ=255)
                                                                                           gCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCg                                                       >  1:816236/1‑149 (MQ=255)
                                                                                            cACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGa                                                      >  4:353690/1‑149 (MQ=255)
                                                                                             aCCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGAt                                                     <  3:454375/149‑1 (MQ=255)
                                                                                              ccTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGCTTTAGCCCTTGCGCTGGCAGCATTGTGCCTGCGGCGAtt                                                    <  2:761559/149‑1 (MQ=255)
                                                                                               cTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATtg                                                   >  4:529367/1‑149 (MQ=255)
                                                                                                     cGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAAt                                             >  3:415214/1‑149 (MQ=255)
                                                                                                           gcAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATg                                       >  1:604026/1‑149 (MQ=255)
                                                                                                             aGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATGGa                                     <  1:174251/149‑1 (MQ=255)
                                                                                                             aGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATGGa                                     <  4:32522/149‑1 (MQ=255)
                                                                                                             aGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATGGa                                     <  3:485545/149‑1 (MQ=255)
                                                                                                              gcgcCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGTCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATGGAc                                    <  1:692046/149‑1 (MQ=255)
                                                                                                              gcgcCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCTTTGTGGCTGGGGCGATTGTGCAATCGTATGGAc                                    >  4:669091/1‑149 (MQ=255)
                                                                                                                gcCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATGGACTa                                  >  2:913231/1‑149 (MQ=255)
                                                                                                                   tCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTa                                                                                       <  4:722112/93‑1 (MQ=255)
                                                                                                                   tCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTa                                                                                       >  3:722112/1‑93 (MQ=255)
                                                                                                                             ttCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATGGACTACGCCGCCGGAttt                     >  1:821581/1‑149 (MQ=255)
                                                                                                                              tCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATGGACTACGCCGCCGGATTTa                    <  4:359199/149‑1 (MQ=255)
                                                                                                                                   gcCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATGGACTACGCCGCCGGATTTATCttt               >  4:306673/1‑149 (MQ=255)
                                                                                                                                    ccTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATGGACTACGCCGCCGGATTTATCTTTa              >  2:150664/1‑149 (MQ=255)
                                                                                                                                     cTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGcc                                                                     >  3:435013/1‑93 (MQ=255)
                                                                                                                                     cTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGcc                                                                     <  4:435013/93‑1 (MQ=255)
                                                                                                                                     cTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATGGACTACGCCGCCGGATTTATCTTTAt             >  3:377857/1‑149 (MQ=255)
                                                                                                                                     cTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATGGACTACGCCGCCGGATTTATCTTTAt             <  3:529367/149‑1 (MQ=255)
                                                                                                                                       gAAACTGTAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATGGACTACGCCGCCGGATTTa                    >  2:499791/1‑140 (MQ=255)
                                                                                                                                       gAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATGGACTACGCCGCCGGATTTa                    <  1:499791/140‑1 (MQ=255)
                                                                                                                                         aaCTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATGGACTACGCCGCCGGATTTATCTTTATGCGt         >  4:771664/1‑149 (MQ=255)
                                                                                                                                           cTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATGGACTACGCCGCCGGATTTATCTTTATGCGTCg       >  3:786049/1‑149 (MQ=255)
                                                                                                                                            tGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAg                                                                                      >  4:190001/1‑69 (MQ=255)
                                                                                                                                            tGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAg                                                                                      <  3:190001/69‑1 (MQ=255)
                                                                                                                                            tGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATGGACTACGCCGCCGGATTTATCTTTATGCGTCGa      >  3:736507/1‑149 (MQ=255)
                                                                                                                                            tGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATGGACTACGCCGCCGGATTTATCTTTATGCGTCGa      >  4:870216/1‑149 (MQ=255)
                                                                                                                                                aTTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATGGACTACGCCGCCGGATTTATCTTTATGCGTCGACTgg  >  2:280275/1‑149 (MQ=255)
                                                                                                                                                aTTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATGGACTACGCCGCCGGATTTATCTTTATGCGTCGACTgg  <  2:421178/149‑1 (MQ=255)
                                                                                                                                                aTTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATGGACTACGCCGCCGGATTTATCTTTATGCGTCGACTgg  >  1:338385/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
GCAATACCGCGAAAGCGAAGCCTGGCAGCGTGATGCGGCATTCTGGGCAGAACAGCGTCGTCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATGGACTACGCCGCCGGATTTATCTTTATGCGTCGACTGG  >  NC_000913/614702‑614994

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