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.

CCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGCC  >  NC_000913/2648314‑2648612
                                                                                                                                                      |                                                                                                                                                     
cccAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAACGCAGCAGACGCTCATTACCCCGGTTAATg                                                                                                                                                        <  1:965041/150‑1 (MQ=255)
  cAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑Gc                                                                                                                                                     <  2:219112/150‑1 (MQ=255)
  cAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑Gc                                                                                                                                                     >  1:631247/1‑150 (MQ=255)
     gggAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGGTTTGAGATTTTCCGCGTACGGGATCGACATCATGCGCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTc                                                                                                                                                  <  1:1737585/149‑1 (MQ=255)
     cGGAGCCTTTTGCTGTCGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTAACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTc                                                                                                                                                  <  2:504575/150‑1 (MQ=255)
     cGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTc                                                                                                                                                  <  1:408502/150‑1 (MQ=255)
     cGGAGCCTGTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGC‑CATTAACCGGGGTAATGAGCGTc                                                                                                                                                  <  2:1547177/150‑1 (MQ=255)
     cGGAGACTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGAATGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTc                                                                                                                                                  <  1:1730562/150‑1 (MQ=255)
       gAGCCTTTTGCTGACTGGCCAACACCAGCGCAGCGTTAGACTGTACTGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCTTGCCCGGATCTTGTAAATAGCGCAGCAGACTCTCATTACCCCGGTTAATG‑GCGTCTg                                                                                                                                                <  1:1826861/150‑1 (MQ=255)
       gAGCCTTTTGCTGACGGGCCAACCCCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCTTACGGGATCGACATCATGCCCGGATCTTGTAAATCGCGCAGCAGACGC‑CATTAACCGGGGTAATGAGCGTCTg                                                                                                                                                <  1:1640781/150‑1 (MQ=255)
       gAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGC‑CATTAACCGGGGTAATGAGCGTCTg                                                                                                                                                <  2:1087621/150‑1 (MQ=255)
       gAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAA‑TTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGC‑CATTAACCGGGGTAATGAGCGTCTgc                                                                                                                                               <  2:1014076/150‑2 (MQ=255)
        aGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTTTCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGC‑CATTAACCGGGGTAATGAGCGTCTgc                                                                                                                                               <  2:1327199/150‑2 (MQ=255)
         gCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGGACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTc                                                                                                                                              <  1:920771/150‑1 (MQ=255)
          ccTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGTCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCg                                                                                                                                             <  2:384846/150‑1 (MQ=255)
          ccTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGC‑CATTAACCGGGGTAATGAGCGTCTgctg                                                                                                                                             <  2:1145668/150‑4 (MQ=255)
          ccTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACGGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCg                                                                                                                                             <  1:952442/150‑1 (MQ=255)
           cTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCgg                                                                                                                                            <  1:1524155/150‑1 (MQ=255)
           cTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCgg                                                                                                                                            <  2:1519107/150‑1 (MQ=255)
                gCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGc                                                                                                                                       >  1:1267310/1‑150 (MQ=255)
                gCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGc                                                                                                                                       <  2:430117/150‑1 (MQ=255)
                 cTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGACCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCt                                                                                                                                      >  1:759890/1‑150 (MQ=255)
                  tGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTg                                                                                                                                     <  2:426541/150‑1 (MQ=255)
                     cGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTaa                                                                                                                                  <  1:1456552/150‑1 (MQ=255)
                       ggCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAAcc                                                                                                                                >  2:1558315/1‑150 (MQ=255)
                         ccAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCt                                                                                                                              >  1:864914/1‑150 (MQ=255)
                            acacCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCAGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGtt                                                                                                                           >  2:1473652/1‑150 (MQ=255)
                             cacCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTc                                                                                                                          <  2:1276680/150‑1 (MQ=255)
                             cacCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTc                                                                                                                          >  1:1532577/1‑150 (MQ=255)
                             cacCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTc                                                                                                                          >  2:889567/1‑150 (MQ=255)
                               ccAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGc                                                                                                                        >  2:25769/1‑150 (MQ=255)
                                 agcgcagcgTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCt                                                                                                                      >  2:1291567/1‑150 (MQ=255)
                                     cagcgTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTgcgc                                                                                                                  <  1:1530592/150‑1 (MQ=255)
                                       gcgTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCgg                                                                                                                >  2:1716551/1‑150 (MQ=255)
                                        cgTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGa                                                                                                               >  2:27218/1‑150 (MQ=255)
                                         gTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGAc                                                                                                              >  1:1685853/1‑150 (MQ=255)
                                            aGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAg                                                                                                           <  1:1753317/150‑1 (MQ=255)
                                            aGACTGTACGGCGAATTTACTGGCTTTGAGATTATCAGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTAACCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCATGCGCGGACCAg                                                                                                           >  1:1645024/1‑150 (MQ=255)
                                            aGACTGTACGGCGAATTTACTGGATTTGAGATTATCCGCGGACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAg                                                                                                           <  2:1156444/150‑1 (MQ=255)
                                             gACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAgg                                                                                                          >  2:1153940/1‑150 (MQ=255)
                                             gACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAgg                                                                                                          >  1:1624239/1‑150 (MQ=255)
                                             gACTGTACGGCGAATTTACTGGCTTTGAGATTATCAGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGAAACCCATTTCGCCTCAGCCGACAAgg                                                                                                          >  1:1624240/1‑150 (MQ=255)
                                                   aCGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATcc                                                                                                    <  2:969472/150‑1 (MQ=255)
                                                    cGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCa                                                                                                   >  1:763385/1‑150 (MQ=255)
                                                    cGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCa                                                                                                   <  2:1650840/150‑1 (MQ=255)
                                                     ggCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCAt                                                                                                  >  2:426203/1‑150 (MQ=255)
                                                     ggCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTAGGCAAGCTGTAACACTGTTCGCCTGCGCGGACCAGGAAATCCAt                                                                                                  >  1:770721/1‑150 (MQ=255)
                                                              aCTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGc                                                                                         >  1:521626/1‑150 (MQ=255)
                                                                tGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCg                                                                                       <  2:933927/150‑1 (MQ=255)
                                                                 ggCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGt                                                                                      >  2:1422647/1‑150 (MQ=255)
                                                                 ggCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGt                                                                                      >  1:1759308/1‑150 (MQ=255)
                                                                      tGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGcc                                                                                 <  2:759197/150‑1 (MQ=255)
                                                                         gaTTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGCCGCTCATTACCCCGGTTTATG‑GCGTCGGTCGGCACGCTGTAACCCTGTTCGCGTGCGCGGACCAGGACAGCCATCACGTAAGCCGTCAGCCAGt                                                                              >  1:1099671/1‑150 (MQ=255)
                                                                           ttATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTAc                                                                            <  1:1436704/150‑1 (MQ=255)
                                                                            tATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGCTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTAct                                                                           >  2:1517124/1‑150 (MQ=255)
                                                                             aTCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTActc                                                                          >  1:1727925/1‑150 (MQ=255)
                                                                              tCCGCGGACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTActct                                                                         <  2:873295/150‑1 (MQ=255)
                                                                               ccGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCtt                                                                        <  2:43086/150‑1 (MQ=255)
                                                                                   gTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCAGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTc                                                                    <  2:1581783/150‑1 (MQ=255)
                                                                                    tACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCa                                                                   >  1:1686082/1‑150 (MQ=255)
                                                                                     aCGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCAc                                                                  >  1:850837/1‑150 (MQ=255)
                                                                                     aCGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCAc                                                                  <  1:590465/150‑1 (MQ=255)
                                                                                      cGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCAcc                                                                 <  2:396356/150‑1 (MQ=255)
                                                                                        ggATCGACATCATGCCCGGATCTTGTTAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCAccct                                                               >  1:91003/1‑148 (MQ=255)
                                                                                        ggATCGACATCATGCCCGGATCTTGTAAATAGCGCCGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGt                                                               >  1:86758/1‑150 (MQ=255)
                                                                                           tCGACATCATGCCCGGATCTTGTAAATATCGCTGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACTCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGtttt                                                            <  1:788386/150‑1 (MQ=255)
                                                                                             gACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTa                                                          >  2:195465/1‑150 (MQ=255)
                                                                                              aCATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCTCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTAt                                                         <  1:1133488/150‑1 (MQ=255)
                                                                                               catcatGCCCGGATCTTGTAAATAGCTCAGCAGACGCTCATTACCACGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACAAGGAAATCCATCACGTAAGCCGTAAGCCAGTACTCTTGGTCACCTTTTTTAta                                                        >  1:116064/1‑149 (MQ=255)
                                                                                                 tcatGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATccc                                                      <  1:1473652/150‑1 (MQ=255)
                                                                                                 tcatGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATccc                                                      >  2:1614697/1‑150 (MQ=255)
                                                                                                      cccGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAgc                                                 >  1:1260832/1‑150 (MQ=255)
                                                                                                      cccGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAgc                                                 >  1:370819/1‑150 (MQ=255)
                                                                                                        cGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAgcgc                                               <  2:242966/150‑1 (MQ=255)
                                                                                                        cGGATCTTGTAAATAGCGCAGCAGACGCTCATTACACCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCAGCGCGGACCAGGAAATCCATCACGTAAGCCGTCATCCAGTACTCTTCGTCAACGTTTTTATACCACAgcgc                                               >  2:773963/1‑150 (MQ=255)
                                                                                                          gATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCaa                                             <  1:1422647/150‑1 (MQ=255)
                                                                                                              ttGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAgcc                                         <  1:1310647/150‑1 (MQ=255)
                                                                                                                  aaaTAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGCCCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCGTCGTCACCGGTTTTATCCCACAGCGCAAAGccgccg                                     >  1:1233483/1‑150 (MQ=255)
                                                                                                                     tAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTa                                  >  2:370275/1‑150 (MQ=255)
                                                                                                                      aGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTAt                                 >  2:338696/1‑150 (MQ=255)
                                                                                                                       gcgcAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATc                                >  1:118671/1‑150 (MQ=255)
                                                                                                                        cgcAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCa                               >  1:825637/1‑150 (MQ=255)
                                                                                                                        cgcAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATAACGTAAGCCGTCCGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATcc                               >  2:1494083/1‑149 (MQ=255)
                                                                                                                          cgacagACGC‑CATTAACCGGGGTAATGAGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACg                             >  1:1137971/4‑150 (MQ=255)
                                                                                                                            gcagACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGtt                           <  1:1558315/150‑1 (MQ=255)
                                                                                                                              agACGCTCATTACCCCGGGTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGGCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGc                         <  1:1648959/150‑1 (MQ=255)
                                                                                                                              agACGC‑CATTAACCGGGGTAATGAGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGc                         >  2:587608/1‑150 (MQ=255)
                                                                                                                                gacgccATTAACCGGGGTAATGAGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCa                        >  2:1250678/6‑150 (MQ=255)
                                                                                                                                 cGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCAtt                      >  2:1405929/1‑150 (MQ=255)
                                                                                                                                    tCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTgc                   <  2:1267310/150‑1 (MQ=255)
                                                                                                                                     cATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTgca                  <  1:535469/150‑1 (MQ=255)
                                                                                                                                     cATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTgca                  <  2:584312/150‑1 (MQ=255)
                                                                                                                                        tACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTgcagca               <  1:25769/150‑1 (MQ=255)
                                                                                                                                        tACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTgcagca               <  2:710672/150‑1 (MQ=255)
                                                                                                                                            ccGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCATACg           >  2:782376/1‑150 (MQ=255)
                                                                                                                                              ggTTACTG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACggg         <  2:1562255/150‑1 (MQ=255)
                                                                                                                                              ggTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACggg         <  1:1716551/150‑1 (MQ=255)
                                                                                                                                                   aTG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATg    >  2:783138/1‑150 (MQ=255)
                                                                                                                                                     g‑gCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGcc  >  2:1434924/1‑150 (MQ=255)
                                                                                                                                                     g‑gCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGGAAATGcc  >  2:1537621/1‑150 (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 < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑