Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,650,497 A→G 7.1% G597G (GGT→GGC) ‡ 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,4970AG 7.1% 324.0 / 10.9 98G597G (GGT→GGC) ‡yfhMbacterial alpha2‑macroglobulin colonization factor ECAM; anti‑host protease defense factor; periplasmic inner membrane‑anchored lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (47/44);  new base G (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 = 7.98e-01

GTTACTGCCTTCATCAACAATGGGCTGTTTCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCGCGCCCTTTATCATCCAGCGTCAACTGAACTTCATCCAGCGTGCGGGAAAGATTT  >  NC_000913/2650353‑2650638
                                                                                                                                                |                                                                                                                                             
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                                                                                                                                          <  3:1037662/140‑1 (MQ=255)
          tcatcaACAATGGGCTGTTTCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGcc                                                                                                                                          >  4:1037662/1‑140 (MQ=255)
          tcatcaACAATGGGCTGTTTCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGAtt                                                                                                                                     >  2:852712/1‑145 (MQ=255)
          tcatcaACAATGGGCTGTTTCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGAtt                                                                                                                                     <  1:852712/145‑1 (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                                                                                                                         <  4:555306/130‑1 (MQ=255)
                                     cTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTAc                                                                                                                         >  3:555306/1‑130 (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                                                                                                                           >  1:480903/1‑121 (MQ=255)
                                            tACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACt                                                                                                                           <  2:480903/121‑1 (MQ=255)
                                            tACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTAc                                                                                                                         >  3:281566/1‑123 (MQ=255)
                                            tACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTAc                                                                                                                         <  4:281566/123‑1 (MQ=255)
                                               gATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAAc                                                                                                                            <  4:30952/117‑1 (MQ=255)
                                               gATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAAc                                                                                                                            >  3:30952/1‑117 (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                                                                                                  <  2:573412/118‑1 (MQ=255)
                                                                        gAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATggg                                                                                                  >  1:573412/1‑118 (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                                                                                                                                  <  2:669336/72‑1 (MQ=255)
                                                                                      tCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCa                                                                                                                                  >  1:669336/1‑72 (MQ=255)
                                                                                        cccGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCa                                                                                                                                  >  2:6708/1‑70 (MQ=255)
                                                                                        cccGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCa                                                                                                                                  <  1:6708/70‑1 (MQ=255)
                                                                                         ccGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGt                                                                     >  1:1031561/1‑130 (MQ=255)
                                                                                         ccGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGt                                                                     <  2:1031561/130‑1 (MQ=255)
                                                                                                tGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTcgccgc                                                       >  4:806586/1‑137 (MQ=255)
                                                                                                tGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTcgccgc                                                       <  3:806586/137‑1 (MQ=255)
                                                                                                   aTCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGttt                                                                                        <  4:667491/101‑1 (MQ=255)
                                                                                                   aTCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGttt                                                                                        >  3:667491/1‑101 (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                                                                                         >  4:63690/1‑74 (MQ=255)
                                                                                                                             cgcggcgcgTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGtt                                                                                         <  3:63690/74‑1 (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                                                                <  4:1421642/86‑1 (MQ=255)
                                                                                                                                        atcGGGCGGTCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGaaa                                                                >  3:1421642/3‑88 (MQ=255)
                                                                                                                                        aCCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTc                                                            <  4:1350134/92‑1 (MQ=255)
                                                                                                                                        aCCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTc                                                            >  3:1350134/1‑92 (MQ=255)
                                                                                                                                        aCCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCGCGCCCTTTATCATCCAGCGTCAACTGAACTTCATCCAGCGTGCGGGAAAGAtt   <  1:872391/149‑1 (MQ=255)
                                                                                                                                               gACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCGCGCCCTTTATCATCCAGCGTCAACTGAACTTCATCCAGCGTGCGGGAAAGAttt  <  2:507961/143‑1 (MQ=255)
                                                                                                                                               gACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCGCGCCCTTTATCATCCAGCGTCAACTGAACTTCATCCAGCGTGCGGGAAAGAttt  >  1:507961/1‑143 (MQ=255)
                                                                                                                                                |                                                                                                                                             
GTTACTGCCTTCATCAACAATGGGCTGTTTCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCGCGCCCTTTATCATCCAGCGTCAACTGAACTTCATCCAGCGTGCGGGAAAGATTT  >  NC_000913/2650353‑2650638

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: