Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,650,498 C→T 7.1% G597D (GGT→GAT) ‡ 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,650,4980CT 7.1% 314.8 / 10.9 98G597D (GGT→GAT) ‡yfhMbacterial alpha2‑macroglobulin colonization factor ECAM; anti‑host protease defense factor; periplasmic inner membrane‑anchored lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (47/44);  new base T (3/4);  total (50/48)
Fisher's exact test for biased strand distribution p-value = 7.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

GTTACTGCCTTCATCAACAATGGGCTGTTTCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCGCGCCCTTTATCATCCAGCGTCAACTGAACTTCATCCAGCGTGCGGGAAAGATTTTCGGCAGC  >  NC_000913/2650353‑2650646
                                                                                                                                                 |                                                                                                                                                    
gTTACTGCCTTCATCAACAATGGGCTGTTTCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGc                                                                                                                                                   >  2:1142739/1‑149 (MQ=255)
  tACTGCCTTCATCAACAATGGGCTGTTTCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGccc                                                                                                                                                 <  1:1142739/149‑1 (MQ=255)
     tGCCTTCATCAACAATGGGCTGTTTCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCATCGCGGCGCGTCACCGGGCGACCGCCCGAt                                                                                                                                              >  3:91841/1‑149 (MQ=255)
      gCCTTCATCAACAATGGGCTGTTTCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGGTCGCCCggtg                                                                                                                                             <  2:843933/149‑4 (MQ=255)
       ccTTCATCAACAATGGGCTGTTTCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATGGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTc                                                                                                                                            >  4:1198257/1‑149 (MQ=255)
         ttCATCAACAATGGGCTGTTTCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCa                                                                                                                                          >  1:71171/1‑149 (MQ=255)
          tcatcaACAATGGGCTGTTTCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGcc                                                                                                                                                  >  4:1037662/1‑140 (MQ=255)
          tcatcaACAATGGGCTGTTTCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGcc                                                                                                                                                  <  3:1037662/140‑1 (MQ=255)
          tcatcaACAATGGGCTGTTTCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGAtt                                                                                                                                             <  1:852712/145‑1 (MQ=255)
          tcatcaACAATGGGCTGTTTCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGAtt                                                                                                                                             >  2:852712/1‑145 (MQ=255)
            atcaACAATGGGCTGTTTCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCa                                                                                                                                       >  2:1311126/1‑149 (MQ=255)
             tcaACAATGGGCTGTTTCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGGTCGCCCggt                                                                                                                                              <  4:1249704/141‑3 (MQ=255)
             tcaACAATGGGCTGTTTCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGGTCGCCCggt                                                                                                                                              >  3:1249704/1‑139 (MQ=255)
               aacaaTGGGCTGTTTCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAAc                                                                                                                                    >  4:413100/1‑149 (MQ=255)
               aacaaTGGGCTGTTTCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAAc                                                                                                                                    >  3:94581/1‑149 (MQ=255)
                 caaTGGGCTGTTTCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTa                                                                                                                                  >  4:562207/1‑149 (MQ=255)
                              ccccGTGCTGTCAGTACGATAATCGTAAACCGATATCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATaa                                                                                                                     <  3:1198257/147‑1 (MQ=255)
                              cACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCATCGCGGCGCGTCACCGGGCGACCGCCCGATTGCAGCAAACTACCCTGGAAAATaa                                                                                                                     <  4:91841/149‑1 (MQ=255)
                               aCCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAAc                                                                                                                    <  1:1311126/149‑1 (MQ=255)
                                ccGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGaaaa                                                                                                                        >  3:1125150/1‑144 (MQ=255)
                                ccGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGaaaa                                                                                                                        <  4:1125150/144‑1 (MQ=255)
                                     cTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTAc                                                                                                                                 >  3:555306/1‑130 (MQ=255)
                                     cTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTAc                                                                                                                                 <  4:555306/130‑1 (MQ=255)
                                     cTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTaa                                                                                                              >  1:569234/1‑149 (MQ=255)
                                      tGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGaaaa                                                                                                                        <  4:1014962/138‑1 (MQ=255)
                                      tGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGaaaa                                                                                                                        >  3:1014962/1‑138 (MQ=255)
                                        tCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATgg                                                                                                           <  1:723111/149‑1 (MQ=255)
                                        tCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATgg                                                                                                           >  2:358088/1‑149 (MQ=255)
                                        tCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATgg                                                                                                           >  2:723111/1‑149 (MQ=255)
                                        tCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATgg                                                                                                           <  2:569234/149‑1 (MQ=255)
                                         cAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATggg                                                                                                          <  3:562207/149‑1 (MQ=255)
                                          aGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATgggg                                                                                                         >  2:1435629/1‑149 (MQ=255)
                                            tACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACt                                                                                                                                   <  2:480903/121‑1 (MQ=255)
                                            tACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACt                                                                                                                                   >  1:480903/1‑121 (MQ=255)
                                            tACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTAc                                                                                                                                 >  3:281566/1‑123 (MQ=255)
                                            tACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTAc                                                                                                                                 <  4:281566/123‑1 (MQ=255)
                                               gATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAAc                                                                                                                                    >  3:30952/1‑117 (MQ=255)
                                               gATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAAc                                                                                                                                    <  4:30952/117‑1 (MQ=255)
                                                 tAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGt                                                                                                  <  2:71171/149‑1 (MQ=255)
                                                   aTCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGttt                                                                                                <  3:413100/149‑1 (MQ=255)
                                                   aTCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGttt                                                                                                <  1:358088/149‑1 (MQ=255)
                                                     cGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTttcc                                                                                              <  4:94581/149‑1 (MQ=255)
                                                       tAAACCGATTTCGAGGCGAACTGCGGATGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTttcctt                                                                                            >  2:792451/1‑149 (MQ=255)
                                                        aaaCCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTttc                                                                                               <  3:48209/145‑1 (MQ=255)
                                                        aaaCCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTttc                                                                                               >  4:48209/1‑145 (MQ=255)
                                                         aaCCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCGGCa                                                                                                                                       >  4:280511/1‑104 (MQ=255)
                                                         aaCCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCGGCa                                                                                                                                       <  3:280511/104‑1 (MQ=255)
                                                         aaCCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTttccttcc                                                                                          <  1:1192924/149‑1 (MQ=255)
                                                          aCCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCa                                                                                         <  1:1435629/149‑1 (MQ=255)
                                                           ccGATTTCGAGGCGAACTGCGGATGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCAc                                                                                        <  1:792451/149‑1 (MQ=255)
                                                           ccGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAAc                                                                                                                    <  4:1348064/121‑1 (MQ=255)
                                                           ccGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAAc                                                                                                                    >  3:1348064/1‑121 (MQ=255)
                                                                        gAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATggg                                                                                                          >  1:573412/1‑118 (MQ=255)
                                                                        gAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATggg                                                                                                          <  2:573412/118‑1 (MQ=255)
                                                                               ggacggaTCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCa                                                                                                                                          <  4:1014427/79‑1 (MQ=255)
                                                                               ggacggaTCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCa                                                                                                                                          >  3:1014427/1‑79 (MQ=255)
                                                                                gacggaTCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGGTCGCCCg                                                                                                                                                <  1:662561/72‑1 (MQ=255)
                                                                                gacggaTCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGGTCGCCCg                                                                                                                                                >  2:662561/1‑72 (MQ=255)
                                                                                   ggaTCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCa                                                                                                                                          >  1:790579/1‑75 (MQ=255)
                                                                                   ggaTCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCAGATTCCa                                                                                                                                          <  2:790579/75‑1 (MQ=255)
                                                                                    gaTCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTcgccgc                                                               >  4:529786/1‑149 (MQ=255)
                                                                                      tCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCa                                                                                                                                          >  1:669336/1‑72 (MQ=255)
                                                                                      tCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCa                                                                                                                                          <  2:669336/72‑1 (MQ=255)
                                                                                        cccGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCa                                                                                                                                          <  1:6708/70‑1 (MQ=255)
                                                                                        cccGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCa                                                                                                                                          >  2:6708/1‑70 (MQ=255)
                                                                                         ccGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGt                                                                             <  2:1031561/130‑1 (MQ=255)
                                                                                         ccGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGt                                                                             >  1:1031561/1‑130 (MQ=255)
                                                                                                tGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTcgccgc                                                               <  3:806586/137‑1 (MQ=255)
                                                                                                tGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTcgccgc                                                               >  4:806586/1‑137 (MQ=255)
                                                                                                   aTCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGttt                                                                                                >  3:667491/1‑101 (MQ=255)
                                                                                                   aTCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGttt                                                                                                <  4:667491/101‑1 (MQ=255)
                                                                                                   aTCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCGCGCCCTTTATCATCCa                                                <  3:529786/149‑1 (MQ=255)
                                                                                                       gCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCGCGCCCTTTATCATCCAGCGt                                            >  3:961331/1‑149 (MQ=255)
                                                                                                         aGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTttccttcc                                                                                          <  2:455461/101‑1 (MQ=255)
                                                                                                         aGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTttccttcc                                                                                          >  1:455461/1‑101 (MQ=255)
                                                                                                              aGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCGCGCCCTTTATCATCCAGCGTCAACTGa                                     >  2:872391/1‑149 (MQ=255)
                                                                                                               gATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTttcct                                                                                             >  3:896500/1‑92 (MQ=255)
                                                                                                               gATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTttcct                                                                                             <  4:896500/92‑1 (MQ=255)
                                                                                                                     cTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAAc                                                                                                                    >  3:864237/1‑63 (MQ=255)
                                                                                                                     cTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAAc                                                                                                                    <  4:864237/63‑1 (MQ=255)
                                                                                                                           aGCGCGGGGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCGCGCCCTTTATCATCCAGCGTCAACTGAACTTCATCCAGCg                        >  2:803259/1‑149 (MQ=255)
                                                                                                                            gcgcggcgcgTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGc                                                                                                    <  1:392448/72‑1 (MQ=255)
                                                                                                                            gcgcggcgcgTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAAGGc                                                                                                    >  2:392448/1‑72 (MQ=255)
                                                                                                                             cgcggcgcgTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGtt                                                                                                 <  3:63690/74‑1 (MQ=255)
                                                                                                                             cgcggcgcgTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGtt                                                                                                 >  4:63690/1‑74 (MQ=255)
                                                                                                                                   gcgTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTTGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCttt                                                                                 <  1:594930/84‑1 (MQ=255)
                                                                                                                                   gcgTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCttt                                                                                 >  2:594930/1‑84 (MQ=255)
                                                                                                                                    cgTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCGCGCCCTTTATCATCCAGCGTCAACTGAACTTCATCCAGCGTGCGGGaaa               <  4:961331/149‑1 (MQ=255)
                                                                                                                                     gTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCGCGCCCTTTATCATCCAGCGTCAACTGAACTTCATCCAGCGTGCGGGAAAg              >  1:542808/1‑149 (MQ=255)
                                                                                                                                      tCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCGCGCCCTTTATCATCCAGCGTCAACTGAACTTCATCCAGCGTGCGGGAAAGa             >  4:485886/1‑149 (MQ=255)
                                                                                                                                       cACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCGCGCCCTTTATCATCCAGCGTCAACTGAACTTCATCCAGCGTGCGGGAAAGAt            >  2:865953/1‑149 (MQ=255)
                                                                                                                                        atcGGGCGGTCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGaaa                                                                        >  3:1421642/3‑88 (MQ=255)
                                                                                                                                        atcGGGCGGTCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGaaa                                                                        <  4:1421642/86‑1 (MQ=255)
                                                                                                                                        aCCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTc                                                                    >  3:1350134/1‑92 (MQ=255)
                                                                                                                                        aCCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTc                                                                    <  4:1350134/92‑1 (MQ=255)
                                                                                                                                        aCCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCGCGCCCTTTATCATCCAGCGTCAACTGAACTTCATCCAGCGTGCGGGAAAGAtt           <  1:872391/149‑1 (MQ=255)
                                                                                                                                               gACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCGCGCCCTTTATCATCCAGCGTCAACTGAACTTCATCCAGCGTGCGGGAAAGAttt          <  2:507961/143‑1 (MQ=255)
                                                                                                                                               gACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCGCGCCCTTTATCATCCAGCGTCAACTGAACTTCATCCAGCGTGCGGGAAAGAttt          >  1:507961/1‑143 (MQ=255)
                                                                                                                                                 ccGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCGCGCCCTTTATCATCCAGCGTCAACTGAACTTCATCCAGCGTGCGGGAAAGATTTTCGGCAGc  <  1:803259/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GTTACTGCCTTCATCAACAATGGGCTGTTTCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCGCGCCCTTTATCATCCAGCGTCAACTGAACTTCATCCAGCGTGCGGGAAAGATTTTCGGCAGC  >  NC_000913/2650353‑2650646

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: