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.

ATCACGTTTACAGGATTATGCCGCCTGCGTAAAGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGTTCATCAGCTTTATT  >  NC_000913/949602‑949890
                                                                                                                                                |                                                                                                                                                 
atcaCGTTTACAGGATTATGCCGCCTGCGTAAAGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCa                                                                                                                                              <  2:368967/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)
           aGGATTATGCCGCCTGCGTAAAGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATgg                                                                                                                                   >  1:633250/1‑149 (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                                                                                                                         <  2:633250/149‑1 (MQ=255)
                     cgccTGCGTAAAGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAAt                                                                                                                         >  4:277760/1‑149 (MQ=255)
                     cgccTGCGTAAAGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGGATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAAt                                                                                                                         >  2:521910/1‑149 (MQ=255)
                      gccTGCGTAAAGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATCGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATa                                                                                                                        >  4:656494/1‑149 (MQ=255)
                        cTGCGTAAAGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACa                                                                                                                      <  2:221172/149‑1 (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)
                                  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                                                               >  2:1876/1‑149 (MQ=255)
                                                                               gTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGaa                                                               <  2:106009/149‑1 (MQ=255)
                                                                               gTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGaa                                                               >  1:276514/1‑149 (MQ=255)
                                                                               gTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCAGaa                                                               >  4:377479/1‑149 (MQ=255)
                                                                                tatttattTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAg                                                              >  1:31944/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)
                                                                                      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:655744/1‑149 (MQ=255)
                                                                                                    aTCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAat                                          <  3:500866/149‑1 (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)
                                                                                                                aTTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGcc                                                                   >  4:158883/1‑112 (MQ=255)
                                                                                                                aTTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGAc                                                   <  1:271105/128‑1 (MQ=255)
                                                                                                                aTTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGt                              >  3:297094/1‑149 (MQ=255)
                                                                                                                  tCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTa                            <  2:527611/149‑1 (MQ=255)
                                                                                                                        aCCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAAt                                                                                       <  1:257986/84‑1 (MQ=255)
                                                                                                                        aCCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAAt                                                                                       >  2:257986/1‑84 (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              <  3:304080/149‑1 (MQ=255)
                                                                                                                                tatTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGTtc              <  2:751699/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                                                                                 >  1:28483/1‑75 (MQ=255)
                                                                                                                                       gATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAAc                                                                                 <  2:28483/75‑1 (MQ=255)
                                                                                                                                       gATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGTTCATCAGCt       <  4:297094/149‑1 (MQ=255)
                                                                                                                                         tGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGTTCATCAGCttt     >  3:238401/1‑149 (MQ=255)
                                                                                                                                           cGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGTTCATCAGCTTTAt   <  3:638239/149‑1 (MQ=255)
                                                                                                                                            gATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGTTCATCAGCTTTAtt  >  4:190622/1‑149 (MQ=255)
                                                                                                                                             aTCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGtt                             <  3:197138/121‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
ATCACGTTTACAGGATTATGCCGCCTGCGTAAAGTGCGCTCCAGAACTTAACGTGGAGGTAAAATTATGCAGTCTGAACGTATTTATTTGGTATGGGCCCATCCTCGTCATGATTCATTGACCGCACATATTGCTGATGCGATCCATCAGCGGGCAATGGAGCGGAAAATACAGGTGACGGAACTCGATTTATATCGGCGTAATTTCAACCCAGTGATGACGCCGGAAGATGAACCAGACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGTTCATCAGCTTTATT  >  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: