Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 949,745:1 +G 5.7% coding (78/591 nt) ycaK → putative NAD(P)H‑dependent oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913949,7451.G 5.7% 478.7 / 22.6 106coding (78/591 nt)ycaKputative NAD(P)H‑dependent oxidoreductase
Reads supporting (aligned to +/- strand):  ref base . (50/50);  new base G (3/3);  total (53/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 = 7.75e-01

Only 100 of 106 total aligned reads displayed.

ATCACGTTTACAGGATTATGCCGCCTGCGTAAAGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGTTCATCAGCTTTATT  >  NC_000913/949602‑949890
                                                                                                                                                |                                                                                                                                                 
atcaCGTTTACAGGATTATGCCGCCTGCGTAAAGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCa                                                                                                                                              <  2:368967/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)
           aGGATTATGCCGCCTGCGTAAAGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATgg                                                                                                                                   >  1:633250/1‑149 (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                                                                                                                         <  2:633250/149‑1 (MQ=255)
                     cgccTGCGTAAAGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAAt                                                                                                                         >  4:277760/1‑149 (MQ=255)
                     cgccTGCGTAAAGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGGATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAAt                                                                                                                         >  2:521910/1‑149 (MQ=255)
                      gccTGCGTAAAGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATCGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATa                                                                                                                        >  4:656494/1‑149 (MQ=255)
                        cTGCGTAAAGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACa                                                                                                                      <  2:221172/149‑1 (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)
                                  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                                                               >  2:1876/1‑149 (MQ=255)
                                                                               gTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGaa                                                               <  2:106009/149‑1 (MQ=255)
                                                                               gTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGaa                                                               >  1:276514/1‑149 (MQ=255)
                                                                               gTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCAGaa                                                               >  4:377479/1‑149 (MQ=255)
                                                                                tatttattTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAg                                                              >  1:31944/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)
                                                                                      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:655744/1‑149 (MQ=255)
                                                                                                    aTCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAat                                          <  3:500866/149‑1 (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‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGcc                                                                   >  4:158883/1‑112 (MQ=255)
                                                                                                                aTTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGAc                                                   <  1:271105/128‑1 (MQ=255)
                                                                                                                aTTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGt                              >  3:297094/1‑149 (MQ=255)
                                                                                                                  tCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTa                            <  2:527611/149‑1 (MQ=255)
                                                                                                                        aCCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAAt                                                                                       <  1:257986/84‑1 (MQ=255)
                                                                                                                        aCCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAAt                                                                                       >  2:257986/1‑84 (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              <  3:304080/149‑1 (MQ=255)
                                                                                                                                tatTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGTtc              <  2:751699/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                                                                                 >  1:28483/1‑75 (MQ=255)
                                                                                                                                       gATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAAc                                                                                 <  2:28483/75‑1 (MQ=255)
                                                                                                                                       gATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGTTCATCAGCt       <  4:297094/149‑1 (MQ=255)
                                                                                                                                         tGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGTTCATCAGCttt     >  3:238401/1‑149 (MQ=255)
                                                                                                                                           cGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGTTCATCAGCTTTAt   <  3:638239/149‑1 (MQ=255)
                                                                                                                                            gATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGTTCATCAGCTTTAtt  >  4:190622/1‑149 (MQ=255)
                                                                                                                                             aTC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGtt                             <  3:197138/121‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
ATCACGTTTACAGGATTATGCCGCCTGCGTAAAGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATC‑CATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGTTCATCAGCTTTATT  >  NC_000913/949602‑949890

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: ‑