Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,648,452 C→A 8.2% G1279V (GGT→GTT)  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,4520CA 8.2% 381.4 / 16.2 136G1279V (GGT→GTT) yfhMbacterial alpha2‑macroglobulin colonization factor ECAM; anti‑host protease defense factor; periplasmic inner membrane‑anchored lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (64/60);  new base A (3/9);  total (67/69)
Fisher's exact test for biased strand distribution p-value = 1.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 142 total aligned reads displayed.

CACGCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGC  >  NC_000913/2648303‑2648598
                                                                                                                                                     |                                                                                                                                                   
cACGCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGATTTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGCTCTTGTAAATAGCGCAGCAGACGCTCAATAc                                                                                                                                                     >  2:206137/1‑150 (MQ=255)
 aCGCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTAcc                                                                                                                                                    <  1:1808295/150‑1 (MQ=255)
 aCGCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTAcc                                                                                                                                                    <  2:545374/150‑1 (MQ=255)
   gCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTCCGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCTTGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTAcccc                                                                                                                                                  >  1:924695/1‑150 (MQ=255)
   gCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATAATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTAcccc                                                                                                                                                  <  1:845450/150‑1 (MQ=255)
         cACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTaa                                                                                                                                            >  2:1303663/1‑150 (MQ=255)
           cccAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATg                                                                                                                                          <  1:1748862/150‑1 (MQ=255)
                cGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTc                                                                                                                                    >  1:229139/1‑150 (MQ=255)
                 ggAGCCTTTTGCTGACGGGCCAACACCAGCGCCGCGTAAGACTGTACAGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTTTAAATAGCGCAGCAGACGCCATTAACCGGGGTAATGAGCGTCt                                                                                                                                   <  2:294593/150‑1 (MQ=255)
                 ggAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATTGCGCAGCCGACGCTCAGTACCCCCGTTAATGGCGTCt                                                                                                                                   >  2:1012658/1‑150 (MQ=255)
                 ggAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCCATTAACCGGGGTAATGAGCGTCt                                                                                                                                   <  2:1542946/150‑1 (MQ=255)
                  gAGCCTTTTGCTGACGGGCCAACCCAAGCGCAGCGTAAGCCTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTCCCCCGGTTAATGGCGACTg                                                                                                                                  >  2:150113/1‑150 (MQ=255)
                  gAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTg                                                                                                                                  >  1:510026/1‑150 (MQ=255)
                  gAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTg                                                                                                                                  <  1:1260073/150‑1 (MQ=255)
                  gAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCACAGCAGACGCTCATTACCCCGGTTAATGGCGTCTg                                                                                                                                  <  2:1316204/150‑1 (MQ=255)
                   aGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTCCTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCCATTAACCGGGTTAATGAGCGTCTgc                                                                                                                                 <  2:987202/150‑2 (MQ=255)
                   aGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCCATTAACCGGGGTAATGAGCGTCTgc                                                                                                                                 <  2:1001433/150‑2 (MQ=255)
                     ccTTTTGCTGCCGGGCCAACACCAGCGCAGCGTAAGACTGGACGGCGAATATACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCg                                                                                                                               <  2:1215312/150‑1 (MQ=255)
                     ccTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCg                                                                                                                               <  1:632943/150‑1 (MQ=255)
                      cTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTTCGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCgg                                                                                                                              >  1:357983/1‑150 (MQ=255)
                      cTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCCATTAACCGGGGTAATGAGCGTCTgctgc                                                                                                                              <  1:883705/150‑5 (MQ=38)
                      cTTTTGCGGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAAGTTACTGGCTTTGAGATTATCCGCGTACGGGAGCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGGCgg                                                                                                                              <  1:291960/150‑1 (MQ=255)
                       ttttGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGc                                                                                                                             >  1:1461997/1‑150 (MQ=255)
                        tttGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTGTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCa                                                                                                                            <  2:1013549/150‑1 (MQ=255)
                              gACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGt                                                                                                                      >  2:1572610/1‑150 (MQ=255)
                               ccGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTa                                                                                                                     <  1:1031671/149‑1 (MQ=255)
                                  ggacAACACCCGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGAGCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAAcc                                                                                                                  <  1:1261924/147‑1 (MQ=255)
                                     cAACCCCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTg                                                                                                               <  2:96057/150‑1 (MQ=255)
                                        cacCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTCCCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTc                                                                                                            >  1:791771/1‑150 (MQ=255)
                                        cacCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTATGAGATTATCCGCGTACGGGATCGCCATCATGACCGGCACTTGTAAATAGCTCAGCAGACGCTCATTCCCCCGGTAAATGGCGTCTGTCGGCACGCTCTACCCCTCTTc                                                                                                            >  2:1157559/1‑150 (MQ=255)
                                          ccAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATTGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGc                                                                                                          >  2:158934/1‑150 (MQ=255)
                                          ccAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGAAGACGCTCATTACCCCGGTTAAAGGCGTCTGTCGGCACGCTGTAACCCTGTTCGc                                                                                                          >  2:1777391/1‑150 (MQ=255)
                                          ccAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGACCGGATCTTGTAAATAGCGCAGCAGACGCTAATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGAAACCCTGTTCGc                                                                                                          >  2:1770035/1‑150 (MQ=255)
                                           cagcgcagcgTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGcc                                                                                                         <  1:1321514/150‑1 (MQ=255)
                                                         aCTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGa                                                                                           >  1:561236/1‑150 (MQ=255)
                                                          cTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGaa                                                                                          >  2:490772/1‑150 (MQ=255)
                                                            gTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAAt                                                                                        <  1:29280/150‑1 (MQ=255)
                                                            gTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAAt                                                                                        <  1:420780/150‑1 (MQ=255)
                                                              aCGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATcc                                                                                      <  2:841290/150‑1 (MQ=255)
                                                               cGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCTCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCa                                                                                     >  1:215694/1‑150 (MQ=255)
                                                               cGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCa                                                                                     >  1:1173906/1‑150 (MQ=255)
                                                                ggCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGTCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCAt                                                                                    <  1:71686/150‑1 (MQ=255)
                                                                 gCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATc                                                                                   >  2:10936/1‑150 (MQ=255)
                                                                 gCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATc                                                                                   >  1:1088028/1‑150 (MQ=255)
                                                                  cGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCCTGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGACGGCACGCTGTTACCCTGTTCGCCTGCGCGGACCAGGATATCCATCa                                                                                  >  2:91363/1‑150 (MQ=255)
                                                                  cGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTCTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCa                                                                                  <  2:194293/150‑1 (MQ=255)
                                                                   gAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGCCATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACTCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCAc                                                                                 >  1:1500564/1‑150 (MQ=255)
                                                                             gCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTc                                                                       <  2:486647/150‑1 (MQ=255)
                                                                               tttGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACTCTGTAACCCTGTTCGCCTGCGCGGACCAGGAATTCCATCACGTAAGCCGTCAg                                                                     <  1:1566700/150‑1 (MQ=255)
                                                                               tttGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCCCGTAAGCCGTCAg                                                                     >  2:1654139/1‑150 (MQ=255)
                                                                               tttGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAg                                                                     <  2:876835/150‑1 (MQ=255)
                                                                               tttGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAg                                                                     <  2:158621/150‑1 (MQ=255)
                                                                                 tGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGcc                                                                   >  2:1822710/1‑150 (MQ=255)
                                                                                 tGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGcc                                                                   >  1:690485/1‑150 (MQ=255)
                                                                                  gagaTTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCa                                                                  <  2:696780/150‑1 (MQ=255)
                                                                                    gaTTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGt                                                                >  2:1042825/1‑150 (MQ=255)
                                                                                      ttATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTAc                                                              >  2:1693183/1‑150 (MQ=255)
                                                                                        aTCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTActc                                                            <  1:999280/150‑1 (MQ=255)
                                                                                         tCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTActct                                                           >  2:1831154/1‑150 (MQ=255)
                                                                                          ccGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCtt                                                          >  2:1834853/1‑150 (MQ=255)
                                                                                             cgTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGt                                                       >  1:741812/1‑150 (MQ=255)
                                                                                              gTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTc                                                      <  1:717063/150‑1 (MQ=255)
                                                                                              gTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTc                                                      <  2:510026/150‑1 (MQ=255)
                                                                                                 cGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCAcc                                                   >  2:1385744/1‑150 (MQ=255)
                                                                                                  gggATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCg                                                  <  1:1303663/150‑1 (MQ=255)
                                                                                                   ggATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGt                                                 >  2:243060/1‑150 (MQ=255)
                                                                                                   ggATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGt                                                 >  1:1589234/1‑150 (MQ=255)
                                                                                                   ggATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATACATAACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGt                                                 >  1:1155263/1‑150 (MQ=255)
                                                                                                     aTCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGttt                                               >  2:41370/1‑150 (MQ=255)
                                                                                                          aatcatGCGCGGTTGTTGTAATTAGCGCTGCAGACGCTCATTACCGCGTTTTTTGGCGTCTGTCGGCACGCTGTAACTCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATc                                          <  1:150113/149‑1 (MQ=255)
                                                                                                           atcatGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATcc                                         <  2:1049625/150‑1 (MQ=255)
                                                                                                            tcatGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATccc                                        <  1:1570303/150‑1 (MQ=255)
                                                                                                             catGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCca                                       >  1:1459312/1‑150 (MQ=255)
                                                                                                               tGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCcaca                                     <  1:1529193/150‑1 (MQ=255)
                                                                                                                gCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGACATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCcaccg                                    >  1:1516300/1‑148 (MQ=255)
                                                                                                                 cccGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAgc                                   <  1:889924/150‑1 (MQ=255)
                                                                                                                  ccGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAgcg                                  <  1:1733185/150‑1 (MQ=255)
                                                                                                                        cTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAgc                            >  1:195218/1‑150 (MQ=255)
                                                                                                                          tGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAgccg                          >  2:1129717/1‑150 (MQ=255)
                                                                                                                               aTAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGtt                     <  2:357983/150‑1 (MQ=255)
                                                                                                                               aTAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCGGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGtt                     <  2:924695/150‑1 (MQ=255)
                                                                                                                                tAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTa                    <  2:1088028/150‑1 (MQ=255)
                                                                                                                                tAGCCCAGCCTACGCTGCTTAACCCGGTTAATGGCGTCTTTCGGCACCCTTTAACCCTTTTCGCCTGCGCGGACCAGGAAATACATCACGTAAGCCTTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTa                    <  2:1797735/150‑1 (MQ=255)
                                                                                                                                  gcgcAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATc                  >  1:921413/1‑150 (MQ=255)
                                                                                                                                    gcagcaGACGCTCATTACCCCGGTTTATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAACTCCATCACGTAAGCCGTCCGCCAGTACTCTTCGTCCGCGTTTTTATCCCCCAGCGCAAAGCCGCCGTTATctc                >  2:282508/1‑148 (MQ=255)
                                                                                                                                    gcagcaGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCAc                >  2:316575/1‑150 (MQ=255)
                                                                                                                                     cagcagACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACg               >  2:690300/1‑150 (MQ=255)
                                                                                                                                      gacagACGCCATTAACCGGGGTAATGAGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGt              >  2:1635422/3‑150 (MQ=255)
                                                                                                                                         agACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCa          <  1:1775034/150‑1 (MQ=255)
                                                                                                                                          gACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCCCAGCGCAAAGCCGCCGTTATCACGTTGCa          >  2:239987/1‑150 (MQ=255)
                                                                                                                                           gacgccATTAACCGGGGTAATGAGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTAGCACGTTGCa          >  1:477173/6‑150 (MQ=255)
                                                                                                                                           gacgccATTAACCGGGGTAATGAGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGGTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCa          >  2:1710899/6‑150 (MQ=255)
                                                                                                                                           aCGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCAt         <  2:897723/150‑1 (MQ=255)
                                                                                                                                           aCGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCAt         <  1:1265715/150‑1 (MQ=255)
                                                                                                                                            cGCTCATTACCCCGGTTAATGGCGTCTGTCGGCATGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTCAGCCTTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCAtt        <  2:1065161/150‑1 (MQ=255)
                                                                                                                                               tCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTgc     >  1:1507312/1‑150 (MQ=255)
                                                                                                                                               tCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTgc     <  2:637888/150‑1 (MQ=255)
                                                                                                                                               tCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTgc     <  2:1291573/150‑1 (MQ=255)
                                                                                                                                               tCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTgc     <  1:1054087/150‑1 (MQ=255)
                                                                                                                                                cATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCCGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTgca    >  1:796571/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CACGCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGC  >  NC_000913/2648303‑2648598

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: