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% 309.1 / NA 94pseudogene (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 (46/48);  total (46/48)

ATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATCTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTAGTTTGTCTGCTCCTCACGATCCC  >  NC_000913/2453145‑2453398
                                                                                                                        |                                                                                                                                     
aTAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTAc                                                                                                                                    >  1:1062574/1‑124 (MQ=255)
aTAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTAc                                                                                                                                    <  2:1062574/124‑1 (MQ=255)
aTAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGcc                                                                                                                     >  2:84191/1‑139 (MQ=255)
aTAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGcc                                                                                                                     >  2:69645/1‑139 (MQ=255)
   gTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCgg                                                                                                                  <  1:1022066/139‑1 (MQ=255)
   gTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCgg                                                                                                                  <  2:1301717/139‑1 (MQ=255)
   gTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCgg                                                                                                                  <  1:84191/139‑1 (MQ=255)
   gTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCgg                                                                                                                  <  1:69645/139‑1 (MQ=255)
    tGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTAc                                                                                                                                    <  2:3882/120‑1 (MQ=255)
    tGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTAc                                                                                                                                    >  1:3882/1‑120 (MQ=255)
    tGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCGAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGt                                                                                                                 >  1:1054843/1‑139 (MQ=255)
       gTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGt                                                                                                              >  2:1131389/1‑139 (MQ=255)
           cgTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATa                                                                                                          >  1:1412607/1‑139 (MQ=255)
              cATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGccc                                                                                                                    >  2:256622/1‑126 (MQ=255)
              cATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGccc                                                                                                                    <  1:256622/126‑1 (MQ=255)
               aTCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCg                                                                                                      <  2:1248222/139‑1 (MQ=255)
               aTCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCg                                                                                                      <  2:153729/139‑1 (MQ=255)
                  gACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAg                                                                                                   <  2:543230/139‑1 (MQ=255)
                       ggCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGGCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTa                                                                                              >  1:841249/1‑139 (MQ=255)
                        gCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAg                                                                                             <  1:1131389/139‑1 (MQ=255)
                                  ccGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAggg                                                                                   >  2:1141155/1‑139 (MQ=255)
                                   cGACCTGGCTGCTGTCCATCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCATGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGa                                                                                  >  2:918329/1‑139 (MQ=255)
                                      ccTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGt                                                                                                              >  2:253813/1‑108 (MQ=255)
                                      ccTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGt                                                                                                              <  1:253813/108‑1 (MQ=255)
                                         ggCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCt                                                                            <  1:1362778/139‑1 (MQ=255)
                                         ggCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCt                                                                            >  1:447589/1‑139 (MQ=255)
                                         ggCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCt                                                                            <  2:584733/139‑1 (MQ=255)
                                            tgctgTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCat                                                                         <  1:133207/139‑1 (MQ=255)
                                            tgctgTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCat                                                                         <  1:495571/139‑1 (MQ=255)
                                                    gtccCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCgg                                                                                                                  <  1:134778/90‑1 (MQ=255)
                                                    gtccCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCgg                                                                                                                  >  2:134778/1‑90 (MQ=255)
                                                    gtccCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCgg                                                                                                                  <  1:159597/90‑1 (MQ=255)
                                                    gtccCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCgg                                                                                                                  >  2:159597/1‑90 (MQ=255)
                                                    gtccCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAg                                                                                             <  2:1366230/111‑1 (MQ=255)
                                                    gtccCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAg                                                                                             >  1:1366230/1‑111 (MQ=255)
                                                      cccACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCg                                                                                                                   >  1:1071906/1‑87 (MQ=255)
                                                      cccACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCg                                                                                                                   <  2:1071906/87‑1 (MQ=255)
                                                          cTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAAtgctgc                                                           >  2:695820/1‑139 (MQ=255)
                                                           tGCCATGGTTGCCGTTATAGCGAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTGGTCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAGCGCTCATATTACGAATGCTGcc                                                          <  1:1291468/139‑1 (MQ=255)
                                                           tGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCt                                                                            >  2:317133/1‑121 (MQ=255)
                                                           tGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCt                                                                            <  1:317133/121‑1 (MQ=255)
                                                           tGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGcc                                                          <  2:339012/139‑1 (MQ=255)
                                                             ccATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCa                                                        <  2:602995/139‑1 (MQ=255)
                                                                tAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTAc                                                                                                                                    >  2:1455947/1‑60 (MQ=255)
                                                                tAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTAc                                                                                                                                    <  1:1455947/60‑1 (MQ=255)
                                                                   ttGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGa                                                  <  2:841249/139‑1 (MQ=255)
                                                                     gCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGcc                                                                                                                               <  2:340491/60‑1 (MQ=255)
                                                                     gCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGcc                                                                                                                               >  1:340491/1‑60 (MQ=255)
                                                                     gCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAggg                                                                                   >  1:864503/1‑104 (MQ=255)
                                                                     gCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAggg                                                                                   <  2:864503/104‑1 (MQ=255)
                                                                      ccGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAAt                                               >  1:217715/1‑139 (MQ=255)
                                                                       cGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATa                                              >  1:956920/1‑139 (MQ=255)
                                                                       cGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATa                                              <  1:918329/139‑1 (MQ=255)
                                                                                 aaCGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACg                                                                              >  2:1178301/1‑97 (MQ=255)
                                                                                 aaCGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACg                                                                              <  1:1178301/97‑1 (MQ=255)
                                                                                   cGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTCCATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTGGCCAGTCAGGGAAACGCTCATGTTACGAATGCTGCCCATATTGAACTAGTGGGAGAGAat                                  >  2:1256393/1‑139 (MQ=255)
                                                                                   cGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAg                                                                                             >  2:943602/1‑80 (MQ=255)
                                                                                   cGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAg                                                                                             <  1:943602/80‑1 (MQ=255)
                                                                                   cGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCat                                                       >  2:913077/1‑118 (MQ=255)
                                                                                   cGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCat                                                       <  1:913077/118‑1 (MQ=255)
                                                                                   cGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGagag                                     >  1:306684/1‑136 (MQ=255)
                                                                                   cGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGagag                                     <  2:306684/136‑1 (MQ=255)
                                                                                    ggTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCatg                                 <  2:33103/139‑1 (MQ=255)
                                                                                    ggTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCatg                                 <  1:695820/139‑1 (MQ=255)
                                                                                           ttgtCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAAGAGTGGGAGAGCATGAtgttgt                          >  2:862633/1‑139 (MQ=255)
                                                                                            tgtCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTa                         >  1:1322041/1‑139 (MQ=255)
                                                                                             gtCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGaaa                                                                                <  1:902599/83‑1 (MQ=255)
                                                                                             gtCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATACGCCCGGTTGTCATACTCGTAGCGGTAGCCATTCAGGGaaa                                                                                >  2:902599/1‑83 (MQ=255)
                                                                                              tCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTAGt                       >  2:504598/1‑139 (MQ=255)
                                                                                              tCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTAGt                       >  1:815486/1‑139 (MQ=255)
                                                                                               cGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTGTAGttt                     >  2:717181/1‑139 (MQ=255)
                                                                                                gCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGcc                                                                                                                     >  1:1214233/1‑43 (MQ=255)
                                                                                                gCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGcc                                                                                                                     <  2:1214233/43‑1 (MQ=255)
                                                                                                gCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGcc                                                          <  1:217838/102‑1 (MQ=255)
                                                                                                gCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGcc                                                          >  2:217838/1‑102 (MQ=255)
                                                                                                gCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTAGttt                     >  1:1175543/1‑139 (MQ=255)
                                                                                                gCCCCACGGCATACTGAGCGAAATTTACATGCACTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTAGttt                     >  1:872712/1‑139 (MQ=255)
                                                                                                     aCGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCt                                                                            <  2:605968/79‑1 (MQ=255)
                                                                                                     aCGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCt                                                                            >  1:605968/1‑79 (MQ=255)
                                                                                                     aCGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTAGTTTGTCTg                >  2:551532/1‑139 (MQ=255)
                                                                                                     aCGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTAGTTTGTCTg                >  1:481406/1‑139 (MQ=255)
                                                                                                          aTACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTAGTTTGTCTGctcct           <  2:1054843/139‑1 (MQ=255)
                                                                                                                aGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGt                                                                                        >  2:1297795/1‑56 (MQ=255)
                                                                                                                aGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGt                                                                                        <  1:1297795/56‑1 (MQ=255)
                                                                                                                aGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGcc                                                          <  1:309710/86‑1 (MQ=255)
                                                                                                                aGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGcc                                                          >  2:309710/1‑86 (MQ=255)
                                                                                                                aGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATATTGGGAGAGCATGATGCTGTAGTTTGTCTGCTCCTCACGAt     <  2:897500/139‑1 (MQ=255)
                                                                                                                aGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTCGTTTGTCTGCTCCTCACGAt     <  1:862633/139‑1 (MQ=255)
                                                                                                                aGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTAGTTTGTCTGCTCCTCACGAt     <  2:481406/139‑1 (MQ=255)
                                                                                                                aGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTAGTTTGTCTGCTCCTCACGAt     <  2:1175543/139‑1 (MQ=255)
                                                                                                                 gCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTAGTTTGTCTGCTCCTCACGATc    >  1:764389/1‑139 (MQ=255)
                                                                                                                   gAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTAGTTTGTCTGCTCCTCACGATccc  <  2:956920/139‑1 (MQ=255)
                                                                                                                     aaTTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGcc                                                          >  2:709619/1‑81 (MQ=255)
                                                                                                                     aaTTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGcc                                                          <  1:709619/81‑1 (MQ=255)
                                                                                                                        |                                                                                                                                     
ATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATCTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTAGTTTGTCTGCTCCTCACGATCCC  >  NC_000913/2453145‑2453398

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