Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 949,745:1 +G 18.5% coding (78/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,7451.G18.5% 427.7 / 119.2 146coding (78/591 nt)ycaKputative NAD(P)H‑dependent oxidoreductase
Reads supporting (aligned to +/- strand):  ref base . (59/60);  new base G (14/13);  total (73/73)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.75e-01

Only 100 of 151 total aligned reads displayed.

AACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGG  >  NC_000913/949646‑949844
                                                                                                    |                                                                                                   
 aCTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATcc                                                                                                    <  1:2538674/100‑1 (MQ=255)
 aCTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATcc                                                                                                    <  1:622470/100‑1 (MQ=255)
  cTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCa                                                                                                   <  2:666100/100‑1 (MQ=255)
    tAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATc                                                                                                 >  2:3249427/1‑100 (MQ=255)
    tAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATc                                                                                                 >  1:2944610/1‑100 (MQ=255)
    tAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATc                                                                                                 <  2:2161763/100‑1 (MQ=255)
       cGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGc                                                                                              >  2:4347574/1‑100 (MQ=255)
        gTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATTCGATCCATCAGCg                                                                                             <  2:2960657/100‑1 (MQ=255)
         tGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCgg                                                                                            <  2:152477/100‑1 (MQ=255)
         tGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGGATCGCATCAGcaa                                                                                            <  2:172457/100‑3 (MQ=255)
          ggAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGGATCGCATCAGcaat                                                                                           <  2:2418098/100‑4 (MQ=255)
          ggAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGGATCGCATCAGcaat                                                                                           <  2:220289/100‑4 (MQ=255)
          ggAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGGATCGCATCAGcaat                                                                                           <  1:1551302/100‑4 (MQ=255)
          ggAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGGATCGCATCAGcaat                                                                                           <  2:1251708/100‑4 (MQ=255)
          ggAGGTAAAATTATGCAGTCTGAACGTATTTAGTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGGATCGCATCAGcaa                                                                                            <  1:435218/99‑3 (MQ=255)
           gAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGc                                                                                          <  1:164874/100‑1 (MQ=255)
           gAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGGATCGCATCAGcaata                                                                                          <  2:84980/100‑5 (MQ=255)
           gAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCCGATGCGATCCATCAGCGGGc                                                                                          <  1:2698114/100‑1 (MQ=255)
             ggTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCaa                                                                                        >  1:1021626/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)
                   aTTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGc                                                                                  >  2:1874101/1‑100 (MQ=255)
                    ttATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCg                                                                                 <  2:2424042/100‑1 (MQ=255)
                     tATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCgg                                                                                <  2:3799071/100‑1 (MQ=255)
                     tATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCgg                                                                                <  1:3249410/100‑1 (MQ=255)
                       tGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGaa                                                                              <  2:1782595/100‑1 (MQ=255)
                       tGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGaa                                                                              >  2:2857471/1‑100 (MQ=255)
                        gCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGaaa                                                                             >  2:2835805/1‑100 (MQ=255)
                        gCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGaa                                                                              >  1:4386955/1‑99 (MQ=255)
                         cAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGaaaa                                                                            <  2:437248/100‑1 (MQ=255)
                            tCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATAc                                                                         <  1:3086750/100‑1 (MQ=255)
                              tGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAg                                                                       <  1:2017603/100‑1 (MQ=255)
                              tGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAg                                                                       >  2:1732032/1‑100 (MQ=255)
                               gAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAgg                                                                      <  1:4003054/100‑1 (MQ=255)
                                  cGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTg                                                                    >  1:3403317/1‑99 (MQ=255)
                                   gTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGAc                                                                  >  1:1293655/1‑100 (MQ=255)
                                      tttatttGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGa                                                               <  1:4162414/100‑1 (MQ=255)
                                      tttatttGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGa                                                               >  2:1314386/1‑100 (MQ=255)
                                      tttatttGGTATGGGCCCATCCTCGTCATGATTCAATTACCGCACAGATTGCTGATGCGATCCACCAGCGGGCAATGCAGCGGAAAGTACAGGTGACGGa                                                               <  1:1214296/100‑1 (MQ=255)
                                        tatttGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAAc                                                             <  1:436539/100‑1 (MQ=255)
                                          tttGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTc                                                           >  2:2937072/1‑100 (MQ=255)
                                            tGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGa                                                         >  2:1419340/1‑100 (MQ=255)
                                             ggTATGGTCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTc                                                           <  1:3463257/97‑1 (MQ=255)
                                             ggTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGAt                                                        <  1:3018055/100‑1 (MQ=255)
                                             ggTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGAAATGGAGCGGAAAATACAGGTGACGGAACTCGAt                                                        >  1:639435/1‑100 (MQ=255)
                                               tATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGAttt                                                      >  2:3822128/1‑100 (MQ=255)
                                               tATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGAttt                                                      >  1:1018075/1‑100 (MQ=255)
                                                aTGGTCCCCTCCTCGTCATGATTCATTGACCGCAAATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGAttt                                                      <  1:2014297/99‑1 (MQ=255)
                                                aTGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGAttt                                                      <  2:1076001/99‑1 (MQ=255)
                                                aTGGGCCCATCCTCGTCATGATTCATTGACCGAACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTta                                                     <  2:648114/100‑1 (MQ=255)
                                                   ggCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTatat                                                  >  1:2139814/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)
                                                               tCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGCTTTATATCGGCGTAATTTc                                      <  1:4213709/100‑1 (MQ=255)
                                                               tCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTc                                      <  2:4402952/100‑1 (MQ=255)
                                                                   gATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAAcc                                  >  1:355188/1‑100 (MQ=255)
                                                                    aTTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAAccc                                 <  1:2592957/100‑1 (MQ=255)
                                                                     ttCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCa                                >  2:2690721/1‑100 (MQ=255)
                                                                       cATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAg                               >  1:2631438/1‑99 (MQ=255)
                                                                       cATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGt                              >  1:762517/1‑100 (MQ=255)
                                                                        aTTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGtg                             >  1:2134870/1‑100 (MQ=255)
                                                                             ccGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGAc                        >  2:4089808/1‑100 (MQ=255)
                                                                              cGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACg                       <  2:267604/100‑1 (MQ=255)
                                                                               gCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGc                      >  1:1727301/1‑100 (MQ=255)
                                                                               gCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGc                      >  2:1449515/1‑100 (MQ=255)
                                                                               ccattgcccgCTGATGGATCGCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGc                      >  2:2261658/10‑100 (MQ=38)
                                                                                cacaTATTGCTCATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACg                       <  2:3610661/98‑1 (MQ=255)
                                                                                 attgcccgCTGATGGATCGCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCg                    >  2:4318481/8‑100 (MQ=38)
                                                                                  caTATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCgg                   >  1:3465025/1‑100 (MQ=255)
                                                                                   atatTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGa                  <  2:684456/100‑1 (MQ=255)
                                                                                    gcccgCTGATGGATCGCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGaa                 >  2:1177074/5‑100 (MQ=255)
                                                                                    gcccgCTGATGGATCGCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGaa                 >  1:797032/5‑100 (MQ=255)
                                                                                     cccgCTGATGGATCGCATCAGCGGGCNATGGNGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGNCGGAAg                >  2:3381702/4‑100 (MQ=255)
                                                                                     cccgCTGATGGATCGCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAg                >  2:1385744/4‑100 (MQ=255)
                                                                                     cccgCTGATGGATCGCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAg                >  1:3347674/4‑100 (MQ=255)
                                                                                     cccgCTGATGGATCGCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAg                >  1:4317150/4‑100 (MQ=255)
                                                                                     atTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAg                >  2:2500223/1‑100 (MQ=255)
                                                                                      ttGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGaa                 >  1:942093/1‑98 (MQ=255)
                                                                                      ccgCTGATGGATCGCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGa               >  2:1776784/3‑100 (MQ=255)
                                                                                       tGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGAt              <  2:3112959/100‑1 (MQ=255)
                                                                                       tGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGAATTATATCGGCGTAATTTCAAAACAGTGATGACGACGCAAGAt              >  1:3097955/1‑100 (MQ=255)
                                                                                       cgCTGATGGATCGCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGAt              >  1:3533108/2‑100 (MQ=255)
                                                                                         cTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGa            >  1:3198193/1‑100 (MQ=255)
                                                                                          tGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGa            <  1:1714255/99‑1 (MQ=255)
                                                                                               cGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGa      >  1:2249584/1‑100 (MQ=255)
                                                                                               cGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGa      >  1:135558/1‑100 (MQ=255)
                                                                                               cGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGAATTGTATAGGCGTAATTTCAACACAGGGATGACGCTGGAAGATGTACCa        >  1:835753/1‑98 (MQ=255)
                                                                                                gATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGAc     >  1:1144947/1‑100 (MQ=255)
                                                                                                gATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGAc     <  2:512847/100‑1 (MQ=255)
                                                                                                 aTCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGa      >  1:2520708/1‑98 (MQ=255)
                                                                                                 aTCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACt    >  1:3150237/1‑100 (MQ=255)
                                                                                                 aTCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACt    <  2:2542357/100‑1 (MQ=255)
                                                                                                  tCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTg   <  1:1333796/100‑1 (MQ=255)
                                                                                                   ccATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTgg  >  2:3863792/1‑100 (MQ=255)
                                                                                                    |                                                                                                   
AACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGG  >  NC_000913/949646‑949844

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