Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,648,458 T→G 9.6% N1277T (AAC→ACC)  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,4580TG 9.6% 309.2 / 20.9 123N1277T (AAC→ACC) yfhMbacterial alpha2‑macroglobulin colonization factor ECAM; anti‑host protease defense factor; periplasmic inner membrane‑anchored lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (56/54);  new base G (3/9);  total (59/64)
Fisher's exact test for biased strand distribution p-value = 1.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.08e-01

Only 100 of 124 total aligned reads displayed.

GCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGG  >  NC_000913/2648309‑2648605
                                                                                                                                                     |                                                                                                                                                    
gcgcACCCTGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGt                                                                                                                                                      <  2:684315/150‑1 (MQ=255)
  gcACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTa                                                                                                                                                    <  1:1003851/150‑1 (MQ=255)
  gcACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTa                                                                                                                                                    >  1:672450/1‑150 (MQ=255)
   cACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTaa                                                                                                                                                   <  1:884163/150‑1 (MQ=255)
     cccAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATg                                                                                                                                                 <  1:906559/150‑1 (MQ=255)
     cccAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAACGCAGCAGACGCTCATTACCCCGGTTAATg                                                                                                                                                 <  1:965041/150‑1 (MQ=255)
       cAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGc                                                                                                                                              <  2:219112/150‑1 (MQ=255)
          gggAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGGTTTGAGATTTTCCGCGTACGGGATCGACATCATGCGCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTc                                                                                                                                           <  1:1737585/149‑1 (MQ=255)
          cGGAGCCTTTTGCTGTCGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTAACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTc                                                                                                                                           <  2:504575/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)
              gCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCCATTAACCGGGGTAATGAGCGTCTgct                                                                                                                                       <  1:1786271/150‑3 (MQ=255)
               ccTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGTCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCg                                                                                                                                      <  2:384846/150‑1 (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)
                        gACGGGCCAACCCCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGt                                                                                                                             >  2:1537270/1‑150 (MQ=255)
                          cGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTaa                                                                                                                           <  1:1456552/150‑1 (MQ=255)
                            ggCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAAcc                                                                                                                         >  2:1558315/1‑150 (MQ=255)
                              ccAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCt                                                                                                                       >  1:819041/1‑150 (MQ=255)
                                 acacCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCAGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGtt                                                                                                                    >  2:1473652/1‑150 (MQ=255)
                                  cacCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTc                                                                                                                   <  2:1276680/150‑1 (MQ=255)
                                  cacCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTc                                                                                                                   >  2:889567/1‑150 (MQ=255)
                                  cacCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTc                                                                                                                   >  1:1532577/1‑150 (MQ=255)
                                    ccAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGc                                                                                                                 >  2:25769/1‑150 (MQ=255)
                                      agcgcagcgTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCt                                                                                                               >  2:1291567/1‑150 (MQ=255)
                                         gcagcgTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGCATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTgcg                                                                                                            >  1:147578/1‑150 (MQ=255)
                                          cagcgTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTgcgc                                                                                                           <  1:1530592/150‑1 (MQ=255)
                                           agcgTAAGACTGTACTGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGcgcg                                                                                                          <  2:751793/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)
                                                 aGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAg                                                                                                    <  1:1648229/150‑1 (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)
                                                           gCGAATTTACTGGCTTTGAGATTATCCGCGTCCGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATc                                                                                          <  1:1727638/150‑1 (MQ=255)
                                                                tttACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTa                                                                                     >  1:680488/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)
                                                                           tGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGcc                                                                          <  2:759197/150‑1 (MQ=255)
                                                                             agaTTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAg                                                                        >  2:1471773/1‑150 (MQ=255)
                                                                              gaTTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGCCGCTCATTACCCCGGTTTATGGCGTCGGTCGGCACGCTGTAACCCTGTTCGCGTGCGCGGACCAGGACAGCCATCACGTAAGCCGTCAGCCAGt                                                                       >  1:1099671/1‑150 (MQ=255)
                                                                                ttATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTAc                                                                     <  1:1436704/150‑1 (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:590465/150‑1 (MQ=255)
                                                                                          aCGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCAc                                                           >  1:850837/1‑150 (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)
                                                                                                tCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGtttt                                                     >  2:1161397/1‑150 (MQ=255)
                                                                                                  gACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTa                                                   >  2:195465/1‑150 (MQ=255)
                                                                                                   aCATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCTCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTAt                                                  <  1:1133488/150‑1 (MQ=255)
                                                                                                    catcatGCCCGGATCTTGTAAATAGCTCAGCAGACGCTCATTACCACGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACAAGGAAATCCATCACGTAAGCCGTAAGCCAGTACTCTTGGTCACCTTTTTTAta                                                 >  1:116064/1‑149 (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)
                                                                                                                     gTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAgccgc                                >  1:161622/1‑150 (MQ=255)
                                                                                                                       aaaTAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGCCCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCGTCGTCACCGGTTTTATCCCACAGCGCAAAGccgccg                              >  1:1233483/1‑150 (MQ=255)
                                                                                                                       aaaTAGCCCAGCAGACGCGCATTACCCCTGATAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCGCCGCGGACCAGGAAAACCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTGTTATCCCACAGCGCAAAGccgccg                              <  1:1396221/150‑1 (MQ=255)
                                                                                                                          tAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTa                           >  2:370275/1‑150 (MQ=255)
                                                                                                                          tAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTa                           <  1:1209703/150‑1 (MQ=255)
                                                                                                                           aGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTAt                          >  2:338696/1‑150 (MQ=255)
                                                                                                                             cgcAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACATCGCAAAGCCGCCGTTATCa                        >  2:1314434/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)
                                                                                                                                   agACGCTCATTACCCCGGGTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGGCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGc                  <  1:1648959/150‑1 (MQ=255)
                                                                                                                                   agACGCCATTAACCGGGGTAATGAGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGc                  >  2:587608/1‑150 (MQ=255)
                                                                                                                                     gacgccATTAACCGGGGTAATGAGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCa                 >  2:1250678/6‑150 (MQ=255)
                                                                                                                                          cATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTgca           <  2:584312/150‑1 (MQ=255)
                                                                                                                                          cATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTgca           <  1:535469/150‑1 (MQ=255)
                                                                                                                                            ttACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTgcagc         >  1:495545/1‑150 (MQ=255)
                                                                                                                                             tACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTgcagca        <  1:25769/150‑1 (MQ=255)
                                                                                                                                                 ccGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCATACg    >  2:782376/1‑150 (MQ=255)
                                                                                                                                                  cGGTTAATTGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGGCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACgg   <  1:981528/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACGGG  >  NC_000913/2648309‑2648605

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: