Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,547,319 C→A 100% 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,3190CA96.8% 329.1 / 0.2 105R742S (CGT→AGT) fimDfimbrial usher outer membrane porin protein; FimCD chaperone‑usher
Reads supporting (aligned to +/- strand):  ref base C (3/1);  new base A (47/54);  total (50/55)
Fisher's exact test for biased strand distribution p-value = 3.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 105 total aligned reads displayed.

CTGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACG  >  NC_000913/4547223‑4547415
                                                                                                |                                                                                                
cTGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTg                                                                                               >  2:1779058/1‑100 (MQ=255)
  ggggCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGt                                                                                             <  1:1939961/100‑1 (MQ=255)
  ggggCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGt                                                                                             <  2:2174127/100‑1 (MQ=255)
  ggggCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGt                                                                                             <  1:938013/100‑1 (MQ=255)
   gggCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGtt                                                                                            >  1:13355/1‑100 (MQ=255)
    ggCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTa                                                                                           <  1:460348/100‑1 (MQ=255)
     gCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGGAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTAt                                                                                          <  1:1413648/100‑1 (MQ=255)
      cAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATg                                                                                         >  2:1530467/1‑100 (MQ=255)
        gCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGcc                                                                                       >  1:1881082/1‑100 (MQ=255)
          cGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGt                                                                                     >  1:1745068/1‑100 (MQ=255)
           gTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGtg                                                                                    <  2:1731130/100‑1 (MQ=255)
           gTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGtg                                                                                    >  2:243164/1‑100 (MQ=255)
            ttAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGtgc                                                                                   <  2:2397050/100‑1 (MQ=255)
                atgATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGtgctgc                                                                                <  1:2075231/97‑1 (MQ=255)
                aCGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGtgctgca                                                                               <  2:2563447/100‑2 (MQ=255)
                aCGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGcc                                                                               <  1:1759407/100‑1 (MQ=255)
                  gATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCtt                                                                             >  2:1030463/1‑100 (MQ=255)
                   aTACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTa                                                                            <  1:996157/100‑1 (MQ=255)
                   aTACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTa                                                                            >  1:2004798/1‑100 (MQ=255)
                   aTACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTa                                                                            <  1:1602584/100‑1 (MQ=255)
                    tACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTAt                                                                           >  1:1768948/1‑100 (MQ=255)
                     aCGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATg                                                                          <  2:673228/100‑1 (MQ=255)
                     aCGGTGGTGCTTGTTAAAGCGCCTGGCCCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATg                                                                          <  2:646087/100‑1 (MQ=255)
                      cGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGc                                                                         >  1:2002188/1‑100 (MQ=255)
                      cGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGc                                                                         >  1:2211747/1‑100 (MQ=255)
                       ggtggtGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGcc                                                                        <  1:641035/100‑1 (MQ=255)
                       ggtggtGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGcc                                                                        >  2:64905/1‑100 (MQ=255)
                        gtggtgCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCa                                                                       >  2:882892/1‑100 (MQ=255)
                           gtgCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTg                                                                    <  2:3031792/100‑1 (MQ=255)
                           gtgCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTg                                                                    <  1:243167/100‑1 (MQ=255)
                            tgCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGa                                                                   >  2:2986447/1‑100 (MQ=255)
                             gCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGaa                                                                  <  2:1920605/100‑1 (MQ=255)
                             gCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGaa                                                                  <  2:915021/100‑1 (MQ=255)
                             gCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGaa                                                                  <  1:2472819/100‑1 (MQ=255)
                             gCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGaa                                                                  >  1:1987706/1‑100 (MQ=255)
                              cTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCACTGAata                                                                <  2:2470484/100‑1 (MQ=255)
                               ttGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAata                                                                <  2:432010/100‑1 (MQ=255)
                               ttGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAat                                                                 <  1:1167270/99‑1 (MQ=255)
                                tGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAatat                                                               >  1:181676/1‑100 (MQ=255)
                                ngTTAAAGCGCCTGGCGCAAAAAATGCAAAAGTCGAAAACAAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCAGTGCTGCCTTATGCCACTGAatat                                                               >  2:1791522/2‑100 (MQ=255)
                                  ttAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCg                                                             <  2:2779234/100‑1 (MQ=255)
                                   tAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCgg                                                            <  1:280662/100‑1 (MQ=255)
                                   tAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCgg                                                            <  1:1855189/100‑1 (MQ=255)
                                     aaGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGa                                                          >  1:1607881/1‑100 (MQ=255)
                                     aaGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGa                                                          >  1:1608794/1‑100 (MQ=255)
                                      aGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGaa                                                         <  2:1738555/100‑1 (MQ=255)
                                      aGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGaa                                                         >  1:2914289/1‑100 (MQ=255)
                                      aGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGaa                                                         >  2:1073617/1‑100 (MQ=255)
                                        cgcCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGaaaa                                                       <  2:2764669/100‑1 (MQ=255)
                                          ccTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATa                                                     >  1:2011572/1‑100 (MQ=255)
                                          ccTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATa                                                     <  2:2651411/100‑1 (MQ=255)
                                            tGGCGTAAAAGCTGCCAACGTCGAAAAACAGACGGGGGTGCGTATCGAATGGAGTGGTTATGCCGTGCTGCCTTATGGCAGTGTATATCTGGAAAATagg                                                   <  2:556175/100‑2 (MQ=255)
                                             ggCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATagag                                                  >  2:2033769/1‑100 (MQ=255)
                                             ggCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATagag                                                  <  2:1685554/100‑1 (MQ=255)
                                               cgcAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATACCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTg                                                >  1:2399606/1‑100 (MQ=255)
                                                 cAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGgc                                              >  2:1693313/1‑100 (MQ=255)
                                                  aaaaGATTGAAAAGTAGAAAACGAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCGGTGGTGCCTTATGCCACTGAATATCGGGAGAATAGAGTg                                                <  2:127100/97‑1 (MQ=255)
                                                     aGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGAGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGgcgc                                            >  2:2061651/1‑98 (MQ=255)
                                                     aGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACAGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATGTCGGGAAAATAGAGTGGCGCTg                                          >  2:650125/1‑100 (MQ=255)
                                                        tGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGAt                                       >  2:2487291/1‑100 (MQ=255)
                                                        tGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGAt                                       >  1:1433633/1‑100 (MQ=255)
                                                         gCAAAAGTCGAAAACCAGACGGGTGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATa                                      <  1:2847707/100‑1 (MQ=255)
                                                         gCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATa                                      >  1:657501/1‑100 (MQ=255)
                                                         gCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATa                                      <  2:2613033/100‑1 (MQ=255)
                                                           aataGTTGAAAAAAAGACGAGAGTTTGTATTGACTGGAGTGGTTATGCCGTGCTGACTTATGACACTGTATTTTTGGTAAATAGAGTGGCGATGGa                                        <  2:552649/93‑1 (MQ=255)
                                                           aaaaGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATAcc                                    <  1:2016162/100‑1 (MQ=255)
                                                           aaaaGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATAcc                                    <  2:2157761/100‑1 (MQ=255)
                                                            aaaGTCGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCa                                   <  2:1361859/100‑1 (MQ=255)
                                                                 cGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATAcc                              >  2:649254/1‑100 (MQ=255)
                                                                 cGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATAcc                              >  2:2830227/1‑100 (MQ=255)
                                                                 cGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATAcc                              >  2:838952/1‑100 (MQ=255)
                                                                 cGAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTAGCGCTGGATACCAATAcc                              >  1:1967979/1‑100 (MQ=255)
                                                                  gAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATAccc                             >  1:1418162/1‑100 (MQ=255)
                                                                  gAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATAccc                             >  1:576824/1‑100 (MQ=255)
                                                                  gAAAACCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATAccc                             >  1:2874016/1‑100 (MQ=255)
                                                                  gAAAACCAGACGGGGGTGCGTACCGACTGGAGTAGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGCAAATAGAGTGGCGCTGGATACCAATAccc                             >  2:1823376/1‑100 (MQ=255)
                                                                    aaaCCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTg                           <  1:431668/100‑1 (MQ=255)
                                                                     aaCCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTgg                          <  1:1311658/100‑1 (MQ=255)
                                                                      actacccGGGGGTGCGTACCGACAGGAATGGTTATGCCGTGCTGCCTTATGCCACTGAATTTCGGGAAAATAAAGTGGCGCTGGATACCAATACCCTGGc                         >  2:1122441/7‑100 (MQ=255)
                                                                      aCCAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGc                         <  2:1654329/100‑1 (MQ=255)
                                                                        cAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTg                       <  2:2166397/100‑1 (MQ=255)
                                                                        cAGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTg                       <  2:2621950/100‑1 (MQ=255)
                                                                         aGACGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCGTGGCTGa                      >  2:632431/1‑100 (MQ=255)
                                                                           aCGGGGGTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATa                    <  1:2525246/100‑1 (MQ=255)
                                                                             gggggTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATaa                   <  2:1733026/99‑1 (MQ=255)
                                                                              ggggTGCGTACCAACTGGAGAGGTGATGACGAGCTGCCTTATGAAACTGAATATCGGGAAAATAGACTGAAACTGGATACCAAAAccat                            >  2:2520668/1‑87 (MQ=255)
                                                                              ggggTGCGTACAGACTGGCGTGGTTATGCCGTGCTGCCTTATGACACTGAATATCGGGAAAATAGAGTCGCGCTGGATACCAATAACATGGCTGATAACg                 >  2:2381957/1‑100 (MQ=255)
                                                                                ggTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTc               <  1:1782156/100‑1 (MQ=255)
                                                                                ggTGCGTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTc               >  2:2359004/1‑100 (MQ=255)
                                                                                     gTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGAttt          >  1:1132213/1‑100 (MQ=255)
                                                                                     gTACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGAttt          <  1:1132931/100‑1 (MQ=255)
                                                                                     gTACCGACTGGAGTGGTTATGCCATGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGAttt          >  1:1880219/1‑100 (MQ=255)
                                                                                      tACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTa         <  1:64905/100‑1 (MQ=255)
                                                                                      tACCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTa         <  2:1924996/100‑1 (MQ=255)
                                                                                       aCCGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAAGATCAGGAAAATAGAGTGGCGCTGAATACCAATACCCTGGCTGATAACGTCGATTTAg        >  2:560061/1‑100 (MQ=255)
                                                                                        ccGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGa       <  1:2830259/100‑1 (MQ=255)
                                                                                        ccGACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGa       <  2:165908/100‑1 (MQ=255)
                                                                                          gACTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTGAGATa     <  2:425829/100‑1 (MQ=255)
                                                                                           aCTGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAATGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATaa    >  1:908026/1‑100 (MQ=255)
                                                                                             tGGAGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAAcg  >  1:1672845/1‑100 (MQ=255)
                                                                                                |                                                                                                
CTGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACG  >  NC_000913/4547223‑4547415

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