Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,648,458 T→G 15.6% N1277T (AAC→ACC)  yfhM ← bacterial alpha2‑macroglobulin colonization factor ECAM; anti‑host protease defense factor; periplasmic inner membrane‑anchored lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,648,4580TG15.6% 492.1 / 79.6 187N1277T (AAC→ACC) yfhMbacterial alpha2‑macroglobulin colonization factor ECAM; anti‑host protease defense factor; periplasmic inner membrane‑anchored lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (83/74);  new base G (10/19);  total (94/93)
Fisher's exact test for biased strand distribution p-value = 1.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.77e-01

Only 100 of 188 total aligned reads displayed.

CTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTA  >  NC_000913/2648361‑2648556
                                                                                                 |                                                                                                   
cTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTa                                                                                                   >  1:2969866/1‑100 (MQ=255)
     cGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑Gc                                                                                             <  1:145550/100‑1 (MQ=255)
      ggCGAATTTACTGGCTTTGAGATTAGCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCg                                                                                            <  2:448167/100‑1 (MQ=255)
      ggCGAATTTAATGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGC‑CATTAACCGGGGTAATgagcg                                                                                            <  2:2499623/100‑5 (MQ=37)
       gCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGACCTTGTAAATAGCGCAACAGACGCTCATTACCCCGGGTAATG‑GCGt                                                                                           >  1:1270192/1‑100 (MQ=255)
       gCGAATTTACTGGCTTTGAGATTAGCCGCGTACGGGATCGACATCATGTCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATg‑g                                                                                              <  2:654952/97‑1 (MQ=255)
        cGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTc                                                                                          <  2:2149613/100‑1 (MQ=255)
        cGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGC‑CATTAACCGGGGTAATGAGCGTc                                                                                          <  1:1668422/100‑1 (MQ=37)
        cGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGC‑CATTAACCGGGGTAATGAGCGTc                                                                                          <  2:2307202/100‑1 (MQ=37)
         gAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGC‑CATTAACCGGGGTAATGAGCGTCt                                                                                         <  1:2272630/100‑1 (MQ=37)
          aaTTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGC‑CATTAACCGGGGTAATGAGCGTCt                                                                                         <  1:3506287/99‑1 (MQ=37)
           aTTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGC‑CATTAACCGGGGTAATGAGCGTCTgc                                                                                       <  1:1001500/100‑2 (MQ=25)
           aTTTACTGGCTTTGAGATTATCCGCGTACGGGAGCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGC‑CATTAACCGGGGTAATGAGCGTCTgc                                                                                       <  2:3544763/100‑2 (MQ=25)
           aTTTACTGGCATTGAGAATATCCGCGTACGAGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGC‑CATTAACCGGGGTAATGAGCGTCt                                                                                         <  1:563012/98‑1 (MQ=25)
            tttACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAAAAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTc                                                                                      <  2:3257472/100‑1 (MQ=255)
             ttACTGGCTTTGCGATTATCCGCGTACGGGATAGACATCATGCCCGGTTCTTGTAAATAGCGCAGCAGACGC‑CATTAACCGGGGTAATGAGCGTCTgc                                                                                       <  2:1268617/98‑2 (MQ=21)
             ttACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCg                                                                                     >  2:3262996/1‑100 (MQ=255)
             ttACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGC‑CATTAACCGGGGTAATGAGCGTCTgctg                                                                                     <  1:3290210/100‑4 (MQ=25)
              tACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCgg                                                                                    >  2:156482/1‑100 (MQ=255)
              tACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCgg                                                                                    <  2:2026758/100‑1 (MQ=255)
              tACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCgg                                                                                    <  2:1912296/100‑1 (MQ=255)
              tACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGC‑CATTAACCGGGGTAATGAGCGTCTgctgc                                                                                    <  1:390810/100‑5 (MQ=25)
               tccgatcttGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGc                                                                                   <  2:3144901/93‑1 (MQ=255)
               aCTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGc                                                                                   <  1:2898908/100‑1 (MQ=255)
               aCTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGC‑CATTAACCGGGGTAATGAGCGTCTgctgcg                                                                                   <  1:2642431/100‑6 (MQ=21)
                cTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCa                                                                                  <  2:2719997/100‑1 (MQ=255)
                 tGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGC‑CATTAACCGGGGTAATGAGCGTCTgctgcgct                                                                                 <  1:1089653/100‑8 (MQ=21)
                  ggCTTTGAGATTATCCGCGTACGGGATCGACACCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACg                                                                                >  1:3618305/1‑100 (MQ=255)
                    cTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCt                                                                              >  2:124777/1‑100 (MQ=255)
                    cTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCt                                                                              >  2:820016/1‑100 (MQ=255)
                     tttGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTg                                                                             <  2:3072585/100‑1 (MQ=255)
                      ttGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCt                                                                              <  1:759308/98‑1 (MQ=255)
                      ttGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTg                                                                             <  2:2620293/99‑1 (MQ=255)
                       tGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCagatcgga                                                                           >  1:3144925/1‑93 (MQ=255)
                       tGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTa                                                                           <  2:1397743/100‑1 (MQ=255)
                       tGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTa                                                                           >  1:2212192/1‑100 (MQ=255)
                           aTTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAAccc                                                                       <  1:3178198/100‑1 (MQ=255)
                           aTTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAAccc                                                                       <  2:1630675/100‑1 (MQ=255)
                           aTTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAAccc                                                                       <  2:2844611/100‑1 (MQ=255)
                           aTTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAAccc                                                                       <  1:1405126/100‑1 (MQ=255)
                            ttATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCt                                                                      >  2:162379/1‑100 (MQ=255)
                             tATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTg                                                                     <  1:3074922/100‑1 (MQ=255)
                             tATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTg                                                                     <  1:1692409/100‑1 (MQ=255)
                              aTCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGt                                                                    >  1:45659/1‑100 (MQ=255)
                                ccGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTc                                                                  >  1:3128907/1‑100 (MQ=255)
                                ccGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTc                                                                  >  1:3005289/1‑100 (MQ=255)
                                   cgTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGcc                                                               >  1:2074808/1‑100 (MQ=255)
                                   cgTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGcc                                                               <  1:2925764/100‑1 (MQ=255)
                                    ggaCGGGATCGACGTCATGCCCGGATCTTGTAAATAGCGCAGAAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCt                                                              >  1:2587539/3‑100 (MQ=255)
                                    gTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCt                                                              >  2:1044856/1‑100 (MQ=255)
                                     tACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTg                                                             <  2:1495233/100‑1 (MQ=255)
                                       cGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTAATTACACCGGTTACTG‑GCGTCTGTCGGCACGATGTAACCCTGTTCGCATgc                                                            >  1:3187770/1‑99 (MQ=255)
                                        gggATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTgcgc                                                          >  1:3445568/1‑100 (MQ=255)
                                         ggATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGcgcg                                                         >  2:543169/1‑100 (MQ=255)
                                          gATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTATAACCCTGTTCGCCTGCGCgg                                                        >  1:3402354/1‑100 (MQ=255)
                                            tCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGAc                                                      <  2:3362145/100‑1 (MQ=255)
                                              gACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCa                                                    >  2:1191418/1‑100 (MQ=255)
                                               aCATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCCCGCTGTAACCCTGTTCGCCTGCGCGGACCAg                                                   >  2:239833/1‑100 (MQ=255)
                                                catcatGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAgg                                                  >  2:1426619/1‑100 (MQ=255)
                                                catcatGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAgg                                                  >  2:1823145/1‑100 (MQ=255)
                                                 atcatGCCTGGATCTTGTAAATAGCGCAGCAGACGCTCCTTACCCCTGTTAATG‑GCCTCTGTCGGCACGCTGTAGCCCTGTTCGCCGGcgcg                                                         <  2:1101521/92‑1 (MQ=255)
                                                  tcatGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGaa                                                >  1:3434281/1‑100 (MQ=255)
                                                  tcatGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGaa                                                >  1:2741196/1‑100 (MQ=255)
                                                   catGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGaaa                                               <  1:3631635/100‑1 (MQ=255)
                                                     tGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAAt                                              >  2:3005918/1‑99 (MQ=255)
                                                      gCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATc                                             >  2:370905/1‑99 (MQ=255)
                                                       cccGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCa                                           >  2:512428/1‑100 (MQ=255)
                                                        ccGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCAt                                          <  1:78138/100‑1 (MQ=255)
                                                        ccGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCAt                                          <  1:478044/100‑1 (MQ=255)
                                                        ccGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCAt                                          <  1:3658521/100‑1 (MQ=255)
                                                          ggATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCa                                        <  1:1015118/100‑1 (MQ=255)
                                                              cTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCTGGAAATCCATCACGTa                                    >  1:3043722/1‑100 (MQ=255)
                                                                tGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAg                                  <  2:474581/100‑1 (MQ=255)
                                                                tGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAg                                  >  1:177124/1‑100 (MQ=255)
                                                                  tAAATAGCGCAGCAGACGCTCATTACTCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGc                                 <  1:717727/99‑1 (MQ=255)
                                                                  tAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGcc                                >  1:1346881/1‑100 (MQ=255)
                                                                   aaaTAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCg                               <  1:946091/100‑1 (MQ=255)
                                                                    aaTAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGt                              >  2:3171219/1‑100 (MQ=255)
                                                                     aTAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTc                             <  2:881525/100‑1 (MQ=255)
                                                                        gcgcAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGc                          >  1:3124298/1‑100 (MQ=255)
                                                                         cgcAGCAGACGCTCATTACCCCGGTTCATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGcc                         <  2:3428559/100‑1 (MQ=255)
                                                                         cgcAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGcc                         <  1:2712131/100‑1 (MQ=255)
                                                                          ccgacagACGC‑CATTAACCGGGGTAATGAGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCa                        >  2:2421201/5‑100 (MQ=25)
                                                                           cagcagACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAg                       <  1:2064213/100‑1 (MQ=255)
                                                                            gacagACGC‑CATTAACCGGGGTAATGAGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGt                      >  1:2439886/3‑100 (MQ=25)
                                                                            gacagACGC‑CATTAACCGGGGTAATGAGCGTCTGTCGGCACGCTGTAACCCTGTTAGCCTGCGCGGACAAGGAAATCCATCACGTAAGCCGTCAGCCCGt                      >  1:3646376/3‑100 (MQ=21)
                                                                            gacagACGC‑CATTAACCGGGGCAATGAGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCAGCGCGGCCCAGGAAATCCATCACGTAGGCCGTCAGCCAGt                      >  2:2272308/3‑100 (MQ=21)
                                                                             gcagACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTa                     >  1:2305664/1‑100 (MQ=255)
                                                                             gcagACGCTCATTACCCCGGATAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTGCGCCACCGCGGACCAGAGCATCCATCACGTAAGCCGTCCGcc                         >  2:1340832/1‑96 (MQ=255)
                                                                              cagACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTAc                    <  1:2019335/100‑1 (MQ=255)
                                                                               agACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTAct                   >  1:626640/1‑100 (MQ=255)
                                                                                gACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTActc                  >  2:369008/1‑100 (MQ=255)
                                                                                  cGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCtt                >  1:3773859/1‑100 (MQ=255)
                                                                                     tCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGt             <  2:2370960/100‑1 (MQ=255)
                                                                                      cATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTc            <  1:1896892/100‑1 (MQ=255)
                                                                                       aTTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCa           >  2:1233176/1‑100 (MQ=255)
                                                                                       aTTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCa           <  1:1132231/100‑1 (MQ=255)
                                                                                          aCCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCg        >  1:862779/1‑100 (MQ=255)
                                                                                            cccGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGtt      >  1:1814577/1‑100 (MQ=255)
                                                                                                gTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTa  >  2:3699786/1‑100 (MQ=255)
                                                                                                 |                                                                                                   
CTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATG‑GCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTA  >  NC_000913/2648361‑2648556

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑