Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 613,549 T→G 5.6% R271R (CGT→CGG fes → enterobactin/ferrienterobactin esterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913613,5490TG 5.6% 364.8 / 5.7 131R271R (CGT→CGGfesenterobactin/ferrienterobactin esterase
Reads supporting (aligned to +/- strand):  ref base T (73/48);  new base G (4/5);  total (78/53)
Fisher's exact test for biased strand distribution p-value = 4.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 132 total aligned reads displayed.

GCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAAGAGGGCGTGTT  >  NC_000913/613402‑613698
                                                                                                                                                   |                                                                                                                                                     
gCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGc                                                                                                                                                     <  2:1505009/150‑1 (MQ=255)
gCAACTTCCTCCCGCCGGGTATGTGTTGATCGACTCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTTGTGAAAGCTATTGCCCCTTTTAGCGATCGTGc                                                                                                                                                     <  1:1329955/150‑1 (MQ=255)
 cAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGcc                                                                                                                                                    <  1:979737/150‑1 (MQ=255)
  aaCTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACAGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCg                                                                                                                                                   >  1:89888/1‑150 (MQ=255)
     ttCCTCCCGGCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATc                                                                                                                                                <  1:658695/150‑1 (MQ=255)
     ttCCTCCCGCCGTGTATGTGTTGATCTACGCTATCGACACCACGCACCTCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTTAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATc                                                                                                                                                <  1:827549/150‑1 (MQ=255)
     ttCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATc                                                                                                                                                >  1:972190/1‑150 (MQ=255)
     ttCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATc                                                                                                                                                <  1:1488832/150‑1 (MQ=255)
        ctccCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCACGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCa                                                                                                                                             >  2:249039/1‑150 (MQ=255)
         tccCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCAc                                                                                                                                            >  1:585688/1‑150 (MQ=255)
           ccgccgTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCg                                                                                                                                          >  2:486994/1‑150 (MQ=255)
            cgccgTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGGCACGATCGctaaaa                                                                                                                                         <  2:548204/150‑6 (MQ=255)
             gccgTGTATGTGTTGATCGACGCTATCTACTCCACGCTCCGCGCCCACTAACTGCCGTGTAATTCTTATTTCTGGCTCGCAGTACAGCATGAGTTTTTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTg                                                                                                                                        <  1:1146240/150‑1 (MQ=255)
               cgTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGt                                                                                                                                      <  1:14296/150‑1 (MQ=255)
                 tgtATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTg                                                                                                                                    <  1:1319791/150‑1 (MQ=255)
                        gtTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCa                                                                                                                             >  2:681220/1‑150 (MQ=255)
                          tGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCaga                                                                                                                           >  1:652870/1‑150 (MQ=255)
                           gATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGGGGTTGCCGGGCagag                                                                                                                          >  2:1395037/1‑150 (MQ=255)
                               gACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGtttt                                                                                                                      <  1:226110/150‑1 (MQ=255)
                                 cGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTtgg                                                                                                                    <  2:240581/150‑1 (MQ=255)
                                 cGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTtgg                                                                                                                    >  2:494782/1‑150 (MQ=255)
                                    tATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTtggtgg                                                                                                                 >  1:1651365/1‑150 (MQ=255)
                                    aaTCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTtggtgg                                                                                                                 <  2:1622470/149‑1 (MQ=255)
                                     aTCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTggg                                                                                                                >  2:921256/1‑150 (MQ=255)
                                     aTCGACACCACGCACCGCGCCCACGAACTGCCGTGGAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGGGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTggg                                                                                                                <  1:678121/150‑1 (MQ=255)
                                      tCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGc                                                                                                               >  1:45267/1‑150 (MQ=255)
                                      tCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGc                                                                                                               >  2:1864839/1‑150 (MQ=255)
                                      tCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGGGGTTGCCGGGCAGAGTTTTGGTGGGc                                                                                                               >  2:893130/1‑150 (MQ=255)
                                         acacCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCttt                                                                                                            >  2:1621113/1‑150 (MQ=255)
                                               cGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCgcgc                                                                                                      >  2:690256/1‑150 (MQ=255)
                                                  aCCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCGCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGt                                                                                                   <  1:1249867/150‑1 (MQ=255)
                                                     gcgcCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATg                                                                                                >  1:571756/1‑150 (MQ=255)
                                                       gcCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGcc                                                                                              <  2:1109411/150‑1 (MQ=255)
                                                        cccACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCg                                                                                             >  2:698327/1‑150 (MQ=255)
                                                         ccACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGGTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCgg                                                                                            <  2:1044003/150‑1 (MQ=255)
                                                             gAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTg                                                                                        <  1:1726226/150‑1 (MQ=255)
                                                             gAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTg                                                                                        <  2:124222/150‑1 (MQ=255)
                                                                cTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCCCCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTATCCGCGCTGTATGCCGGACTGCAc                                                                                     >  2:1020332/1‑150 (MQ=255)
                                                                cTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCATATCGCAACGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGACGTACTGCAc                                                                                     >  2:311440/1‑150 (MQ=255)
                                                                     gtgtAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGAGCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGcc                                                                                <  1:1062406/150‑1 (MQ=255)
                                                                      tgtAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGTGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCt                                                                               >  1:341127/1‑150 (MQ=255)
                                                                         aaTGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGaa                                                                            <  2:1538014/150‑1 (MQ=255)
                                                                             cGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCt                                                                        >  1:1690469/1‑150 (MQ=255)
                                                                                   tCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCt                                                                  >  1:693503/1‑150 (MQ=255)
                                                                                    cTGGCTAGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTACGCGCTGTATGACGGACTGCACTGGCCTGATCGCTTTGGCtg                                                                 >  1:727202/1‑150 (MQ=255)
                                                                                      ggCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCtgtg                                                               <  2:1372690/150‑1 (MQ=255)
                                                                                       gCTCGCAGAACAGGAAGAGTTATTACCCCTGGTGAAAGCTATTGGCCCTTTTAGCGATCGGGCCGATCGCACCGTGGTTGCCGGGCTGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTgtgt                                                              <  2:1720695/150‑1 (MQ=255)
                                                                                        cTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTa                                                             >  1:362936/1‑150 (MQ=255)
                                                                                        cTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTa                                                             <  1:1435022/150‑1 (MQ=255)
                                                                                        cTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGAAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTa                                                             >  1:368344/1‑150 (MQ=255)
                                                                                          cGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTAtt                                                           <  2:640852/150‑1 (MQ=255)
                                                                                          cGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTAtt                                                           >  2:1566101/1‑150 (MQ=255)
                                                                                                  aGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAgg                                                   >  2:459800/1‑149 (MQ=255)
                                                                                                  aGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGt                                                   >  2:1190696/1‑150 (MQ=255)
                                                                                                  aGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGACAGt                                                   >  1:1030801/1‑150 (MQ=255)
                                                                                                   gCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTc                                                  <  1:1613604/150‑1 (MQ=255)
                                                                                                   gCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTc                                                  >  1:434016/1‑150 (MQ=255)
                                                                                                   gCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTc                                                  >  1:498114/1‑150 (MQ=255)
                                                                                                   gCAAGAGTTATTACCCCTGGTGAAAGCTATTGCACCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGACGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGCACGATTTGGCTGTGTATTAAGCCAGTc                                                  >  1:1336320/1‑150 (MQ=255)
                                                                                                     aaGCGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAg                                                <  1:884742/150‑1 (MQ=255)
                                                                                                      agagTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAgg                                               >  1:1493193/1‑150 (MQ=255)
                                                                                                        agTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGAt                                             >  2:771876/1‑150 (MQ=255)
                                                                                                          ttattaCCGCTGGTGTAAGCTATTGCTCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTTCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCg                                           <  2:301676/150‑1 (MQ=255)
                                                                                                          ttattaCCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGGGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCg                                           <  1:362178/150‑1 (MQ=255)
                                                                                                             ttaCCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTACGCCAGGCAGGATCGTAc                                        >  2:1062320/1‑150 (MQ=255)
                                                                                                              taCCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACt                                       <  2:567490/150‑1 (MQ=255)
                                                                                                               aCCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACtg                                      >  1:198230/1‑150 (MQ=255)
                                                                                                               aCCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGAAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACtg                                      >  1:1776485/1‑150 (MQ=255)
                                                                                                                ccccTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACtgg                                     <  1:1908338/150‑1 (MQ=255)
                                                                                                                ccccTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACtgg                                     >  1:1067443/1‑150 (MQ=255)
                                                                                                                ccccTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACtgg                                     >  1:1280853/1‑150 (MQ=255)
                                                                                                                ccccTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACtgg                                     >  1:28019/1‑150 (MQ=255)
                                                                                                                ccccTGGTGAAAGCTATTGCCCCTTTTAGCGATCGGGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACtgg                                     <  2:73899/150‑1 (MQ=255)
                                                                                                                 cccTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACtggt                                    >  2:1363141/1‑150 (MQ=255)
                                                                                                                     ggTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGcc                                <  1:249039/150‑1 (MQ=255)
                                                                                                                     ggTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCAGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGcc                                >  2:1008246/1‑150 (MQ=255)
                                                                                                                        taaaGCTATTGCCCCTTTTAGCGATCGGGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCa                             <  2:545182/149‑1 (MQ=255)
                                                                                                                        gAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCa                             >  2:1112745/1‑150 (MQ=255)
                                                                                                                          aaGCTATTGCCACTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATc                           >  2:53277/1‑150 (MQ=255)
                                                                                                                           aGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCg                          >  1:1505317/1‑150 (MQ=255)
                                                                                                                             cTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTATATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCggg                        >  1:5775/1‑150 (MQ=255)
                                                                                                                              tATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCgggg                       >  1:1425210/1‑150 (MQ=255)
                                                                                                                               aTTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTGGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCGGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGgggc                      >  1:1113593/1‑150 (MQ=255)
                                                                                                                                ttGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGgggcg                     <  2:254965/150‑1 (MQ=255)
                                                                                                                                ttGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGGGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGgggcg                     <  1:1136601/150‑1 (MQ=255)
                                                                                                                                 tGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTTCACTGGCATGAAAGATATGGCTGTTGATTAAGCCAGTCAGGAACGTAATGGTGGCCGCATCGggtcgg                    >  2:28934/1‑150 (MQ=255)
                                                                                                                                  gCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGgggcggg                   >  2:612851/1‑150 (MQ=255)
                                                                                                                                  gCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGgggcggg                   >  2:223003/1‑150 (MQ=255)
                                                                                                                                  gCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGgggcggg                   >  1:1521185/1‑150 (MQ=255)
                                                                                                                                    cccTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGgca                 >  2:110401/1‑150 (MQ=255)
                                                                                                                                    cccTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGgca                 >  1:865781/1‑150 (MQ=255)
                                                                                                                                        tttAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGCCTGCACTGGCCTGCACGCTTTGGCTGTGTAGTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCaa             >  1:289565/1‑150 (MQ=255)
                                                                                                                                        cggtgCGATCGGCACGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCaa             >  2:1260389/5‑150 (MQ=255)
                                                                                                                                        cggtgCGATCGGCACGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGAAGTGGTGGCCGCATCGGGGCGGgcagcat             >  1:1256611/5‑149 (MQ=255)
                                                                                                                                           aGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAagag          >  2:876983/1‑150 (MQ=255)
                                                                                                                                               aTCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAAGAGGGCg      <  2:1906438/150‑1 (MQ=255)
                                                                                                                                                 cGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAAGAGGGCgtg    <  1:1190696/150‑1 (MQ=255)
                                                                                                                                                  gTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAAGAGGGCgtgt   >  2:1860174/1‑150 (MQ=255)
                                                                                                                                                  gTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTCGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAAGAGGGCgtgt   <  1:997705/150‑1 (MQ=255)
                                                                                                                                                   tGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCGGGATGCCGGACTGCACTGTGCTGAACGCTATGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAAGAGGGCgggct  <  2:854212/150‑5 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAAGAGGGCGTGTT  >  NC_000913/613402‑613698

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