Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,049,989 A→C 5.7% noncoding (91/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,9890AC 5.7% 502.0 / 13.9 145noncoding (91/109 nt)spfSpot 42 sRNA antisense regulator of galK translation, Hfq‑dependent
Reads supporting (aligned to +/- strand):  ref base A (76/59);  new base C (4/5);  total (80/65)
Fisher's exact test for biased strand distribution p-value = 5.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 146 total aligned reads displayed.

TGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTAT  >  NC_000913/4049840‑4050136
                                                                                                                                                     |                                                                                                                                                     
tgtaaaaaaaTTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTaa                                                                                                                                                       >  2:277724/4‑150 (MQ=255)
   aaaaaaaTTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGa                                                                                                                                                    >  1:1409754/1‑150 (MQ=255)
    aaaaaaTTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGAc                                                                                                                                                   >  1:308886/1‑150 (MQ=255)
      aaaaTTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCAGCCCCAGTCAGTAAAGACTg                                                                                                                                                 >  2:1867945/1‑150 (MQ=255)
       aaaTTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTgg                                                                                                                                                >  1:1839832/1‑150 (MQ=255)
        aaTTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCATTACTGACTggg                                                                                                                                               <  2:1567694/150‑1 (MQ=255)
         aTTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTgggg                                                                                                                                              <  2:89058/150‑1 (MQ=255)
              aaaaaGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGttt                                                                                                                                         >  1:881857/1‑150 (MQ=255)
                aaaGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCCCGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCATTACTGACTGGGGCggctaa                                                                                                                                       <  2:1082604/150‑6 (MQ=255)
                 aaGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTGAGCCGCCCCAGTCAGTAATGACTGGGGCGtttttt                                                                                                                                      <  2:1537470/150‑1 (MQ=255)
                    tGCTTTCTTAACTTAAAAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGATGTAAGATTTTCTATATTTCAGACCTTTTACTTCACTTAATCGTATTTGGCTTAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTAtt                                                                                                                                  <  1:193051/150‑1 (MQ=255)
                    tGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTTGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTAtt                                                                                                                                  <  1:895527/150‑1 (MQ=255)
                     gCTTTCAGAACTGAACAAAAAAGAGTACAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCGATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTg                                                                                                                                 <  2:1712931/150‑1 (MQ=255)
                     gCTATCTGAACTGAACAAAAAAGAGTAAACTAAGTCGCGTAGGGAACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGAATTTTATTg                                                                                                                                 <  2:550384/150‑1 (MQ=255)
                      cTTTCTGAACTGACCAAAAAAGAGGAAAGTTAGTCGCGTAGGGTACAGAGGTAAGCTGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTGGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTgg                                                                                                                                <  1:2713231/150‑1 (MQ=255)
                       tttCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTATGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTggg                                                                                                                               <  1:2478073/150‑1 (MQ=255)
                       tttCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTggg                                                                                                                               <  1:2469976/150‑1 (MQ=255)
                       tttCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTggg                                                                                                                               >  1:2259348/1‑150 (MQ=255)
                       tttCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTggg                                                                                                                               <  1:1367053/150‑1 (MQ=255)
                         tCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCg                                                                                                                             >  1:2682166/1‑150 (MQ=255)
                          ctgaactgaaCAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGCGCGTTTTTTATTGGGCga                                                                                                                            >  1:1297342/1‑150 (MQ=255)
                             aactgaacAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCgaaag                                                                                                                         >  2:2332634/1‑150 (MQ=255)
                             aactgaacAAAAAAGAGAAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCAGCCCCACACAGTAATGACTGGGGCGTTTTTTATTTTGCgaaat                                                                                                                         >  1:583238/1‑149 (MQ=255)
                                tgaacAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCgaaagaaa                                                                                                                      >  2:11575/1‑150 (MQ=255)
                                tgaacAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCgaaagaaa                                                                                                                      <  1:42139/150‑1 (MQ=255)
                                tgaacAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCgaaagaaa                                                                                                                      <  1:42164/150‑1 (MQ=255)
                                      aaaaaGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATcc                                                                                                                <  2:1880035/150‑1 (MQ=255)
                                       aaaaGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCg                                                                                                               >  2:1275582/1‑150 (MQ=255)
                                       aaaaGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCg                                                                                                               >  1:166142/1‑150 (MQ=255)
                                        aaaGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATATGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGAATGGGGCGTTTTTTAATTGGCGAAAGAAAAGATACGt                                                                                                              >  2:250534/1‑150 (MQ=255)
                                           gagTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAAt                                                                                                           <  2:2308600/150‑1 (MQ=255)
                                            agTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTAATTCACGTAATCGGATTTGTCTGAATATTTTAGCCGCCACAGTCAGTAATGACTGTGGCGTTTTTTATTGGGAGAAAGAAAAGATCAGTAAtt                                                                                                          >  2:1905968/1‑149 (MQ=255)
                                             gTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGc                                                                                                         >  2:1023412/1‑150 (MQ=255)
                                                aaGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTg                                                                                                      >  2:749398/1‑150 (MQ=255)
                                                         gcgTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTAt                                                                                             <  2:891389/150‑1 (MQ=255)
                                                          cgTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATg                                                                                            <  1:1018219/150‑1 (MQ=255)
                                                            tAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTCTTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGtt                                                                                          <  2:1814534/150‑1 (MQ=255)
                                                             aGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGACGATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGGAATGCCTGCTGCGCGATGttt                                                                                         >  2:1860496/1‑150 (MQ=255)
                                                             aGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCATTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAAAGACTGGGGCGTTTTATATAGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGttt                                                                                         >  2:1715694/1‑150 (MQ=255)
                                                                gTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATAAGTAATGCCTGATGCGCTATGTTTAta                                                                                      >  2:98165/1‑149 (MQ=255)
                                                                 tACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCa                                                                                     >  2:1371729/1‑150 (MQ=255)
                                                                   cAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAgg                                                                                   <  1:624882/150‑1 (MQ=255)
                                                                        gTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGTTACGTAATGCCTGATGCGCTATGTTTATCATGCCAAc                                                                             >  2:247947/1‑150 (MQ=255)
                                                                        gTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCCAc                                                                             >  1:2419501/1‑150 (MQ=255)
                                                                        gTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAAc                                                                             >  1:2419488/1‑150 (MQ=255)
                                                                        gTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAAc                                                                             <  2:1409754/150‑1 (MQ=255)
                                                                         tAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACg                                                                            <  1:1626250/150‑1 (MQ=255)
                                                                           aGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGt                                                                          >  1:1876904/1‑150 (MQ=255)
                                                                               gTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGaa                                                                      <  2:881719/150‑1 (MQ=255)
                                                                                 tCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAAtt                                                                    >  2:2673289/1‑150 (MQ=255)
                                                                                  cTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCTGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTg                                                                   <  2:221187/150‑1 (MQ=255)
                                                                                   tATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGt                                                                  >  2:2435810/1‑150 (MQ=255)
                                                                                   tATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATAATGCCAACGGTAGAATTGt                                                                  >  1:752227/1‑150 (MQ=255)
                                                                                   tATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATAGTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTGTTAATGTGCGAAAGAAAAGATAAGTAATTCCTGATGCCCGATGTTTATCAGGCCAACGGGATAATTGt                                                                  >  1:17834/1‑150 (MQ=255)
                                                                                    aTCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTa                                                                 >  1:1984039/1‑150 (MQ=255)
                                                                                     tCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTaa                                                                <  1:2297454/150‑1 (MQ=255)
                                                                                     tCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTaa                                                                >  2:2300460/1‑150 (MQ=255)
                                                                                      cTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAAt                                                               >  2:1161208/1‑150 (MQ=255)
                                                                                       tttCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATc                                                              <  1:2411200/150‑1 (MQ=255)
                                                                                       tttCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATGTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCGGATGCGCTATGTTTATAAGGCCAACGGTAGAATTGTAATc                                                              >  1:1886964/1‑150 (MQ=255)
                                                                                         tCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGGAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTCGAATTGTAATCTa                                                            >  1:2261704/1‑150 (MQ=255)
                                                                                           aGACCTTTTAGTTCAGGTACGGGGATTTGGGTGATTTTTTTTGCCGCCCCCGTCAGTATTGCGTGCGGGGTTTTTTTTGGGGCTAATGAATAGTTCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAAGGTTTGAATtgtattgtatg                                                          <  2:880561/150‑2 (MQ=255)
                                                                                           aGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTACTGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTAtt                                                          >  2:851726/1‑150 (MQ=255)
                                                                                           aGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTAtt                                                          >  2:1475891/1‑150 (MQ=255)
                                                                                             aCCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGa                                                        >  2:406051/1‑150 (MQ=255)
                                                                                               cTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAAt                                                      >  2:1704278/1‑150 (MQ=255)
                                                                                               cTTTTAATTCACGTAATAGGATTTGGCTGAATATGTTAGCCGCCCCAGTCAGTAATTACTGGGGCGTTTTTTATTGGGCGAAAGACAAGATCCGTAATGCCTGATGCGCGATGGTTATCAGTCCAACGGTAGAATTTTAATCTATTGAAt                                                      >  2:1742973/1‑150 (MQ=255)
                                                                                                ttttACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAAtt                                                     >  1:78047/1‑150 (MQ=255)
                                                                                                    aCTTCACGTAATCTGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACg                                                 <  2:493769/150‑1 (MQ=255)
                                                                                                    aCTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACg                                                 <  2:493841/150‑1 (MQ=255)
                                                                                                     cTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAACAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATAGTAATCTATTGAAATTACgg                                                >  2:603792/1‑150 (MQ=255)
                                                                                                      ttCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACggg                                               >  2:2461920/1‑150 (MQ=255)
                                                                                                      ttCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACggg                                               >  2:2461860/1‑150 (MQ=255)
                                                                                                      ttCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACggg                                               >  2:2430382/1‑150 (MQ=255)
                                                                                                       tCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCACTCACTAATGACTGGGGCGAATTTTATTGGGCGAAAGAAAAGATCCCTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGc                                              <  1:1014544/150‑1 (MQ=255)
                                                                                                         aCGTACTCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTACTGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCg                                            <  1:2052395/150‑1 (MQ=255)
                                                                                                          cGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCgg                                           <  1:809049/150‑1 (MQ=255)
                                                                                                            tAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGACAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGAt                                         >  2:1878156/1‑150 (MQ=255)
                                                                                                              aTCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTCTTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTGTATCAGGCCAACGGTAGAATTGTAATCTAGTGAAGTTACGGGCCGGAGAc                                       >  2:1699237/1‑150 (MQ=255)
                                                                                                               tCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACg                                      >  2:1766134/1‑150 (MQ=255)
                                                                                                                cGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGc                                     <  1:2332634/150‑1 (MQ=255)
                                                                                                                 ggATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTTTTTGGGGAAAGAAAAGATACCTAATTCCTGAAGCGCAAAGTTTATAAGGCCCAACGTAGAATTTTTATATATTTTAATTAAGGGCCGGATaaccc                                    >  1:2082698/1‑146 (MQ=255)
                                                                                                                       ggCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATcc                              >  2:1457198/1‑150 (MQ=255)
                                                                                                                        gCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCg                             <  1:767175/150‑1 (MQ=255)
                                                                                                                        gCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCg                             <  1:767283/150‑1 (MQ=255)
                                                                                                                         cTGAATATTTTAGCCGCCCCAGTCAGTTATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCgg                            <  1:1636744/150‑1 (MQ=255)
                                                                                                                              tatTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACaa                       <  1:1194876/150‑1 (MQ=255)
                                                                                                                              tatTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTGACGGGCCGGATACGCCACATCCGGCACaa                       <  2:1936404/150‑1 (MQ=255)
                                                                                                                              tatTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAAGTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACaa                       >  2:1485651/1‑150 (MQ=255)
                                                                                                                                   tAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCAtt                  <  1:2149378/150‑1 (MQ=255)
                                                                                                                                     gccgccCCAGTCAGTAATGACTGGGGCGTTTTTTATTTGGCGGAAGAAAACATCCGTAATTCCTGATTCGCGATGTTTATAAGGCCAAACGGAGGATTTTAAACTAATGTATTTAAGGGCCGGATACGCCCCAAACGGCACAACCAATaa                >  1:1928647/1‑150 (MQ=255)
                                                                                                                                     gccgccCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTaa                >  2:378204/1‑150 (MQ=255)
                                                                                                                                      ccgccCCAGTCAGTAATGACTGGGGCGTTTTTTAATGGGGCCAAAGAAAGAAGCCGAAATCCATGAGCCCCGATGTTTACAGGCCAAACGGACGAATGTAAAACAATTGAATTAACGGCCCGCAAACCCCCCACCCGCACAACCAATAAg               >  2:2276850/1‑150 (MQ=255)
                                                                                                                                      ccgccCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAg               <  1:749398/150‑1 (MQ=255)
                                                                                                                                         ccccAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCa            >  2:659995/1‑150 (MQ=255)
                                                                                                                                         ccccAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACTGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCa            <  1:744049/150‑1 (MQ=255)
                                                                                                                                         ccccAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCa            >  2:1576565/1‑150 (MQ=255)
                                                                                                                                               tcagtAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGaaaa      <  1:283556/150‑1 (MQ=255)
                                                                                                                                                   tAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAattat  <  1:1275582/150‑1 (MQ=255)
                                                                                                                                                   tAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATAAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAattat  >  2:977107/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTAT  >  NC_000913/4049840‑4050136

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: