Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 949,741 Δ1 bp 5.6% 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. 5.6% 375.9 / 15.8 107coding (74/591 nt)ycaKputative NAD(P)H‑dependent oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (50/50);  new base . (3/3);  total (54/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

Only 100 of 110 total aligned reads displayed.

TCCTGAATCATCACGTTTACAGGATTATGCCGCCTGCGTAAAGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGTTCATCAGCTTTAT  >  NC_000913/949593‑949889
                                                                                                                                                    |                                                                                                                                                     
tCCTGAATCATCACGTTTACAGGATTATGCCGCCTGCGTAAAGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGc                                                                                                                                                       <  3:341081/149‑1 (MQ=255)
 ccTGAATCATCACGTTTACAGGATTATGCCGCCTGCGTAAAGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCg                                                                                                                                                      >  2:485345/1‑149 (MQ=255)
  cTGAATCATCACGTTTACAGGATTATGCCGCCTGCGTAAAGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGa                                                                                                                                                     >  2:88554/1‑149 (MQ=255)
   tGAATCATCACGTTTACAGGATTATGCCGCCTGCGTAAAGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGa                                                                                                                                                     <  2:275186/148‑1 (MQ=255)
   tGAATCATCACGTTTACAGGATTATGCCGCCTGCGTAAAGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGa                                                                                                                                                     >  1:275186/1‑148 (MQ=255)
   tGAATCATCACGTTTACAGGATTATGCCGCCTGCGTAAAGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGAt                                                                                                                                                    <  2:301973/149‑1 (MQ=255)
         atcaCGTTTACAGGATTATGCCGCCTGCGTAAAGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCa                                                                                                                                             >  4:268869/1‑149 (MQ=255)
          tcaCGTTTACAGGATTATGCCGCCTGCGTAAAGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATTCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAg                                                                                                                                            >  3:447001/1‑149 (MQ=255)
          tcaCGTTTACAGGATTATGCCGCCTGCGTAAAGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAg                                                                                                                                            >  2:336528/1‑149 (MQ=255)
            aCGTTTACAGGATTATGCCGCCTGCGTAAAGTGCGCTCCAGATCTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCg                                                                                                                                          <  1:203006/149‑1 (MQ=255)
              gTTTACAGGATTATGCCGCCTGCGTAAAGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATG‑GATCGCATCAGcaat                                                                                                                                        <  1:597381/149‑4 (MQ=255)
                   cAGGATTATGCCGCCTGCGTAAAGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATg                                                                                                                                   <  1:501874/149‑1 (MQ=255)
                      gATTATGCCGCCTGCGTAAAGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAg                                                                                                                                >  4:378187/1‑149 (MQ=255)
                       aTTATGCCGCCTGCGTAAAGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGc                                                                                                                               <  3:589250/149‑1 (MQ=255)
                       aTTATGCCGCCTGCGTAAAGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGc                                                                                                                               <  1:485345/149‑1 (MQ=255)
                         tATGCCGCCTGCGTAAAGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCgg                                                                                                                             >  3:173039/1‑149 (MQ=255)
                              cgccTGCGTAAAGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAAt                                                                                                                        >  2:412642/1‑149 (MQ=255)
                              cgccTGCGTAAAGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAAt                                                                                                                        >  4:277760/1‑149 (MQ=255)
                              cgccTGCGTAAAGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAAt                                                                                                                        <  2:633250/149‑1 (MQ=255)
                              cgccTGCGTAAAGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGGATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAAt                                                                                                                        >  2:521910/1‑149 (MQ=255)
                               gccTGCGTAAAGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATCGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATa                                                                                                                       >  4:656494/1‑149 (MQ=255)
                                  tGCGTAAAGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAg                                                                                                                    <  4:148558/149‑1 (MQ=255)
                                   gCGTAAAGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAgg                                                                                                                   >  1:751699/1‑149 (MQ=255)
                                       aaaGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGAc                                                                                                               <  4:185607/149‑1 (MQ=255)
                                       aaaGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGAc                                                                                                               >  4:500866/1‑149 (MQ=255)
                                        aaGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGGATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACg                                                                                                              <  1:521910/149‑1 (MQ=255)
                                           tGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGaa                                                                                                           <  4:173039/149‑1 (MQ=255)
                                                tCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGa                                                                                                      <  1:816413/149‑1 (MQ=255)
                                                 ccAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGAt                                                                                                     <  1:88554/149‑1 (MQ=255)
                                                    gAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTta                                                                                                  >  1:192430/1‑149 (MQ=255)
                                                     aaCTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTtat                                                                                                 <  1:353860/149‑1 (MQ=255)
                                                      aCTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTtata                                                                                                <  2:753547/149‑1 (MQ=255)
                                                        ttAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATc                                                                                              >  3:550511/1‑149 (MQ=255)
                                                           aCGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGc                                                                                           <  3:378187/149‑1 (MQ=255)
                                                               ggAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTaa                                                                                       <  2:293502/149‑1 (MQ=255)
                                                                 aGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAAtt                                                                                     <  1:412642/149‑1 (MQ=255)
                                                                  ggTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAAttt                                                                                    >  2:388304/1‑149 (MQ=255)
                                                                   gTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTc                                                                                   <  3:591145/149‑1 (MQ=255)
                                                                     aaaaTTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTtat                                                                                                 >  2:313631/1‑133 (MQ=255)
                                                                     aaaaTTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTtat                                                                                                 <  1:313631/133‑1 (MQ=255)
                                                                     aaaaTTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCaa                                                                                 <  2:192430/149‑1 (MQ=255)
                                                                       aaTTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCaa                                                                                 <  3:30682/147‑1 (MQ=255)
                                                                       aaTTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCaa                                                                                 <  3:66832/147‑1 (MQ=255)
                                                                       aaTTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCaa                                                                                 >  4:30682/1‑147 (MQ=255)
                                                                       aaTTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCaa                                                                                 >  4:66832/1‑147 (MQ=255)
                                                                             gCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGtga                                                                         >  1:110625/1‑149 (MQ=255)
                                                                                     aaCGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCg                                                                 >  1:106009/1‑149 (MQ=255)
                                                                                        gTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGaa                                                              >  1:276514/1‑149 (MQ=255)
                                                                                        gTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGaa                                                              <  2:106009/149‑1 (MQ=255)
                                                                                        gTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGaa                                                              >  2:1876/1‑149 (MQ=255)
                                                                                          atttatttGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGa                                                            >  3:193704/1‑149 (MQ=255)
                                                                                           tttatttGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGAt                                                           >  4:278712/1‑149 (MQ=255)
                                                                                           tttatttGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGAt                                                           <  1:738402/149‑1 (MQ=255)
                                                                                              atttGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGaa                                                        <  2:781309/149‑1 (MQ=255)
                                                                                               tttGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATAc                                                                                                                      <  3:775693/86‑1 (MQ=255)
                                                                                               tttGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAAc                                                       <  4:167483/149‑1 (MQ=255)
                                                                                                    tATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGAc                                                  >  1:445807/1‑149 (MQ=255)
                                                                                                         gCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGgaa                                             <  1:810905/149‑1 (MQ=255)
                                                                                                         gCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGgaa                                             <  2:110625/149‑1 (MQ=255)
                                                                                                           ccATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGAtt                                                                                                    >  4:533174/1‑92 (MQ=255)
                                                                                                             aTCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAat                                         <  3:500866/149‑1 (MQ=255)
                                                                                                             aTCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAat                                         >  3:655744/1‑149 (MQ=255)
                                                                                                               ccTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAatat                                       <  1:336528/149‑1 (MQ=255)
                                                                                                                cTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATg                                      >  2:267330/1‑149 (MQ=255)
                                                                                                                cTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATg                                      >  4:53220/1‑149 (MQ=255)
                                                                                                                cTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCAGAAGATGAACCAGACTGGAAGAATATg                                      <  3:377479/149‑1 (MQ=255)
                                                                                                                 tCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAata                                        >  2:372629/1‑146 (MQ=255)
                                                                                                                 tCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAata                                        <  1:372629/146‑1 (MQ=255)
                                                                                                                 tCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCACACTGGAAGAATATgg                                     >  2:235584/1‑149 (MQ=255)
                                                                                                                  cGTCATGATTCATTGACCGCACATATTGATGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGa                                    <  4:193704/149‑1 (MQ=255)
                                                                                                                         aTTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGcc                                                                  <  3:158883/112‑1 (MQ=255)
                                                                                                                         aTTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGAc                                                  <  1:271105/128‑1 (MQ=255)
                                                                                                                         aTTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGAc                                                  >  2:271105/1‑128 (MQ=255)
                                                                                                                         aTTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGt                             >  3:297094/1‑149 (MQ=255)
                                                                                                                           tCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTa                           <  2:527611/149‑1 (MQ=255)
                                                                                                                                 aCCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAAt                                                                                      >  2:257986/1‑84 (MQ=255)
                                                                                                                                 aCCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAAt                                                                                      <  1:257986/84‑1 (MQ=255)
                                                                                                                                    gCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGaa                                                                                                           <  3:391955/60‑1 (MQ=255)
                                                                                                                                    gCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGaa                                                                                                           >  4:391955/1‑60 (MQ=255)
                                                                                                                                    gCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCCACCCAGTGATGACGCCGGAAGATGAAc                                                       <  3:533595/112‑1 (MQ=255)
                                                                                                                                    gCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAAc                                                       >  4:533595/1‑112 (MQ=255)
                                                                                                                                    gCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTATTCTCCaga                  >  1:693025/1‑149 (MQ=255)
                                                                                                                                    gCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTATTCTCCaga                  >  3:19361/1‑149 (MQ=255)
                                                                                                                                       tatatTGCTGATGTGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAGTTTCAACCCAGTGACGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGt               >  4:517716/2‑149 (MQ=255)
                                                                                                                                        atatTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCTATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGtt              >  2:619933/1‑149 (MQ=255)
                                                                                                                                         tatTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGTtc             <  2:751699/149‑1 (MQ=255)
                                                                                                                                         tatTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGTtc             <  3:304080/149‑1 (MQ=255)
                                                                                                                                          cccgCTGATG‑GATCGCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAatat                                       >  4:146192/4‑122 (MQ=255)
                                                                                                                                          cccgCTGATG‑GATCGCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAatat                                       <  3:146192/119‑1 (MQ=255)
                                                                                                                                           ttGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGTtcat           >  2:425822/1‑149 (MQ=255)
                                                                                                                                           ccgCTGATG‑GATCGCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGa                                                                                                      >  4:194884/3‑58 (MQ=25)
                                                                                                                                           ccgCTGATG‑GATCGCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGa                                                                                                      <  3:194884/56‑1 (MQ=25)
                                                                                                                                           ccgCTGATG‑GATCGCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGTtcat           >  1:620485/3‑149 (MQ=255)
                                                                                                                                               tGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGTTCATCAGc       <  2:586840/149‑1 (MQ=255)
                                                                                                                                               tGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGTTCATCAGc       >  4:638239/1‑149 (MQ=255)
                                                                                                                                                gATGCGATC‑CATTAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGTTCATCAGCt      <  3:371758/149‑1 (MQ=255)
                                                                                                                                                gATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAAc                                                                                <  2:28483/75‑1 (MQ=255)
                                                                                                                                                gATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAAc                                                                                >  1:28483/1‑75 (MQ=255)
                                                                                                                                                gATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGTTCATCAGCt      <  4:297094/149‑1 (MQ=255)
                                                                                                                                                    cGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGTTCATCAGCTTTAt  <  3:638239/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TCCTGAATCATCACGTTTACAGGATTATGCCGCCTGCGTAAAGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGTTCATCAGCTTTAT  >  NC_000913/949593‑949889

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