Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 613,550 G→C 6.1% A272P (GCC→CCC) ‡ fes → enterobactin/ferrienterobactin esterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913613,5500GC 6.1% 355.1 / 12.8 115A272P (GCC→CCC) ‡fesenterobactin/ferrienterobactin esterase
Reads supporting (aligned to +/- strand):  ref base G (58/49);  new base C (3/4);  total (62/53)
Fisher's exact test for biased strand distribution p-value = 7.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.03e-01

Only 100 of 115 total aligned reads displayed.

GCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAAGAGGGCGTGT  >  NC_000913/613402‑613697
                                                                                                                                                    |                                                                                                                                                   
gCAACTTCCTCCCGCCGTGTATGTGTTGATCTACGATATCGACACAACAAAAAACAAAAAAGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGc                                                                                                                                                    >  2:44989/1‑150 (MQ=255)
gCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGCCACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGc                                                                                                                                                    >  1:1632979/1‑150 (MQ=255)
gCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGc                                                                                                                                                    >  1:1012157/1‑150 (MQ=255)
gCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGc                                                                                                                                                    <  1:846263/150‑1 (MQ=255)
     ttCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATc                                                                                                                                               <  2:52185/150‑1 (MQ=255)
      tcctccCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCg                                                                                                                                              >  1:978748/1‑150 (MQ=255)
      tcctccCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCg                                                                                                                                              <  2:1421460/150‑1 (MQ=255)
        ctccCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGGCACGATCGct                                                                                                                                            <  1:206417/150‑2 (MQ=255)
        ctccCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGGCACGATCGct                                                                                                                                            <  1:848794/150‑2 (MQ=255)
         tccCGCCGTGTATGTGTTGATCGACGATATCGCCAACACACAAAACACCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCAc                                                                                                                                           >  2:761987/1‑150 (MQ=255)
          cccGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCAcc                                                                                                                                          >  1:747691/1‑150 (MQ=255)
          cccGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCAcc                                                                                                                                          >  1:1439802/1‑150 (MQ=255)
          cccGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCAcc                                                                                                                                          >  2:1453721/1‑150 (MQ=255)
              ccgTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTgg                                                                                                                                      <  2:595588/150‑1 (MQ=255)
              ccgTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGGCACGATCGctaaaagg                                                                                                                                      <  1:1404616/150‑8 (MQ=255)
              ccgTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGGCACGATCGctaaaagg                                                                                                                                      <  1:1741790/150‑8 (MQ=255)
                     tgtgTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCgg                                                                                                                               <  1:1440897/150‑1 (MQ=255)
                     tgtgTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCgg                                                                                                                               >  1:949430/1‑150 (MQ=255)
                       tgtTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGc                                                                                                                             >  1:789026/1‑150 (MQ=255)
                       tgtTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGc                                                                                                                             >  2:903768/1‑150 (MQ=255)
                       tgtTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGc                                                                                                                             >  1:151360/1‑150 (MQ=255)
                       tgtTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGc                                                                                                                             >  2:695141/1‑150 (MQ=255)
                            aTCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGt                                                                                                                        <  1:1385750/150‑1 (MQ=255)
                              cgccGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCTATCGCACAGTGGTTGCCGGGCAGACttt                                                                                                                      >  2:1605120/4‑150 (MQ=255)
                                  gCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTtggt                                                                                                                  <  2:1741409/150‑1 (MQ=255)
                                  gCTATCGACACCACGCACCGCGCCCACGAACTGCCGGGTAAGGCTGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTtggt                                                                                                                  <  1:1687453/150‑1 (MQ=255)
                                    tATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCGCCGTGGTTGCCGGGCAGAGTTTtggtgg                                                                                                                >  1:145632/1‑150 (MQ=255)
                                     aTCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTggg                                                                                                               >  2:1345385/1‑150 (MQ=255)
                                     aTCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTggg                                                                                                               >  2:1161686/1‑150 (MQ=255)
                                      tCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGc                                                                                                              >  1:1154888/1‑150 (MQ=255)
                                       cGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCt                                                                                                             >  2:874068/1‑150 (MQ=255)
                                        gACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTATTGGTGGGCtt                                                                                                            >  2:10936/1‑150 (MQ=255)
                                             cacGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCcgc                                                                                                       >  2:1391400/1‑150 (MQ=255)
                                               cGCACCGCGCCCACGCACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACGCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCgcgc                                                                                                     <  2:1451649/150‑1 (MQ=255)
                                                gCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCt                                                                                                    <  2:789026/150‑1 (MQ=255)
                                                gCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCt                                                                                                    <  1:1282863/150‑1 (MQ=255)
                                                gCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTAACCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCt                                                                                                    <  1:1303754/150‑1 (MQ=255)
                                                 cACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTg                                                                                                   <  2:1213203/150‑1 (MQ=255)
                                                   ccGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTa                                                                                                 >  1:212352/1‑150 (MQ=255)
                                                   ccGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTa                                                                                                 <  1:374022/150‑1 (MQ=255)
                                                    cgcgCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTAt                                                                                                >  2:859376/1‑150 (MQ=255)
                                                    cgcgCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTAt                                                                                                >  2:1307543/1‑150 (MQ=255)
                                                       gcCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAGGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGcc                                                                                             >  1:303342/1‑150 (MQ=255)
                                                           aCGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCAGCGTGGTTGCCGGGCAGAGTTTTTGTGGGCTTTCCGCGCTGTATGCCGGAc                                                                                         <  2:251086/150‑1 (MQ=255)
                                                             gAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGTCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTg                                                                                       <  1:1298087/150‑1 (MQ=255)
                                                                cTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTCGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCAc                                                                                    >  2:1339120/1‑150 (MQ=255)
                                                                cTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCAc                                                                                    >  2:1035643/1‑150 (MQ=255)
                                                                     gtgtAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGACAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGcc                                                                               >  1:631981/1‑150 (MQ=255)
                                                                       gtAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTg                                                                             <  1:8349/150‑1 (MQ=255)
                                                                       gtAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGAGCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTg                                                                             <  1:1233437/150‑1 (MQ=255)
                                                                        tAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGa                                                                            <  2:743156/150‑1 (MQ=255)
                                                                          aTGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGGTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAAc                                                                          <  1:191498/150‑1 (MQ=255)
                                                                           tGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTGCGCGCTGTATGCCGGACTGCACTGGCCTGAACg                                                                         <  2:309460/150‑1 (MQ=255)
                                                                           tGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACg                                                                         >  2:836258/1‑150 (MQ=255)
                                                                            gCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGc                                                                        >  1:1403573/1‑150 (MQ=255)
                                                                                aTTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGGGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTg                                                                    <  1:1339120/150‑1 (MQ=255)
                                                                                aTTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTg                                                                    <  1:1035643/150‑1 (MQ=255)
                                                                                 tttCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCTGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTgg                                                                   >  2:715782/1‑150 (MQ=255)
                                                                                 tttCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTgg                                                                   <  1:1345385/150‑1 (MQ=255)
                                                                                  ttCTGGCTCGCAGTACAGCAAGAGTTATTCCCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCTCCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGc                                                                  <  1:138756/150‑1 (MQ=255)
                                                                                  ttCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGc                                                                  <  1:1161686/150‑1 (MQ=255)
                                                                                      ggCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCTCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCtgtg                                                              <  2:932951/150‑1 (MQ=255)
                                                                                       gCTCGCCGTACAGCAAGAGTTATTACCCCTGGTGATAGCTATTGCCCTTTTTAGCGAGCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTgtgt                                                             <  1:839245/150‑1 (MQ=255)
                                                                                       gCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTgtgt                                                             <  2:843516/150‑1 (MQ=255)
                                                                                       gCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTgtgt                                                             >  2:1352905/1‑150 (MQ=255)
                                                                                        cTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGCTCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTa                                                            >  1:317854/1‑150 (MQ=255)
                                                                                         tCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTAt                                                           <  1:761987/150‑1 (MQ=255)
                                                                                          cGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTAtt                                                          <  1:464685/150‑1 (MQ=255)
                                                                                             aGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAg                                                       <  1:1060582/150‑1 (MQ=255)
                                                                                                   gCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTc                                                 >  1:1546870/1‑150 (MQ=255)
                                                                                                    cAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCa                                                <  2:1332709/150‑1 (MQ=255)
                                                                                                        agTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCTGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGAt                                            <  1:856415/150‑1 (MQ=255)
                                                                                                          ttattaCCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCg                                          >  2:1533984/1‑150 (MQ=255)
                                                                                                           tattaCCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCAGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGt                                         >  2:1272696/1‑150 (MQ=255)
                                                                                                               aCCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCTCACCGTGGTTGCCGGGAAGAGTTTTGGTGGGCTTTCCGCGCTGTATGACGGACTGCACTGGCCTGAACGCTTGGGCGGTGGATTTAGCCTGTCCGGCTCGTACtg                                     >  1:287145/1‑150 (MQ=255)
                                                                                                                ccccTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACtgg                                    >  2:461855/1‑150 (MQ=255)
                                                                                                                ccccTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACtgg                                    >  1:1113351/1‑150 (MQ=255)
                                                                                                                  ccTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACtggtg                                  <  1:1307543/150‑1 (MQ=255)
                                                                                                                   cTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACtggtgg                                 >  2:1492271/1‑150 (MQ=255)
                                                                                                                    ttgTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTGGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGc                                <  1:1117126/148‑1 (MQ=255)
                                                                                                                    gggTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGc                                <  1:636799/149‑1 (MQ=255)
                                                                                                                      gTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCg                              >  2:47525/1‑150 (MQ=255)
                                                                                                                        gAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCa                            >  2:1617957/1‑150 (MQ=255)
                                                                                                                                 tGCCCCTTTACGCGATCGTGCCGATCGCACCGTGGTGGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGgggcgg                   <  2:918109/150‑1 (MQ=255)
                                                                                                                                  gCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCgggg                      >  1:1783452/1‑146 (MQ=255)
                                                                                                                                  gCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCgggg                      <  2:1783452/146‑1 (MQ=255)
                                                                                                                                  gCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGgggcggg                  >  2:406426/1‑150 (MQ=255)
                                                                                                                                    cccTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGgca                <  1:542747/150‑1 (MQ=255)
                                                                                                                                      cTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGgcagc              >  2:850364/1‑150 (MQ=255)
                                                                                                                                      cTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGgcagc              <  2:1113351/150‑1 (MQ=255)
                                                                                                                                        cggtgCGATCGGCACGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCaa            >  1:1822609/5‑150 (MQ=255)
                                                                                                                                         ggtgCGATCGGCACGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAag           >  1:357239/4‑150 (MQ=255)
                                                                                                                                         ggtgCGATCGACACGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAag           >  1:749804/4‑150 (MQ=255)
                                                                                                                                          tAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGCACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGTGGCCGCATCGGGGCGGGCAGCAagag         >  1:30011/1‑150 (MQ=255)
                                                                                                                                           aGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAagag         >  2:1329620/1‑150 (MQ=255)
                                                                                                                                            gCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAAGAgg        >  2:379926/1‑150 (MQ=255)
                                                                                                                                             cGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGGTGTGGGGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGGGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAAGAggg       <  1:724850/150‑1 (MQ=255)
                                                                                                                                             cGATCGTGCCGAGCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCGTTCCGCGCGGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAAGAggg       <  2:1142238/150‑1 (MQ=255)
                                                                                                                                                 cGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGCTCGTACTGGTGGCCGCATCGGGGCGGGCAGCAAGAGGGCgtg   >  1:384373/1‑150 (MQ=255)
                                                                                                                                                  gTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAAGAGGGCgtgt  <  2:1439802/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAAGAGGGCGTGT  >  NC_000913/613402‑613697

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: