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.

AACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGG  >  NC_000913/949646‑949844
                                                                                                    |                                                                                                   
 aCTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATc‑c                                                                                                    <  1:2538674/100‑1 (MQ=255)
 aCTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATc‑c                                                                                                    <  1:622470/100‑1 (MQ=255)
  cTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑Ca                                                                                                   <  2:666100/100‑1 (MQ=255)
    tAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATc                                                                                                 >  2:3249427/1‑100 (MQ=255)
    tAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATc                                                                                                 >  1:2944610/1‑100 (MQ=255)
    tAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATc                                                                                                 <  2:2161763/100‑1 (MQ=255)
       cGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGc                                                                                              >  2:4347574/1‑100 (MQ=255)
        gTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATTCGATC‑CATCAGCg                                                                                             <  2:2960657/100‑1 (MQ=255)
         tGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCgg                                                                                            <  2:152477/100‑1 (MQ=255)
         tGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATG‑GATCGCATCAGcaa                                                                                            <  2:172457/100‑3 (MQ=255)
          ggAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATG‑GATCGCATCAGcaat                                                                                           <  2:2418098/100‑4 (MQ=255)
          ggAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATG‑GATCGCATCAGcaat                                                                                           <  2:220289/100‑4 (MQ=255)
          ggAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATG‑GATCGCATCAGcaat                                                                                           <  1:1551302/100‑4 (MQ=255)
          ggAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATG‑GATCGCATCAGcaat                                                                                           <  2:1251708/100‑4 (MQ=255)
          ggAGGTAAAATTATGCAGTCTGAACGTATTTAGTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATG‑GATCGCATCAGcaa                                                                                            <  1:435218/99‑3 (MQ=255)
           gAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGc                                                                                          <  1:164874/100‑1 (MQ=255)
           gAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATG‑GATCGCATCAGcaata                                                                                          <  2:84980/100‑5 (MQ=255)
           gAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCCGATGCGATC‑CATCAGCGGGc                                                                                          <  1:2698114/100‑1 (MQ=255)
             ggTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCaa                                                                                        >  1:1021626/1‑100 (MQ=255)
              gTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAAt                                                                                       <  2:303947/100‑1 (MQ=255)
                aaaaTTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATgg                                                                                     <  1:30681/100‑1 (MQ=255)
                aaaaTTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATA‑CATCAGCGGGCAATgg                                                                                     >  1:2577616/1‑100 (MQ=255)
                aaaaTTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATCGCTGATG‑GATCGCATCAGcaatatgtgc                                                                                     <  2:2102026/100‑10 (MQ=37)
                   aTTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGc                                                                                  >  2:1874101/1‑100 (MQ=255)
                    ttATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCg                                                                                 <  2:2424042/100‑1 (MQ=255)
                     tATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCgg                                                                                <  2:3799071/100‑1 (MQ=255)
                     tATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCgg                                                                                <  1:3249410/100‑1 (MQ=255)
                       tGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGaa                                                                              <  2:1782595/100‑1 (MQ=255)
                       tGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGaa                                                                              >  2:2857471/1‑100 (MQ=255)
                        gCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGaaa                                                                             >  2:2835805/1‑100 (MQ=255)
                        gCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGaa                                                                              >  1:4386955/1‑99 (MQ=255)
                         cAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGaaaa                                                                            <  2:437248/100‑1 (MQ=255)
                            tCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATAc                                                                         <  1:3086750/100‑1 (MQ=255)
                              tGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAg                                                                       <  1:2017603/100‑1 (MQ=255)
                              tGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAg                                                                       >  2:1732032/1‑100 (MQ=255)
                               gAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAgg                                                                      <  1:4003054/100‑1 (MQ=255)
                                  cGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTg                                                                    >  1:3403317/1‑99 (MQ=255)
                                   gTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGAc                                                                  >  1:1293655/1‑100 (MQ=255)
                                      tttatttGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGa                                                               <  1:4162414/100‑1 (MQ=255)
                                      tttatttGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGa                                                               >  2:1314386/1‑100 (MQ=255)
                                      tttatttGGTATGGGCCCATCCTCGTCATGATTCAATTACCGCACAGATTGCTGATGCGATC‑CACCAGCGGGCAATGCAGCGGAAAGTACAGGTGACGGa                                                               <  1:1214296/100‑1 (MQ=255)
                                        tatttGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAAc                                                             <  1:436539/100‑1 (MQ=255)
                                          tttGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTc                                                           >  2:2937072/1‑100 (MQ=255)
                                            tGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGa                                                         >  2:1419340/1‑100 (MQ=255)
                                             ggTATGGTCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTc                                                           <  1:3463257/97‑1 (MQ=255)
                                             ggTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGAt                                                        <  1:3018055/100‑1 (MQ=255)
                                             ggTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGAAATGGAGCGGAAAATACAGGTGACGGAACTCGAt                                                        >  1:639435/1‑100 (MQ=255)
                                               tATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGAttt                                                      >  2:3822128/1‑100 (MQ=255)
                                               tATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGAttt                                                      >  1:1018075/1‑100 (MQ=255)
                                                aTGGTCCCCTCCTCGTCATGATTCATTGACCGCAAATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGAttt                                                      <  1:2014297/99‑1 (MQ=255)
                                                aTGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGAttt                                                      <  2:1076001/99‑1 (MQ=255)
                                                aTGGGCCCATCCTCGTCATGATTCATTGACCGAACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTta                                                     <  2:648114/100‑1 (MQ=255)
                                                   ggCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTatat                                                  >  1:2139814/1‑100 (MQ=255)
                                                    gCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATc                                                 >  2:2020385/1‑100 (MQ=255)
                                                        aTCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCg                                             <  2:1371622/100‑1 (MQ=255)
                                                            tCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAAt                                         <  2:4217064/100‑1 (MQ=255)
                                                            tCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGAAATGGAGCGGAAAATACAGGTTACGGAACTCGATTTATAgcggcg                                             <  1:138558/96‑1 (MQ=255)
                                                             cGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAAtt                                        >  1:512845/1‑100 (MQ=255)
                                                             cGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAAt                                         >  1:2291445/1‑99 (MQ=255)
                                                               tCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGCTTTATATCGGCGTAATTTc                                      <  1:4213709/100‑1 (MQ=255)
                                                               tCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTc                                      <  2:4402952/100‑1 (MQ=255)
                                                                   gATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAAcc                                  >  1:355188/1‑100 (MQ=255)
                                                                    aTTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAAccc                                 <  1:2592957/100‑1 (MQ=255)
                                                                     ttCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCa                                >  2:2690721/1‑100 (MQ=255)
                                                                       cATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAg                               >  1:2631438/1‑99 (MQ=255)
                                                                       cATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGt                              >  1:762517/1‑100 (MQ=255)
                                                                        aTTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGtg                             >  1:2134870/1‑100 (MQ=255)
                                                                             ccGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGAc                        >  2:4089808/1‑100 (MQ=255)
                                                                              cGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACg                       <  2:267604/100‑1 (MQ=255)
                                                                               gCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGc                      >  1:1727301/1‑100 (MQ=255)
                                                                               gCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGc                      >  2:1449515/1‑100 (MQ=255)
                                                                               ccattgcccgCTGATG‑GATCGCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGc                      >  2:2261658/10‑100 (MQ=38)
                                                                                cacaTATTGCTCATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACg                       <  2:3610661/98‑1 (MQ=255)
                                                                                 attgcccgCTGATG‑GATCGCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCg                    >  2:4318481/8‑100 (MQ=38)
                                                                                  caTATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCgg                   >  1:3465025/1‑100 (MQ=255)
                                                                                   atatTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGa                  <  2:684456/100‑1 (MQ=255)
                                                                                    gcccgCTGATG‑GATCGCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGaa                 >  2:1177074/5‑100 (MQ=255)
                                                                                    gcccgCTGATG‑GATCGCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGaa                 >  1:797032/5‑100 (MQ=255)
                                                                                     cccgCTGATG‑GATCGCATCAGCGGGCNATGGNGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGNCGGAAg                >  2:3381702/4‑100 (MQ=255)
                                                                                     cccgCTGATG‑GATCGCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAg                >  2:1385744/4‑100 (MQ=255)
                                                                                     cccgCTGATG‑GATCGCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAg                >  1:3347674/4‑100 (MQ=255)
                                                                                     cccgCTGATG‑GATCGCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAg                >  1:4317150/4‑100 (MQ=255)
                                                                                     atTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAg                >  2:2500223/1‑100 (MQ=255)
                                                                                      ttGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGaa                 >  1:942093/1‑98 (MQ=255)
                                                                                      ccgCTGATG‑GATCGCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGa               >  2:1776784/3‑100 (MQ=255)
                                                                                       tGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGAt              <  2:3112959/100‑1 (MQ=255)
                                                                                       tGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGAATTATATCGGCGTAATTTCAAAACAGTGATGACGACGCAAGAt              >  1:3097955/1‑100 (MQ=255)
                                                                                       cgCTGATG‑GATCGCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGAt              >  1:3533108/2‑100 (MQ=255)
                                                                                         cTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGa            >  1:3198193/1‑100 (MQ=255)
                                                                                          tGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGa            <  1:1714255/99‑1 (MQ=255)
                                                                                               cGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGa      >  1:2249584/1‑100 (MQ=255)
                                                                                               cGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGa      >  1:135558/1‑100 (MQ=255)
                                                                                               cGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGAATTGTATAGGCGTAATTTCAACACAGGGATGACGCTGGAAGATGTACCa        >  1:835753/1‑98 (MQ=255)
                                                                                                gATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGAc     >  1:1144947/1‑100 (MQ=255)
                                                                                                gATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGAc     <  2:512847/100‑1 (MQ=255)
                                                                                                 aTC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGa      >  1:2520708/1‑98 (MQ=255)
                                                                                                 aTC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACt    >  1:3150237/1‑100 (MQ=255)
                                                                                                 aTC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACt    <  2:2542357/100‑1 (MQ=255)
                                                                                                  tC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTg   <  1:1333796/100‑1 (MQ=255)
                                                                                                   c‑cATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTgg  >  2:3863792/1‑100 (MQ=255)
                                                                                                    |                                                                                                   
AACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGG  >  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: ‑