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.

TCCTGAATCATCACGTTTACAGGATTATGCCGCCTGCGTAAAGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGTTCATCAGCTTTAT  >  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)
         atcaCGTTTACAGGATTATGCCGCCTGCGTAAAGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCa                                                                                                                                             >  4:268869/1‑149 (MQ=255)
          tcaCGTTTACAGGATTATGCCGCCTGCGTAAAGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATTCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAg                                                                                                                                            >  3:447001/1‑149 (MQ=255)
          tcaCGTTTACAGGATTATGCCGCCTGCGTAAAGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAg                                                                                                                                            >  2:336528/1‑149 (MQ=255)
            aCGTTTACAGGATTATGCCGCCTGCGTAAAGTGCGCTCCAGATCTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCg                                                                                                                                          <  1:203006/149‑1 (MQ=255)
              gTTTACAGGATTATGCCGCCTGCGTAAAGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGGATCGCATCAGcaat                                                                                                                                        <  1:597381/149‑4 (MQ=255)
                   cAGGATTATGCCGCCTGCGTAAAGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATg                                                                                                                                   <  1:501874/149‑1 (MQ=255)
                      gATTATGCCGCCTGCGTAAAGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAg                                                                                                                                >  4:378187/1‑149 (MQ=255)
                       aTTATGCCGCCTGCGTAAAGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGc                                                                                                                               <  3:589250/149‑1 (MQ=255)
                       aTTATGCCGCCTGCGTAAAGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGc                                                                                                                               <  1:485345/149‑1 (MQ=255)
                         tATGCCGCCTGCGTAAAGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCgg                                                                                                                             >  3:173039/1‑149 (MQ=255)
                              cgccTGCGTAAAGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAAt                                                                                                                        >  2:412642/1‑149 (MQ=255)
                              cgccTGCGTAAAGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAAt                                                                                                                        >  4:277760/1‑149 (MQ=255)
                              cgccTGCGTAAAGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAAt                                                                                                                        <  2:633250/149‑1 (MQ=255)
                              cgccTGCGTAAAGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGGATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAAt                                                                                                                        >  2:521910/1‑149 (MQ=255)
                               gccTGCGTAAAGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATCGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATa                                                                                                                       >  4:656494/1‑149 (MQ=255)
                                  tGCGTAAAGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAg                                                                                                                    <  4:148558/149‑1 (MQ=255)
                                   gCGTAAAGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAgg                                                                                                                   >  1:751699/1‑149 (MQ=255)
                                       aaaGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGAc                                                                                                               <  4:185607/149‑1 (MQ=255)
                                       aaaGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGAc                                                                                                               >  4:500866/1‑149 (MQ=255)
                                        aaGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGGATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACg                                                                                                              <  1:521910/149‑1 (MQ=255)
                                           tGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGaa                                                                                                           <  4:173039/149‑1 (MQ=255)
                                                tCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGa                                                                                                      <  1:816413/149‑1 (MQ=255)
                                                 ccAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGAt                                                                                                     <  1:88554/149‑1 (MQ=255)
                                                    gAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTta                                                                                                  >  1:192430/1‑149 (MQ=255)
                                                     aaCTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTtat                                                                                                 <  1:353860/149‑1 (MQ=255)
                                                      aCTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTtata                                                                                                <  2:753547/149‑1 (MQ=255)
                                                        ttAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATc                                                                                              >  3:550511/1‑149 (MQ=255)
                                                           aCGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGc                                                                                           <  3:378187/149‑1 (MQ=255)
                                                               ggAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTaa                                                                                       <  2:293502/149‑1 (MQ=255)
                                                                 aGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAAtt                                                                                     <  1:412642/149‑1 (MQ=255)
                                                                  ggTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAAttt                                                                                    >  2:388304/1‑149 (MQ=255)
                                                                   gTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTc                                                                                   <  3:591145/149‑1 (MQ=255)
                                                                     aaaaTTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTtat                                                                                                 >  2:313631/1‑133 (MQ=255)
                                                                     aaaaTTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTtat                                                                                                 <  1:313631/133‑1 (MQ=255)
                                                                     aaaaTTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCaa                                                                                 <  2:192430/149‑1 (MQ=255)
                                                                       aaTTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCaa                                                                                 <  3:30682/147‑1 (MQ=255)
                                                                       aaTTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCaa                                                                                 <  3:66832/147‑1 (MQ=255)
                                                                       aaTTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCaa                                                                                 >  4:30682/1‑147 (MQ=255)
                                                                       aaTTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCaa                                                                                 >  4:66832/1‑147 (MQ=255)
                                                                             gCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGtga                                                                         >  1:110625/1‑149 (MQ=255)
                                                                                     aaCGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCg                                                                 >  1:106009/1‑149 (MQ=255)
                                                                                        gTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGaa                                                              >  1:276514/1‑149 (MQ=255)
                                                                                        gTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGaa                                                              <  2:106009/149‑1 (MQ=255)
                                                                                        gTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGaa                                                              >  2:1876/1‑149 (MQ=255)
                                                                                          atttatttGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGa                                                            >  3:193704/1‑149 (MQ=255)
                                                                                           tttatttGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGAt                                                           >  4:278712/1‑149 (MQ=255)
                                                                                           tttatttGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGAt                                                           <  1:738402/149‑1 (MQ=255)
                                                                                              atttGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGaa                                                        <  2:781309/149‑1 (MQ=255)
                                                                                               tttGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATAc                                                                                                                      <  3:775693/86‑1 (MQ=255)
                                                                                               tttGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAAc                                                       <  4:167483/149‑1 (MQ=255)
                                                                                                    tATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGAc                                                  >  1:445807/1‑149 (MQ=255)
                                                                                                         gCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGgaa                                             <  1:810905/149‑1 (MQ=255)
                                                                                                         gCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGgaa                                             <  2:110625/149‑1 (MQ=255)
                                                                                                           ccATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGAtt                                                                                                    >  4:533174/1‑92 (MQ=255)
                                                                                                             aTCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAat                                         <  3:500866/149‑1 (MQ=255)
                                                                                                             aTCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAat                                         >  3:655744/1‑149 (MQ=255)
                                                                                                               ccTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAatat                                       <  1:336528/149‑1 (MQ=255)
                                                                                                                cTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATg                                      >  2:267330/1‑149 (MQ=255)
                                                                                                                cTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATg                                      >  4:53220/1‑149 (MQ=255)
                                                                                                                cTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCAGAAGATGAACCAGACTGGAAGAATATg                                      <  3:377479/149‑1 (MQ=255)
                                                                                                                 tCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAata                                        >  2:372629/1‑146 (MQ=255)
                                                                                                                 tCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAata                                        <  1:372629/146‑1 (MQ=255)
                                                                                                                 tCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCACACTGGAAGAATATgg                                     >  2:235584/1‑149 (MQ=255)
                                                                                                                  cGTCATGATTCATTGACCGCACATATTGATGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGa                                    <  4:193704/149‑1 (MQ=255)
                                                                                                                         aTTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGcc                                                                  <  3:158883/112‑1 (MQ=255)
                                                                                                                         aTTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGAc                                                  <  1:271105/128‑1 (MQ=255)
                                                                                                                         aTTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGAc                                                  >  2:271105/1‑128 (MQ=255)
                                                                                                                         aTTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGt                             >  3:297094/1‑149 (MQ=255)
                                                                                                                           tCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTa                           <  2:527611/149‑1 (MQ=255)
                                                                                                                                 aCCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAAt                                                                                      >  2:257986/1‑84 (MQ=255)
                                                                                                                                 aCCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAAt                                                                                      <  1:257986/84‑1 (MQ=255)
                                                                                                                                    gCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGaa                                                                                                           <  3:391955/60‑1 (MQ=255)
                                                                                                                                    gCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGaa                                                                                                           >  4:391955/1‑60 (MQ=255)
                                                                                                                                    gCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCCACCCAGTGATGACGCCGGAAGATGAAc                                                       <  3:533595/112‑1 (MQ=255)
                                                                                                                                    gCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAAc                                                       >  4:533595/1‑112 (MQ=255)
                                                                                                                                    gCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTATTCTCCaga                  >  1:693025/1‑149 (MQ=255)
                                                                                                                                    gCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTATTCTCCaga                  >  3:19361/1‑149 (MQ=255)
                                                                                                                                       tatatTGCTGATGTGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAGTTTCAACCCAGTGACGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGt               >  4:517716/2‑149 (MQ=255)
                                                                                                                                        atatTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCTATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGtt              >  2:619933/1‑149 (MQ=255)
                                                                                                                                         tatTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGTtc             <  2:751699/149‑1 (MQ=255)
                                                                                                                                         tatTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGTtc             <  3:304080/149‑1 (MQ=255)
                                                                                                                                          cccgCTGATGGATCGCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAatat                                       >  4:146192/4‑122 (MQ=255)
                                                                                                                                          cccgCTGATGGATCGCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAatat                                       <  3:146192/119‑1 (MQ=255)
                                                                                                                                           ttGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGTtcat           >  2:425822/1‑149 (MQ=255)
                                                                                                                                           ccgCTGATGGATCGCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGa                                                                                                      >  4:194884/3‑58 (MQ=25)
                                                                                                                                           ccgCTGATGGATCGCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGa                                                                                                      <  3:194884/56‑1 (MQ=25)
                                                                                                                                           ccgCTGATGGATCGCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGTtcat           >  1:620485/3‑149 (MQ=255)
                                                                                                                                               tGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGTTCATCAGc       <  2:586840/149‑1 (MQ=255)
                                                                                                                                               tGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGTTCATCAGc       >  4:638239/1‑149 (MQ=255)
                                                                                                                                                gATGCGATCCATTAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGTTCATCAGCt      <  3:371758/149‑1 (MQ=255)
                                                                                                                                                gATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAAc                                                                                <  2:28483/75‑1 (MQ=255)
                                                                                                                                                gATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAAc                                                                                >  1:28483/1‑75 (MQ=255)
                                                                                                                                                gATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGTTCATCAGCt      <  4:297094/149‑1 (MQ=255)
                                                                                                                                                    cGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGTTCATCAGCTTTAt  <  3:638239/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TCCTGAATCATCACGTTTACAGGATTATGCCGCCTGCGTAAAGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGTTCATCAGCTTTAT  >  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: