Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,049,986 G→T 6.1% 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,9860GT 6.1% 251.4 / 3.8 74noncoding (88/109 nt)spfSpot 42 sRNA antisense regulator of galK translation, Hfq‑dependent
Reads supporting (aligned to +/- strand):  ref base G (31/38);  new base T (2/3);  total (33/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.68e-01

GTGGGGAGTGGCGAAAACTGGGATCAGGCGCACTAAGATTCGCCTGAACATGCCTTTTTTCGTAAGTAAGCAACATAAGCTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGC  >  NC_000913/4049717‑4050247
                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                     
gTGGGGTGTGGCGAAAACTGGGATCAGGCGCACTAAGATTCGCCTGAACATGCCTTTTTTCGTAAGTAATCAACATAAGCTGTCACGTTTTGTGATAGCTATTAGAAATTCCTATGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTCTAAGACCTGTTACTTCACGTAATCGGATTTGGCTGATTCTTTTTGCCGCACCAGTCAGTAATg                                                                                                                                                                                                                                                                  >  2:462595/1‑275 (MQ=255)
      aGTGGCGAAAACTGGGATCAGGCGCACTAAGATTCGCCTGAACATGCCTTTTTTCGTAAGTAAGCAACATAAGCTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTggg                                                                                                                                                                                                                                                            <  1:104555/275‑1 (MQ=255)
       gTGGCGAAAACTGGGATCAGGCGCACTAAGATTCGCCTGAACATGCCTTTTTTCGTAAGTAAGCAACATAAGCTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTggg                                                                                                                                                                                                                                                            <  4:460906/274‑1 (MQ=255)
       gTGGCGAAAACTGGGATCAGGCGCACTAAGATTCGCCTGAACATGCCTTTTTTCGTAAGTAAGCAACATAAGCTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTggg                                                                                                                                                                                                                                                            <  4:386106/274‑1 (MQ=255)
                          ggCGCACTAAGATTCGCCTGAACATGCCTTTTTTCGTAAGTAAGCAACATAAGCTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGGGTTTTTTATTGGGCgaaa                                                                                                                                                                                                                                        >  4:353470/1‑275 (MQ=255)
                                cTAAGATCCGCCTGAACATGCCTTTTTTCGTAAGTAAGCAACATAAGCTGTCACGTTTTCTGCTGGCTCTTAGGAATTCCTATGCAACAAATGAAAAACATTTACAAAAAGTGCCTTCTGAACTGAACACAAACGAGCAAAGATAGTCGCGCAGGGTAACGAGGGAAGATGTTCTATCTTTCCGACCTTTTACGTAACGTAATCGGATTTGGCTGACTTTTTTTGCCCCCCCATTCAGTAATGATGGGGGCGTTTTTAATGTGCCGaatgacaag                                                                                                                                                                                                                                  >  4:498565/1‑275 (MQ=255)
                                        cGCCTGAACATGCCTTTTTTCGTAAGTAAGCAACATAAGCTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTaa                                                                                                                                                                                                                          >  2:314530/1‑275 (MQ=255)
                                              aaCATGCCTTTTTTCGTAAGTAAGCAACATAAGCTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTg                                                                                                                                                                                                                    <  2:515080/275‑1 (MQ=255)
                                                    ccTTTTTTCGTAAGTAAGCAACATAAGCTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATgcg                                                                                                                                                                                                               <  4:412568/274‑1 (MQ=255)
                                                        ttttCGTAAGTAAGCAACATAAGCTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATg                                                                                                                                                                                                          >  2:441790/1‑275 (MQ=255)
                                                             gtaagtaaGCAACATAAGCTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTAt                                                                                                                                                                                                     <  4:314529/275‑1 (MQ=255)
                                                               aagtaagCAACATAAGCTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCCCGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCa                                                                                                                                                                                                   >  4:452117/1‑275 (MQ=255)
                                                                         cATAAGCAGTCACGTTTTGTGATCGCTATTAGAAATTCCTATGCAACAACTGAAAAAAAATTCCCAAAAGTCCTTTCTGAACTGACCAAAAAACAGTAAAGTTAGTCGCGTAGGGTACAGAGTTAAGATGTTCTATCTTTCATACCTTTTACTTCACCTAATCGGATTTCGCTGAATATTTTATCCGCCACAGTCAGTAATGACTGGGCCGTTTTTCATTGGGCGAATGAACTGACCGGTCATGCTTGATGCCCTATGTTTATCAGGCCAAcagg                                                                                                                                                                                         >  4:182179/1‑272 (MQ=255)
                                                                                  tCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTa                                                                                                                                                                                <  1:8795/275‑1 (MQ=255)
                                                                                  tCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTa                                                                                                                                                                                <  2:428895/275‑1 (MQ=255)
                                                                                   cACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTCGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTaa                                                                                                                                                                               >  2:149489/1‑275 (MQ=255)
                                                                                      gTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATc                                                                                                                                                                             >  2:401253/1‑274 (MQ=255)
                                                                                         ttGTGATGGCTATTAGAAATTCCTATGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGTAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACAGTAGAATTGTAATCTAtt                                                                                                                                                                         >  2:277463/1‑275 (MQ=255)
                                                                                          tgtgATGGCTATTAGAAATTCCTATGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTg                                                                                                                                                                        >  2:119712/1‑275 (MQ=255)
                                                                                                     ttAGAAATTCCTATGAAAAAAATGAAAAAAAAATACAAAAAGGGGTTTTTAAACTGAACAAAAAAGAGTAAAGTTAGTCGCTTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGGAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGc                                                                                                                                                             <  4:26185/275‑1 (MQ=255)
                                                                                                           aTTCCTATGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCGGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGAGAAAGAAAAGATCCGTAAAGCATGAGGCGCTATGTTTATGAGGCAAAGGGTAGATTTGAAAACATTGGAATTTAGGGGCCGGaa                                                                                                                                                        >  4:209823/1‑273 (MQ=255)
                                                                                                               cTAGGCAAAAACTTAAAAAAAATTAAAAAAAAGGGTTTTAGAACTGAACAAAAAAGGGGAAAGGTAGTCGCGAAGGAGACAGAGGTTAGATTTTCTAGCATTGAGACCTGTTGCTTCACGTATTCGGATTTGGCTGAAAATTTTAGCCGCCAGAGTCAGTAATGTCTGGGGCATTTTTTGTTGGGCGAAAGAAAAGATCCGTAATGCCTGGTACTCTATGTTTATCAGGCCAACGGTAGAATGGTAATCTGTTGAATTTACGGGACGGATACGcc                                                                                                                                                   <  4:462594/275‑1 (MQ=255)
                                                                                                                tATGCAAAAAATGAAAAAAAATTAAAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGcc                                                                                                                                                   <  4:441789/274‑1 (MQ=255)
                                                                                                                   gCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATa                                                                                                                                                       >  4:395689/1‑267 (MQ=255)
                                                                                                                   gCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACAt                                                                                                                                               <  2:353471/275‑1 (MQ=255)
                                                                                                                       caacTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGAc                                                                                                                                                                                                                                                                >  4:555486/1‑158 (MQ=255)
                                                                                                                       caacTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGcc                                                                                                                                                   <  2:395690/267‑1 (MQ=255)
                                                                                                                          atGAAAAAAAATAACAAAAAGTGCGGTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATTTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCATTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGTCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGcac                                                                                                                                        <  4:413510/274‑1 (MQ=255)
                                                                                                                           tGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCTTAGGGTACAGATGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGcaca                                                                                                                                       <  2:209824/275‑1 (MQ=255)
                                                                                                                           tGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTggg                                                                                                                                                                                                                                                            <  2:555487/158‑1 (MQ=255)
                                                                                                                                                          ctgaacAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTcc                                                                                                        <  2:315252/275‑1 (MQ=255)
                                                                                                                                                             aacAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTg                                                                                                                                                                                                                                                              >  2:258281/1‑122 (MQ=255)
                                                                                                                                                             aacAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGa                                                                                                                                                                                                                   >  3:160726/1‑165 (MQ=255)
                                                                                                                                                             aacAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAAt                                                                                                                                                                                    >  4:374503/1‑196 (MQ=255)
                                                                                                                                                               cAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGt                                                                                                   <  2:452118/275‑1 (MQ=255)
                                                                                                                                                                 aaaaaGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGc                                                                                                                                                                                                                                                          <  4:258280/122‑1 (MQ=255)
                                                                                                                                                                 aaaaaGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATgcg                                                                                                                                                                                                               <  1:160726/165‑1 (MQ=255)
                                                                                                                                                                 aaaaaGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTa                                                                                                                                                                                <  2:374504/196‑1 (MQ=255)
                                                                                                                                                                    aaGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACt                                                                                                                                                                                                                                                               <  6:186421/114‑1 (MQ=255)
                                                                                                                                                                    aaGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACt                                                                                                                                                                                                                                                               >  5:186421/1‑114 (MQ=255)
                                                                                                                                                                       agTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTa                                                                                           <  2:417942/275‑1 (MQ=255)
                                                                                                                                                                                 gTCGCTTAGGGTACAGAGGTAAGATGGTCTATCTTTAAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCAt                                                                                 <  2:133199/275‑1 (MQ=255)
                                                                                                                                                                                       tAGGGTACAAGGGTAAGATGTTCTAACTTTCAGACCTTTTATTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCt                                                                           <  4:119711/275‑1 (MQ=255)
                                                                                                                                                                                        aGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTa                                                                                                                                                                                                                                                                     <  6:451784/88‑1 (MQ=255)
                                                                                                                                                                                        aGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTa                                                                                                                                                                                                                                                                     >  5:451784/1‑88 (MQ=255)
                                                                                                                                                                                            aaaggagGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAAc                                                                       <  2:119993/270‑1 (MQ=255)
                                                                                                                                                                                                 agGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGt                                                                  >  4:56640/1‑274 (MQ=255)
                                                                                                                                                                                                    tAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCATTactgactg                                                                                                                                                                                                                                                              <  8:574273/83‑1 (MQ=255)
                                                                                                                                                                                                    tAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCATTactgactg                                                                                                                                                                                                                                                              >  7:574273/1‑83 (MQ=255)
                                                                                                                                                                                                    tAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCt                                                                                     >  4:450463/1‑252 (MQ=255)
                                                                                                                                                                                                    tAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATc                                                               >  2:45111/1‑274 (MQ=255)
                                                                                                                                                                                                      aGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCa                                                                                                                                   >  4:291965/1‑204 (MQ=255)
                                                                                                                                                                                                        aTGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCAt                                                                                 <  2:450464/252‑1 (MQ=255)
                                                                                                                                                                                                          gTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTaa                                                                                                                               <  2:291966/204‑1 (MQ=255)
                                                                                                                                                                                                           ttCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGACGTCCTGCGTTTCTTCTCCCGGCGGCATACCGCTAAACCAGGTATCCAgctgtc                                                       >  4:119992/1‑270 (MQ=255)
                                                                                                                                                                                                             cTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTttct                                                                                               >  2:120044/1‑233 (MQ=255)
                                                                                                                                                                                                             cTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGc                                                     <  4:45110/275‑1 (MQ=255)
                                                                                                                                                                                                             cTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGc                                                     <  2:510606/275‑1 (MQ=255)
                                                                                                                                                                                                               aTCTTTCAGACCTTTTACTTCACGTAATCGGAGTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTGATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGt                                                                                                   <  2:283020/227‑1 (MQ=255)
                                                                                                                                                                                                                 cTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATgcg                                                                                                                                                                                                               <  6:624632/117‑1 (MQ=255)
                                                                                                                                                                                                                 cTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATgcg                                                                                                                                                                                                               >  5:624632/1‑117 (MQ=255)
                                                                                                                                                                                                                 cTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTa                                                                                           <  4:120043/233‑1 (MQ=255)
                                                                                                                                                                                                                            tttACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGt                                                                                                                                                                                         >  3:75342/1‑128 (MQ=255)
                                                                                                                                                                                                                             ttACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCt                                                                                     >  2:128923/1‑227 (MQ=255)
                                                                                                                                                                                                                                cTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGaa                                                                                                                                                                                     <  1:75342/128‑1 (MQ=255)
                                                                                                                                                                                                                                 ttCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCAt                                                                                 <  4:128922/227‑1 (MQ=255)
                                                                                                                                                                                                                                    aCGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGttt                              >  2:276591/1‑275 (MQ=255)
                                                                                                                                                                                                                                     cGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCATTactgactg                                                                                                                                                                                                                                                              >  8:1009921/1‑50 (MQ=37)
                                                                                                                                                                                                                                     cGTAATCGGATTTGGCTGAATATGTTAGCCGCCCCAGTCATTactgactg                                                                                                                                                                                                                                                              <  7:1009921/50‑1 (MQ=25)
                                                                                                                                                                                                                                                           ttttAGCCGCCACAGTCAGTAATCACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCa        <  4:149488/274‑1 (MQ=255)
                                                                                                                                                                                                                                                                 ccgccCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACTAAGAAACCGTTTCAACCTGc  >  2:472668/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                 ccgacCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAGATGATCCGTAATGCCTGATGCGCTATGTTTCGCAAGCCAACCGTAGAATTGTAATCTATTGAATTTACGGGCCGGGTACGCCACATCCGGCACAAGCATTa                                                                                                                                <  6:299024/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                      ccAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAAt                                                                                                                                                                                    >  5:666846/1‑91 (MQ=255)
                                                                                                                                                                                                                                                                      ccAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAAt                                                                                                                                                                                    <  6:666846/91‑1 (MQ=255)
                                                                                                                                                                                                                                                                       cagtcattACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTg                                                                                                                                                                        <  6:118890/102‑1 (MQ=255)
                                                                                                                                                                                                                                                                       cagtcattACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTg                                                                                                                                                                        >  5:118890/1‑102 (MQ=255)
                                                                                                                                                                                                                                                                       cagtcagtAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGt                                                                                                                                                                                                                            <  6:1052519/50‑1 (MQ=255)
                                                                                                                                                                                                                                                                       cagtcagtAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGt                                                                                                                                                                                                                            >  5:1052519/1‑50 (MQ=255)
                                                                                                                                                                                                                                                                       cagtcagtAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACg                                                                                                                                                     <  6:577459/121‑1 (MQ=255)
                                                                                                                                                                                                                                                                       cagtcagtAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACg                                                                                                                                                     <  6:629191/121‑1 (MQ=255)
                                                                                                                                                                                                                                                                       cagtcagtAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACg                                                                                                                                                     >  5:629191/1‑121 (MQ=255)
                                                                                                                                                                                                                                                                       cagtcagtAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACg                                                                                                                                                     >  5:577459/1‑121 (MQ=255)
                                                                                                                                                                                                                                                                             gtAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGaaa                                                                                                                      <  6:339710/146‑1 (MQ=255)
                                                                                                                                                                                                                                                                             gtAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGaaa                                                                                                                      >  5:339710/1‑146 (MQ=255)
                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                     
GTGGGGAGTGGCGAAAACTGGGATCAGGCGCACTAAGATTCGCCTGAACATGCCTTTTTTCGTAAGTAAGCAACATAAGCTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGC  >  NC_000913/4049717‑4050247

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