Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 614,854 A→T 5.1% D233V (GAC→GTC)  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,8540AT 5.1% 268.0 / 4.1 78D233V (GAC→GTC) entFenterobactin synthase multienzyme complex component, ATP‑dependent
Reads supporting (aligned to +/- strand):  ref base A (40/34);  new base T (2/2);  total (42/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 = 9.15e-01

CCGCGAAAGCGAAGCCTGGCAGCGTGATGCGGCATTCTGGGCAGAACAGCGTCGTCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATGGACTACGCCGCCGGATTTATCTTTATGCGTCGACTGG  >  NC_000913/614708‑614994
                                                                                                                                                  |                                                                                                                                            
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                                                                                                                                           >  3:713756/1‑121 (MQ=255)
                             cGGCATTCTGGGCAGAACAGCGTCGTCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACgg                                                                                                                                           <  4:713756/121‑1 (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                                                                                     >  4:68008/1‑149 (MQ=255)
                                                       cAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGc                                                                                     <  3:97992/149‑1 (MQ=255)
                                                        aaCTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGcc                                                                                    >  2:692046/1‑149 (MQ=255)
                                                        aaCTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGcc                                                                                    >  3:703924/1‑149 (MQ=255)
                                                           tGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCtt                                                                                 >  4:429211/1‑149 (MQ=255)
                                                           tGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCtt                                                                                 <  4:363666/149‑1 (MQ=255)
                                                           tGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCtt                                                                                 <  2:289187/149‑1 (MQ=255)
                                                           tGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCtt                                                                                 >  3:311761/1‑149 (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                                     <  3:485545/149‑1 (MQ=255)
                                                                                                       aGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATGGa                                     <  4:32522/149‑1 (MQ=255)
                                                                                                       aGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATGGa                                     <  1:174251/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                                                                                       >  3:722112/1‑93 (MQ=255)
                                                                                                             tCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTa                                                                                       <  4:722112/93‑1 (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                                                                                      <  3:190001/69‑1 (MQ=255)
                                                                                                                                      tGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAg                                                                                      >  4:190001/1‑69 (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)
                                                                                                                                                  |                                                                                                                                            
CCGCGAAAGCGAAGCCTGGCAGCGTGATGCGGCATTCTGGGCAGAACAGCGTCGTCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATGGACTACGCCGCCGGATTTATCTTTATGCGTCGACTGG  >  NC_000913/614708‑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: ‑