Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,547,319 C→A R742S (CGT→AGT)  fimD → fimbrial usher outer membrane porin protein; FimCD chaperone‑usher

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,547,3190CA99.1% 377.8 / ‑4.0 111R742S (CGT→AGT) fimDfimbrial usher outer membrane porin protein; FimCD chaperone‑usher
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (47/63);  total (48/63)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

Only 100 of 111 total aligned reads displayed.

CTGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAAC  >  NC_000913/4547223‑4547414
                                                                                                |                                                                                               
cTGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTg                                                                                              <  1:1405884/100‑1 (MQ=255)
  ggggCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGt                                                                                            >  2:3385379/1‑100 (MQ=255)
    ggCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTa                                                                                          >  1:1657779/1‑100 (MQ=255)
     gCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTAt                                                                                         >  2:1993689/1‑100 (MQ=255)
        gCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGcc                                                                                      <  2:2661719/100‑1 (MQ=255)
         ccGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCg                                                                                     <  2:935453/100‑1 (MQ=255)
          cGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGt                                                                                    >  2:38943/1‑100 (MQ=255)
          cGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGt                                                                                    <  1:683270/100‑1 (MQ=255)
           gTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGtg                                                                                   <  1:187782/100‑1 (MQ=255)
           gTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGtg                                                                                   >  1:361069/1‑100 (MQ=255)
            ttAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGtgc                                                                                  <  2:323052/100‑1 (MQ=255)
            ttAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGtgc                                                                                  <  2:1271670/100‑1 (MQ=255)
             tAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGtgct                                                                                 <  1:2518912/100‑1 (MQ=255)
             tAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGtgc                                                                                  <  1:2697866/99‑1 (MQ=255)
             tAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGtgc                                                                                  <  1:3170584/99‑1 (MQ=255)
                aCGATACGGTGGTGCATGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGAGTACAGACTGGAGTGGTTATGCCGTGCTGcc                                                                              >  2:971412/1‑100 (MQ=255)
                 cgatctggtggtGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCt                                                                             <  1:443337/94‑1 (MQ=255)
                 cGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCt                                                                             <  1:3409114/100‑1 (MQ=255)
                  gATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCtt                                                                            >  1:888710/1‑100 (MQ=255)
                   aTACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTa                                                                           >  2:1691791/1‑100 (MQ=255)
                    tACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTAt                                                                          <  2:1846327/100‑1 (MQ=255)
                    tACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTAt                                                                          <  2:1040760/100‑1 (MQ=255)
                     aCGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATg                                                                         <  2:1221555/100‑1 (MQ=255)
                     aCGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATg                                                                         <  2:1612205/100‑1 (MQ=255)
                     aCGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATg                                                                         <  2:2189673/100‑1 (MQ=255)
                     aCGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATg                                                                         <  2:37227/100‑1 (MQ=255)
                      cGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGc                                                                        >  1:952213/1‑100 (MQ=255)
                       ggtggtGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCtagatcg                                                                       >  2:443342/1‑94 (MQ=255)
                       ggtggtGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGcc                                                                       <  2:3105135/100‑1 (MQ=255)
                       ggtggtGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGcc                                                                       <  1:1355986/100‑1 (MQ=255)
                       ggtggtGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGcc                                                                       >  1:1391794/1‑100 (MQ=255)
                         tggtgCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAGAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGc                                                                        <  1:1083876/97‑1 (MQ=255)
                         tggtgCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCAc                                                                     <  2:1073540/100‑1 (MQ=255)
                          ggtgCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACt                                                                    >  1:990150/1‑100 (MQ=255)
                          ggtgCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACt                                                                    <  1:808227/100‑1 (MQ=255)
                            tgCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGa                                                                  >  2:2259565/1‑100 (MQ=255)
                            tgCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGa                                                                  <  2:1755/100‑1 (MQ=255)
                             gCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGaa                                                                 >  2:2858649/1‑100 (MQ=255)
                              cTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGaa                                                                 >  1:2719635/1‑99 (MQ=255)
                              cTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAat                                                                >  2:3219256/1‑100 (MQ=255)
                                tGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAatat                                                              >  1:2916273/1‑100 (MQ=255)
                                 gTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATc                                                             >  1:2294654/1‑100 (MQ=255)
                                   tAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCgg                                                           <  2:2488035/100‑1 (MQ=255)
                                   tAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCgg                                                           <  1:1108766/100‑1 (MQ=255)
                                   tAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCgg                                                           >  1:202849/1‑100 (MQ=255)
                                    aaaGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCggg                                                          <  2:2175852/100‑1 (MQ=255)
                                    aaaGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCggg                                                          <  1:2548488/100‑1 (MQ=255)
                                    aaaGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCggg                                                          >  2:563486/1‑100 (MQ=255)
                                    aaaGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCgg                                                           <  2:2798020/99‑1 (MQ=255)
                                      aGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGaa                                                        <  2:3304275/100‑1 (MQ=255)
                                       gcgcCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGaaa                                                       <  1:2602062/100‑1 (MQ=255)
                                        cgcCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGaaaa                                                      <  2:1248049/100‑1 (MQ=255)
                                        cgcCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGaaaa                                                      <  2:1621870/100‑1 (MQ=255)
                                              gcgcAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGt                                                >  2:77919/1‑100 (MQ=255)
                                                gcAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTgg                                              <  2:1657802/100‑1 (MQ=255)
                                                  aaaaGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGgcg                                            >  1:663730/1‑100 (MQ=255)
                                                  aaaaGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGgcg                                            <  1:775896/100‑1 (MQ=255)
                                                     aGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTg                                         <  1:3433289/100‑1 (MQ=255)
                                                      gATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTgg                                        >  1:2015402/1‑100 (MQ=255)
                                                      gATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTgg                                        <  1:1660136/100‑1 (MQ=255)
                                                       aTGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTGATGCCACTGAATATCGGGAAAATAGAGTGGCGCt                                          >  1:2395763/1‑97 (MQ=255)
                                                        tGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGAt                                      >  1:1672560/1‑100 (MQ=255)
                                                        tGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGAt                                      >  1:2062112/1‑100 (MQ=255)
                                                        tGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCAGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGAt                                      >  1:2225913/1‑100 (MQ=255)
                                                          cAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATAc                                    >  2:2508925/1‑100 (MQ=255)
                                                          cAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATAc                                    >  1:2707703/1‑100 (MQ=255)
                                                          cAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATAc                                    <  1:1076513/100‑1 (MQ=255)
                                                           aaaaGTCGAAACCCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATAcc                                   <  1:2506768/100‑1 (MQ=255)
                                                             aaGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCaa                                 <  2:887815/100‑1 (MQ=255)
                                                               gTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATa                               <  2:61097/100‑1 (MQ=255)
                                                                tCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATAc                              >  1:346592/1‑100 (MQ=255)
                                                                tCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATAc                              <  1:559778/100‑1 (MQ=255)
                                                                tCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATAc                              <  1:2845814/100‑1 (MQ=255)
                                                                  gAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATAccc                            >  2:1846936/1‑100 (MQ=255)
                                                                  gAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATAccc                            >  2:680480/1‑100 (MQ=255)
                                                                     aaCCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTgg                         <  2:886113/100‑1 (MQ=255)
                                                                       ccAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCt                       <  2:1426517/100‑1 (MQ=255)
                                                                       ccAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCt                       >  2:2564707/1‑100 (MQ=255)
                                                                       ccAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCt                       >  2:1586019/1‑100 (MQ=255)
                                                                        cAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTg                      <  2:202852/100‑1 (MQ=255)
                                                                        cAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTg                      <  2:1990357/100‑1 (MQ=255)
                                                                         aGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGa                     >  1:442485/1‑100 (MQ=255)
                                                                         aGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGa                     >  2:1276231/1‑100 (MQ=255)
                                                                         aGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGa                     <  1:2169868/100‑1 (MQ=255)
                                                                         aGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGa                     <  1:3077730/100‑1 (MQ=255)
                                                                          gACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGAt                    <  1:2332717/100‑1 (MQ=255)
                                                                            cGGGGGTTCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATaa                  <  2:2015431/100‑1 (MQ=255)
                                                                            cGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATaa                  >  1:2549955/1‑100 (MQ=255)
                                                                               gggTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGt               >  2:475814/1‑100 (MQ=255)
                                                                                ggTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTc              <  2:2885019/100‑1 (MQ=255)
                                                                                 gTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCg             <  1:2213366/100‑1 (MQ=255)
                                                                                  tGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGa            <  1:475809/100‑1 (MQ=255)
                                                                                     gTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGAttt         >  2:2657866/1‑100 (MQ=255)
                                                                                     gTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGAttt         <  1:3398885/100‑1 (MQ=255)
                                                                                     gTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGAttt         >  2:827473/1‑100 (MQ=255)
                                                                                     gTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGAttt         <  2:1164418/100‑1 (MQ=255)
                                                                                       aCCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAg       >  2:25903/1‑100 (MQ=255)
                                                                                       aCCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAg       >  1:71397/1‑100 (MQ=255)
                                                                                        ccGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGa      <  1:2940011/100‑1 (MQ=255)
                                                                                            cTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAAc  <  1:77918/100‑1 (MQ=255)
                                                                                                |                                                                                               
CTGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAAC  >  NC_000913/4547223‑4547414

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