Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,648,463:1 +A 8.1% coding (3825/4962 nt) yfhM ← bacterial alpha2‑macroglobulin colonization factor ECAM; anti‑host protease defense factor; periplasmic inner membrane‑anchored lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,648,4631.A 8.1% 502.2 / 43.6 135coding (3825/4962 nt)yfhMbacterial alpha2‑macroglobulin colonization factor ECAM; anti‑host protease defense factor; periplasmic inner membrane‑anchored lipoprotein
Reads supporting (aligned to +/- strand):  ref base . (62/62);  new base A (3/8);  total (65/70)
Fisher's exact test for biased strand distribution p-value = 2.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.26e-02

Only 100 of 140 total aligned reads displayed.

CCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCC  >  NC_000913/2648314‑2648612
                                                                                                                                                      |                                                                                                                                                     
cccAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATg                                                                                                                                                        <  1:1748862/150‑1 (MQ=255)
cccAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATg                                                                                                                                                        >  2:1613227/1‑150 (MQ=255)
   aGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCg                                                                                                                                                    <  2:820049/150‑1 (MQ=255)
    gCGGAGCCTTTTGCTGACGGGCCAACACCCGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGt                                                                                                                                                   <  1:1495793/150‑1 (MQ=255)
     cGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTc                                                                                                                                                  >  1:229139/1‑150 (MQ=255)
      ggAGCCTTTTGCTGACGGGCCAACACCAGCGCCGCGTAAGACTGTACAGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTTTAAATAGCGCAGCAGACGC‑CATTAACCGGGGTAATGAGCGTCt                                                                                                                                                 <  2:294593/150‑1 (MQ=255)
      ggAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATTGCGCAGCCGACGCTCAGTACCCCCGTTAATG‑GCGTCt                                                                                                                                                 >  2:1012658/1‑150 (MQ=255)
       gAGCCTTTTGCTGACGGGCCAACCCAAGCGCAGCGTAAGCCTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTCCCCCGGTTAATG‑GCGACTg                                                                                                                                                >  2:150113/1‑150 (MQ=255)
       gAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTg                                                                                                                                                <  1:1260073/150‑1 (MQ=255)
       gAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTg                                                                                                                                                >  1:510026/1‑150 (MQ=255)
        aGCCTTTTGCTGCCGGGCCCACACCAGCGCAGCGTAAGACTGTACGGCGACTTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGC‑CATTAACCGGGGTAATGAGCGTCTgc                                                                                                                                               <  2:31985/150‑2 (MQ=255)
        aGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGt                                                                                                                                               <  1:1559701/150‑1 (MQ=255)
        aGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGC‑CATTAACCGGGGTAATGAGCGTCTgc                                                                                                                                               <  2:1001433/150‑2 (MQ=255)
         gCCTTTTGCTGACGGGCCAAACCCATCGAAGCGTAATACTGTCCGTCTAATTTACTGGCTATGAGATTATCCGCGTACGGTATCGACATCATTCCCGGATCTTGTAAATAGCGCAGAAGACGC‑CATTAACCGGGGTAATGAGCGTCTgct                                                                                                                                              <  2:1304224/150‑3 (MQ=255)
          ccTTTTGCTGCCGGGCCAACACCAGCGCAGCGTAAGACTGGACGGCGAATATACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCg                                                                                                                                             <  2:1215312/150‑1 (MQ=255)
          ccTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCg                                                                                                                                             <  1:632943/150‑1 (MQ=255)
           cTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTTCGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCgg                                                                                                                                            >  1:357983/1‑150 (MQ=255)
           cTTTTGCGGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAAGTTACTGGCTTTGAGATTATCCGCGTACGGGAGCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGGCgg                                                                                                                                            <  1:291960/150‑1 (MQ=255)
            ttttGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGc                                                                                                                                           >  1:1461997/1‑150 (MQ=255)
             tttGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTGTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCa                                                                                                                                          <  2:1013549/150‑1 (MQ=255)
                   gACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGt                                                                                                                                    >  2:1572610/1‑150 (MQ=255)
                    ccGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTa                                                                                                                                   <  1:1031671/149‑1 (MQ=255)
                    aCGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTa                                                                                                                                   <  2:1809818/150‑1 (MQ=255)
                       ggacAACACCCGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGAGCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAAcc                                                                                                                                <  1:1261924/147‑1 (MQ=255)
                         ccAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGT‑ACCCTg                                                                                                                             >  2:1775034/1‑150 (MQ=255)
                          cAACCCCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTg                                                                                                                             <  2:96057/150‑1 (MQ=255)
                             cacCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTCCCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTc                                                                                                                          >  1:791771/1‑150 (MQ=255)
                             cacCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTc                                                                                                                          >  2:141009/1‑150 (MQ=255)
                             cacCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTATGAGATTATCCGCGTACGGGATCGCCATCATGACCGGCACTTGTAAATAGCTCAGCAGACGCTCATTCCCCCGGTAAATG‑GCGTCTGTCGGCACGCTCTACCCCTCTTc                                                                                                                          >  2:1157559/1‑150 (MQ=255)
                               ccAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATTGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGc                                                                                                                        >  2:158934/1‑150 (MQ=255)
                               ccAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGAAGACGCTCATTACCCCGGTTAAAG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGc                                                                                                                        >  2:1777391/1‑150 (MQ=255)
                               ccAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGACCGGATCTTGTAAATAGCGCAGCAGACGCTAATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGAAACCCTGTTCGc                                                                                                                        >  2:1770035/1‑150 (MQ=255)
                                cagcgcagcgTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGcc                                                                                                                       <  1:1321514/150‑1 (MQ=255)
                                            aGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAg                                                                                                           <  1:976860/150‑1 (MQ=255)
                                            aGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAg                                                                                                           <  1:1821668/150‑1 (MQ=255)
                                              aCTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGa                                                                                                         >  1:561236/1‑150 (MQ=255)
                                                 gTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAAt                                                                                                      <  1:420780/150‑1 (MQ=255)
                                                 gTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAAt                                                                                                      <  1:29280/150‑1 (MQ=255)
                                                    cGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCTCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCa                                                                                                   >  1:215694/1‑150 (MQ=255)
                                                    cGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCa                                                                                                   >  1:1173906/1‑150 (MQ=255)
                                                     ggCGAAGTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCTGATAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGCCCAGGAAATCCAt                                                                                                  >  2:899020/1‑150 (MQ=255)
                                                      gCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATc                                                                                                 >  1:1088028/1‑150 (MQ=255)
                                                      gCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATc                                                                                                 >  2:10936/1‑150 (MQ=255)
                                                       cGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATcc                                                                                                >  2:462337/1‑149 (MQ=255)
                                                             tACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAg                                                                                          <  2:68345/150‑1 (MQ=255)
                                                                  gCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTc                                                                                     <  2:486647/150‑1 (MQ=255)
                                                                   cTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCa                                                                                    >  1:47423/1‑150 (MQ=255)
                                                                    tttGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAg                                                                                   <  2:876835/150‑1 (MQ=255)
                                                                      tGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGcc                                                                                 >  1:690485/1‑150 (MQ=255)
                                                                      tGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGcc                                                                                 >  2:1822710/1‑150 (MQ=255)
                                                                       gagaTTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCa                                                                                <  2:696780/150‑1 (MQ=255)
                                                                         gaTTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGt                                                                              >  2:1042825/1‑150 (MQ=255)
                                                                         gaTTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAAAAGCGCAGCAGACGCTCATTACACCGGGTAATG‑GCGTCTCTCGGCACGCTGAAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGAAAGCCGTCAGCCAGt                                                                              >  1:889615/1‑150 (MQ=255)
                                                                           ttATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTAc                                                                            >  2:1693183/1‑150 (MQ=255)
                                                                             aTCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTActc                                                                          <  1:999280/150‑1 (MQ=255)
                                                                               ccGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCtt                                                                        >  2:1834853/1‑150 (MQ=255)
                                                                                  cgTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGt                                                                     >  1:741812/1‑150 (MQ=255)
                                                                                   gTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTc                                                                    <  2:510026/150‑1 (MQ=255)
                                                                                   gTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTc                                                                    <  1:717063/150‑1 (MQ=255)
                                                                                      cGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCAcc                                                                 >  2:1385744/1‑150 (MQ=255)
                                                                                       gggATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCg                                                                <  1:1303663/150‑1 (MQ=255)
                                                                                        ggATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGt                                                               >  1:1589234/1‑150 (MQ=255)
                                                                                        ggATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATACATAACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGt                                                               >  1:1155263/1‑150 (MQ=255)
                                                                                               aatcatGCGCGGTTGTTGTAATTAGCGCTGCAGACGCTCATTACCGCGTTTTTTG‑GCGTCTGTCGGCACGCTGTAACTCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATc                                                        <  1:150113/149‑1 (MQ=255)
                                                                                                atcatGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATcc                                                       <  2:1049625/150‑1 (MQ=255)
                                                                                                 tcatGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATccc                                                      <  1:1570303/150‑1 (MQ=255)
                                                                                                  catGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCca                                                     >  1:1459312/1‑150 (MQ=255)
                                                                                                    tGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCcaca                                                   <  1:1529193/150‑1 (MQ=255)
                                                                                                    tGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCcaca                                                   <  1:820766/150‑1 (MQ=255)
                                                                                                     gCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGACATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCcaccg                                                  >  1:1516300/1‑148 (MQ=255)
                                                                                                      cccGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAgc                                                 <  1:889924/150‑1 (MQ=255)
                                                                                                       ccGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAgcg                                                <  1:1733185/150‑1 (MQ=255)
                                                                                                       ccGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAgcc                                                >  2:158292/1‑149 (MQ=255)
                                                                                                               tGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAgccg                                        >  2:1129717/1‑150 (MQ=255)
                                                                                                                gTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAgccgc                                       >  2:515581/1‑150 (MQ=255)
                                                                                                                   aaTAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGt                                    >  2:1700880/1‑150 (MQ=255)
                                                                                                                    aTAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGtt                                   <  2:357983/150‑1 (MQ=255)
                                                                                                                    aTAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCGGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGtt                                   <  2:924695/150‑1 (MQ=255)
                                                                                                                     tAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTa                                  <  2:1088028/150‑1 (MQ=255)
                                                                                                                       gcgcAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATc                                >  1:1758027/1‑150 (MQ=255)
                                                                                                                         gcagcaGACGCTCATTACCCCGGTTTATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAACTCCATCACGTAAGCCGTCCGCCAGTACTCTTCGTCCGCGTTTTTATCCCCCAGCGCAAAGCCGCCGTTATctc                              >  2:282508/1‑148 (MQ=255)
                                                                                                                           gacagACGC‑CATTAACCGGGGTAATGAGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGt                            >  2:1635422/3‑150 (MQ=255)
                                                                                                                            gcagACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGGTATCACGtt                           >  2:188206/1‑150 (MQ=255)
                                                                                                                              agACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGc                         <  1:703246/150‑1 (MQ=255)
                                                                                                                              agACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTATTCTACACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGc                         >  2:245014/1‑150 (MQ=255)
                                                                                                                               gACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCCCAGCGCAAAGCCGCCGTTATCACGTTGCa                        >  2:239987/1‑150 (MQ=255)
                                                                                                                                aCGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCAt                       <  2:897723/150‑1 (MQ=255)
                                                                                                                                aCGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCAt                       <  1:1265715/150‑1 (MQ=255)
                                                                                                                                 cGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCATGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTCAGCCTTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCAtt                      <  2:1065161/150‑1 (MQ=255)
                                                                                                                                    tCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTgc                   <  1:1054087/150‑1 (MQ=255)
                                                                                                                                    tCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTgc                   <  2:637888/150‑1 (MQ=255)
                                                                                                                                    tCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTgc                   >  1:1507312/1‑150 (MQ=255)
                                                                                                                                    tCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTgc                   <  2:1291573/150‑1 (MQ=255)
                                                                                                                                     cATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCCGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTgca                  >  1:796571/1‑150 (MQ=255)
                                                                                                                                       ttACCCCGTTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTgcagc                <  2:229139/150‑1 (MQ=255)
                                                                                                                                            ccGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACg           <  2:47423/150‑1 (MQ=255)
                                                                                                                                            ccGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACg           <  2:1854133/150‑1 (MQ=255)
                                                                                                                                             gggTAAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACgg          <  1:248219/149‑1 (MQ=255)
                                                                                                                                                     g‑gCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGcc  <  1:1837693/150‑1 (MQ=255)
                                                                                                                                                     g‑gCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGCCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGcc  <  2:544833/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
CCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCC  >  NC_000913/2648314‑2648612

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