Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,679,402 T→A 6.1% intergenic (+31/‑155) ydgH → / → ydgI DUF1471 family periplasmic protein/putative arginine/ornithine antiporter transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,679,4020TA 6.1% 351.7 / 10.3 98intergenic (+31/‑155)ydgH/ydgIDUF1471 family periplasmic protein/putative arginine/ornithine antiporter transporter
Reads supporting (aligned to +/- strand):  ref base T (50/42);  new base A (2/4);  total (52/46)
Fisher's exact test for biased strand distribution p-value = 4.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.41e-01

CTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCACGTTTTAATCACT  >  NC_000913/1679257‑1679548
                                                                                                                                                 |                                                                                                                                                  
ctcCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATgcc                                                                                                                                                 <  2:341428/149‑1 (MQ=255)
    tATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCATCAAGCCGCCt                                                                                                                                             <  2:235950/149‑1 (MQ=255)
     aTCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCATCAAGCCGCCtt                                                                                                                                            <  2:718456/149‑1 (MQ=255)
      tCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTCAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCttt                                                                                                                                           >  4:352767/1‑149 (MQ=255)
        aGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCg                                                                                                                                         >  2:123012/1‑149 (MQ=255)
        aGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCATCAAGCCGCCttata                                                                                                                                         <  3:522701/149‑4 (MQ=255)
           ttGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCAt                                                                                                                                      >  4:846994/1‑149 (MQ=255)
           ttGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCATCAAGCCGCCttatagga                                                                                                                                      <  3:343436/149‑7 (MQ=255)
            tGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCAtt                                                                                                                                     <  1:378475/149‑1 (MQ=255)
               gAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCAttttt                                                                                                                                  >  3:667317/1‑149 (MQ=255)
                aaaCGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTc                                                                                                                                          <  1:102095/140‑1 (MQ=255)
                aaaCGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTc                                                                                                                                          >  2:663916/1‑140 (MQ=255)
                aaaCGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTc                                                                                                                                          <  1:663916/140‑1 (MQ=255)
                aaaCGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTc                                                                                                                                          >  2:102095/1‑140 (MQ=255)
                 aaCGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAg                                                                                                                                <  1:123012/149‑1 (MQ=255)
                 aaCGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAg                                                                                                                                >  4:397028/1‑149 (MQ=255)
                       ccGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACaa                                                                                                                          >  1:277723/1‑149 (MQ=255)
                       ccGCGAAGAAAAGTGCTAATTATTACCCCATCTCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACaa                                                                                                                          <  2:299561/149‑1 (MQ=255)
                        cgcgAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACaaa                                                                                                                         <  1:435052/149‑1 (MQ=255)
                          cgAAGAAAGGTGCTAAGTATTACCACATCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAAtt                                                                                                                       >  3:793389/1‑149 (MQ=255)
                               aaaGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTc                                                                                                                     <  3:789367/146‑1 (MQ=255)
                               aaaGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTc                                                                                                                     >  4:789367/1‑146 (MQ=255)
                                      cTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCa                                                                                                           >  4:387329/1‑149 (MQ=255)
                                          gTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATtgc                                                                                                       <  2:214376/149‑1 (MQ=255)
                                           tATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATtgca                                                                                                      >  1:714679/1‑149 (MQ=255)
                                            aTTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATtgcat                                                                                                     <  3:793004/149‑1 (MQ=255)
                                             ttACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATtgcatg                                                                                                    >  1:785438/1‑149 (MQ=255)
                                             ttACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATtgc                                                                                                       >  3:611817/1‑146 (MQ=255)
                                             ttACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATtgc                                                                                                       <  4:611817/146‑1 (MQ=255)
                                              tACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATtgcatgc                                                                                                   >  3:134212/1‑149 (MQ=255)
                                                ccACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAAttt                                                                                                                      >  4:460984/1‑128 (MQ=255)
                                                ccACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAAttt                                                                                                                      <  3:460984/128‑1 (MQ=255)
                                                 cacaTCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGa                                                                                                >  3:667140/1‑149 (MQ=255)
                                                    aTCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGacac                                                                                             >  3:430073/1‑149 (MQ=255)
                                                     tCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATGAGGCGGCTTGATGCCGCCTTTCGCATTTTTAg                                                                                                                                >  4:548838/1‑113 (MQ=255)
                                                     tCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAg                                                                                                                                <  3:548838/113‑1 (MQ=255)
                                                     tCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAg                                                                                                                                >  2:464279/1‑113 (MQ=255)
                                                     tCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAg                                                                                                                                <  1:464279/113‑1 (MQ=255)
                                                         ccGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACAtctc                                                                                        >  2:466510/1‑149 (MQ=255)
                                                          cGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATctct                                                                                       >  3:373476/1‑149 (MQ=255)
                                                           gCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCtt                                                                                      >  3:392769/1‑149 (MQ=255)
                                                            ccAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTc                                                                                     <  4:373476/149‑1 (MQ=255)
                                                              aGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGt                                                                                   <  1:128925/149‑1 (MQ=255)
                                                                    aGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATc                                                                             >  1:904819/1‑149 (MQ=255)
                                                                      gAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATccc                                                                           <  3:846994/149‑1 (MQ=255)
                                                                      gAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATccc                                                                           <  4:3022/149‑1 (MQ=255)
                                                                                  aaCCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTgg                                                               >  3:317095/1‑149 (MQ=255)
                                                                                     cTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCaa                                                            >  3:12258/1‑149 (MQ=255)
                                                                                      tGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCaaa                                                           <  2:319517/149‑1 (MQ=255)
                                                                                      tGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCaaa                                                           >  3:48195/1‑149 (MQ=255)
                                                                                                 gcAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTc                                                                                     >  3:807424/1‑112 (MQ=255)
                                                                                                 gcAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTc                                                                                     <  4:807424/112‑1 (MQ=255)
                                                                                                  cAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTAt                                               <  2:904819/149‑1 (MQ=255)
                                                                                                  cAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTAt                                               <  4:793389/149‑1 (MQ=255)
                                                                                                        tGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGttt                                         <  2:312603/149‑1 (MQ=255)
                                                                                                        tGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGttt                                         >  1:312603/1‑149 (MQ=255)
                                                                                                          tataAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTg                                       >  2:266759/1‑149 (MQ=255)
                                                                                                           ataaataaCTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGc                                      >  1:792802/1‑149 (MQ=255)
                                                                                                           ataaataaCTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGc                                      >  3:736095/1‑149 (MQ=255)
                                                                                                             aaataaCTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGccc                                    >  4:334798/1‑149 (MQ=255)
                                                                                                                taaCTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTaa                                 <  1:858624/149‑1 (MQ=255)
                                                                                                                 aaCTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAat                                >  3:865084/1‑149 (MQ=255)
                                                                                                                    tATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAAtatt                             >  4:846996/1‑149 (MQ=255)
                                                                                                                     aTTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAAtattt                            >  3:355189/1‑149 (MQ=255)
                                                                                                                       tCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTGGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAAtatttat                          <  3:352767/149‑1 (MQ=255)
                                                                                                                       tCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAAtatttat                          <  3:429511/149‑1 (MQ=255)
                                                                                                                        cAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAAtatttatt                         <  2:792802/149‑1 (MQ=255)
                                                                                                                          gTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTAttct                       >  3:22961/1‑149 (MQ=255)
                                                                                                                           ttACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCa                                                      >  1:753730/1‑117 (MQ=255)
                                                                                                                           ttACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCa                                                      <  2:753730/117‑1 (MQ=255)
                                                                                                                           ttACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTg                                       <  1:205702/132‑1 (MQ=255)
                                                                                                                           ttACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTg                                       >  2:205702/1‑132 (MQ=255)
                                                                                                                               tCCTATAAGGTGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATACAGACAGATCTCTTCGTAAAATCCCGCGCTTCAGTGGCACATTCCATTTTTATGCGTTTTGCCCTAATATTTATTCCTCTTa                  >  4:600230/1‑149 (MQ=255)
                                                                                                                               tCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTa                  >  3:113441/1‑149 (MQ=255)
                                                                                                                               tCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTa                  >  3:255075/1‑149 (MQ=255)
                                                                                                                                ccTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTAt                 <  4:134212/149‑1 (MQ=255)
                                                                                                                                ccTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTAt                 <  4:736095/149‑1 (MQ=255)
                                                                                                                                ccTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTAt                 <  3:718986/149‑1 (MQ=255)
                                                                                                                                 cTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATc                >  1:899096/1‑149 (MQ=255)
                                                                                                                                  tataAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCa               <  2:899096/149‑1 (MQ=255)
                                                                                                                                  tataAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCa               <  2:617060/149‑1 (MQ=255)
                                                                                                                                  tataAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGCGGCAAACTCCATTTTTATGCGTTTTGCCCTAAtattta                           >  1:709728/1‑137 (MQ=255)
                                                                                                                                  tataAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGCGGCAAACTCCATTTTTATGCGTTTTGCCCTAAtattta                           <  2:709728/137‑1 (MQ=255)
                                                                                                                                   gaaaGGCGGCATCAAGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCAc              >  3:319335/3‑149 (MQ=255)
                                                                                                                                   ataAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCAc              <  4:113441/149‑1 (MQ=255)
                                                                                                                                    aaaGGCGGCATCAAGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCACg             >  1:597898/2‑149 (MQ=255)
                                                                                                                                     aaGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCACGt            >  1:931450/1‑149 (MQ=255)
                                                                                                                                      aGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCACGtt           <  1:910131/149‑1 (MQ=255)
                                                                                                                                       ggcggcTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCACGttt          <  4:667140/149‑1 (MQ=255)
                                                                                                                                        gcggcTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCACGtttt         >  3:671304/1‑149 (MQ=255)
                                                                                                                                          ggcTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCACGTTTTaa       <  4:430073/149‑1 (MQ=255)
                                                                                                                                          ggcTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCACGTTTTaa       <  1:853064/149‑1 (MQ=255)
                                                                                                                                          ggcTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCACGTTTTaa       >  4:32207/1‑149 (MQ=255)
                                                                                                                                            cTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCACGTTTTAATc     >  2:5656/1‑149 (MQ=255)
                                                                                                                                             ttGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCACGTTTTAATCa    >  4:455202/1‑149 (MQ=255)
                                                                                                                                             ttGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCACGTTTTAATCa    <  4:392769/149‑1 (MQ=255)
                                                                                                                                              tGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCACGTTTTAATCAc   <  2:277723/149‑1 (MQ=255)
                                                                                                                                               gATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCACGTTTTAATCACt  >  3:208060/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCACGTTTTAATCACT  >  NC_000913/1679257‑1679548

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