Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,127,938 T→A 7.7% intergenic (‑44/+140) 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,9380TA 7.7% 273.7 / 11.1 78intergenic (‑44/+140)yiiX/metJputative lipid binding hydrolase, DUF830 family protein/transcriptional repressor, S‑adenosylmethionine‑binding
Reads supporting (aligned to +/- strand):  ref base T (29/43);  new base A (3/3);  total (32/46)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

TGCGATGAGCGAGAGATCTGAAAGATGATGTCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGT  >  NC_000913/4127791‑4128082
                                                                                                                                                   |                                                                                                                                                
tGCGATGAGCGAGAGATCTGAAAGATGATGTCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTg                                                                                                                                                 <  1:717797/149‑1 (MQ=255)
   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                                                                                                                 <  4:466087/149‑1 (MQ=255)
                                gCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGc                                                                                                                 <  2:255670/149‑1 (MQ=255)
                                gCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGc                                                                                                                 <  4:190203/149‑1 (MQ=255)
                                 ccGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGcc                                                                                                                <  1:674593/149‑1 (MQ=255)
                                 ccGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGcc                                                                                                                >  4:441099/1‑149 (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)
                                                                                                                                                   |                                                                                                                                                
TGCGATGAGCGAGAGATCTGAAAGATGATGTCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGCAGGCTGAGGATCAGCAGCCTGTTTTTCATAGTTAAACGTCCATGTATAAAAAGCGGTGGGTCGCAGACAACGTGCTCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGT  >  NC_000913/4127791‑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: ‑