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% 332.5 / NA 104pseudogene (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 (49/55);  total (49/55)

Only 100 of 104 total aligned reads displayed.

CCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATCTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTAGTTTGTCTGCTCCTCACGAT  >  NC_000913/2453131‑2453395
                                                                                                                                      |                                                                                                                                  
ccGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACa                                                                                                                                >  1:139193/1‑139 (MQ=255)
                aGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCg                                                                                                                >  2:1007428/1‑139 (MQ=255)
                aGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCg                                                                                                                >  1:1233793/1‑139 (MQ=255)
                 gTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCgg                                                                                                               >  1:1342288/1‑139 (MQ=255)
                    aGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTg                                                                                                            >  2:554040/1‑139 (MQ=255)
                       cacgTCAATGACACGGCTAGAATAACCGACCTGGCGGCAGTCCGTGCCACTGCCATAGGTGCCGTTATAGCGAACGGTACGGTTGTCGCCCAACGGCATACTGAGCGAAAGTTACATGCCCTTATCCGCCCGGTTGTCa                                                                                                         <  2:1243311/137‑1 (MQ=255)
                             aTCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCg                                                                                                   <  1:1007428/139‑1 (MQ=255)
                             aTCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCg                                                                                                   <  1:1127718/139‑1 (MQ=255)
                              tCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGt                                                                                                  >  1:1195108/1‑139 (MQ=255)
                                      gCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAg                                                                                          <  2:375710/139‑1 (MQ=255)
                                                 cGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGccc                                                                                                                 <  1:337603/105‑1 (MQ=255)
                                                 cGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGa                                                                               <  1:342193/139‑1 (MQ=255)
                                                  gACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGaa                                                                              >  2:442192/1‑139 (MQ=255)
                                                   aCCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGaaa                                                                             <  2:1370231/139‑1 (MQ=255)
                                                    ccTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGa                                                                               >  1:201030/1‑136 (MQ=255)
                                                    ccTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGa                                                                               <  2:201030/136‑1 (MQ=255)
                                                     cTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACg                                                                           <  2:438105/139‑1 (MQ=255)
                                                     cTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACg                                                                           <  2:1159673/139‑1 (MQ=255)
                                                       ggCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCt                                                                         <  1:1154596/139‑1 (MQ=255)
                                                        gctgctGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTc                                                                        >  2:250428/1‑139 (MQ=255)
                                                              gtccgtccCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCg                                                                                                   <  2:637488/106‑1 (MQ=255)
                                                              gtccgtccCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCg                                                                                                   >  1:637488/1‑106 (MQ=255)
                                                               tccgtccCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCAt                                                                                                        >  1:330037/1‑100 (MQ=255)
                                                               tccgtccCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCAt                                                                                                        <  2:330037/100‑1 (MQ=255)
                                                                  gtccCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCt                                                                         <  2:1302936/128‑1 (MQ=255)
                                                                  gtccCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCt                                                                         >  1:1302936/1‑128 (MQ=255)
                                                                  gtccCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGaa                                                              >  1:42402/1‑139 (MQ=255)
                                                                   tccCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAAt                                                             <  2:1233793/139‑1 (MQ=255)
                                                                      cACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACa                                                                                                                                >  1:1406972/1‑69 (MQ=255)
                                                                      cACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACa                                                                                                                                <  2:1406972/69‑1 (MQ=255)
                                                                      cACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGaaa                                                                             >  2:1185891/1‑120 (MQ=255)
                                                                      cACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGaaa                                                                             <  1:1185891/120‑1 (MQ=255)
                                                                      cACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAAtgct                                                          <  1:442192/139‑1 (MQ=255)
                                                                      cACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAAtgct                                                          >  1:757584/1‑139 (MQ=255)
                                                                         tGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGcc                                                       <  2:1342288/139‑1 (MQ=255)
                                                                         tGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGcc                                                       <  1:1112406/139‑1 (MQ=255)
                                                                         tGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGcc                                                       <  1:250428/139‑1 (MQ=255)
                                                                         tGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGcc                                                       <  2:139193/139‑1 (MQ=255)
                                                                           ccATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCa                                                     >  2:913758/1‑139 (MQ=255)
                                                                              tAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGcc                                                                                                                            <  1:1376385/65‑1 (MQ=255)
                                                                              tAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGcc                                                                                                                            >  2:1376385/1‑65 (MQ=255)
                                                                                   gCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGaaa                                             <  2:757584/139‑1 (MQ=255)
                                                                                   gCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGaaa                                             >  2:410235/1‑139 (MQ=255)
                                                                                     cGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTAc                                                                                                                                 >  2:120967/1‑53 (MQ=255)
                                                                                     cGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTAc                                                                                                                                 <  1:120967/53‑1 (MQ=255)
                                                                                         ataGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTgg                                       >  2:788330/1‑139 (MQ=255)
                                                                                                aCGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAg                                                                                          >  2:1042016/1‑81 (MQ=255)
                                                                                                aCGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAg                                                                                          <  1:1042016/81‑1 (MQ=255)
                                                                                                aCGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCa                                <  1:788330/139‑1 (MQ=255)
                                                                                                 cGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGc                                                                          <  2:456728/96‑1 (MQ=255)
                                                                                                 cGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGc                                                                          >  1:456728/1‑96 (MQ=255)
                                                                                                  ggTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCatg                              <  1:913758/139‑1 (MQ=255)
                                                                                                  ggTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCatg                              <  1:410235/139‑1 (MQ=255)
                                                                                                      cTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGga                                     >  1:932445/1‑128 (MQ=255)
                                                                                                       tgttgtCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACg                                                                           >  1:425175/1‑89 (MQ=255)
                                                                                                       tgttgtCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACg                                                                           <  2:425175/89‑1 (MQ=255)
                                                                                                       tgttgtCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGAtgtt                         <  2:42402/139‑1 (MQ=255)
                                                                                                            tCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGaaa                                                                             <  2:751722/82‑1 (MQ=255)
                                                                                                            tCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGaaa                                                                             >  1:751722/1‑82 (MQ=255)
                                                                                                            tCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTAGt                    >  1:574982/1‑139 (MQ=255)
                                                                                                             cGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTa                                                                                                 <  2:1408264/61‑1 (MQ=255)
                                                                                                             cGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTa                                                                                                 >  1:1408264/1‑61 (MQ=255)
                                                                                                              gCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAg                                                                                          <  2:1089486/67‑1 (MQ=255)
                                                                                                              gCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAg                                                                                          >  1:1089486/1‑67 (MQ=255)
                                                                                                              gCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCatg                              <  2:1288203/127‑1 (MQ=255)
                                                                                                              gCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCatg                              >  1:1288203/1‑127 (MQ=255)
                                                                                                                cccACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTAGTTTGt                <  2:574982/139‑1 (MQ=255)
                                                                                                                   aCGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGa                                                                               <  1:250131/73‑1 (MQ=255)
                                                                                                                   aCGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGa                                                                               >  2:250131/1‑73 (MQ=255)
                                                                                                                   aCGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGa                                               >  1:1002461/1‑105 (MQ=255)
                                                                                                                            tGAGCGAAATTTACATGCCCTTATCCGTCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTAGTTTGTCTGCTCCTCACg    <  2:808928/139‑1 (MQ=255)
                                                                                                                            tGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGt                                                                                     <  2:477448/58‑1 (MQ=255)
                                                                                                                            tGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGt                                                                                     >  1:477448/1‑58 (MQ=255)
                                                                                                                             gAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTAGTTTGTCTGCTCCTCACGa   <  2:1195108/139‑1 (MQ=255)
                                                                                                                              aGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCt                                                                         <  2:34145/68‑1 (MQ=255)
                                                                                                                              aGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCt                                                                         >  1:34145/1‑68 (MQ=255)
                                                                                                                              aGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAAtg                                                            <  2:819540/81‑1 (MQ=255)
                                                                                                                              aGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAAtg                                                            >  1:819540/1‑81 (MQ=255)
                                                                                                                              aGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGa                                               >  2:269699/1‑94 (MQ=255)
                                                                                                                              aGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGAtgttgt                       >  1:832403/1‑118 (MQ=255)
                                                                                                                              aGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTAGTTTGTCTGCTCCTCACGa   <  1:1237879/138‑1 (MQ=255)
                                                                                                                              aGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTAGTTTGTCTGCTCCTCACGa   >  2:1237879/1‑138 (MQ=255)
                                                                                                                              aGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTAGTTTGTCTGCTCCTCACGAt  >  2:1371373/1‑139 (MQ=255)
                                                                                                                              aGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTAGTTTGTCTGCTCCTCACGAt  >  1:179537/1‑139 (MQ=255)
                                                                                                                              aGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTAGTTTGTCTGCTCCTCACGAt  >  1:148986/1‑139 (MQ=255)
                                                                                                                              aGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGGAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGAttttgt                       <  2:832403/118‑1 (MQ=255)
                                                                                                                               gCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCt                                                                         >  1:36625/1‑67 (MQ=255)
                                                                                                                               gCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCt                                                                         <  2:36625/67‑1 (MQ=255)
                                                                                                                               gCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAAt                                                             <  1:631284/79‑1 (MQ=255)
                                                                                                                               gCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAAt                                                             >  2:631284/1‑79 (MQ=255)
                                                                                                                               gCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCa                                >  1:629982/1‑108 (MQ=255)
                                                                                                                               gCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCa                                <  2:629982/108‑1 (MQ=255)
                                                                                                                                cGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGa                                                                               >  2:92065/1‑60 (MQ=255)
                                                                                                                                cGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGa                                                                               <  1:92065/60‑1 (MQ=255)
                                                                                                                                cGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTAGTTTGTCTGCTCCTCAc     <  2:459339/134‑1 (MQ=255)
                                                                                                                                cGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTAGTTTGTCTGCTCCTCAc     >  1:459339/1‑134 (MQ=255)
                                                                                                                                cGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTAGTTTGTCTGCTCCTCAc     <  1:1010723/134‑1 (MQ=255)
                                                                                                                                cGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCAGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTAGTTTGTCTGCTCCTCAc     >  2:1010723/1‑134 (MQ=255)
                                                                                                                                 gAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGATGCCCata                                                   >  2:876250/1‑87 (MQ=255)
                                                                                                                                 gAAATTTACACGCCCTTATCCGCCCGGTTGTCATAATCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCata                                                   <  1:876250/87‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                  
CCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATCTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTAGTTTGTCTGCTCCTCACGAT  >  NC_000913/2453131‑2453395

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