Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,548,617 A→C L111L (CTA→CTC)  fimG → minor component of type 1 fimbriae

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,548,6170AC99.3% 516.3 / NA 148L111L (CTA→CTC) fimGminor component of type 1 fimbriae
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (68/79);  minor base T (0/1);  total (68/80)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 148 total aligned reads displayed.

GGGAACGTCGAGGGTCACTGCCAGCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTACAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAGGTCAGA  >  NC_000913/4548521‑4548713
                                                                                                |                                                                                                
 ggAACGTCGAGGGTCACTGCCAGCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAgg                                                                                              >  1:1985342/1‑100 (MQ=255)
  gAACGTCGAGGGTCACTGCCAGCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGa                                                                                             <  2:667798/100‑1 (MQ=255)
   aaCGTCGAGGGTCACTGCCAGCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGAt                                                                                            <  2:1216130/100‑1 (MQ=255)
      gTCGAGGGTCACTGCCAGCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGAc                                                                                         <  2:3759772/100‑1 (MQ=255)
       tCGAGGGTCACTGCCAGCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACa                                                                                        <  1:1521692/100‑1 (MQ=255)
          aGGGTCACTGCCAGCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTg                                                                                     <  2:2876883/100‑1 (MQ=255)
           gggTCACTGCCAGCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTgg                                                                                    <  1:1796412/100‑1 (MQ=255)
            ggTCACTGCCAGCTTCAGCGGGGCAGCCGACAGTACCGGATGTTATAAAGACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACTGTGGc                                                                                   >  1:3868780/1‑100 (MQ=255)
             gTCACTGCCAGCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCa                                                                                  <  2:3481470/100‑1 (MQ=255)
             gTCACTGCCAGCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCa                                                                                  <  1:2324620/100‑1 (MQ=255)
              tCACTGCCAGCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCaa                                                                                 <  2:1536778/100‑1 (MQ=255)
              tCACTGCCAGCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCaa                                                                                 >  2:2651446/1‑100 (MQ=255)
               cACTGCCAGCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAGCCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAac                                                                                <  2:1310781/100‑1 (MQ=255)
               cACTGCCAGCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAac                                                                                <  1:178401/100‑1 (MQ=255)
                aCTGCCAGCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAaca                                                                               >  2:2723696/1‑100 (MQ=255)
                aCTGCCAGCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAaca                                                                               >  2:3016809/1‑100 (MQ=255)
                 cTGCCAGCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAacac                                                                              >  2:3124766/1‑100 (MQ=255)
                  tGCCAGCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAAcaca                                                                             <  2:1813901/100‑1 (MQ=255)
                   gCCAGCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACAt                                                                            <  2:2487474/100‑1 (MQ=255)
                   gCCAGCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACAt                                                                            >  1:632243/1‑100 (MQ=255)
                   gCCAGCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCGCCAGGATGACAGTGGCAACACAt                                                                            <  1:1153918/100‑1 (MQ=255)
                     cAGCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTg                                                                          <  1:2784998/100‑1 (MQ=255)
                     cAGCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTg                                                                          <  1:759838/100‑1 (MQ=255)
                     cAGCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTg                                                                          <  1:3410636/100‑1 (MQ=255)
                      aGCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGa                                                                         <  1:1411809/100‑1 (MQ=255)
                      aGCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGa                                                                         <  1:1411840/100‑1 (MQ=255)
                       gCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGaa                                                                        <  2:226681/100‑1 (MQ=255)
                       gCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGaa                                                                        <  2:719127/100‑1 (MQ=255)
                        cTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGCGGCAACACATGGAAt                                                                       >  2:501817/1‑100 (MQ=255)
                        cTTCAACGGGGCAGCCGACAGTACCGGATATTATAAAGACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAAt                                                                       <  2:2626486/100‑1 (MQ=255)
                          tCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATAc                                                                     <  2:2018939/100‑1 (MQ=255)
                            aGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGTCCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTg                                                                   >  2:527570/1‑100 (MQ=255)
                            aGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTg                                                                   >  1:2848700/1‑100 (MQ=255)
                            aGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTg                                                                   >  2:1937279/1‑100 (MQ=255)
                             gCGGGGCAGCCGCCAGGACCGGATATTATAAAAACCAGGGGACCGCGAAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTgg                                                                  <  2:1991727/100‑1 (MQ=255)
                             gCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAATACATTGAATACTgg                                                                  <  2:2308012/100‑1 (MQ=255)
                              cGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGgc                                                                 >  2:3713132/1‑100 (MQ=255)
                              cGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGgc                                                                 >  1:1221112/1‑100 (MQ=255)
                              cGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGgc                                                                 <  1:2110604/100‑1 (MQ=255)
                               ggggCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGgcg                                                                >  2:2418101/1‑100 (MQ=255)
                               ggggCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGgcg                                                                <  2:900415/100‑1 (MQ=255)
                                gggCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGgcgc                                                               >  1:2558694/1‑100 (MQ=255)
                                gggCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCACCAGGATGACAGTGGCAACACATTGAATACTGgcgc                                                               >  1:3840936/1‑100 (MQ=255)
                                     gCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCa                                                          >  2:3147331/1‑100 (MQ=255)
                                        gACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCaaaa                                                       <  1:1875488/100‑1 (MQ=255)
                                        gACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCaaaa                                                       >  2:3537093/1‑100 (MQ=255)
                                            gTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGt                                                   >  1:2859885/1‑100 (MQ=255)
                                              aCCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTc                                                 <  1:3180037/100‑1 (MQ=255)
                                               ccGGATAGTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGTTGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCa                                                >  2:2500423/1‑100 (MQ=255)
                                                  gATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGt                                             <  1:3087848/100‑1 (MQ=255)
                                                   atatTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTg                                            <  2:1368335/100‑1 (MQ=255)
                                                    tattatAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTgg                                           <  1:952458/100‑1 (MQ=255)
                                                      ttatAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTggct                                         >  2:450567/1‑98 (MQ=255)
                                                      ttatAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGgat                                         <  1:590734/100‑1 (MQ=255)
                                                         taAAAACCAGGGGACCGCGCAAAACATCCAGTGAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGgatgat                                      <  2:2884540/100‑1 (MQ=255)
                                                             aaCCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCt                                  >  1:2106293/1‑100 (MQ=255)
                                                             aaCCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCt                                  <  2:2427214/100‑1 (MQ=255)
                                                             aaCCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCt                                  >  1:661702/1‑100 (MQ=255)
                                                              aataGGGGACCGCGCAAAACATCCGGTTAGAGCTCCACGCTGACAGTGGAAAAACATAGAATACTGGCGCAAACAAAACAGATCAGGTGGACAATTccgc                                 >  2:3185590/4‑98 (MQ=255)
                                                              aCCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTc                                 <  2:823501/100‑1 (MQ=255)
                                                               ccAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTca                                <  1:3692091/100‑1 (MQ=255)
                                                                cAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTcac                               <  1:3953871/100‑1 (MQ=255)
                                                                cAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTcac                               >  2:3100624/1‑100 (MQ=255)
                                                                 aGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTcaca                              >  1:2629321/1‑100 (MQ=255)
                                                                  ggggACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACaa                             >  2:2668615/1‑100 (MQ=255)
                                                                   gggACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAAt                            <  2:3601114/100‑1 (MQ=255)
                                                                   gggACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAAt                            <  1:1866481/100‑1 (MQ=255)
                                                                     gACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCa                          >  1:2550029/1‑100 (MQ=255)
                                                                      aCCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAg                         <  2:2080912/100‑1 (MQ=255)
                                                                      aCCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAg                         >  1:1189538/1‑100 (MQ=255)
                                                                         gcgcAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCCCAATCAgcgc                      >  2:2044527/1‑100 (MQ=255)
                                                                         gcgcAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAgcgc                      <  1:273819/100‑1 (MQ=255)
                                                                         gcgcAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAgcgc                      >  1:3592550/1‑100 (MQ=255)
                                                                          cgcAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCa                     <  2:292290/100‑1 (MQ=255)
                                                                           gcAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCAc                    <  1:2692436/100‑1 (MQ=255)
                                                                           gcAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCAc                    >  1:953088/1‑100 (MQ=255)
                                                                           gcAAAACATCCAGTTAGAGCTCCAGGAGGACAGTGGCAACACATTGAATACTGCCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGAAc                    <  2:875484/100‑1 (MQ=255)
                                                                            cAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACt                   >  2:1449160/1‑100 (MQ=255)
                                                                            cAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACt                   >  2:2995699/1‑100 (MQ=255)
                                                                             aaaaCATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACtt                  >  2:1718259/1‑100 (MQ=255)
                                                                              aaaCATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTc                 >  1:1929950/1‑100 (MQ=255)
                                                                              aaaCATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTc                 <  2:407109/100‑1 (MQ=255)
                                                                               aaCATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTcc                <  1:697832/100‑1 (MQ=255)
                                                                                aCATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTccc               <  1:2959210/100‑1 (MQ=255)
                                                                                 cATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCg              <  2:3571506/100‑1 (MQ=255)
                                                                                  aTCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGt             <  2:1295012/100‑1 (MQ=255)
                                                                                   tCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGtt            >  2:3607433/1‑100 (MQ=255)
                                                                                     cAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTAc          <  2:2062761/100‑1 (MQ=255)
                                                                                     cAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTAc          <  1:2329787/100‑1 (MQ=255)
                                                                                      aGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACa         >  1:3997570/1‑100 (MQ=255)
                                                                                       gTTCGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAg        >  1:3142391/1‑100 (MQ=255)
                                                                                       gTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAg        <  2:1890382/100‑1 (MQ=255)
                                                                                        ttAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAATAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAgg       <  2:1188085/100‑1 (MQ=255)
                                                                                         tAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAGGt      <  1:51502/100‑1 (MQ=255)
                                                                                         tAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAGGt      >  1:3373024/1‑100 (MQ=255)
                                                                                          agagCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAGGtn     <  1:3239078/100‑2 (MQ=255)
                                                                                            agCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAGGTCag   <  2:477584/100‑1 (MQ=255)
                                                                                             gCTCCAGGATGATAGTGGCAACACATTGAATCCTGGCGCAACCAAAAGAGTACAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAGGTCaga  >  2:3563099/1‑100 (MQ=255)
                                                                                             gCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAGGTCaga  >  2:1874037/1‑100 (MQ=255)
                                                                                             gCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAGGTCaga  >  1:1543040/1‑100 (MQ=255)
                                                                                                |                                                                                                
GGGAACGTCGAGGGTCACTGCCAGCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTACAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAGGTCAGA  >  NC_000913/4548521‑4548713

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: ‑