Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 949,741 Δ1 bp 19.3% coding (74/591 nt) ycaK → putative NAD(P)H‑dependent oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913949,7410C.19.3% 397.6 / 99.4 140coding (74/591 nt)ycaKputative NAD(P)H‑dependent oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (54/59);  new base . (14/13);  total (68/72)
Fisher's exact test for biased strand distribution p-value = 8.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.21e-01

Only 100 of 148 total aligned reads displayed.

CCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGA  >  NC_000913/949642‑949840
                                                                                                   |                                                                                                    
ccAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGc                                                                                                      <  2:105356/100‑1 (MQ=255)
ccAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGc                                                                                                      >  2:3193133/1‑100 (MQ=255)
ccAGAACGTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGAACGCACATATTGCTGATGc                                                                                                      >  1:818459/1‑100 (MQ=255)
    aaCTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCGCGTCATCATTCATTGACCGCACATATTGCTGATGCGATc                                                                                                  >  1:3091130/1‑100 (MQ=255)
      cTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCa                                                                                               <  2:666100/100‑1 (MQ=255)
        tAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATc                                                                                             >  1:2944610/1‑100 (MQ=255)
        tAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATc                                                                                             <  2:2161763/100‑1 (MQ=255)
        tAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATc                                                                                             >  2:3249427/1‑100 (MQ=255)
           cGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGc                                                                                          >  2:4347574/1‑100 (MQ=255)
           cGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGGATCGCATCAGc                                                                                          <  1:1283972/100‑1 (MQ=255)
             tGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCgg                                                                                        <  2:152477/100‑1 (MQ=255)
             tGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGGATCGCATCAGcaa                                                                                        <  2:172457/100‑3 (MQ=255)
              ggAGGTAAAATTATGCAGTTTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCggg                                                                                       <  2:4043809/100‑1 (MQ=255)
              ggAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGGATCGCATCAGcaat                                                                                       <  1:1551302/100‑4 (MQ=255)
              ggAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGGATCGCATCAGcaat                                                                                       <  2:2418098/100‑4 (MQ=255)
              ggAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGGATCGCATCAGcaat                                                                                       <  2:220289/100‑4 (MQ=255)
              ggAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGGATCGCATCAGcaat                                                                                       <  2:1251708/100‑4 (MQ=255)
              ggAGGTAAAATTATGCAGTCTGAACGTATTTAGTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGGATCGCATCAGcaa                                                                                        <  1:435218/99‑3 (MQ=255)
               gAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCgggn                                                                                      <  1:2189764/100‑2 (MQ=255)
               gAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGc                                                                                      <  1:164874/100‑1 (MQ=255)
               gAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGGATCGCATCAGcaata                                                                                      <  2:28036/100‑5 (MQ=255)
               gAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGGATCGCATCAGcaata                                                                                      <  2:84980/100‑5 (MQ=255)
                 ggTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCaa                                                                                    >  1:1021626/1‑100 (MQ=255)
                 ggTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCAATCCATCAGCGGGCaa                                                                                    >  2:3183579/1‑100 (MQ=255)
                  gTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAAt                                                                                   >  1:2804595/1‑100 (MQ=255)
                  gTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAAt                                                                                   <  2:303947/100‑1 (MQ=255)
                    aaaaTTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATgg                                                                                 <  1:30681/100‑1 (MQ=255)
                    aaaaTTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATACATCAGCGGGCAATgg                                                                                 >  1:2577616/1‑100 (MQ=255)
                    aaaaTTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATCGCTGATGGATCGCATCAGcaatatgtgc                                                                                 <  2:2102026/100‑10 (MQ=37)
                     aaaTTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGGATCGCATCAGcaatatgtgcg                                                                                <  2:1101737/100‑11 (MQ=37)
                       aTTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGc                                                                              >  2:1874101/1‑100 (MQ=255)
                        ttATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCg                                                                             <  2:2424042/100‑1 (MQ=255)
                         tATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCgg                                                                            <  1:3249410/100‑1 (MQ=255)
                         tATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCgg                                                                            <  2:3799071/100‑1 (MQ=255)
                           tGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGaa                                                                          <  2:1782595/100‑1 (MQ=255)
                           tGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGaa                                                                          >  2:436541/1‑100 (MQ=255)
                            gCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGaaa                                                                         >  2:2835805/1‑100 (MQ=255)
                            gCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGaa                                                                          >  1:4386955/1‑99 (MQ=255)
                             cAGTCTGAACGTATTTGTTTGGTATGGGCCCATCCGCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGaaaa                                                                        <  1:2569511/100‑1 (MQ=255)
                                  tGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAg                                                                   >  2:1732032/1‑100 (MQ=255)
                                  tGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAg                                                                   <  1:2017603/100‑1 (MQ=255)
                                   gAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAgg                                                                  <  1:4003054/100‑1 (MQ=255)
                                   gAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAgg                                                                  >  2:817921/1‑100 (MQ=255)
                                   gAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAgg                                                                  <  2:2645767/100‑1 (MQ=255)
                                   gAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAg                                                                   <  1:904358/99‑1 (MQ=255)
                                      cGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTg                                                                >  1:3403317/1‑99 (MQ=255)
                                       gTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGAc                                                              >  1:1293655/1‑100 (MQ=255)
                                          tttatttGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGa                                                           <  1:3502591/100‑1 (MQ=255)
                                          tttatttGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGa                                                           <  1:4162414/100‑1 (MQ=255)
                                          tttatttGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGa                                                           >  2:1314386/1‑100 (MQ=255)
                                            tatttGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAAc                                                         >  2:4238740/1‑100 (MQ=255)
                                            tatttGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAAc                                                         <  1:436539/100‑1 (MQ=255)
                                              tttGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTc                                                       >  2:2937072/1‑100 (MQ=255)
                                              tttGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTc                                                       <  1:328919/100‑1 (MQ=255)
                                                tGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGa                                                     >  2:1419340/1‑100 (MQ=255)
                                                 ggTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGAAATGGAGCGGAAAATACAGGTGACGGAACTCGAt                                                    >  1:639435/1‑100 (MQ=255)
                                                  gTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATCTTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGAtt                                                   >  1:3155189/1‑100 (MQ=255)
                                                   tATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGAttt                                                  >  1:1018075/1‑100 (MQ=255)
                                                   tATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGAttt                                                  >  2:3822128/1‑100 (MQ=255)
                                                    aTGGTCCCCTCCTCGTCATGATTCATTGACCGCAAATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGAttt                                                  <  1:2014297/99‑1 (MQ=255)
                                                    aTGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGAttt                                                  <  2:1076001/99‑1 (MQ=255)
                                                       ggCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTatat                                              >  1:2139814/1‑100 (MQ=255)
                                                       ggCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTatat                                              >  1:3445214/1‑100 (MQ=255)
                                                        gCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATc                                             >  2:2020385/1‑100 (MQ=255)
                                                            aTCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCg                                         <  2:1371622/100‑1 (MQ=255)
                                                                tCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAAt                                     <  2:4217064/100‑1 (MQ=255)
                                                                tCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGAAATGGAGCGGAAAATACAGGTTACGGAACTCGATTTATAgcggcg                                         <  1:138558/96‑1 (MQ=255)
                                                                 cGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAAtt                                    >  1:512845/1‑100 (MQ=255)
                                                                 cGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAAt                                     >  1:2291445/1‑99 (MQ=255)
                                                                  cacATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGTAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAAttt                                   <  1:2442274/98‑1 (MQ=255)
                                                                   tCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGCTTTATATCGGCGTAATTTc                                  <  1:4213709/100‑1 (MQ=255)
                                                                   tCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTc                                  <  2:4402952/100‑1 (MQ=255)
                                                                      tGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAAc                               >  1:3063288/1‑100 (MQ=255)
                                                                      tGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAAc                               <  1:706061/100‑1 (MQ=255)
                                                                       gATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAAcc                              >  1:355188/1‑100 (MQ=255)
                                                                        aTTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAAccc                             <  2:482504/100‑1 (MQ=255)
                                                                         ttCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCa                            >  2:2690721/1‑100 (MQ=255)
                                                                           cATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAg                           >  1:2631438/1‑99 (MQ=255)
                                                                              tGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGtgat                       >  2:2923591/1‑100 (MQ=255)
                                                                                aCCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGtgatta                     <  2:2987250/100‑3 (MQ=255)
                                                                                   gCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGc                  >  2:1449515/1‑100 (MQ=255)
                                                                                   gCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGc                  >  1:1727301/1‑100 (MQ=255)
                                                                                   ccattgcccgCTGATGGATCGCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGc                  >  2:2261658/10‑100 (MQ=38)
                                                                                    cacaTATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGcc                 >  1:1120603/1‑100 (MQ=255)
                                                                                    cacaTATTGCTCATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACg                   <  2:3610661/98‑1 (MQ=255)
                                                                                      caTATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCgg               >  1:3465025/1‑100 (MQ=255)
                                                                                      caTATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCgg               <  1:260915/100‑1 (MQ=255)
                                                                                        tatTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGa              >  1:3583298/1‑99 (MQ=255)
                                                                                        gcccgCTGATGGATCGCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGaa             >  1:797032/5‑100 (MQ=255)
                                                                                         cccgCTGATGGATCGCATCAGCGGGCNATGGNGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGNCGGAAg            >  2:3381702/4‑100 (MQ=255)
                                                                                         cccgCTGATGGATCGCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAg            >  1:4317150/4‑100 (MQ=255)
                                                                                         cccgCTGATGGATCGCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAg            >  2:1262751/4‑100 (MQ=255)
                                                                                         cccgCTGATGGATCGCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAg            >  2:6532/4‑100 (MQ=255)
                                                                                         cccgCTGATGGATCGCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAg            >  2:1385744/4‑100 (MQ=255)
                                                                                         atTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAg            >  2:2500223/1‑100 (MQ=255)
                                                                                          ccgCTGATGGATCGCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGa           >  2:1776784/3‑100 (MQ=255)
                                                                                           tGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGAATTATATCGGCGTAATTTCAAAACAGTGATGACGACGCAAGAt          >  1:3097955/1‑100 (MQ=255)
                                                                                              tGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGa        <  1:1714255/99‑1 (MQ=255)
                                                                                                   cGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGa  >  1:135558/1‑100 (MQ=255)
                                                                                                   cGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGAATTGTATAGGCGTAATTTCAACACAGGGATGACGCTGGAAGATGTACCa    >  1:835753/1‑98 (MQ=255)
                                                                                                   |                                                                                                    
CCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGA  >  NC_000913/949642‑949840

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