Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,127,941 T→A 7.8% intergenic (‑47/+137) yiiX ← / ← metJ putative lipid binding hydrolase, DUF830 family protein/transcriptional repressor, S‑adenosylmethionine‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,127,9410TA 7.8% 261.0 / 11.6 77intergenic (‑47/+137)yiiX/metJputative lipid binding hydrolase, DUF830 family protein/transcriptional repressor, S‑adenosylmethionine‑binding
Reads supporting (aligned to +/- strand):  ref base T (29/42);  new base A (3/3);  total (32/45)
Fisher's exact test for biased strand distribution p-value = 6.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.54e-01

GATGAGCGAGAGATCTGAAAGATGATGTCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGT  >  NC_000913/4127794‑4128082
                                                                                                                                                   |                                                                                                                                             
gATGAGCGAGAGATCTGAAAGATGATGTCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTc                                                                                                                                              >  3:585724/1‑149 (MQ=255)
 aTGAGCGAGAGATCTGAAAGATGATGTCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCg                                                                                                                                             >  2:133065/1‑149 (MQ=255)
     gCGAGAGATCTGAAAGATGATGTCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCgttgt                                                                                                                                         <  2:428878/149‑1 (MQ=255)
             tCTGAAAGATGATGTCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCgtt                                                                                                                                           <  2:15130/139‑1 (MQ=255)
             tCTGAAAGATGATGTCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCgtt                                                                                                                                           >  1:15130/1‑139 (MQ=255)
               tGAAAGATGATGTCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGa                                                                                                                               >  3:628603/1‑149 (MQ=255)
                gAAAGATGATGTCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGAt                                                                                                                              >  3:572691/1‑149 (MQ=255)
                 aaaGATGATGTCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGAta                                                                                                                             >  4:382029/1‑148 (MQ=255)
                  aaGATGATGTCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATgc                                                                                                                            >  2:446118/1‑149 (MQ=255)
                     atgatgTCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATgcggc                                                                                                                         <  1:576523/149‑1 (MQ=255)
                     atgatgTCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATgcggc                                                                                                                         <  3:392217/149‑1 (MQ=255)
                          gTCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATgcggc                                                                                                                         >  4:224596/1‑144 (MQ=255)
                          gTCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATgcggc                                                                                                                         <  3:224596/144‑1 (MQ=255)
                          gTCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGaa                                                                                                                    <  1:195095/149‑1 (MQ=255)
                           tCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAAc                                                                                                                   <  1:446118/149‑1 (MQ=255)
                             gCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGc                                                                                                                 <  2:255670/149‑1 (MQ=255)
                             gCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGc                                                                                                                 <  4:466087/149‑1 (MQ=255)
                             gCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGc                                                                                                                 <  4:190203/149‑1 (MQ=255)
                              ccGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGcc                                                                                                                >  4:441099/1‑149 (MQ=255)
                              ccGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGcc                                                                                                                <  1:674593/149‑1 (MQ=255)
                              ccGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGcc                                                                                                                <  4:572691/149‑1 (MQ=255)
                              ccGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGcc                                                                                                                <  3:441099/149‑1 (MQ=255)
                                    tgtgGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATcc                                                                                                          <  1:133065/149‑1 (MQ=255)
                                        gCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATAAGGCGTTCACGCCGCATCCGGca                                                                                                      <  3:382029/149‑2 (MQ=255)
                                             cAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACaa                                                                                                 <  4:218445/149‑1 (MQ=255)
                                             cAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACaa                                                                                                 <  2:350215/149‑1 (MQ=255)
                                              aGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAg                                                                                                <  4:3415/149‑1 (MQ=255)
                                               ggCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATgcggc                                                                                                                         <  3:441392/123‑1 (MQ=255)
                                               ggCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATgcggc                                                                                                                         >  4:441392/1‑123 (MQ=255)
                                                 cAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTc                                                                                             <  1:288631/149‑1 (MQ=255)
                                                 cAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTc                                                                                             <  2:147608/149‑1 (MQ=255)
                                                   aaGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGt                                                                                           <  2:329847/149‑1 (MQ=255)
                                                        aGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCaaa                                                                                      >  2:44816/1‑149 (MQ=255)
                                                         ggTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAAt                                                                                     <  1:598014/149‑1 (MQ=255)
                                                          gTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAAtt                                                                                    >  2:317029/1‑149 (MQ=255)
                                                          gTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAAtt                                                                                    >  4:617486/1‑149 (MQ=255)
                                                          gTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAAtt                                                                                    <  3:617486/149‑1 (MQ=255)
                                                             cAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCaa                                                                                 >  1:62795/1‑149 (MQ=255)
                                                                  aCCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAAt                                                                            >  2:185319/1‑149 (MQ=255)
                                                                    cagcagGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTg                                                                          >  3:700218/1‑149 (MQ=255)
                                                                       cagGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGCCCTACATGTTCGTGCAAATTCAATAATTTGCaa                                                                       <  4:362739/149‑1 (MQ=255)
                                                                        agGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAat                                                                      <  1:44816/149‑1 (MQ=255)
                                                                          gCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCaaa                                                                                      >  4:636924/1‑131 (MQ=255)
                                                                          gCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCaaa                                                                                      <  3:636924/131‑1 (MQ=255)
                                                                                gATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCgttg                                                                                                                                          <  2:574801/73‑1 (MQ=255)
                                                                                gATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCgttg                                                                                                                                          >  1:574801/1‑73 (MQ=255)
                                                                                gATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCaaa                                                                                      <  2:645409/125‑1 (MQ=255)
                                                                                gATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCaaa                                                                                      >  1:645409/1‑125 (MQ=255)
                                                                                       agcCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGa                                                       <  4:628603/149‑1 (MQ=255)
                                                                                       agcCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGa                                                       <  4:700218/149‑1 (MQ=255)
                                                                                         ccTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATa                                                     >  4:398956/1‑149 (MQ=255)
                                                                                           tGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAg                                                   <  1:317029/149‑1 (MQ=255)
                                                                                             tttttCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCg                                                 >  1:401537/1‑149 (MQ=255)
                                                                                               tttCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTa                                               <  4:585724/149‑1 (MQ=255)
                                                                                                     aGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGa                                                       >  3:656660/1‑135 (MQ=255)
                                                                                                     aGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGa                                                       <  4:656660/135‑1 (MQ=255)
                                                                                                            cGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGa                                  <  2:401537/149‑1 (MQ=255)
                                                                                                               ccATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTc                               <  3:784408/149‑1 (MQ=255)
                                                                                                               ccATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTc                               >  4:784408/1‑149 (MQ=255)
                                                                                                                       aaaaaGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATCGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCcgc                       >  4:739456/1‑149 (MQ=255)
                                                                                                                         aaaGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCgcgc                     >  4:788228/1‑149 (MQ=255)
                                                                                                                          aaGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGcgcc                    <  3:788228/149‑1 (MQ=255)
                                                                                                                                   gTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCttttt           <  2:62795/149‑1 (MQ=255)
                                                                                                                                   gTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCttttt           >  1:83046/1‑149 (MQ=255)
                                                                                                                                    gcataaaCAACGAGCACGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCcgc                       <  2:608612/130‑1 (MQ=255)
                                                                                                                                    gcataaaCAACGAGCACGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCcgc                       <  2:275682/130‑1 (MQ=255)
                                                                                                                                    gcataaaCAACGAGCACGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCcgc                       <  2:102573/130‑1 (MQ=255)
                                                                                                                                    gcataaaCAACGAGCACGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCcgc                       >  1:608612/7‑136 (MQ=21)
                                                                                                                                    gcataaaCAACGAGCACGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCcgc                       >  1:275682/7‑136 (MQ=255)
                                                                                                                                    gcataaaCAACGAGCACGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCcgc                       >  1:102573/7‑136 (MQ=21)
                                                                                                                                         gACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCttt     >  3:55059/1‑149 (MQ=255)
                                                                                                                                          aCAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTa    <  2:83046/149‑1 (MQ=255)
                                                                                                                                           cAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCt       <  2:489057/145‑1 (MQ=255)
                                                                                                                                           cAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCt       >  1:489057/1‑145 (MQ=255)
                                                                                                                                            aaCGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCt       <  2:244866/144‑1 (MQ=255)
                                                                                                                                            aaCGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCt       >  1:244866/1‑144 (MQ=255)
                                                                                                                                            aaCGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGt  <  3:398956/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
GATGAGCGAGAGATCTGAAAGATGATGTCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGT  >  NC_000913/4127794‑4128082

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: