Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 613,551 C→A 6.3% A272D (GCC→GAC) ‡ fes → enterobactin/ferrienterobactin esterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913613,5510CA 6.3% 315.9 / 8.9 111A272D (GCC→GAC) ‡fesenterobactin/ferrienterobactin esterase
Reads supporting (aligned to +/- strand):  ref base C (56/47);  new base A (3/4);  total (59/52)
Fisher's exact test for biased strand distribution p-value = 7.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.54e-01

Only 100 of 116 total aligned reads displayed.

GCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAAGAGGGCGTGTTAC  >  NC_000913/613402‑613700
                                                                                                                                                     |                                                                                                                                                     
gCAACTTCCTCCCGCCGTGTATGTGTTGATCTACGATATCGACACAACAAAAAACAAAAAAGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGc                                                                                                                                                       >  2:44989/1‑150 (MQ=255)
gCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGCCACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGc                                                                                                                                                       >  1:1632979/1‑150 (MQ=255)
gCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGc                                                                                                                                                       >  1:1012157/1‑150 (MQ=255)
     ttCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATc                                                                                                                                                  <  2:52185/150‑1 (MQ=255)
      tcctccCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCg                                                                                                                                                 >  1:978748/1‑150 (MQ=255)
       cctccCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGc                                                                                                                                                >  2:374022/1‑150 (MQ=255)
        ctccCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGGCACGATCGct                                                                                                                                               <  1:206417/150‑2 (MQ=255)
        ctccCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGGCACGATCGct                                                                                                                                               <  1:848794/150‑2 (MQ=255)
         tccCGCCGTGTATTTGTTGATCGACGCTATCTACCCCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTGTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCAc                                                                                                                                              <  1:1380260/150‑1 (MQ=255)
         tccCGCCGTGTATGTGTTGATCGACGATATCGCCAACACACAAAACACCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCAc                                                                                                                                              >  2:761987/1‑150 (MQ=255)
          cccGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCAcc                                                                                                                                             >  2:1453721/1‑150 (MQ=255)
          cccGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCAcc                                                                                                                                             >  1:1439802/1‑150 (MQ=255)
              ccgTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTgg                                                                                                                                         <  2:595588/150‑1 (MQ=255)
              ccgTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGGCACGATCGctaaaagg                                                                                                                                         <  1:1404616/150‑8 (MQ=255)
                    aTGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGACGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATAACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGCTCGTGCCGATCGCACCGTGGTTGCCg                                                                                                                                   >  1:312764/1‑150 (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)
                            aTCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGt                                                                                                                           <  1:1385750/150‑1 (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:1161686/1‑150 (MQ=255)
                                     aTCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTggg                                                                                                                  >  2:1345385/1‑150 (MQ=255)
                                      tCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGc                                                                                                                 >  1:1154888/1‑150 (MQ=255)
                                       cGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCt                                                                                                                >  2:874068/1‑150 (MQ=255)
                                       cGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCGCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCt                                                                                                                <  1:578905/150‑1 (MQ=255)
                                        gACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTATTGGTGGGCtt                                                                                                               >  2:10936/1‑150 (MQ=255)
                                             cacGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCcgc                                                                                                          >  2:1361057/1‑150 (MQ=255)
                                             cacGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCcgc                                                                                                          >  2:1391400/1‑150 (MQ=255)
                                             cacGCACCGCGCCCACGAACTGACGTGTAATGCGGATTTCAGGCTCGCAGTACAGCAAGAGTTATTCCCCCTGTAGAAAGCTAATGCCCCTTTTAGCGATCGTTCCTATCTCACCGAGCTTGCCGGGCAGAGTTTTGGTGTGCTTTCcgc                                                                                                          >  2:1240058/1‑150 (MQ=255)
                                               cGCACCGCGCCCACGCACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACGCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCgcgc                                                                                                        <  2:1451649/150‑1 (MQ=255)
                                                gCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCt                                                                                                       <  1:1282863/150‑1 (MQ=255)
                                                gCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCt                                                                                                       <  2:789026/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:1307543/1‑150 (MQ=255)
                                                    cgcgCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTAt                                                                                                   >  2:859376/1‑150 (MQ=255)
                                                      cgcCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGc                                                                                                 <  1:1803873/150‑1 (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)
                                                                                  ttCTGGCTCGCAGTACAGCAAGAGTTATTCCCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCTCCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGc                                                                     <  1:138756/150‑1 (MQ=255)
                                                                                  ttCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGc                                                                     <  1:1161686/150‑1 (MQ=255)
                                                                                    cTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCtg                                                                   <  1:1453721/150‑1 (MQ=255)
                                                                                      ggCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCTCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCtgtg                                                                 <  2:932951/150‑1 (MQ=255)
                                                                                       gCTCGCCGTACAGCAAGAGTTATTACCCCTGGTGATAGCTATTGCCCTTTTTAGCGAGCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTgtgt                                                                <  1:839245/150‑1 (MQ=255)
                                                                                       gCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTgtgt                                                                >  2:1352905/1‑150 (MQ=255)
                                                                                       gCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTgtgt                                                                <  2:843516/150‑1 (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)
                                                                                                           tattaCCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGt                                            <  2:64014/150‑1 (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)
                                                                                                                              tATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCgggg                         >  2:1029234/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)
                                                                                                                                    cccTTTTAGCGATCGTGCCGATCGCCCCGTGGATGCCGGGCAGTGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTTTCTATTAAGCCAGTCAGGATCGTACTCGTGGCCTCATCGgtgcgggcc                   >  2:1577911/1‑149 (MQ=255)
                                                                                                                                    cccTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGgca                   <  1:542747/150‑1 (MQ=255)
                                                                                                                                      cTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGgcagc                 >  2:850364/1‑150 (MQ=255)
                                                                                                                                        cggtgCGATCGGCACGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCaa               >  1:1822609/5‑150 (MQ=255)
                                                                                                                                         ggtgCGATCGGCACGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAag              >  1:357239/4‑150 (MQ=255)
                                                                                                                                         ggtgCGATCGACACGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAag              >  1:749804/4‑150 (MQ=255)
                                                                                                                                          tAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAaga             >  1:809253/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)
                                                                                                                                                  gTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAAGAGGGCgtgt     <  2:1439802/150‑1 (MQ=255)
                                                                                                                                                     ccGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGGCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAAGAGGGCGTGTTAc  >  1:346051/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACTGGCCTGAACGCTTTGGCTGTGTATTAAGCCAGTCAGGATCGTACTGGTGGCCGCATCGGGGCGGGCAGCAAGAGGGCGTGTTAC  >  NC_000913/613402‑613700

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: