Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 659,205 A→T 6.0% intergenic (+55/+46) tatE → / ← lipA TatABCE protein translocation system subunit/lipoate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913659,2050AT 6.0% 295.1 / 6.4 84intergenic (+55/+46)tatE/lipATatABCE protein translocation system subunit/lipoate synthase
Reads supporting (aligned to +/- strand):  ref base A (40/38);  new base T (2/3);  total (42/42)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GGGTTCAAGAAGGCGATGAATGATGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATGGGTAAAGCCCATCGCCAGCGCTTCGGC  >  NC_000913/659058‑659352
                                                                                                                                                   |                                                                                                                                                   
gggTTCAAGAAGGCGATGAATGATGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCaa                                                                                                                                                    <  3:775434/149‑1 (MQ=255)
   ttCAAGAAGGCGATGAATGATGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCaaaaa                                                                                                                                                 <  4:765425/149‑1 (MQ=255)
   ttCAAGAAGGCGATGAATGATGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCaaaaa                                                                                                                                                 <  4:602344/149‑1 (MQ=255)
     cAAGAAGGCGATGAATGATGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGc                                                                                                                                               <  3:831411/149‑1 (MQ=255)
       agaagGCGATGAATGATGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCgg                                                                                                                                             <  3:425669/149‑1 (MQ=255)
       agaagGCGATGAATGATGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCg                                                                                                                                              >  1:492607/1‑148 (MQ=255)
       agaagGCGATGAATGATGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGAGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCgg                                                                                                                                             <  4:145904/149‑1 (MQ=255)
           ggCGATGAATGATGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGttt                                                                                                                                         <  3:784933/149‑1 (MQ=255)
           ggCGATGAATGATGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGtt                                                                                                                                          >  1:136190/1‑148 (MQ=255)
           ggCGATGAATGATGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGtt                                                                                                                                          >  2:480651/1‑148 (MQ=255)
           ggCGATGAATGATGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGtt                                                                                                                                          <  2:136190/148‑1 (MQ=255)
            gCGATGAATGATGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTTTTGAAGCGGGtttt                                                                                                                                        <  4:852118/149‑1 (MQ=255)
            gCGATGAATGATGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGttt                                                                                                                                         <  1:506007/148‑1 (MQ=255)
             cGATGAATGATGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTTACGTGGCGAGCAGGACGCTCCCTCTATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGtttt                                                                                                                                        <  2:492607/148‑1 (MQ=255)
             cGATGAATGATGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATATTGTTCGATACACAATTCCCGCTTCCGAAAGCGGGtttt                                                                                                                                        <  2:469108/148‑1 (MQ=255)
                tgaatgaTGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTa                                                                                                                                     >  2:403526/1‑148 (MQ=255)
                  aatgaTGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATca                                                                                                                                  <  4:742143/149‑1 (MQ=255)
                  aatgaTGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATc                                                                                                                                   <  2:184850/148‑1 (MQ=255)
                    tgatgaCGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATcaga                                                                                                                                >  3:851552/1‑149 (MQ=255)
                       tgaCGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATcagaca                                                                                                                              >  1:282058/1‑148 (MQ=255)
                        gaCGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATcagacaga                                                                                                                            >  4:794943/1‑149 (MQ=255)
                             tgctgcGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATgtaa                                                                                                                       <  4:824414/149‑1 (MQ=255)
                             tgctgcGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATgta                                                                                                                        <  1:403526/148‑1 (MQ=255)
                                    gcgAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAatt                                                                                                                 >  1:628167/1‑148 (MQ=255)
                                       aaaaaaGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATATAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAAttatta                                                                                                             <  3:855050/149‑1 (MQ=255)
                                       aaaaaaGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAAttatta                                                                                                             <  3:794943/149‑1 (MQ=255)
                                          aaaGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACa                                                                                                           <  2:104868/148‑1 (MQ=255)
                                             ggCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGAt                                                                                                       <  4:851552/149‑1 (MQ=255)
                                                gcAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTTTTGAAGCGGGttt                                                                                                                                         <  4:26815/112‑1 (MQ=255)
                                                gcAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTTTTGAAGCGGGttt                                                                                                                                         >  3:26815/1‑112 (MQ=255)
                                                       ttGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACt                                                                                              >  1:844349/1‑148 (MQ=255)
                                                          aTCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATc                                                                                                                                   <  1:400201/108‑1 (MQ=255)
                                                          aTCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATc                                                                                                                                   >  2:400201/1‑108 (MQ=255)
                                                          aTCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCa                                                                                          >  3:99857/1‑149 (MQ=255)
                                                            cTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAat                                                                                                                  >  2:547929/1‑123 (MQ=255)
                                                            cTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAat                                                                                                                  <  1:547929/123‑1 (MQ=255)
                                                            cTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCAt                                                                                         <  2:628167/148‑1 (MQ=255)
                                                             ttCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATcc                                                                                       >  3:663398/1‑149 (MQ=255)
                                                              tCAGGCTGAAAAGCTCTCTCATGAAGAGTGACGTGGCGAGCAGGACGCTCCGTCAATATGTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTCTCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCGTccc                                                                                      >  4:762253/1‑149 (MQ=255)
                                                               cAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCt                                                                                     >  4:627075/1‑149 (MQ=255)
                                                                              tctcATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCg                                                                      <  4:699088/149‑1 (MQ=255)
                                                                               ctcATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCg                                                                      >  2:361179/1‑148 (MQ=255)
                                                                                 cATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCg                                                                                 <  1:124973/135‑1 (MQ=255)
                                                                                 cATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCg                                                                                 >  2:124973/1‑135 (MQ=255)
                                                                                 cATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGc                                                                    >  1:209795/1‑148 (MQ=255)
                                                                                  ctAAAGAGTGACGTGGGGCGCAGGCGTCTCCCTCAATATCTTGTTCGTTAGAAAAAACCCGCTTAAAAAAGCGGGTTTTTTATCAGACAGATGTACGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAATTCGGCGt                                                                  >  4:50887/2‑149 (MQ=255)
                                                                                  aTAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCg                                                                   <  2:172475/148‑1 (MQ=255)
                                                                                   tAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGt                                                                  >  1:431651/1‑148 (MQ=255)
                                                                                    aaaGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTgg                                                                <  4:727505/149‑1 (MQ=255)
                                                                                    aaaGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTgg                                                                >  3:727505/1‑149 (MQ=255)
                                                                                     aaGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTgg                                                                >  2:883572/1‑148 (MQ=255)
                                                                                          tGACGTGGCGAGCAGGACGCTTCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCACCCCTTTCGCCTGCAAATCGGCGTGGTaagaa                                                          <  3:557072/149‑1 (MQ=255)
                                                                                             cGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAagaaga                                                        >  2:355853/1‑148 (MQ=255)
                                                                                               tGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGc                                                      >  1:341757/1‑148 (MQ=255)
                                                                                                 gCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGt                                                                  >  3:133494/1‑134 (MQ=255)
                                                                                                 gCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGt                                                                  <  4:133494/134‑1 (MQ=255)
                                                                                                  cGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAAttatta                                                                                                             <  1:793609/90‑1 (MQ=255)
                                                                                                  cGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAAttatta                                                                                                             >  2:793609/1‑90 (MQ=255)
                                                                                                     gCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACaa                                                <  2:282058/148‑1 (MQ=255)
                                                                                                          aCGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTg                                                                 <  1:655043/126‑1 (MQ=255)
                                                                                                          aCGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTg                                                                 >  2:655043/1‑126 (MQ=255)
                                                                                                          aCGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGa                                           >  1:199212/1‑148 (MQ=255)
                                                                                                           cGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGAc                                          <  2:199212/148‑1 (MQ=255)
                                                                                                             cTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCg                                        <  2:431651/148‑1 (MQ=255)
                                                                                                               cccTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCa                                      <  2:341757/148‑1 (MQ=255)
                                                                                                               cccTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCa                                      >  1:807800/1‑148 (MQ=255)
                                                                                                                 cTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGt                                                                  <  3:679074/118‑1 (MQ=255)
                                                                                                                 cTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGt                                                                  >  4:679074/1‑118 (MQ=255)
                                                                                                                     atatCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATgcagca                               <  3:304028/149‑1 (MQ=255)
                                                                                                                      tatCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGcc                                                                               <  4:143739/100‑1 (MQ=255)
                                                                                                                      tatCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGcc                                                                               >  3:143739/1‑100 (MQ=255)
                                                                                                                      tatCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTaag                                                            >  3:65553/1‑119 (MQ=255)
                                                                                                                      tatCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTaag                                                            <  4:65553/119‑1 (MQ=255)
                                                                                                                                gATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATGGGTAAAGcc                    >  3:349276/1‑149 (MQ=255)
                                                                                                                                 aTACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATGGGTAAAGcc                    >  1:233176/1‑148 (MQ=255)
                                                                                                                                  tACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATGGGTAAAGccc                   >  2:548880/1‑148 (MQ=255)
                                                                                                                                   aCAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATGGGTaaa                       >  2:870265/1‑143 (MQ=255)
                                                                                                                                   aCAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATGGGTaaa                       <  1:870265/143‑1 (MQ=255)
                                                                                                                                   aCAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATGGGTAAAGCCCAt                 >  3:207622/1‑149 (MQ=255)
                                                                                                                                    taaaaaaCCCGCTTTTTGAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATgg                            >  3:755833/2‑137 (MQ=255)
                                                                                                                                    taaaaaaCCCGCTTTTTGAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATgg                            <  4:755833/136‑1 (MQ=255)
                                                                                                                                      aaaaaCCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATGGGTAAAGCCCATCg               <  1:310040/148‑1 (MQ=255)
                                                                                                                                      aaaaaCCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATGGGTAAAGCCCATCGc              >  4:192875/1‑149 (MQ=255)
                                                                                                                                            ccGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATGGGTAAAGCCCATCGCCAgcg         <  2:129810/148‑1 (MQ=255)
                                                                                                                                                ttCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATGGGTAAAGCCCATCGCCAGCGCTTCg    >  4:271963/1‑149 (MQ=255)
                                                                                                                                                 tCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATGGGTAAAGCCCATCGCCAGCGCTTCg    >  2:685851/1‑148 (MQ=255)
                                                                                                                                                  cAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATGGGTAAAGCCCATCGCCAGCGCTTCGGc  <  3:271963/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GGGTTCAAGAAGGCGATGAATGATGACGATGCTGCGGCGAAAAAAGGCGCAGACGTTGATCTTCAGGCTGAAAAGCTCTCTCATAAAGAGTGACGTGGCGAGCAGGACGCTCCCTCAATATCTTGTTCGATACAAAAAACCCGCTTCAAAAAGCGGGTTTTTTATCAGACAGATGTAAGTAATTATTACAGGATTACTTAACTTCCATCCCTTTCGCCTGCAAATCGGCGTGGTAAGAAGAGCGGACAAACGGACCGCATGCAGCATGGGTAAAGCCCATCGCCAGCGCTTCGGC  >  NC_000913/659058‑659352

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: