Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 949,745:1 +G 16.7% 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.G16.7% 392.2 / 87.3 126coding (78/591 nt)ycaKputative NAD(P)H‑dependent oxidoreductase
Reads supporting (aligned to +/- strand):  ref base . (56/49);  new base G (10/11);  total (66/60)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.58e-01

Only 100 of 129 total aligned reads displayed.

ACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGG  >  NC_000913/949647‑949844
                                                                                                   |                                                                                                   
aCTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATc‑c                                                                                                    >  2:2524583/1‑100 (MQ=255)
 cTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATc‑c                                                                                                    <  2:1649620/99‑1 (MQ=255)
   tAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATc                                                                                                 <  2:1176270/100‑1 (MQ=255)
    aaCGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCa                                                                                                <  2:2270093/100‑1 (MQ=255)
      cGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATG‑GATCGCATCAGc                                                                                              <  2:2741809/100‑1 (MQ=255)
        tGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATG‑GATCGCATCAGcaa                                                                                            <  2:2888570/100‑3 (MQ=255)
        tGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATG‑GATCGCATCAGcaa                                                                                            <  2:249774/100‑3 (MQ=255)
        tGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATG‑GATCGCATCAGcaa                                                                                            <  2:1018510/100‑3 (MQ=255)
        tGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGGCATGATTCATTGACCGCACATATTGCTGATG‑GATCGCATCAGcaa                                                                                            <  1:785513/100‑3 (MQ=255)
         ggtggTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCggg                                                                                           <  1:976031/97‑1 (MQ=255)
         ggAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATG‑GATCGCATCAGcaat                                                                                           <  2:590954/100‑4 (MQ=255)
          gAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATG‑GATCGCATCAGcaa                                                                                            <  1:785434/98‑3 (MQ=255)
           aGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATG‑GATCGCATCAGcaatat                                                                                         <  2:1478261/100‑6 (MQ=255)
            ggTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCaa                                                                                        >  2:1755202/1‑100 (MQ=255)
            ggTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCaa                                                                                        >  1:1511391/1‑100 (MQ=255)
             gTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAAt                                                                                       >  1:104619/1‑100 (MQ=255)
              tAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATg                                                                                      >  1:2680648/1‑100 (MQ=255)
               aaaaTTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATgg                                                                                     <  2:2516058/100‑1 (MQ=255)
                  aTTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGc                                                                                  <  2:1661651/100‑1 (MQ=255)
                   ttATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCg                                                                                 <  2:236321/100‑1 (MQ=255)
                   ttATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCg                                                                                 <  1:1111708/100‑1 (MQ=255)
                   ttATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCg                                                                                 <  1:491517/100‑1 (MQ=255)
                   ttATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCg                                                                                 <  1:297055/100‑1 (MQ=255)
                     aTGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCgg                                                                                >  1:1018672/1‑99 (MQ=255)
                      tGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGaa                                                                              >  1:1569020/1‑100 (MQ=255)
                        cAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGaaaa                                                                            <  1:498986/100‑1 (MQ=255)
                         aGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAAt                                                                           <  2:1484219/100‑1 (MQ=255)
                         aGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATGGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAAt                                                                           <  1:192647/100‑1 (MQ=255)
                           tCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATAc                                                                         >  2:2031923/1‑100 (MQ=255)
                            cTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACa                                                                        <  2:1387610/100‑1 (MQ=255)
                              gAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAgg                                                                      <  2:1350230/100‑1 (MQ=255)
                              gAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAgg                                                                      <  2:2549773/100‑1 (MQ=255)
                              gAACGTATTTATTTGGTATGGACCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAgg                                                                      <  2:1625565/100‑1 (MQ=255)
                                   tatttattTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACg                                                                 >  2:1851579/1‑100 (MQ=255)
                                     tggatttGGTATGGGCCTATCGTCGTGGTTATTGATTTGCCGCACAAATTGCTGATGCGTTC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGa                                                               <  1:1631039/97‑1 (MQ=255)
                                       tatttGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAAc                                                             <  2:1028803/100‑1 (MQ=255)
                                       tatttGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAAc                                                             <  1:2697049/100‑1 (MQ=255)
                                       tatttGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAAc                                                             <  1:1407411/100‑1 (MQ=255)
                                        atttGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACt                                                            <  2:2726678/100‑1 (MQ=255)
                                         tttGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTc                                                           >  2:257171/1‑100 (MQ=255)
                                           tGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACCGAACTCGa                                                         >  1:1426563/1‑100 (MQ=255)
                                              tATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGAttt                                                      <  2:1081660/100‑1 (MQ=255)
                                               aTGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTta                                                     <  2:2822571/100‑1 (MQ=255)
                                               aTGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTta                                                     <  2:486887/100‑1 (MQ=255)
                                                 gggCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTtata                                                   >  1:1163373/1‑100 (MQ=255)
                                                  ggCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTatat                                                  >  1:1258677/1‑100 (MQ=255)
                                                  ggCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTatat                                                  >  1:21369/1‑100 (MQ=255)
                                                  ggCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTatat                                                  >  1:2885456/1‑100 (MQ=255)
                                                   gCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATc                                                 >  2:1683201/1‑100 (MQ=255)
                                                     ccATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCgg                                               <  2:2496807/100‑1 (MQ=255)
                                                     ccATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCgg                                               <  2:1282718/100‑1 (MQ=255)
                                                      cATCCTCGACATGATTCATTGACGGCACATATTGCTGCTGGTATC‑GATCTGCGGGCTATGACGCGGAAAATACAGGTGCTGGACGTCGGTTtata                                                   >  1:1723818/1‑95 (MQ=255)
                                                       aTCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCg                                             <  2:1890437/100‑1 (MQ=255)
                                                       aTCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCg                                             <  2:2221990/100‑1 (MQ=255)
                                                         ccTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTa                                           >  1:2599296/1‑100 (MQ=255)
                                                            cGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAAtt                                        <  1:2827327/100‑1 (MQ=255)
                                                               cATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCa                                     >  2:2164093/1‑100 (MQ=255)
                                                                 tgtttCATTGACCGCAACTATTGCTTACCCTTTC‑CCTCAGCCAGCCATGTAGTTGAAATTACACCTGACGCAACTCGATTTATATCGGCGTAATGTCaa                                    <  1:642514/96‑1 (MQ=255)
                                                                  gATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTaaaac                                  >  1:1423097/1‑98 (MQ=255)
                                                                  gATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAAcc                                  >  2:2475730/1‑100 (MQ=255)
                                                                   aTTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAAccc                                 <  2:1251538/100‑1 (MQ=255)
                                                                      cATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGt                              <  2:700787/100‑1 (MQ=255)
                                                                       aTTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCCATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGtg                             >  1:2659052/1‑100 (MQ=255)
                                                                       aTTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGtg                             >  1:223660/1‑100 (MQ=255)
                                                                       aTTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGtg                             >  2:2297724/1‑100 (MQ=255)
                                                                       aTTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGt                              >  1:602604/1‑99 (MQ=255)
                                                                        ttGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGtga                            <  2:1807317/100‑1 (MQ=255)
                                                                          gACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGtgatg                          >  2:969903/1‑100 (MQ=255)
                                                                          gACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGtgatg                          >  1:998507/1‑100 (MQ=255)
                                                                          gACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGtgatg                          <  1:2105671/100‑1 (MQ=255)
                                                                           aCCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGtgatga                         >  2:1936673/1‑100 (MQ=255)
                                                                           aCCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGtgatga                         >  2:1349009/1‑100 (MQ=255)
                                                                             cGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACg                       >  1:530629/1‑100 (MQ=255)
                                                                             cGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAAATCGATTTATATCGGCGTAATTTCAACCCAGtgatgaag                       >  1:1655896/1‑98 (MQ=255)
                                                                               cacaTATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGcc                     >  1:282305/1‑100 (MQ=255)
                                                                                  atatTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGa                  <  1:1619011/100‑1 (MQ=255)
                                                                                   tatTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGaa                 <  1:2228152/100‑1 (MQ=255)
                                                                                   tatTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGaa                 <  1:47580/100‑1 (MQ=255)
                                                                                   gcccgCTGATG‑GATCGCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGaa                 >  2:1535266/5‑100 (MQ=255)
                                                                                    cccgCTGATG‑GATCGCATCAGCGGGCAATGGAGCGGCAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAg                >  2:2342296/4‑100 (MQ=255)
                                                                                    cccgCTGATG‑GATCGCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAg                >  1:2370841/4‑100 (MQ=255)
                                                                                    cccgCTGATG‑GATCGCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAg                >  1:2424998/4‑100 (MQ=255)
                                                                                    cccgCTGATG‑GATCGCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAg                >  2:529436/4‑100 (MQ=255)
                                                                                     ttGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTGATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAg                <  1:1075848/99‑1 (MQ=255)
                                                                                     ccgCTGATG‑GATCGCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGa               >  2:643864/3‑100 (MQ=255)
                                                                                     ccgCTGATG‑GATCGCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGa               >  2:2165407/3‑100 (MQ=255)
                                                                                      cgCTGATG‑GATCGCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGAt              >  1:1488450/2‑100 (MQ=255)
                                                                                       gCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATg             <  2:842531/100‑1 (MQ=255)
                                                                                       gCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATg             <  1:1867297/100‑1 (MQ=255)
                                                                                        cTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGa            >  1:1520497/1‑100 (MQ=255)
                                                                                        cTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGa            >  2:509913/1‑100 (MQ=255)
                                                                                         tGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGaa           >  1:78216/1‑100 (MQ=255)
                                                                                          gATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACACAGTGATGACGCCGGAGGATGAAc          >  1:1003651/1‑100 (MQ=255)
                                                                                           aTGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAAcc         <  1:2855512/100‑1 (MQ=255)
                                                                                           aTGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAAcc         <  2:1714685/100‑1 (MQ=255)
                                                                                              cGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGa      >  1:2677607/1‑100 (MQ=255)
                                                                                              cGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATGTATATCGGCGTAATTTCAACCCAGTGATGACGCCCGAAGATGAACCAGa      <  1:201700/100‑1 (MQ=255)
                                                                                               gATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGAc     >  1:2681550/1‑100 (MQ=255)
                                                                                               gATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGAc     >  2:833739/1‑100 (MQ=255)
                                                                                               gATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGAc     >  1:581615/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
ACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGG  >  NC_000913/949647‑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: ‑