Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 949,741 Δ1 bp 17.2% 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.17.2% 352.0 / 75.6 122coding (74/591 nt)ycaKputative NAD(P)H‑dependent oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (51/49);  new base . (10/11);  total (62/60)
Fisher's exact test for biased strand distribution p-value = 8.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

Only 100 of 130 total aligned reads displayed.

CCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGA  >  NC_000913/949642‑949840
                                                                                                   |                                                                                                    
tccgatcTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGc                                                                                                      <  2:124169/94‑1 (MQ=255)
ccAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGc                                                                                                      <  2:1572359/100‑1 (MQ=255)
ccAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGc                                                                                                      >  2:78509/1‑100 (MQ=255)
ccAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGc                                                                                                      >  1:1348718/1‑100 (MQ=255)
     aCTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATc‑c                                                                                                >  2:2524583/1‑100 (MQ=255)
      cTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATc‑c                                                                                                <  2:1649620/99‑1 (MQ=255)
        tAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGcagatcgga                                                                                              >  1:124168/1‑92 (MQ=255)
        tAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATc                                                                                             >  2:2505571/1‑100 (MQ=255)
        tAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATc                                                                                             <  2:1176270/100‑1 (MQ=255)
         aaCGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCa                                                                                            <  2:2270093/100‑1 (MQ=255)
          aCGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGGATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAg                                                                                           <  1:229424/100‑1 (MQ=255)
           cGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGc                                                                                          >  2:2228166/1‑100 (MQ=255)
           cGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATG‑GATCGCATCAGc                                                                                          <  2:2741809/100‑1 (MQ=255)
             tGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATG‑GATCGCATCAGcaa                                                                                        <  2:1018510/100‑3 (MQ=255)
             tGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATG‑GATCGCATCAGcaa                                                                                        <  2:249774/100‑3 (MQ=255)
             tGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATG‑GATCGCATCAGcaa                                                                                        <  2:2888570/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‑GATCGCATCAGcaata                                                                                      <  2:1334065/100‑5 (MQ=255)
               gAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATG‑GATCGCATCAGcaa                                                                                        <  1:785434/98‑3 (MQ=255)
               gAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATCGACCGCACATATTGCTGATG‑GATCGCATCAGcaata                                                                                      <  2:2056978/100‑5 (MQ=255)
                aGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCa                                                                                     >  1:2686246/1‑100 (MQ=255)
                aGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATG‑GATCGCATCAGcaatat                                                                                     <  2:1478261/100‑6 (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)
                    aaaaTTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATG‑GATCGCATCAGcaatatgtgc                                                                                 <  2:1837456/100‑10 (MQ=38)
                        ttATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCg                                                                             <  2:236321/100‑1 (MQ=255)
                        ttATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCg                                                                             <  1:491517/100‑1 (MQ=255)
                        ttATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCg                                                                             <  1:1111708/100‑1 (MQ=255)
                          aTGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCgg                                                                            >  1:1018672/1‑99 (MQ=255)
                              aGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAAt                                                                       <  2:1484219/100‑1 (MQ=255)
                                tCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATAc                                                                     >  2:2031923/1‑100 (MQ=255)
                                  tGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAg                                                                   <  2:1687266/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)
                                      cGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGa                                                               <  2:1636206/100‑1 (MQ=255)
                                        tatttattTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACg                                                             >  2:1851579/1‑100 (MQ=255)
                                          tggatttGGTATGGGCCTATCGTCGTGGTTATTGATTTGCCGCACAAATTGCTGATGCGTTC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGa                                                           <  1:1631039/97‑1 (MQ=255)
                                           ttatttGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGaa                                                          >  1:1859547/1‑100 (MQ=255)
                                            tatttGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAAc                                                         <  1:2697049/100‑1 (MQ=255)
                                            tatttGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAAc                                                         <  2:1028803/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)
                                                   tATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGAttt                                                  <  2:1081660/100‑1 (MQ=255)
                                                    aTGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTta                                                 <  2:486887/100‑1 (MQ=255)
                                                    aTGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTta                                                 <  2:2822571/100‑1 (MQ=255)
                                                      gggCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTtata                                               >  1:1163373/1‑100 (MQ=255)
                                                       ggCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTatat                                              >  1:2885456/1‑100 (MQ=255)
                                                       ggCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTatat                                              >  1:21369/1‑100 (MQ=255)
                                                        gCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATc                                             >  1:538790/1‑100 (MQ=255)
                                                        gCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATc                                             >  2:1683201/1‑100 (MQ=255)
                                                          ccATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCgg                                           <  2:2496807/100‑1 (MQ=255)
                                                           cATCCTCGACATGATTCATTGACGGCACATATTGCTGCTGGTATC‑GATCTGCGGGCTATGACGCGGAAAATACAGGTGCTGGACGTCGGTTtata                                               >  1:1723818/1‑95 (MQ=255)
                                                            aTCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCg                                         <  2:2221990/100‑1 (MQ=255)
                                                              ccTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTa                                       >  1:2599296/1‑100 (MQ=255)
                                                               cTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTaa                                      >  1:2031808/1‑100 (MQ=255)
                                                                    cATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCa                                 >  2:2164093/1‑100 (MQ=255)
                                                                     aTGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCaa                                <  1:1577431/100‑1 (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‑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:1349009/1‑100 (MQ=255)
                                                                                aCCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGtgatga                     >  2:1936673/1‑100 (MQ=255)
                                                                                 ccGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGAc                    >  1:1504607/1‑100 (MQ=255)
                                                                                  cGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACg                   >  1:530629/1‑100 (MQ=255)
                                                                                   gCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGc                  >  1:2634329/1‑100 (MQ=255)
                                                                                   gCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGc                  >  1:352719/1‑100 (MQ=255)
                                                                                    cacaTATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGcc                 >  1:282305/1‑100 (MQ=255)
                                                                                       atatTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGa              <  1:1619011/100‑1 (MQ=255)
                                                                                        tatTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGaa             <  1:47580/100‑1 (MQ=255)
                                                                                         cccgCTGATG‑GATCGCATCAGCGGGCAATGGAGCGGCAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAg            >  2:2342296/4‑100 (MQ=255)
                                                                                         cccgCTGATG‑GATCGCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAg            >  1:2370841/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)
                                                                                           cgCTGATG‑GATCGCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGAt          >  1:1488450/2‑100 (MQ=255)
                                                                                            gCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATg         <  1:1867297/100‑1 (MQ=255)
                                                                                            gCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATg         >  2:142666/1‑100 (MQ=255)
                                                                                            gCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATg         <  2:842531/100‑1 (MQ=255)
                                                                                            gCTGATG‑GATCGCATCAGCGGGCAATGTAGCGGAAATTACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGtgatgatgacgcaa             >  1:2833899/1‑88 (MQ=255)
                                                                                            gCTGATG‑GATCGCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATg         >  2:2439409/1‑100 (MQ=255)
                                                                                             cTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGa        >  2:509913/1‑100 (MQ=255)
                                                                                             cTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGa        >  1:1520497/1‑100 (MQ=255)
                                                                                               gATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACACAGTGATGACGCCGGAGGATGAAc      >  1:1003651/1‑100 (MQ=255)
                                                                                                aTGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAAcc     <  2:1714685/100‑1 (MQ=255)
                                                                                                aTGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAAcc     <  1:2855512/100‑1 (MQ=255)
                                                                                                   cGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATGTATATCGGCGTAATTTCAACCCAGTGATGACGCCCGAAGATGAACCAGa  <  1:201700/100‑1 (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: ‑