Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,679,398 T→A 6.2% intergenic (+27/‑159) 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,3980TA 6.2% 343.0 / 10.4 97intergenic (+27/‑159)ydgH/ydgIDUF1471 family periplasmic protein/putative arginine/ornithine antiporter transporter
Reads supporting (aligned to +/- strand):  ref base T (49/42);  new base A (2/4);  total (51/46)
Fisher's exact test for biased strand distribution p-value = 4.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.66e-01

CGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCACGTTTTAATCA  >  NC_000913/1679251‑1679546
                                                                                                                                                   |                                                                                                                                                    
cGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCtt                                                                                                                                                     <  2:172759/149‑1 (MQ=255)
  aaGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGa                                                                                                                                                   <  3:561727/149‑1 (MQ=255)
      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:663916/140‑1 (MQ=255)
                      aaaCGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTc                                                                                                                                        >  2:663916/1‑140 (MQ=255)
                      aaaCGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTc                                                                                                                                        <  1:102095/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                                                                                                                   >  4:789367/1‑146 (MQ=255)
                                     aaaGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTc                                                                                                                   <  3:789367/146‑1 (MQ=255)
                                      aaGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGa                                                                                                                                                   <  1:88143/113‑1 (MQ=255)
                                      aaGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGa                                                                                                                                                   >  2:88143/1‑113 (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                                                                                                                    <  3:460984/128‑1 (MQ=255)
                                                      ccACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAAttt                                                                                                                    >  4:460984/1‑128 (MQ=255)
                                                       cacaTCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGa                                                                                              >  3:667140/1‑149 (MQ=255)
                                                          aTCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGacac                                                                                           >  3:430073/1‑149 (MQ=255)
                                                           tCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATGAGGCGGCTTGATGCCGCCTTTCGCATTTTTAg                                                                                                                              >  4:548838/1‑113 (MQ=255)
                                                           tCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAg                                                                                                                              <  1:464279/113‑1 (MQ=255)
                                                           tCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAg                                                                                                                              >  2:464279/1‑113 (MQ=255)
                                                           tCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAg                                                                                                                              <  3:548838/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                                                                                   <  4:807424/112‑1 (MQ=255)
                                                                                                       gcAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTc                                                                                   >  3:807424/1‑112 (MQ=255)
                                                                                                        cAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTAt                                             <  4:793389/149‑1 (MQ=255)
                                                                                                        cAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTAt                                             <  2:904819/149‑1 (MQ=255)
                                                                                                              tGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGttt                                       >  1:312603/1‑149 (MQ=255)
                                                                                                              tGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGttt                                       <  2:312603/149‑1 (MQ=255)
                                                                                                                tataAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTg                                     >  2:266759/1‑149 (MQ=255)
                                                                                                                 ataaataaCTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGc                                    >  3:736095/1‑149 (MQ=255)
                                                                                                                 ataaataaCTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGc                                    >  1:792802/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:736095/149‑1 (MQ=255)
                                                                                                                                      ccTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTAt               <  4:134212/149‑1 (MQ=255)
                                                                                                                                      ccTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTAt               <  3:718986/149‑1 (MQ=255)
                                                                                                                                       cTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATc              >  1:899096/1‑149 (MQ=255)
                                                                                                                                        tataAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCa             <  2:617060/149‑1 (MQ=255)
                                                                                                                                        tataAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCa             <  2:899096/149‑1 (MQ=255)
                                                                                                                                        tataAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGCGGCAAACTCCATTTTTATGCGTTTTGCCCTAAtattta                         <  2:709728/137‑1 (MQ=255)
                                                                                                                                        tataAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGCGGCAAACTCCATTTTTATGCGTTTTGCCCTAAtattta                         >  1:709728/1‑137 (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     >  4:32207/1‑149 (MQ=255)
                                                                                                                                                ggcTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCACGTTTTaa     <  1:853064/149‑1 (MQ=255)
                                                                                                                                                  cTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCACGTTTTAATc   >  2:5656/1‑149 (MQ=255)
                                                                                                                                                   ttGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCACGTTTTAATCa  <  4:392769/149‑1 (MQ=255)
                                                                                                                                                   ttGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCACGTTTTAATCa  >  4:455202/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGCAGACACATCTCTTCGTAAAATCCCGCGCTTCAGTGGCAAACTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCACGTTTTAATCA  >  NC_000913/1679251‑1679546

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: