Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,648,463:1 +A 8.3% 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.3% 440.9 / 38.7 120coding (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 . (57/53);  new base A (3/7);  total (60/60)
Fisher's exact test for biased strand distribution p-value = 3.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.74e-01

Only 100 of 124 total aligned reads displayed.

CCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCC  >  NC_000913/2648314‑2648612
                                                                                                                                                      |                                                                                                                                                     
cccAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAACGCAGCAGACGCTCATTACCCCGGTTAATg                                                                                                                                                        <  1:965041/150‑1 (MQ=255)
  cAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGc                                                                                                                                                     <  2:219112/150‑1 (MQ=255)
  cAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGc                                                                                                                                                     >  1:631247/1‑150 (MQ=255)
     gggAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGGTTTGAGATTTTCCGCGTACGGGATCGACATCATGCGCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTc                                                                                                                                                  <  1:1737585/149‑1 (MQ=255)
     cGGAGCCTTTTGCTGTCGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTAACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTc                                                                                                                                                  <  2:504575/150‑1 (MQ=255)
     cGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTc                                                                                                                                                  <  1:408502/150‑1 (MQ=255)
     cGGAGCCTGTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCCATTAACCGGGGTAATGAGCGTc                                                                                                                                                  <  2:1547177/150‑1 (MQ=255)
     cGGAGACTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGAATGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTc                                                                                                                                                  <  1:1730562/150‑1 (MQ=255)
       gAGCCTTTTGCTGACTGGCCAACACCAGCGCAGCGTTAGACTGTACTGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCTTGCCCGGATCTTGTAAATAGCGCAGCAGACTCTCATTACCCCGGTTAATGGCGTCTg                                                                                                                                                <  1:1826861/150‑1 (MQ=255)
       gAGCCTTTTGCTGACGGGCCAACCCCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCTTACGGGATCGACATCATGCCCGGATCTTGTAAATCGCGCAGCAGACGCCATTAACCGGGGTAATGAGCGTCTg                                                                                                                                                <  1:1640781/150‑1 (MQ=255)
       gAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCCATTAACCGGGGTAATGAGCGTCTg                                                                                                                                                <  2:1087621/150‑1 (MQ=255)
       gAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCCATTAACCGGGGTAATGAGCGTCTgc                                                                                                                                               <  2:1014076/150‑2 (MQ=255)
        aGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTTTCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCCATTAACCGGGGTAATGAGCGTCTgc                                                                                                                                               <  2:1327199/150‑2 (MQ=255)
         gCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGGACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTc                                                                                                                                              <  1:920771/150‑1 (MQ=255)
          ccTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGTCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCg                                                                                                                                             <  2:384846/150‑1 (MQ=255)
          ccTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCCATTAACCGGGGTAATGAGCGTCTgctg                                                                                                                                             <  2:1145668/150‑4 (MQ=255)
          ccTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACGGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCg                                                                                                                                             <  1:952442/150‑1 (MQ=255)
           cTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCgg                                                                                                                                            <  1:1524155/150‑1 (MQ=255)
           cTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCgg                                                                                                                                            <  2:1519107/150‑1 (MQ=255)
                gCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGc                                                                                                                                       >  1:1267310/1‑150 (MQ=255)
                gCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGc                                                                                                                                       <  2:430117/150‑1 (MQ=255)
                 cTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGACCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCt                                                                                                                                      >  1:759890/1‑150 (MQ=255)
                  tGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTg                                                                                                                                     <  2:426541/150‑1 (MQ=255)
                     cGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTaa                                                                                                                                  <  1:1456552/150‑1 (MQ=255)
                       ggCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAAcc                                                                                                                                >  2:1558315/1‑150 (MQ=255)
                         ccAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCt                                                                                                                              >  1:864914/1‑150 (MQ=255)
                            acacCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCAGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGtt                                                                                                                           >  2:1473652/1‑150 (MQ=255)
                             cacCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTc                                                                                                                          <  2:1276680/150‑1 (MQ=255)
                             cacCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTc                                                                                                                          >  1:1532577/1‑150 (MQ=255)
                             cacCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTc                                                                                                                          >  2:889567/1‑150 (MQ=255)
                               ccAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGc                                                                                                                        >  2:25769/1‑150 (MQ=255)
                                 agcgcagcgTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCt                                                                                                                      >  2:1291567/1‑150 (MQ=255)
                                     cagcgTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTgcgc                                                                                                                  <  1:1530592/150‑1 (MQ=255)
                                       gcgTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCgg                                                                                                                >  2:1716551/1‑150 (MQ=255)
                                        cgTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGa                                                                                                               >  2:27218/1‑150 (MQ=255)
                                         gTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGAc                                                                                                              >  1:1685853/1‑150 (MQ=255)
                                            aGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAg                                                                                                           <  1:1753317/150‑1 (MQ=255)
                                            aGACTGTACGGCGAATTTACTGGCTTTGAGATTATCAGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTAACCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCATGCGCGGACCAg                                                                                                           >  1:1645024/1‑150 (MQ=255)
                                            aGACTGTACGGCGAATTTACTGGATTTGAGATTATCCGCGGACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAg                                                                                                           <  2:1156444/150‑1 (MQ=255)
                                             gACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAgg                                                                                                          >  2:1153940/1‑150 (MQ=255)
                                             gACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAgg                                                                                                          >  1:1624239/1‑150 (MQ=255)
                                             gACTGTACGGCGAATTTACTGGCTTTGAGATTATCAGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGAAACCCATTTCGCCTCAGCCGACAAgg                                                                                                          >  1:1624240/1‑150 (MQ=255)
                                                   aCGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATcc                                                                                                    <  2:969472/150‑1 (MQ=255)
                                                    cGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCa                                                                                                   >  1:763385/1‑150 (MQ=255)
                                                    cGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCa                                                                                                   <  2:1650840/150‑1 (MQ=255)
                                                     ggCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCAt                                                                                                  >  2:426203/1‑150 (MQ=255)
                                                     ggCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTAGGCAAGCTGTAACACTGTTCGCCTGCGCGGACCAGGAAATCCAt                                                                                                  >  1:770721/1‑150 (MQ=255)
                                                              aCTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGc                                                                                         >  1:521626/1‑150 (MQ=255)
                                                                tGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCg                                                                                       <  2:933927/150‑1 (MQ=255)
                                                                 ggCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGt                                                                                      >  2:1422647/1‑150 (MQ=255)
                                                                 ggCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGt                                                                                      >  1:1759308/1‑150 (MQ=255)
                                                                      tGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGcc                                                                                 <  2:759197/150‑1 (MQ=255)
                                                                         gaTTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGCCGCTCATTACCCCGGTTTATGGCGTCGGTCGGCACGCTGTAACCCTGTTCGCGTGCGCGGACCAGGACAGCCATCACGTAAGCCGTCAGCCAGt                                                                              >  1:1099671/1‑150 (MQ=255)
                                                                           ttATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTAc                                                                            <  1:1436704/150‑1 (MQ=255)
                                                                            tATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGCTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTAct                                                                           >  2:1517124/1‑150 (MQ=255)
                                                                             aTCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTActc                                                                          >  1:1727925/1‑150 (MQ=255)
                                                                              tCCGCGGACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTActct                                                                         <  2:873295/150‑1 (MQ=255)
                                                                               ccGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCtt                                                                        <  2:43086/150‑1 (MQ=255)
                                                                                   gTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCAGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTc                                                                    <  2:1581783/150‑1 (MQ=255)
                                                                                    tACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCa                                                                   >  1:1686082/1‑150 (MQ=255)
                                                                                     aCGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCAc                                                                  >  1:850837/1‑150 (MQ=255)
                                                                                     aCGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCAc                                                                  <  1:590465/150‑1 (MQ=255)
                                                                                      cGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCAcc                                                                 <  2:396356/150‑1 (MQ=255)
                                                                                        ggATCGACATCATGCCCGGATCTTGTTAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCAccct                                                               >  1:91003/1‑148 (MQ=255)
                                                                                        ggATCGACATCATGCCCGGATCTTGTAAATAGCGCCGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGt                                                               >  1:86758/1‑150 (MQ=255)
                                                                                           tCGACATCATGCCCGGATCTTGTAAATATCGCTGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACTCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGtttt                                                            <  1:788386/150‑1 (MQ=255)
                                                                                             gACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTa                                                          >  2:195465/1‑150 (MQ=255)
                                                                                              aCATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCTCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTAt                                                         <  1:1133488/150‑1 (MQ=255)
                                                                                               catcatGCCCGGATCTTGTAAATAGCTCAGCAGACGCTCATTACCACGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACAAGGAAATCCATCACGTAAGCCGTAAGCCAGTACTCTTGGTCACCTTTTTTAta                                                        >  1:116064/1‑149 (MQ=255)
                                                                                                 tcatGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATccc                                                      <  1:1473652/150‑1 (MQ=255)
                                                                                                 tcatGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATccc                                                      >  2:1614697/1‑150 (MQ=255)
                                                                                                      cccGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAgc                                                 >  1:1260832/1‑150 (MQ=255)
                                                                                                      cccGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAgc                                                 >  1:370819/1‑150 (MQ=255)
                                                                                                        cGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAgcgc                                               <  2:242966/150‑1 (MQ=255)
                                                                                                        cGGATCTTGTAAATAGCGCAGCAGACGCTCATTACACCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCAGCGCGGACCAGGAAATCCATCACGTAAGCCGTCATCCAGTACTCTTCGTCAACGTTTTTATACCACAgcgc                                               >  2:773963/1‑150 (MQ=255)
                                                                                                          gATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCaa                                             <  1:1422647/150‑1 (MQ=255)
                                                                                                              ttGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAgcc                                         <  1:1310647/150‑1 (MQ=255)
                                                                                                                  aaaTAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGCCCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCGTCGTCACCGGTTTTATCCCACAGCGCAAAGccgccg                                     >  1:1233483/1‑150 (MQ=255)
                                                                                                                     tAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTa                                  >  2:370275/1‑150 (MQ=255)
                                                                                                                      aGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTAt                                 >  2:338696/1‑150 (MQ=255)
                                                                                                                       gcgcAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATc                                >  1:118671/1‑150 (MQ=255)
                                                                                                                        cgcAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCa                               >  1:825637/1‑150 (MQ=255)
                                                                                                                        cgcAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATAACGTAAGCCGTCCGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATcc                               >  2:1494083/1‑149 (MQ=255)
                                                                                                                          cgacagACGCCATTAACCGGGGTAATGAGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACg                             >  1:1137971/4‑150 (MQ=255)
                                                                                                                            gcagACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGtt                           <  1:1558315/150‑1 (MQ=255)
                                                                                                                              agACGCTCATTACCCCGGGTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGGCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGc                         <  1:1648959/150‑1 (MQ=255)
                                                                                                                              agACGCCATTAACCGGGGTAATGAGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGc                         >  2:587608/1‑150 (MQ=255)
                                                                                                                                gacgccATTAACCGGGGTAATGAGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCa                        >  2:1250678/6‑150 (MQ=255)
                                                                                                                                 cGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCAtt                      >  2:1405929/1‑150 (MQ=255)
                                                                                                                                    tCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTgc                   <  2:1267310/150‑1 (MQ=255)
                                                                                                                                     cATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTgca                  <  1:535469/150‑1 (MQ=255)
                                                                                                                                     cATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTgca                  <  2:584312/150‑1 (MQ=255)
                                                                                                                                        tACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTgcagca               <  1:25769/150‑1 (MQ=255)
                                                                                                                                        tACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTgcagca               <  2:710672/150‑1 (MQ=255)
                                                                                                                                            ccGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCATACg           >  2:782376/1‑150 (MQ=255)
                                                                                                                                              ggTTACTGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACggg         <  2:1562255/150‑1 (MQ=255)
                                                                                                                                              ggTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACggg         <  1:1716551/150‑1 (MQ=255)
                                                                                                                                                   aTGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATg    >  2:783138/1‑150 (MQ=255)
                                                                                                                                                     ggCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGcc  >  2:1434924/1‑150 (MQ=255)
                                                                                                                                                     ggCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGcc  >  2:1537621/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
CCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCC  >  NC_000913/2648314‑2648612

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