Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,049,986 G→T 10.0% noncoding (88/109 nt) spf → Spot 42 sRNA antisense regulator of galK translation, Hfq‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,049,9860GT10.0% 466.9 / 32.4 141noncoding (88/109 nt)spfSpot 42 sRNA antisense regulator of galK translation, Hfq‑dependent
Reads supporting (aligned to +/- strand):  ref base G (66/60);  new base T (6/8);  total (72/69)
Fisher's exact test for biased strand distribution p-value = 5.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

Only 100 of 146 total aligned reads displayed.

TGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATT  >  NC_000913/4049840‑4050134
                                                                                                                                                  |                                                                                                                                                     
tGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTa                                                                                                                                                     >  2:70527/1‑149 (MQ=255)
   aaaaaaaTTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATg                                                                                                                                                  >  2:934393/1‑149 (MQ=255)
    aaaaaaTTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGa                                                                                                                                                 >  1:878239/1‑149 (MQ=255)
        aaTTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCATTACTGACTgg                                                                                                                                             <  4:327015/149‑1 (MQ=255)
        aaTTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCATTACTGACTgg                                                                                                                                             <  4:52418/149‑1 (MQ=255)
         aTTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCATTACTGACTggg                                                                                                                                            <  4:936831/149‑1 (MQ=255)
          ttACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTgggg                                                                                                                                           <  4:178013/149‑1 (MQ=255)
           tACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCATTACTGACTGGGGc                                                                                                                                          <  1:843275/149‑1 (MQ=255)
           tACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAAATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCATTACTGACTGGGGc                                                                                                                                          <  1:309830/149‑1 (MQ=255)
              aaaaaGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGtt                                                                                                                                       >  3:519448/1‑149 (MQ=255)
              aaaaaGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTAAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGtt                                                                                                                                       >  4:216923/1‑149 (MQ=255)
               aaaaGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGttt                                                                                                                                      <  1:604478/149‑1 (MQ=255)
                  aGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGtttttt                                                                                                                                   <  1:589471/149‑1 (MQ=255)
                  aGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGtttttt                                                                                                                                   <  2:1066891/149‑1 (MQ=255)
                     gCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTAtt                                                                                                                                <  1:70527/149‑1 (MQ=255)
                     gCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTAtt                                                                                                                                <  3:176410/149‑1 (MQ=255)
                      cTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTg                                                                                                                               <  1:1002210/149‑1 (MQ=255)
                       tttCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTgg                                                                                                                              <  2:923254/149‑1 (MQ=255)
                           tgaactgaacAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCga                                                                                                                          >  3:302235/1‑149 (MQ=255)
                             aactgaacAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCgaaa                                                                                                                        <  1:659897/149‑1 (MQ=255)
                               ctgaacAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCgaaaga                                                                                                                      <  4:114975/149‑1 (MQ=255)
                                tgaacAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCgaaagaa                                                                                                                     >  3:782082/1‑149 (MQ=255)
                                 gaacAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCgaaagaaa                                                                                                                    <  4:685607/149‑1 (MQ=255)
                                  aacAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGaaagaaaa                                                                                                                   <  1:1023723/149‑1 (MQ=255)
                                  aacAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGaaagaaaa                                                                                                                   >  1:163229/1‑149 (MQ=255)
                                  aacAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGaaagaaaa                                                                                                                   <  4:191370/149‑1 (MQ=255)
                                    cAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGaaagaaaaga                                                                                                                 >  1:148269/1‑149 (MQ=255)
                                     aaaaaaGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGTTTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGAGTGTTTTTTTTTAGGCGAATGAACAGAt                                                                                                                <  2:300098/149‑1 (MQ=255)
                                     aaaaaaGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGAt                                                                                                                >  3:487789/1‑149 (MQ=255)
                                     aaaaaaGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGAt                                                                                                                <  4:454874/149‑1 (MQ=255)
                                     aaaaaaGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGAt                                                                                                                <  2:271447/149‑1 (MQ=255)
                                       aaaaGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATcc                                                                                                              >  4:426428/1‑149 (MQ=255)
                                       aaaaGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATcc                                                                                                              <  4:931100/149‑1 (MQ=255)
                                              tAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGc                                                                                                       >  2:817017/1‑149 (MQ=255)
                                               aaaGTTAGTCGCGTAGGGTACAGAGGTAGGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGcc                                                                                                      <  1:1020465/149‑1 (MQ=255)
                                                          cgTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTAt                                                                                           <  4:963570/149‑1 (MQ=255)
                                                           gTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCATTACTGACTgggg                                                                                                                                           <  2:999185/100‑1 (MQ=255)
                                                               ggTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTAt                                                                                     <  4:332186/149‑1 (MQ=255)
                                                                gTACAGAGGTAAGATTTTCTATCTTTCAGACCTTTTACATCTCGTAATCTGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGATTGGGGCGTTTTTTATTGGGCGTAAGAAATGATCCGTAATGCCTGATGCGCTATGTTTAt                                                                                     <  4:97016/149‑1 (MQ=255)
                                                                  aCAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATc                                                                                                               <  3:817871/121‑1 (MQ=255)
                                                                  aCAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCa                                                                                   <  1:817017/149‑1 (MQ=255)
                                                                  aCAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCa                                                                                   <  4:488808/149‑1 (MQ=255)
                                                                       ggTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAg                                                                                  >  4:108091/1‑145 (MQ=255)
                                                                       ggTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAg                                                                                  <  3:108091/145‑1 (MQ=255)
                                                                        gTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCaa                                                                             >  2:428825/1‑149 (MQ=255)
                                                                        gTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCaa                                                                             <  2:115774/149‑1 (MQ=255)
                                                                          aaGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACg                                                                           <  3:547172/149‑1 (MQ=255)
                                                                           aGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACgg                                                                          >  2:1055191/1‑149 (MQ=255)
                                                                           aGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACgg                                                                          >  1:943336/1‑149 (MQ=255)
                                                                                ttCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGGAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGaa                                                                     >  2:126509/1‑149 (MQ=255)
                                                                                ttCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGaa                                                                     >  3:798261/1‑149 (MQ=255)
                                                                                ttCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTAGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGaa                                                                     >  3:652259/1‑149 (MQ=255)
                                                                                  cTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAAtt                                                                   <  1:897668/149‑1 (MQ=255)
                                                                                   tATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTg                                                                  >  1:826743/1‑149 (MQ=255)
                                                                                       tttCAGACCTTTTACTTCACGTAATCGGGCTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGAGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCGATGTTTATCAGGCCAACGGAAGAATTGTAAt                                                             >  2:112690/1‑149 (MQ=255)
                                                                                       tttCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAAt                                                             <  2:975977/149‑1 (MQ=255)
                                                                                       tttCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAAt                                                             >  1:881507/1‑149 (MQ=255)
                                                                                       tttCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGACCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAAt                                                             <  3:370489/149‑1 (MQ=255)
                                                                                             aCCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCa                                                                                   >  3:166750/1‑122 (MQ=255)
                                                                                             aCCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCa                                                                                   <  4:166750/122‑1 (MQ=255)
                                                                                               cTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTgg                                                                                                                              >  3:351546/1‑77 (MQ=255)
                                                                                               cTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGaa                                                     >  3:474289/1‑149 (MQ=255)
                                                                                                ttttACTTCACGTAATCGGATTTGGCTGAATATTTTAGTCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACTGTAGAATTGTAATCTATTGAAt                                                    >  4:862960/1‑149 (MQ=255)
                                                                                                ttttACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAg                                                                                  <  1:167110/120‑1 (MQ=255)
                                                                                                ttttACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAg                                                                                  >  2:167110/1‑120 (MQ=255)
                                                                                                ttttACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAAt                                                    >  1:76270/1‑149 (MQ=255)
                                                                                                      ttCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACgg                                              >  4:386485/1‑149 (MQ=255)
                                                                                                      ttCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACggg                                             >  3:825170/1‑149 (MQ=255)
                                                                                                        cACGTAATCGGTTTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGc                                            >  4:410104/1‑149 (MQ=255)
                                                                                                             aaTCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGAt                                       <  1:934393/149‑1 (MQ=255)
                                                                                                              aTCGGATTTGGCTGAATATTTTAGCCGCCCCAACAAGTAATGAGTGGGGCGTTTTTAATAGGCCGAAAGAAAGGTTTCGAAATGCTTAATCTGCTATGTTAATCAGTCAAACGTTAGAATGTAATTCTATGGATTTAACGGCCCGGATa                                      >  4:1008909/1‑149 (MQ=255)
                                                                                                               tCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATGCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATAc                                     >  4:688583/1‑149 (MQ=255)
                                                                                                                        gCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATcc                            >  2:149644/1‑149 (MQ=255)
                                                                                                                          tGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCgg                          <  1:916083/149‑1 (MQ=255)
                                                                                                                          tGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCgg                          >  1:687364/1‑149 (MQ=255)
                                                                                                                          tGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTCTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCgg                          <  3:1035028/149‑1 (MQ=255)
                                                                                                                               atTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACaa                     >  4:788854/1‑149 (MQ=255)
                                                                                                                               atTATAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGCAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACaa                     >  4:894398/1‑149 (MQ=255)
                                                                                                                                   tAGCCGCCCCAGTCAGTAATGACTAGGGCGTTTTTTATTGGGCGAAAGAAAAAATCTGTAATGCCTGCTGCGCTATGTTTATCAGGCAAACGGTAGAATTGTAATCTATTGAATTTACGAGCCGGATACGCCACATCCGGCACAAGCAt                 <  1:959598/149‑1 (MQ=255)
                                                                                                                                     gccgccCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTa               >  1:49217/1‑149 (MQ=255)
                                                                                                                                      ccgccCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTaa              <  2:826743/149‑1 (MQ=255)
                                                                                                                                      ccgccCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTaa              <  3:237742/149‑1 (MQ=255)
                                                                                                                                       cgccCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAg             <  1:966343/149‑1 (MQ=255)
                                                                                                                                       cgccCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAg             <  1:654915/149‑1 (MQ=255)
                                                                                                                                       cgccCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAg             <  3:203113/149‑1 (MQ=255)
                                                                                                                                       cgccCCAGTCAGTAATAACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAg             <  3:894398/149‑1 (MQ=255)
                                                                                                                                        gccCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCTGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAgg            >  4:566663/1‑149 (MQ=255)
                                                                                                                                        gccCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAgg            <  2:380145/149‑1 (MQ=255)
                                                                                                                                          cccAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTaa              >  1:1058698/1‑145 (MQ=255)
                                                                                                                                          cccAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTaa              <  2:1058698/145‑1 (MQ=255)
                                                                                                                                          cccAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCa          >  3:947080/1‑149 (MQ=255)
                                                                                                                                          cccAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGca                        <  2:537045/135‑1 (MQ=255)
                                                                                                                                          cccAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCa          >  4:636170/1‑149 (MQ=255)
                                                                                                                                          cccAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCa          >  3:229053/1‑149 (MQ=255)
                                                                                                                                           ccAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCaa         >  4:165244/1‑149 (MQ=255)
                                                                                                                                           ccAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCaa         >  2:755109/1‑149 (MQ=255)
                                                                                                                                                cagtAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGaaaa    <  4:487789/149‑1 (MQ=255)
                                                                                                                                                cagtAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGaaaa    <  3:426428/149‑1 (MQ=255)
                                                                                                                                                 agtAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAat   <  4:1033915/149‑1 (MQ=255)
                                                                                                                                                  gtAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAatt  >  4:809198/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATT  >  NC_000913/4049840‑4050134

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