Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,547,319 C→A 23.3% 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,3190CA23.3% 439.4 / 127.5 217R742S (CGT→AGT) fimDfimbrial usher outer membrane porin protein; FimCD chaperone‑usher
Reads supporting (aligned to +/- strand):  ref base C (78/87);  new base A (24/27);  total (103/114)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.02e-01

Only 100 of 221 total aligned reads displayed.

CGCTGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGG  >  NC_000913/4547221‑4547418
                                                                                                  |                                                                                                   
 gCTGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGt                                                                                                   >  2:465066/1‑100 (MQ=255)
 gCTGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGt                                                                                                   <  1:1279070/100‑1 (MQ=255)
 gCTGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGt                                                                                                   <  2:2099391/100‑1 (MQ=255)
   tGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTgg                                                                                                 <  2:2925586/100‑1 (MQ=255)
      ggCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTa                                                                                              <  1:3250025/100‑1 (MQ=255)
      ggCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTa                                                                                              >  2:2621404/1‑100 (MQ=255)
      ggCAGCCGGTAAACGATACGGTGGTGCTTGTTAAAGCGCCGGACGCGAAAGATGCAAAAGTCGAAAACCAGAGGAGGGGACGTACGGACTGGCg                                                                                                    >  2:2462885/1‑94 (MQ=255)
        cAGCCGTTAACCGATACGGTAGTGCATTTTACCGCTCCCGTCGCCAAAGATGCAAAAGTCGAAAACCAGACGGCGGGGCATACCGACTGGCGTGGTTATg                                                                                            >  2:181059/1‑100 (MQ=255)
          gCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGTCGTGGTTATGcc                                                                                          <  1:2724343/100‑1 (MQ=255)
          gCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGcc                                                                                          >  2:1757904/1‑100 (MQ=255)
            cGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGcc                                                                                          >  1:900200/1‑98 (MQ=255)
            cGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGt                                                                                        >  1:792939/1‑100 (MQ=255)
            cGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGt                                                                                        >  1:570784/1‑100 (MQ=255)
             gTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGtg                                                                                       >  1:1079812/1‑100 (MQ=255)
              ttAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGtgc                                                                                      >  1:2891226/1‑100 (MQ=255)
              ttAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGtgc                                                                                      >  2:1059803/1‑100 (MQ=255)
              ttAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGtgc                                                                                      <  2:2136822/100‑1 (MQ=255)
               tAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGtgct                                                                                     >  1:2614528/1‑100 (MQ=255)
                aaaCGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGtgctg                                                                                    <  1:1167340/100‑1 (MQ=255)
                  aCGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGcc                                                                                  <  1:3174045/100‑1 (MQ=255)
                      tACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTAt                                                                              <  1:2183012/100‑1 (MQ=255)
                      tACGGTGGGGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGAGGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCt                                                                                 <  2:1850990/97‑1 (MQ=255)
                       aCGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTAt                                                                              >  2:1478164/1‑99 (MQ=255)
                        cGGTGGTGCTTGTTAAAGCACCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGc                                                                            <  1:3983988/100‑1 (MQ=255)
                         ggtggtGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGcc                                                                           <  2:3109994/100‑1 (MQ=255)
                          gtggtgCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCa                                                                          <  1:2516849/100‑1 (MQ=255)
                             gtgCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTg                                                                       <  1:2472236/100‑1 (MQ=255)
                             gtgCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTg                                                                       <  1:3472725/100‑1 (MQ=255)
                             gtgCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTg                                                                       <  1:4046074/100‑1 (MQ=255)
                              tgCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACt                                                                        >  2:474225/1‑98 (MQ=255)
                               gCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGaa                                                                     >  2:1955564/1‑100 (MQ=255)
                               gCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGaa                                                                     >  2:36244/1‑100 (MQ=255)
                               gCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGa                                                                      >  2:1303868/1‑99 (MQ=255)
                               gCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGaa                                                                     >  1:679333/1‑100 (MQ=255)
                               gCTTGTTAAAGCGCCTGGCGCAAAAGATGAAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTg                                                                       <  2:2394971/98‑1 (MQ=255)
                                  tGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAatat                                                                  >  2:3668912/1‑100 (MQ=255)
                                    ttAAAGCGCCTGTCGCAAAATATGCAAAAGTCGTAAACCAGACGGGGGTACGCACCGACTGGCGTGGTTCTGCCGTGCTGCCTTATGCCACTGaacgtcg                                                                >  2:2737867/1‑95 (MQ=255)
                                    ttAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCg                                                                <  1:336780/100‑1 (MQ=255)
                                    ttAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCg                                                                <  2:9535/100‑1 (MQ=255)
                                     tAAAGCGCGTGGCACAAAAGATGCAAAAGTGGAAAACCAGACGTGGGTGGGTACCGACTGTGGTGGTTATTCCGTGCGGTCTTATGCCGCTGAATATCgg                                                               >  2:3445699/1‑100 (MQ=255)
                                     tAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCgg                                                               >  1:382731/1‑100 (MQ=255)
                                     tAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCgg                                                               <  1:164185/100‑1 (MQ=255)
                                     tAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCgg                                                               <  1:3746996/100‑1 (MQ=255)
                                     tAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCgg                                                               <  1:3441231/100‑1 (MQ=255)
                                      aaaGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCggg                                                              >  2:1826258/1‑100 (MQ=255)
                                       aaGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCGACTGAATATCggg                                                              >  1:2663460/1‑99 (MQ=255)
                                       aaGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGa                                                             >  1:3001365/1‑100 (MQ=255)
                                       aaGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTCGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGa                                                             <  2:2218751/100‑1 (MQ=255)
                                       aaGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGa                                                             >  2:3952109/1‑100 (MQ=255)
                                        aGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGaa                                                            >  2:648617/1‑100 (MQ=255)
                                          cgcCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGaaaa                                                          <  1:1972234/100‑1 (MQ=255)
                                          cgcCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGaaa                                                           >  2:3323957/1‑99 (MQ=255)
                                           gcCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCATACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGaaa                                                           >  2:233451/1‑98 (MQ=255)
                                            ccTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATa                                                        <  2:3354600/100‑1 (MQ=255)
                                              tGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATaga                                                      >  2:4027053/1‑100 (MQ=255)
                                               ggCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATagag                                                     >  2:1731993/1‑100 (MQ=255)
                                               ggCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATagag                                                     <  2:4109323/100‑1 (MQ=255)
                                               ggCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATagag                                                     >  1:1930203/1‑100 (MQ=255)
                                                    aaaaGATGAAAAAGTCGAAAACGAGTCGGAGGTGCGTACCGACTGGCGTGGTGATGCCGTGCTGGCTTATGCCACTGAATATCGGGAAAATagag                                                     <  2:2207780/95‑1 (MQ=255)
                                                     aaaGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGgcgc                                               >  1:4131341/1‑100 (MQ=255)
                                                       aGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTg                                             >  2:2353463/1‑100 (MQ=255)
                                                           gCAAAAGTCGAAAACCAGACTGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCAACTGAATATCCGGAAAATAGACTGGTGCTGGAt                                          <  2:2159128/99‑1 (MQ=255)
                                                           gCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATa                                         <  1:606941/100‑1 (MQ=255)
                                                           gCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATa                                         >  2:1229781/1‑100 (MQ=255)
                                                            cAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGCTAc                                        >  2:162614/1‑100 (MQ=255)
                                                            cAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATAc                                        >  2:1420062/1‑100 (MQ=255)
                                                            cAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATAc                                        >  2:2578481/1‑100 (MQ=255)
                                                            cAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATAc                                        <  2:679327/100‑1 (MQ=255)
                                                            cAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATAc                                        >  2:1019523/1‑100 (MQ=255)
                                                             aaaaGTCGAAAACCAGCCGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATAcc                                       <  2:3784514/100‑1 (MQ=255)
                                                             aaaaGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTAGAtcag                                       >  2:2511201/1‑97 (MQ=255)
                                                              aaaGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATAc                                        >  2:549871/1‑98 (MQ=255)
                                                                aGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGagatcgga                                    >  1:1858156/1‑93 (MQ=255)
                                                                aGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAAt                                    <  1:3109358/100‑1 (MQ=255)
                                                                 gTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATa                                   >  2:877020/1‑100 (MQ=255)
                                                                   cGAAAACCAGACGGGGGTGCGTACCGACTGGAGAGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATAcc                                 >  2:2709255/1‑100 (MQ=255)
                                                                    gAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATAccc                                >  1:145064/1‑100 (MQ=255)
                                                                      aaaCCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTg                              <  1:2215649/100‑1 (MQ=255)
                                                                       aaCCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTgg                             >  1:1778112/1‑100 (MQ=255)
                                                                       aaCCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTgg                             <  1:3323977/100‑1 (MQ=255)
                                                                        aCCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGc                            <  1:25643/100‑1 (MQ=255)
                                                                        aCCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGc                            <  2:4105519/100‑1 (MQ=255)
                                                                           aGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGa                         >  1:734927/1‑100 (MQ=255)
                                                                             aCGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATa                       <  2:1106991/100‑1 (MQ=255)
                                                                             aCGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATa                       <  2:32485/100‑1 (MQ=255)
                                                                                ggggTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGTAAATAGCGTGGCGCTGGATACCAATACCCTGGCTGATAACg                    >  2:3458743/1‑100 (MQ=255)
                                                                                  ggTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTc                  <  1:2952567/100‑1 (MQ=255)
                                                                                      cGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGAtt              >  1:468384/1‑100 (MQ=255)
                                                                                        tACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTa            <  2:1517100/100‑1 (MQ=255)
                                                                                        tACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATTGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGAttt             >  1:2415406/1‑99 (MQ=255)
                                                                                         aCCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAg           <  2:1900756/100‑1 (MQ=255)
                                                                                         aCCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAg           <  1:3145158/100‑1 (MQ=255)
                                                                                         aCCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAg           <  2:1458000/100‑1 (MQ=255)
                                                                                          ccGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGa          >  1:1311710/1‑100 (MQ=255)
                                                                                              cTGGAGTGGTTATGCCGTGCTGCCTGATGACACTGAATATCCGGAAAAAAGAGGGGCGTGGGATTGCAATACCCTGGCTGATAACGTCGATTtatata        <  2:1362878/98‑1 (MQ=255)
                                                                                                ggCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAAcgc    >  1:2945821/1‑100 (MQ=255)
                                                                                                ggCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAAcgc    <  1:2353476/100‑1 (MQ=255)
                                                                                                 gCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAAcgc    >  2:2701510/1‑99 (MQ=255)
                                                                                                  cGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCgg  <  1:3886426/100‑1 (MQ=255)
                                                                                                  cGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCgg  >  1:3611972/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
CGCTGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGG  >  NC_000913/4547221‑4547418

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: