Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 374,295 C→T intergenic (+190/‑343) mhpE → / → mhpT 4‑hydroxy‑2‑ketovalerate aldolase/putative transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913374,2950CT99.3% 399.8 / ‑4.1 148intergenic (+190/‑343)mhpE/mhpT4‑hydroxy‑2‑ketovalerate aldolase/putative transport protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (69/78);  minor base A (0/1);  total (69/79)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.54e-01

Only 100 of 159 total aligned reads displayed.

GTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCG  >  NC_000913/374205‑374383
                                                                                          |                                                                                        
gTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGgtgc                                                                                        <  2:2464871/93‑1 (MQ=37)
  aGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCGCGACTGCCGGATGCGTCGTGAACGCCTTATCCTGCCTACGGATGGTGCgg                                                                                      <  2:2964138/93‑1 (MQ=11)
  aGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCgg                                                                                      <  1:3786868/93‑1 (MQ=37)
  aGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCgg                                                                                      <  2:756242/93‑1 (MQ=37)
  aGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCgg                                                                                      >  1:974525/1‑93 (MQ=37)
    gTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGa                                                                                    >  2:2206485/1‑93 (MQ=25)
      cTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAAt                                                                                  <  2:4062008/93‑1 (MQ=21)
       tGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAAtt                                                                                 >  1:178616/1‑93 (MQ=18)
        gATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAAttt                                                                                <  2:797332/93‑1 (MQ=18)
        gATAAGACGCGTTAGCGCCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAAttt                                                                                <  2:3682856/93‑1 (MQ=16)
         aTAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTg                                                                               <  2:675317/93‑1 (MQ=18)
          tAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGt                                                                              >  2:513677/1‑93 (MQ=18)
          tAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGt                                                                              <  1:2676217/93‑1 (MQ=18)
           aaGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTa                                                                             <  1:456245/93‑1 (MQ=18)
           aaGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTa                                                                             <  2:1148412/93‑1 (MQ=18)
            aGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAg                                                                            >  1:3742953/1‑93 (MQ=18)
            aGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAg                                                                            <  1:991916/93‑1 (MQ=18)
             gACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAgg                                                                           >  2:2546051/1‑93 (MQ=18)
               cgcgTTAGCTTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGcc                                                                         <  2:701017/93‑1 (MQ=16)
                gcgTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCt                                                                        >  2:2987374/1‑93 (MQ=18)
                gcgTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGAGCGGGAATTTGTAGGCCt                                                                        <  2:2332013/93‑1 (MQ=18)
                 cgTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTg                                                                       >  2:3138681/1‑93 (MQ=18)
                  gTTATCTTCGCATCAGCCATCTGCTCCCGACTGCCGGATGCGGCTTGAACGCCTTATCCGGCCTAGGGATGGTGCGGGAATTTGTAGGCCTGa                                                                      <  2:416981/93‑1 (MQ=11)
                  gTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGa                                                                      >  2:1016147/1‑93 (MQ=18)
                  gTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGa                                                                      >  2:1482287/1‑93 (MQ=18)
                     aGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATaa                                                                   >  1:39496/1‑93 (MQ=18)
                      gCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATCAg                                                                  >  1:536769/1‑93 (MQ=16)
                      gCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAg                                                                  >  1:2663814/1‑93 (MQ=18)
                      gCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAg                                                                  <  1:4405492/93‑1 (MQ=18)
                       cGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGa                                                                 <  1:2915196/93‑1 (MQ=18)
                       cGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGa                                                                 <  2:13870/93‑1 (MQ=18)
                        gTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGAc                                                                >  1:2706277/1‑93 (MQ=18)
                        gTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGAc                                                                <  1:2128364/93‑1 (MQ=18)
                         tCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGAcg                                                               >  1:947374/1‑93 (MQ=18)
                         gcGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGAcg                                                               <  1:3921165/92‑1 (MQ=18)
                          cGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGAcgc                                                              <  1:4657352/93‑1 (MQ=18)
                          cGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGAcgc                                                              >  2:1086118/1‑93 (MQ=18)
                             aTCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGtt                                                           >  1:2527633/1‑93 (MQ=18)
                             aTCAGACATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGtt                                                           <  2:3755865/93‑1 (MQ=16)
                              tCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTa                                                          <  1:233158/93‑1 (MQ=18)
                              tCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTa                                                          >  1:1748144/1‑93 (MQ=18)
                              tCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTa                                                          <  1:3857651/93‑1 (MQ=18)
                              gcAGGCATCTGGGCGCGACTGCCGGATGCGGCGTGAACGCGTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTa                                                          <  2:1177409/92‑1 (MQ=11)
                               cAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAg                                                         >  2:617956/1‑93 (MQ=18)
                                aGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGgtgc                                                                                        >  2:4635646/1‑61 (MQ=25)
                                aGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGAcg                                                               >  1:3629469/1‑86 (MQ=18)
                                aGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGc                                                        <  1:211662/93‑1 (MQ=18)
                                aGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGc                                                        <  2:4637079/93‑1 (MQ=18)
                                  gCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGt                                                      >  2:1910435/1‑93 (MQ=18)
                                   cATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTc                                                     >  2:66667/1‑93 (MQ=18)
                                     tCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGc                                                   <  2:3599621/93‑1 (MQ=18)
                                         cgcACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCg                                                       >  1:1299355/1‑85 (MQ=18)
                                          gcACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAg                                              >  2:2053069/1‑93 (MQ=18)
                                           cACGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAgg                                             <  2:2779029/93‑1 (MQ=16)
                                            aCGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGc                                            >  2:2866643/1‑93 (MQ=21)
                                               aCTGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATc                                         <  1:2812881/93‑1 (MQ=16)
                                                cTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCt                                        >  1:9999/1‑93 (MQ=18)
                                                 tGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTg                                       >  2:4250286/1‑93 (MQ=21)
                                                 tGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTg                                       <  1:736000/93‑1 (MQ=18)
                                                  gCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTgc                                      <  1:128185/93‑1 (MQ=18)
                                                    cGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTgcgc                                    <  1:148899/93‑1 (MQ=25)
                                                    cGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTgcgc                                    >  1:4640233/1‑93 (MQ=37)
                                                     ggATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCa                                   <  2:1843753/93‑1 (MQ=37)
                                                      gATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGGGCAc                                  >  1:3026855/1‑93 (MQ=255)
                                                      gATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCAc                                  >  2:3699785/1‑93 (MQ=37)
                                                         ggtgcgTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGAc                               <  1:1235663/90‑1 (MQ=255)
                                                         gcggcgTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGAc                               <  2:3623398/93‑1 (MQ=37)
                                                         gcggcgTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGAc                               <  1:1031607/93‑1 (MQ=37)
                                                          cggcgTGAACGCGTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACt                              <  2:294578/93‑1 (MQ=255)
                                                          cggcgTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACt                              <  1:3823878/93‑1 (MQ=37)
                                                           ggcgTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTg                             >  1:42623/1‑93 (MQ=255)
                                                             cgTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGcc                           >  1:1464135/1‑93 (MQ=37)
                                                              gTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCg                          >  2:2175553/1‑93 (MQ=37)
                                                              gTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCg                          >  1:3189483/1‑93 (MQ=37)
                                                              gTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCg                          >  1:2828199/1‑93 (MQ=37)
                                                               tGAGCGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCgg                         >  2:1987407/1‑93 (MQ=255)
                                                               tGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCgg                         >  1:1558281/1‑93 (MQ=25)
                                                                 aaCGCTTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGAt                       <  2:3105624/93‑1 (MQ=16)
                                                                 aaCGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGAt                       >  1:1199013/1‑93 (MQ=18)
                                                                  aCGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATg                      <  1:4588697/93‑1 (MQ=21)
                                                                  aCGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATg                      >  1:2390692/1‑93 (MQ=18)
                                                                    gCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATgcg                    <  2:3332099/93‑1 (MQ=21)
                                                                    gCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATgcg                    <  1:135247/93‑1 (MQ=21)
                                                                      cTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATgcggc                  >  2:476645/1‑93 (MQ=25)
                                                                      cTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATgaggc                  >  2:207612/1‑93 (MQ=255)
                                                                       ttATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGTCGGATgcggcg                 >  1:1512556/1‑93 (MQ=255)
                                                                       ttATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATgcggcg                 >  2:4019812/1‑93 (MQ=255)
                                                                       ttATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATgcggcg                 >  1:245373/1‑93 (MQ=37)
                                                                        tATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGt                <  2:1616609/93‑1 (MQ=37)
                                                                          tCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGa              >  1:610883/1‑93 (MQ=255)
                                                                            cgtgcTACGGATGGTGCGGTAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAAc            <  1:758070/89‑1 (MQ=14)
                                                                               ccTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGcc         <  2:1327591/93‑1 (MQ=255)
                                                                                cTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCt        <  1:4096645/93‑1 (MQ=255)
                                                                                cTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCt        >  1:3036794/1‑93 (MQ=255)
                                                                                 tACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGc                            >  1:1962270/1‑72 (MQ=21)
                                                                                 tACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCtt       <  2:4033254/93‑1 (MQ=255)
                                                                                 tACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCtt       >  1:1299207/1‑93 (MQ=255)
                                                                                   cGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTAt     <  1:2258611/93‑1 (MQ=255)
                                                                                    ggATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATc    <  1:3269647/93‑1 (MQ=255)
                                                                                     gATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATcc   <  1:612354/93‑1 (MQ=255)
                                                                                          |                                                                                        
GTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCG  >  NC_000913/374205‑374383

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