Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,650,498 C→T 6.7% 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 6.7% 333.4 / 10.1 105G597D (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 (44/54);  new base T (2/5);  total (46/59)
Fisher's exact test for biased strand distribution p-value = 4.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.48e-01

Only 100 of 106 total aligned reads displayed.

GCGTTACTGCCTTCATCAACAATGGGCTGTTTCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCGCGCCCTTTATCATCCAGCGTCAACTGAACTTCATCCAGCGTGCGGGAAAGATTTTCGGCA  >  NC_000913/2650351‑2650644
                                                                                                                                                   |                                                                                                                                                  
gcgTTACTGCCTTCATCAACAATGGGCTGTTTCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGAc                                                                                                                                                    >  2:599875/1‑148 (MQ=255)
    tACTGCCTTCATCAACAATGGGCTGTTTCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGccc                                                                                                                                               >  4:613893/1‑149 (MQ=255)
       tGCCTTCATCAACAATGGGCTGTTTCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGGTCGCCCgg                                                                                                                                             <  2:46300/148‑2 (MQ=255)
       tGCCTTCATCAACAATGGGCTGTTTCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGa                                                                                                                                             >  2:890252/1‑148 (MQ=255)
         ccTTCATCAACAATGGGCTGTTTCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGGTCGCCCggtga                                                                                                                                          <  4:124867/149‑5 (MQ=255)
         ccTTCATCAACAATGGGCTGTTTCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGGTCGCCCggtg                                                                                                                                           <  1:960364/148‑4 (MQ=255)
          cTTCATCAACAATGGGCTGTTTCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTcc                                                                                                                                         <  4:832528/149‑1 (MQ=255)
           ttCATCAACAATGGGCTGTTTCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTcc                                                                                                                                         >  1:254043/1‑148 (MQ=255)
           ttCATCAACAATGGGCTGTTTCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTcc                                                                                                                                         <  1:314355/148‑1 (MQ=255)
           ttCATCAACAATGGGCTGTTTCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTcc                                                                                                                                         >  1:923117/1‑148 (MQ=255)
           ttCATCAACAATGGGCTGTTTCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTcc                                                                                                                                         <  1:920065/148‑1 (MQ=255)
            tcatcaACAATGGGCTGTTTCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCa                                                                                                                                        <  1:597679/148‑1 (MQ=255)
              atcaACAATGGGCTGTTTCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCa                                                                                                                                     >  3:26050/1‑149 (MQ=255)
               tcaACAATGGGCTGTTTCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGGTCGCCCggt                                                                                                                                            <  3:345078/141‑3 (MQ=255)
               tcaACAATGGGCTGTTTCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCa                                                                                                                                     >  2:202/1‑148 (MQ=255)
                caacaaTGGGCTGTTTCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCATCGGGCGACCGCCCGATTCCAGCaa                                                                                                                                    >  2:632758/1‑148 (MQ=255)
                     aTGGGCTGTTTCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTAc                                                                                                                               >  1:811937/1‑148 (MQ=255)
                     aTGGGCTGTTTCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTAc                                                                                                                               <  2:105624/148‑1 (MQ=255)
                      tGGGCTGTTTCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTAccc                                                                                                                             <  4:623769/149‑1 (MQ=255)
                       gggCTGTTTCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTAccc                                                                                                                             <  2:510844/148‑1 (MQ=255)
                       gggCTGTTTCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTAccc                                                                                                                             <  2:963532/148‑1 (MQ=255)
                         gCTGTTTCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTgg                                                                                                                          >  4:35354/1‑149 (MQ=255)
                          cTGTTTCACCGTGCTGTCAGTAAGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTgg                                                                                                                          >  2:226919/1‑148 (MQ=255)
                            gTTTCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGaaa                                                                                                                       <  4:223414/149‑1 (MQ=255)
                              ttCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAAt                                                                                                                     <  4:889238/149‑1 (MQ=255)
                                cACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATaa                                                                                                                   <  3:613862/149‑1 (MQ=255)
                                 aCCGTGCAGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGCGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATaa                                                                                                                   >  2:340595/1‑148 (MQ=255)
                                       cTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTaa                                                                                                            <  3:613893/149‑1 (MQ=255)
                                         gTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAAt                                                                                                           <  1:733500/148‑1 (MQ=255)
                                         gTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATg                                                                                                          >  4:132624/1‑149 (MQ=255)
                                         gTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATg                                                                                                          <  3:324349/149‑1 (MQ=255)
                                           cAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATggg                                                                                                        <  4:306033/149‑1 (MQ=255)
                                           cAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATggg                                                                                                        <  3:854230/149‑1 (MQ=255)
                                           cAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATgg                                                                                                         <  2:617630/148‑1 (MQ=255)
                                             gTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGa                                                                                                      >  4:373224/1‑149 (MQ=255)
                                              tACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGa                                                                                                      <  2:819245/148‑1 (MQ=255)
                                              tACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGa                                                                                                      <  1:740539/148‑1 (MQ=255)
                                                 gATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGc                                                                                                  <  3:931879/149‑1 (MQ=255)
                                                  aTAATCGTAAACCGATTTCGAGGAGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCg                                                                                                 >  3:929371/1‑149 (MQ=255)
                                                   tAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTGCCCTGGAAAATAACCTGTAATGGGGAATGCg                                                                                                 <  1:599875/148‑1 (MQ=255)
                                                    aaTCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGtt                                                                                               <  4:102386/149‑1 (MQ=255)
                                                     aTCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGtt                                                                                               >  1:715912/1‑148 (MQ=255)
                                                      tCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGttt                                                                                              <  2:697667/148‑1 (MQ=255)
                                                          aaaCCGATTTCGAGGCGAACTGCTGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTttccttc                                                                                         <  4:760521/149‑1 (MQ=255)
                                                          aaaCCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTttcctt                                                                                          >  1:806460/1‑148 (MQ=255)
                                                          aaaCCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTttcctt                                                                                          >  1:776596/1‑148 (MQ=255)
                                                           aaCCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTttccttcc                                                                                        <  3:35354/149‑1 (MQ=255)
                                                                aTTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTgg                                                                                   <  3:438511/149‑1 (MQ=255)
                                                                aTTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTg                                                                                    <  2:781551/148‑1 (MQ=255)
                                                                     gAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTc                                                                              <  3:564297/149‑1 (MQ=255)
                                                                      aGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTc                                                                              <  1:420325/148‑1 (MQ=255)
                                                                      aGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTc                                                                              <  1:649536/148‑1 (MQ=255)
                                                                      aGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTc                                                                              <  1:933431/148‑1 (MQ=255)
                                                                       ggCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTg                                                                            >  3:689981/1‑149 (MQ=255)
                                                                       ggCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTg                                                                            >  3:540922/1‑149 (MQ=255)
                                                                        gCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGt                                                                           >  4:523533/1‑149 (MQ=255)
                                                                         cGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGt                                                                           <  2:776596/148‑1 (MQ=255)
                                                                         cGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTa                                                                          <  4:289779/149‑1 (MQ=255)
                                                                            aCTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGa                                                                        <  1:107429/148‑1 (MQ=255)
                                                                            aCTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGa                                                                        >  1:441133/1‑148 (MQ=255)
                                                                             cTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGTAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGaa                                                                       <  2:828857/148‑1 (MQ=255)
                                                                             cTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGaa                                                                       <  1:348231/148‑1 (MQ=255)
                                                                             cTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGaa                                                                       <  2:399032/148‑1 (MQ=255)
                                                                                 ggacggaTCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTc                                                                  >  3:798001/1‑149 (MQ=255)
                                                                                  gacggaTCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAAc                                                                     >  4:642606/1‑145 (MQ=255)
                                                                                  gacggaTCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAAc                                                                     <  3:642606/145‑1 (MQ=255)
                                                                                    cggaTCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTcgcc                                                               <  3:915582/149‑1 (MQ=255)
                                                                                     ggaTCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTcgccg                                                              >  4:155270/1‑149 (MQ=255)
                                                                                        tCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGGTCGCCCg                                                                                                                                              >  1:115947/1‑66 (MQ=255)
                                                                                        tCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGGTCGCCCg                                                                                                                                              <  2:115947/66‑1 (MQ=255)
                                                                                        tCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCgcgc                                                           >  4:409123/1‑149 (MQ=255)
                                                                                          cccGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCGCGcc                                                          >  1:786726/1‑148 (MQ=255)
                                                                                           ccGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCGCGccc                                                         >  1:757179/1‑148 (MQ=255)
                                                                                            cGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCGCGCCCt                                                        >  1:352110/1‑148 (MQ=255)
                                                                                             ggCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCGCGCCCtt                                                       <  2:760177/148‑1 (MQ=255)
                                                                                              gCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCGCGCCCttt                                                      <  2:186426/148‑1 (MQ=255)
                                                                                                aaTGCATCGGCGGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTCGAAACTTCGCCGCGCCCTTTatc                                                   <  4:62656/149‑1 (MQ=255)
                                                                                                 aTGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCGCGCCCTTTatca                                                  <  4:114435/149‑1 (MQ=255)
                                                                                                    cATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCGCGCCCTTTatcatc                                                >  2:174660/1‑148 (MQ=255)
                                                                                                     aTCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCGCGCCCTTTATCATcc                                               <  2:352110/148‑1 (MQ=255)
                                                                                                     aTCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCGCGCCCTTTATCATcc                                               >  2:219207/1‑148 (MQ=255)
                                                                                                     aTCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCGCGCCCTTTATCATCCa                                              <  4:674186/149‑1 (MQ=255)
                                                                                                     aTCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCGCGCCCTTTATCATCCa                                              >  3:674186/1‑149 (MQ=255)
                                                                                                         gCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCATCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCGCGCCCTTGATCATCCGGCGt                                          <  4:829032/149‑1 (MQ=255)
                                                                                                             gCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCGCGCCCTTTATCATCCAGCGTCaa                                       >  2:698961/1‑148 (MQ=255)
                                                                                                                 gATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCGCGCCCTTTATCATCCAGCGTCAACTGaa                                  <  3:217608/149‑1 (MQ=255)
                                                                                                                 gATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCGCGCCCTTTATCATCCAGCGTCAACTGa                                   <  2:31658/148‑1 (MQ=255)
                                                                                                                        tGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCGCGCCCTTTATCATCCAGCGTCAACTGAACTTCAt                            >  2:399964/1‑148 (MQ=255)
                                                                                                                        tGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCGCGCCCTTTATCATCCAGCGTCAACTGAACTTCAt                            >  1:542361/1‑148 (MQ=255)
                                                                                                                           tCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCGCGCCCTTTATCATCCAGCGTCAACTGAACTTCATCCAg                        >  3:750194/1‑149 (MQ=255)
                                                                                                                              gcgcggcgcgTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCGCGCCCTTTATCATCCAGCGTCAACTGAACTTCATCCAGCGt                     >  4:602079/1‑149 (MQ=255)
                                                                                                                                gcggcgcgTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCGCGCCCTTTATCATCCAGCGTCAACTGAACTTCATCCAGCGTGc                   <  3:396175/149‑1 (MQ=255)
                                                                                                                                  ggcgcgTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCGCGCCCTTTATCATCCAGCGTCAACTGAACTTCATCCAGCGTGCgg                 <  4:689981/149‑1 (MQ=255)
                                                                                                                                      cgTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCGCGCCCTTTATCATCCAGCGTCAACTGAACTTCATCCAGCGTGCGGGaa              <  2:342575/148‑1 (MQ=255)
                                                                                                                                          aCCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCGCGCCCTTTATCATCCAGCGTCAACTGAACTTCATCCAGCGTGCGGGAAAGAtt         >  4:529674/1‑149 (MQ=255)
                                                                                                                                            cGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCGCGCCCTTTATCATCCAGCGTCAACTGAACTTCATCCAGCGTGCGGGAAAGAtttt       <  3:474441/149‑1 (MQ=255)
                                                                                                                                             gggCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCGCGCCCTTTATCATCCAGCGTCAACTGAACTTCATCCAGCGTGCGGGAAAGAtttt       <  2:441133/148‑1 (MQ=255)
                                                                                                                                             gggCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCGCGCCCTTTATCATCCAGCGTCAACTGAACTTCATCCAGCGTGCGGGAAAGATTTTc      >  3:138141/1‑149 (MQ=255)
                                                                                                                                                cGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCGCGCCCTTTATCATCCAGCGTCAACTGAACTTCATCCAGCGTGCGGGAAAGATTTTCgg    >  1:575490/1‑148 (MQ=255)
                                                                                                                                                  aCCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACATCGCCGCGCCCTTTATCATCCAGCGTCAACTGAACTTCATCCAGCGTGCGGGAAAGATTTTCGGCa  >  1:819505/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCGTTACTGCCTTCATCAACAATGGGCTGTTTCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCGCGCCCTTTATCATCCAGCGTCAACTGAACTTCATCCAGCGTGCGGGAAAGATTTTCGGCA  >  NC_000913/2650351‑2650644

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