Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,453,265 C→T pseudogene (1737/2646 nt) yfcU ← pseudogene, FimD fimbrial export usher family,putative membrane, Not classified, putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,453,2650CT100.0% 331.1 / NA 100pseudogene (1737/2646 nt)yfcUpseudogene, FimD fimbrial export usher family,putative membrane, Not classified, putative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (47/53);  total (47/53)

GCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATCTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTAGTTTGTCTGCTCCTCACGATCC  >  NC_000913/2453130‑2453397
                                                                                                                                       |                                                                                                                                    
gCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTAc                                                                                                                                   <  1:1196684/139‑1 (MQ=255)
             tGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCg                                                                                                                      <  2:1528864/139‑1 (MQ=255)
             tGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCg                                                                                                                      <  2:1123054/139‑1 (MQ=255)
                tAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACa                                                                                                                                  <  1:1887963/124‑1 (MQ=255)
                tAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACa                                                                                                                                  >  2:1887963/1‑124 (MQ=255)
                  gTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGACCGTCCCACTGCCATAGTTGCCGTTATAGCTATCGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCgg                                                                                                                 <  1:337940/139‑1 (MQ=255)
                   tGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCg                                                                                                                  <  1:335834/137‑1 (MQ=255)
                   tGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCg                                                                                                                  >  2:335834/1‑137 (MQ=255)
                      gTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGccc                                                                                                                   <  1:1273830/133‑1 (MQ=255)
                      gTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGccc                                                                                                                   >  2:1273830/1‑133 (MQ=255)
                      gTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGt                                                                                                             >  1:837580/1‑139 (MQ=255)
                       tCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTc                                                                                                            <  1:690331/139‑1 (MQ=255)
                        cgcgTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCa                                                                                                           <  2:38624/139‑1 (MQ=255)
                         gcgTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCAt                                                                                                          >  2:1168585/1‑139 (MQ=255)
                          cgTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGcc                                                                                                                    <  1:1814186/128‑1 (MQ=255)
                          cgTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGcc                                                                                                                    >  2:1814186/1‑128 (MQ=255)
                               tCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGt                                                                                                    >  2:1930639/1‑139 (MQ=255)
                                        cTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGc                                                                                           <  1:921776/139‑1 (MQ=255)
                                         tAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGcc                                                                                          <  2:1548448/139‑1 (MQ=255)
                                           aaaTAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTAc                                                                                                                                   >  1:1920924/1‑96 (MQ=255)
                                           aaaTAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTAc                                                                                                                                   <  2:1920924/96‑1 (MQ=255)
                                                aCCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAg                                                                                            >  1:833023/1‑130 (MQ=255)
                                                aCCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAg                                                                                            <  2:833023/130‑1 (MQ=255)
                                                 ccGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTAt                                                                                                                         <  2:231469/100‑1 (MQ=255)
                                                 ccGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTAt                                                                                                                         >  1:231469/1‑100 (MQ=255)
                                                 ccGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAggg                                                                                  <  2:1976096/139‑1 (MQ=255)
                                                  cGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCgata                                                                                             <  1:1806373/127‑4 (MQ=255)
                                                  cGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCgata                                                                                             >  2:1806373/1‑124 (MQ=255)
                                                        ggCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCt                                                                           <  2:837580/139‑1 (MQ=255)
                                                         gctgctGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCa                                                                                                           >  2:1748365/1‑106 (MQ=255)
                                                         gctgctGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCa                                                                                                           <  1:1748365/106‑1 (MQ=255)
                                                         gctgctGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGcc                                                                                          <  1:692852/123‑1 (MQ=255)
                                                         gctgctGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGcc                                                                                          >  2:692852/1‑123 (MQ=255)
                                                           tgctgTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCat                                                                        <  1:1168585/139‑1 (MQ=255)
                                                                   gtccCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGaa                                                                >  2:1412726/1‑139 (MQ=255)
                                                                   gtccCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGaa                                                                >  1:83279/1‑139 (MQ=255)
                                                                          tGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGcc                                                         <  1:1952979/139‑1 (MQ=255)
                                                                          tGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCCTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGcc                                                         <  1:1452801/139‑1 (MQ=255)
                                                                            ccATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCa                                                       >  2:921347/1‑139 (MQ=255)
                                                                                aGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGcc                                                         <  1:1660195/133‑1 (MQ=255)
                                                                                aGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGcc                                                         >  2:1660195/1‑133 (MQ=255)
                                                                                  ttGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGa                                                 <  1:2012015/139‑1 (MQ=255)
                                                                                    gCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCgg                                                                                                                 <  1:1500819/73‑1 (MQ=255)
                                                                                    gCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCgg                                                                                                                 >  2:1500819/1‑73 (MQ=255)
                                                                                    gCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAAtg                                                              <  2:51914/124‑1 (MQ=255)
                                                                                    gCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAAtg                                                              <  2:1400202/124‑1 (MQ=255)
                                                                                    gCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAAtg                                                              >  1:1400202/1‑124 (MQ=255)
                                                                                    gCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAAtg                                                              >  1:51914/1‑124 (MQ=255)
                                                                                      cGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCt                                                                           <  1:1139669/109‑1 (MQ=255)
                                                                                      cGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCt                                                                           >  2:1139669/1‑109 (MQ=255)
                                                                                      cGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATa                                             >  1:486374/1‑139 (MQ=255)
                                                                                          ataGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAAtgc                                                             >  2:452933/1‑119 (MQ=255)
                                                                                          ataGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAAtgc                                                             <  1:452933/119‑1 (MQ=255)
                                                                                                aaCGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGcc                                                                                                                    <  2:1755811/58‑1 (MQ=255)
                                                                                                aaCGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGcc                                                                                                                    >  1:1755811/1‑58 (MQ=255)
                                                                                                 aCGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATTCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCa                                  >  1:1457437/1‑139 (MQ=255)
                                                                                                     tACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTAc                                                                                                                                   <  2:477789/38‑1 (MQ=255)
                                                                                                     tACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTAc                                                                                                                                   >  1:477789/1‑38 (MQ=255)
                                                                                                        tgttgtCGCCCCACGGCATACTGAGCGAAATTTAc                                                                                                                                   >  1:713953/1‑35 (MQ=255)
                                                                                                        tgttgtCGCCCCACGGCATACTGAGCGAAATTTAc                                                                                                                                   <  2:713953/35‑1 (MQ=255)
                                                                                                            gtCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCa                                                                                                           <  2:1962244/55‑1 (MQ=255)
                                                                                                            gtCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCa                                                                                                           >  1:1962244/1‑55 (MQ=255)
                                                                                                            gtCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTAg                       >  1:252978/1‑139 (MQ=255)
                                                                                                             tCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTAGt                      >  1:320268/1‑139 (MQ=255)
                                                                                                               gCCCCACGGCATACTGAGCGAAATTTACATGCCCtt                                                                                                                           >  1:1918630/1‑36 (MQ=255)
                                                                                                               gCCCCACGGCATACTGAGCGAAATTTACATGCCCtt                                                                                                                           <  2:1918630/36‑1 (MQ=255)
                                                                                                               gCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGa                                                                                 <  1:42361/78‑1 (MQ=255)
                                                                                                               gCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGa                                                                                 >  2:42361/1‑78 (MQ=255)
                                                                                                               gCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACg                                                                             >  2:1609405/1‑82 (MQ=255)
                                                                                                               gCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACg                                                                             <  1:1609405/82‑1 (MQ=255)
                                                                                                               gCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTAGttt                    >  2:1894414/1‑139 (MQ=255)
                                                                                                                 cccACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTAGTTTGt                  <  2:83279/139‑1 (MQ=255)
                                                                                                                 cccACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTAGTTTGt                  <  2:1559112/139‑1 (MQ=255)
                                                                                                                    aCGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTAc                                                                   <  2:1508907/87‑1 (MQ=255)
                                                                                                                    aCGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTAc                                                                   >  1:1508907/1‑87 (MQ=255)
                                                                                                                    aCGGCATACTGAGCGAAATTTACATCCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGcc                                                         <  2:1401139/97‑1 (MQ=255)
                                                                                                                    aCGGCATACTGAGCGAAATTTACATCCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGcc                                                         >  1:1401139/1‑97 (MQ=255)
                                                                                                                            cTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTAGTTTGTCTGCTCCTCAc       <  2:320268/139‑1 (MQ=255)
                                                                                                                            cTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTAGTTTGTCTGCTCCTCAc       <  1:1930639/139‑1 (MQ=255)
                                                                                                                             tGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGa                                                                                 <  2:1816277/64‑1 (MQ=255)
                                                                                                                             tGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGa                                                                                 >  1:1816277/1‑64 (MQ=255)
                                                                                                                             tGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTAGTTTGTCTGCTCCTCACg      <  2:252978/139‑1 (MQ=255)
                                                                                                                              gAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAg                                                                                            >  2:1326514/1‑52 (MQ=255)
                                                                                                                              gAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAg                                                                                            <  1:1326514/52‑1 (MQ=255)
                                                                                                                              gAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTgg                                         >  2:398317/1‑103 (MQ=255)
                                                                                                                              gAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTgg                                         <  1:398317/103‑1 (MQ=255)
                                                                                                                               aGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGc                                                                                           >  2:1771502/1‑52 (MQ=255)
                                                                                                                               aGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGc                                                                                           <  1:1771502/52‑1 (MQ=255)
                                                                                                                               aGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCt                                                                           >  2:982438/1‑68 (MQ=255)
                                                                                                                               aGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCt                                                                           <  1:982438/68‑1 (MQ=255)
                                                                                                                               aGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCa                                                       >  1:1364334/1‑88 (MQ=255)
                                                                                                                               aGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCa                                                       <  2:1364334/88‑1 (MQ=255)
                                                                                                                                gCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTAGTTTGTCTGCTCCTCACGATc   >  1:40088/1‑139 (MQ=255)
                                                                                                                                 cGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTAGTTTGTCTGCTCCTCACGATcc  >  2:786431/1‑139 (MQ=255)
                                                                                                                                  gAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGc                                                                            <  1:1493835/64‑1 (MQ=255)
                                                                                                                                  gAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGc                                                                            >  2:1493835/1‑64 (MQ=255)
                                                                                                                                  gAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAg                                            >  2:1887971/1‑96 (MQ=255)
                                                                                                                                  gAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAg                                            <  1:1887971/96‑1 (MQ=255)
                                                                                                                                    aaTTTACATGCCCTTATCCGCCCGGTTGTCATACTCg                                                                                                     >  2:457899/1‑37 (MQ=255)
                                                                                                                                    aaTTTACATGCCCTTATCCGCCCGGTTGTCATACTCg                                                                                                     <  1:457899/37‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                    
GCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATCTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTAGTTTGTCTGCTCCTCACGATCC  >  NC_000913/2453130‑2453397

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