Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 172,496 A→G 8.8% Q345Q (CAA→CAG fhuB → iron(III) hydroxamate ABC transporter membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913172,4960AG 8.8% 256.0 / 12.8 91Q345Q (CAA→CAGfhuBiron(III) hydroxamate ABC transporter membrane subunit
Reads supporting (aligned to +/- strand):  ref base A (47/36);  new base G (2/6);  total (49/42)
Fisher's exact test for biased strand distribution p-value = 1.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.42e-01

CTCTGGCTGACTCGCGTGTGGATGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTAT  >  NC_000913/172347‑172645
                                                                                                                                                     |                                                                                                                                                     
ctctGGCTGACTCGCGTGTGGATGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCaa                                                                                                                                                       >  1:1405951/1‑150 (MQ=255)
    ggCTGACTCGCGTGTGGATGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACAtg                                                                                                                                                   >  2:1035187/1‑150 (MQ=255)
       tGACTCGCGTGTGGATGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGc                                                                                                                                                >  1:99197/1‑150 (MQ=255)
        gACTCGCGTGTGGATGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCt                                                                                                                                               <  1:1593438/150‑1 (MQ=255)
            cgcgTGTGGATGGAAGTGTCCAGCGGTGCTGTCACTGCGTTGATGGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAGCACATGTTGGCg                                                                                                                                           <  2:587425/150‑1 (MQ=255)
             gcgTGTGGATGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAGCACATGTTGGCGt                                                                                                                                          >  1:823611/1‑150 (MQ=255)
             gcgTGTGGATGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGt                                                                                                                                          <  1:446094/150‑1 (MQ=255)
              cgTGTGGATGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAGCACATGTTGGCGtt                                                                                                                                         <  2:823611/150‑1 (MQ=255)
               gtgtGGATGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAGCACATGTTGGCGttc                                                                                                                                        <  2:359850/150‑2 (MQ=255)
                tgtgGATGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTg                                                                                                                                       >  2:951551/1‑150 (MQ=255)
                tgtgGATGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTg                                                                                                                                       >  1:990660/1‑150 (MQ=255)
                  tggatggaAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAGCACATGTTGGCGttcagc                                                                                                                                     <  2:910091/150‑5 (MQ=255)
                  tggatggaAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAGCACATGTTGGCGttcagc                                                                                                                                     <  2:874835/150‑5 (MQ=255)
                  gggatggaAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAGCACATGTTGGCGttcagc                                                                                                                                     <  2:723005/149‑5 (MQ=255)
                   ggatggaAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGccc                                                                                                                                    >  1:751915/1‑150 (MQ=255)
                           gtgtCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCgcgggc                                                                                                                            >  2:1107134/1‑150 (MQ=255)
                            tgtCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCgcgggcg                                                                                                                           <  2:1252993/150‑1 (MQ=255)
                             gtCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCgcgggcgg                                                                                                                          >  2:164495/1‑150 (MQ=255)
                               ccACCGGTTCGGTCACTGCGGTGATGGGGGCTCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGGAGCATTAGCGCGCCGGATATGAAGGTCAACGTTCGTGTCGCGGCTGAACTCCAACATGTGCTGGCGTTTGGCCTCGCGGgcggcg                                                                                                                        <  2:1301633/150‑1 (MQ=255)
                                cACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGt                                                                                                                       >  2:1118666/1‑150 (MQ=255)
                                cACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGt                                                                                                                       >  2:1145181/1‑150 (MQ=255)
                                  ccGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGtgc                                                                                                                     <  2:130252/150‑1 (MQ=255)
                                  ccGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGtgc                                                                                                                     <  1:545797/150‑1 (MQ=255)
                                    ggTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGtgctg                                                                                                                   <  1:1143256/150‑1 (MQ=255)
                                    ggTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGtgctg                                                                                                                   >  2:988673/1‑150 (MQ=255)
                                       tCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGctgctg                                                                                                                <  2:709167/150‑1 (MQ=255)
                                       tCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGctgctg                                                                                                                >  1:863206/1‑150 (MQ=255)
                                         ggTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGtt                                                                                                              >  2:805413/1‑150 (MQ=255)
                                         ggTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGtt                                                                                                              >  1:443779/1‑150 (MQ=255)
                                          gTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTg                                                                                                             >  2:1077398/1‑150 (MQ=255)
                                           tCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGa                                                                                                            >  1:971481/1‑150 (MQ=255)
                                             aCTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATg                                                                                                          <  2:859113/150‑1 (MQ=255)
                                              cTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCCGGCGGCGTCCTGCTGTTGATgg                                                                                                         >  2:364772/1‑150 (MQ=255)
                                               tGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGc                                                                                                        <  2:1570589/150‑1 (MQ=255)
                                               tGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGc                                                                                                        <  2:472936/150‑1 (MQ=255)
                                                gCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCt                                                                                                       >  1:558225/1‑150 (MQ=255)
                                                 cGTTGATCGGTGCGCCGCTTCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCtg                                                                                                      >  1:1164596/1‑150 (MQ=255)
                                                   ttGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCtgtg                                                                                                    >  2:177223/1‑150 (MQ=255)
                                                    tGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTgtgg                                                                                                   >  1:28013/1‑150 (MQ=255)
                                                     gATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTgtggt                                                                                                  >  1:879280/1‑150 (MQ=255)
                                                        cGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTgtggtggt                                                                                               >  2:488240/1‑150 (MQ=255)
                                                          gTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGtggtggtgg                                                                                             >  2:486993/1‑150 (MQ=255)
                                                             cgccgcTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGgcgc                                                                                          <  2:751915/150‑1 (MQ=255)
                                                               ccgcTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTg                                                                                        <  1:988673/150‑1 (MQ=255)
                                                                  ctgctACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCg                                                                                     <  1:166806/150‑1 (MQ=255)
                                                                  ctgctACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCg                                                                                     <  2:663188/150‑1 (MQ=255)
                                                                    gctACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGtt                                                                                   <  1:810098/150‑1 (MQ=255)
                                                                     ctACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGttt                                                                                  >  2:832445/1‑150 (MQ=255)
                                                                       aCTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTgg                                                                                <  2:99197/150‑1 (MQ=255)
                                                                       aCTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTgg                                                                                <  1:1603908/150‑1 (MQ=255)
                                                                        cTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGt                                                                               >  2:849431/1‑150 (MQ=255)
                                                                         tgtgGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTc                                                                              >  2:1408317/1‑150 (MQ=255)
                                                                           tgGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGt                                                                            >  1:596354/1‑150 (MQ=255)
                                                                            ggCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTg                                                                           >  1:747130/1‑150 (MQ=255)
                                                                            ggCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTg                                                                           >  2:715119/1‑150 (MQ=255)
                                                                              ctaatgctgcGTAAACGCGGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGACAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGATGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGGTTTGTCGTGAt                                                                         >  2:1577955/5‑150 (MQ=255)
                                                                              cTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGAt                                                                         >  2:573242/1‑150 (MQ=255)
                                                                               tGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATg                                                                        >  2:313277/1‑150 (MQ=255)
                                                                               agTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTCGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCGGTGGTGGTGGCGCTGTCGTTTGGTCGTGATg                                                                        >  2:1503250/2‑150 (MQ=255)
                                                                                gTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATgc                                                                       <  1:716461/150‑1 (MQ=255)
                                                                                       cgTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGc                                                                <  2:1164596/150‑1 (MQ=255)
                                                                                          ttACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTgg                                                             <  2:28013/150‑1 (MQ=255)
                                                                                           tACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGGCGTGATGCGCACGGCTGGa                                                            <  1:334227/150‑1 (MQ=255)
                                                                                           tACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCACTGTCGTTTGGTCGTGATGCGCACGGCTGGa                                                            <  2:1431801/150‑1 (MQ=255)
                                                                                             cGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGTACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACg                                                          <  2:730888/150‑1 (MQ=255)
                                                                                             cGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACg                                                          <  1:164495/150‑1 (MQ=255)
                                                                                                 gcaTTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTggg                                                      <  1:612904/150‑1 (MQ=255)
                                                                                                 gcaTTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTggg                                                      <  2:443779/150‑1 (MQ=255)
                                                                                                        cgcgcCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCgg                                               >  1:981267/1‑150 (MQ=255)
                                                                                                              ggATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGtt                                         <  2:1239940/150‑1 (MQ=255)
                                                                                                               gATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTg                                        <  1:832445/150‑1 (MQ=255)
                                                                                                                   tGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCg                                    >  2:875814/1‑150 (MQ=255)
                                                                                                                   tGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGAGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGAGAGCGGGGCGTTGCTCg                                    >  2:1191995/1‑150 (MQ=255)
                                                                                                                    gAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGa                                   >  1:459919/1‑150 (MQ=255)
                                                                                                                     aaGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAg                                  <  2:55498/150‑1 (MQ=255)
                                                                                                                      aGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAgg                                 <  2:971481/150‑1 (MQ=255)
                                                                                                                       ggTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGa                                >  1:639519/1‑150 (MQ=255)
                                                                                                                         tCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGGTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGAtt                              >  1:222282/1‑150 (MQ=255)
                                                                                                                          ccaggttCGTTTCGCGGCTGCTCGGCAAGCTGTGCTGGCGTTTGCGCTCGCGGGCGGGGTGCTGCTGTTGATGGCTGTGGTGGGGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGAttt                             <  2:721759/144‑1 (MQ=255)
                                                                                                                          cAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGAttt                             >  2:859132/1‑150 (MQ=255)
                                                                                                                           aaCGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCGCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGGGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTa                            <  2:1405951/150‑1 (MQ=255)
                                                                                                                              gATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATg                         >  2:256394/1‑150 (MQ=255)
                                                                                                                                  gtgtCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCt                     <  1:572740/150‑1 (MQ=255)
                                                                                                                                    gtCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTGTGGTCGTGATGCGCACGGCTGGACGTGGGCGCGCGGGGCGTTGCTCGAGGATTTAATGCCCTgg                   >  1:1460975/1‑150 (MQ=255)
                                                                                                                                      cgcgGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGgcg                 >  1:1548219/1‑150 (MQ=255)
                                                                                                                                           cTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCa                                                                    >  2:1541653/1‑94 (MQ=255)
                                                                                                                                           cTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCa                                                                    <  1:1541653/94‑1 (MQ=255)
                                                                                                                                             aaaCGCCAGCACATGTTGGCGTTTGCCATCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCcg          >  2:1191394/2‑150 (MQ=255)
                                                                                                                                               aCGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTTATGCGCACGGCTGTACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCcgcg        <  2:1607762/150‑1 (MQ=255)
                                                                                                                                                  ccAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCTTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAAt     >  2:425674/1‑150 (MQ=255)
                                                                                                                                                  ccAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAAt     >  2:1227274/1‑150 (MQ=255)
                                                                                                                                                   cAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAAtt    >  1:787996/1‑150 (MQ=255)
                                                                                                                                                     aCATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTAt  <  2:1072408/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTCTGGCTGACTCGCGTGTGGATGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTAT  >  NC_000913/172347‑172645

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