Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 949,741 Δ1 bp 15.1% 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.15.1% 270.2 / 46.9 93coding (74/591 nt)ycaKputative NAD(P)H‑dependent oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (40/38);  new base . (9/5);  total (50/43)
Fisher's exact test for biased strand distribution p-value = 4.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

CCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGA  >  NC_000913/949642‑949840
                                                                                                   |                                                                                                    
ncAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGc                                                                                                      >  2:1300114/2‑100 (MQ=255)
ccAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGc                                                                                                      <  2:839470/100‑1 (MQ=255)
ccAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGc                                                                                                      >  2:418825/1‑100 (MQ=255)
 cAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGc                                                                                                      <  2:2174410/99‑1 (MQ=255)
  aGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGa                                                                                                    >  1:114170/1‑100 (MQ=255)
   gAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGAt                                                                                                   <  1:2381838/100‑1 (MQ=255)
    aaCTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATc                                                                                                  <  1:2156766/100‑1 (MQ=255)
    aaCTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATc                                                                                                  <  2:1865951/100‑1 (MQ=255)
      cTTAACGTGGAGGTAAAATTATGCAGTCTGATCGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑Ca                                                                                               <  2:1179275/100‑1 (MQ=255)
      cTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑Ca                                                                                               <  2:1962356/100‑1 (MQ=255)
      cTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑Ca                                                                                               >  2:1234556/1‑100 (MQ=255)
        tAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATc                                                                                             <  1:1293156/100‑1 (MQ=255)
          aCGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAg                                                                                           >  1:56347/1‑100 (MQ=255)
           cGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGc                                                                                          >  1:403934/1‑100 (MQ=255)
            gTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATG‑GATCGCATCAGca                                                                                         <  2:495415/100‑2 (MQ=255)
             tGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATG‑GATCGCATCAGcaa                                                                                        <  1:2559255/100‑3 (MQ=255)
             tGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATG‑GATCGCATCAGcaa                                                                                        <  1:2335227/100‑3 (MQ=255)
               gAGGTAAAATTATGCAGTCTGAACGTATTTATGTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGc                                                                                      <  2:225902/100‑1 (MQ=255)
               gAGGTAAAATTATGCAGTCGGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACAGATTGCTGATGCGATC‑CATCAGCGGGc                                                                                      <  2:425246/100‑1 (MQ=255)
                aGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCGTCCTCGTCATGATTCATTGACCGCACATATTGCTGATG‑GATCGCATCAGcaa                                                                                        <  2:1550105/97‑3 (MQ=255)
                aGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCa                                                                                     >  2:1745242/1‑100 (MQ=255)
                 ggTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCaa                                                                                    >  2:2235831/1‑100 (MQ=255)
                   tAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCACTg                                                                                  >  2:2249244/1‑100 (MQ=255)
                   tAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATg                                                                                  >  1:1372474/1‑100 (MQ=255)
                    aaaaTTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATgg                                                                                 <  1:621829/100‑1 (MQ=255)
                    aaaaTTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATG‑GATCGCATCAGcaatatgtgc                                                                                 <  1:1992697/100‑10 (MQ=38)
                     aaaTTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGa                                                                                >  2:321298/1‑100 (MQ=255)
                      aaTTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAg                                                                               <  1:462973/100‑1 (MQ=255)
                          aTGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGGCAAGATTCATTGACCGCACATATTGCAGATGCGTAC‑CATCAGCGGGCAATGGAGCgg                                                                            >  2:1253826/1‑99 (MQ=255)
                             cAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGaaaa                                                                        >  2:393314/1‑100 (MQ=255)
                              agccTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATGGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAAt                                                                       <  2:420434/97‑1 (MQ=255)
                              aGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAAt                                                                       >  1:2115935/1‑100 (MQ=255)
                               gTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATa                                                                      <  1:1040811/100‑1 (MQ=255)
                                tCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATAc                                                                     >  2:2064102/1‑100 (MQ=255)
                                tCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATAc                                                                     >  1:2128449/1‑100 (MQ=255)
                                tCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATAc                                                                     >  1:1504570/1‑100 (MQ=255)
                                tCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATAc                                                                     <  1:418829/100‑1 (MQ=255)
                                 cgatctgTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACa                                                                    <  2:1043518/94‑1 (MQ=255)
                                   gAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAgg                                                                  <  2:216545/100‑1 (MQ=255)
                                      cGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CCTCAGCGGGAAATGGCACGGAAAATaaccgt                                                                 >  1:1108431/1‑93 (MQ=255)
                                       gTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACaagatcg                                                              >  1:1043526/1‑94 (MQ=255)
                                         atttatttGGTATGGGCCCATCCTCGTCATGATTCATTTACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCTGAAAATACAGGTGAc                                                              >  2:2246480/1‑98 (MQ=255)
                                         atttatttGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACgg                                                            >  1:1303892/1‑100 (MQ=255)
                                          tttatttGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGa                                                           >  2:1888975/1‑100 (MQ=255)
                                          tttatttGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGa                                                           >  1:1427814/1‑100 (MQ=255)
                                           ttatttGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGaa                                                          <  1:431607/100‑1 (MQ=255)
                                            tatttGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAAc                                                         >  2:998002/1‑100 (MQ=255)
                                                tGGTATGGGCCCATCCTCGTCATGATTCATTGACCTCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGa                                                     <  1:1910944/100‑1 (MQ=255)
                                                 ggTATGGGCCCATCCTCGTCATGATTCATTTACCGCACATATTGCTGATACGATC‑CATCAGCGCGCAATGAAGCGGAAAATACAGGAGACAGAACTc                                                       >  2:713104/1‑97 (MQ=255)
                                                   tATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGAttt                                                  <  1:1624623/100‑1 (MQ=255)
                                                   tATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGAttt                                                  <  1:1598824/100‑1 (MQ=255)
                                                   gatctgCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGAttt                                                  <  1:853004/95‑1 (MQ=255)
                                                    aTGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTta                                                 <  1:388638/100‑1 (MQ=255)
                                                       ggCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTatat                                              >  2:440791/1‑100 (MQ=255)
                                                        gCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTAGATc                                             >  2:852996/1‑100 (MQ=255)
                                                         cccATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCg                                            >  2:659624/1‑100 (MQ=255)
                                                           cATCCTCGTAATGATTCATTGACCGCACATATTGCTGATGAGATC‑CATCAGCGGGCAATGGAGCGGAAAATACATGTGACGGAACTCGAGTTATATCGGc                                          >  2:408301/1‑100 (MQ=255)
                                                            aTCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCg                                         >  2:840535/1‑100 (MQ=255)
                                                            aTCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCg                                         >  1:1978626/1‑100 (MQ=255)
                                                            aTCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCg                                         <  1:2151602/100‑1 (MQ=255)
                                                             tCCTCGTCATGATTCATTGACAGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCg                                         <  2:2300964/99‑1 (MQ=255)
                                                              ccTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTa                                       >  1:583405/1‑100 (MQ=255)
                                                                   tCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTc                                  <  2:1086093/100‑1 (MQ=255)
                                                                    cATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCa                                 >  2:84536/1‑100 (MQ=255)
                                                                       gATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAAcc                              <  1:2537823/100‑1 (MQ=255)
                                                                           cATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGt                          >  1:802117/1‑100 (MQ=255)
                                                                           cATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGt                          <  1:1734648/100‑1 (MQ=255)
                                                                           cATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATAGCGGCGTAATTTCAACCCAGt                          <  2:547114/100‑1 (MQ=255)
                                                                            aTTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGtg                         >  1:1147182/1‑100 (MQ=255)
                                                                            aTTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGAGTTATATCGGCGTAATTTCAACCCAGtg                         >  2:569180/1‑100 (MQ=255)
                                                                             ttGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGtga                        >  1:609004/1‑100 (MQ=255)
                                                                              tggccGCACATATTGCTGATGCGAGC‑CATCGGCGGGCAATGGAGCGGAAAATAAAGGTGACGGCACGCGATTTATAACGGCGCAATTTCAACCCCGtga                        >  2:2429941/4‑99 (MQ=255)
                                                                              tGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGtgat                       <  1:808973/100‑1 (MQ=255)
                                                                               gACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGtgatg                      >  1:1190506/1‑100 (MQ=255)
                                                                                aCCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGGACTCGATTTATATCGGCGTAATTTCAACCCAGtgatg                      <  2:1370490/99‑1 (MQ=255)
                                                                                  cGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGTCGTAATTTCAACCCAGTGATGACg                   <  2:929707/100‑1 (MQ=255)
                                                                                  cGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACg                   <  2:2459481/100‑1 (MQ=255)
                                                                                   gCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGc                  <  1:706346/100‑1 (MQ=255)
                                                                                   gCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGc                  >  1:292535/1‑100 (MQ=255)
                                                                                   gCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGc                  >  1:2268702/1‑100 (MQ=255)
                                                                                    cattgcccgCTGATG‑GATCGCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGcc                 >  1:1633611/9‑100 (MQ=38)
                                                                                      caTATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCgg               <  2:1840429/100‑1 (MQ=255)
                                                                                        tatTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGGCGGCGGa              >  1:1665437/1‑99 (MQ=255)
                                                                                        gcccgCTGATG‑GATCGCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGaa             >  1:2064094/5‑100 (MQ=255)
                                                                                        gcccgCTGATG‑GATCGCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGaa             >  1:1223303/5‑100 (MQ=255)
                                                                                         cccgCTGATG‑GATCGCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATGGGCGTAATTTCAACCCAGTGATGACGCCGGAAg            >  2:294242/4‑100 (MQ=255)
                                                                                         cccgCTGATG‑GATCGCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAg            >  1:697636/4‑100 (MQ=255)
                                                                                         cccgCTGATG‑GATCGCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAg            >  1:1507335/4‑100 (MQ=255)
                                                                                         cccgCTGATG‑GATCGCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAg            >  1:1295061/4‑100 (MQ=255)
                                                                                          ccgCTGATG‑GATCGCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGa           >  1:2027048/3‑100 (MQ=255)
                                                                                          ccgCTGATG‑GATCGCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGAATTATATCGGAGTAATTTCAACCCAGTGATGACGCAGGAAGa           >  2:2083797/3‑100 (MQ=255)
                                                                                             cTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGa        <  2:1020621/100‑1 (MQ=255)
                                                                                               gATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAAc      <  2:673370/100‑1 (MQ=255)
                                                                                               gATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAAc      >  1:959335/1‑100 (MQ=255)
                                                                                                aTGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAAcc     <  1:2188456/100‑1 (MQ=255)
                                                                                                aTGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAAcc     <  1:590471/100‑1 (MQ=255)
                                                                                                 tGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCa    <  2:528076/100‑1 (MQ=255)
                                                                                                   cGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGa  >  2:1045899/1‑100 (MQ=255)
                                                                                                   cGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGa  >  1:592560/1‑100 (MQ=255)
                                                                                                   |                                                                                                    
CCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGA  >  NC_000913/949642‑949840

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