Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,049,986 G→T 7.2% 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 7.2% 632.7 / 29.3 181noncoding (88/109 nt)spfSpot 42 sRNA antisense regulator of galK translation, Hfq‑dependent
Reads supporting (aligned to +/- strand):  ref base G (81/87);  new base T (9/4);  total (90/91)
Fisher's exact test for biased strand distribution p-value = 1.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.89e-01

Only 100 of 190 total aligned reads displayed.

GTGGGGAGTGGCGAAAACTGGGATCAGGCGCACTAAGATTCGCCTGAACATGCCTTTTTTCGTAAGTAAGCAACATAAGCTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTA  >  NC_000913/4049717‑4050259
                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                 
   gggAGTGGCGAAAACTGGGATCAGGCGCACTAAGATTCGCCTGAACATGCCTTTTTTCGTAAGTAAGCAACATAAGCTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTGAAAACAAATTACAAAAAGTTCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTTGCGTAGGGTACAGCGGTAAGATGTTCAATCTTTCAGACCTTTTACTTCACGTAATAGGATTTGGCTGACTTTTTTAGCCCACCCTGTCAGTAATGAcg                                                                                                                                                                                                                                                                           >  4:261961/1‑274 (MQ=255)
                                                   gCCTTTTTTCGTAAGTAAGCAACATAAGCTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATgc                                                                                                                                                                                                                            >  4:242387/1‑274 (MQ=255)
                                                                   aagCAACATAAGCTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGcc                                                                                                                                                                                                           <  2:242387/275‑1 (MQ=255)
                                                                                   cACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGTGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTaa                                                                                                                                                                                           >  2:151548/1‑275 (MQ=255)
                                                                                                ggCTATTAGAAATTCCTATGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGAc                                                                                                                                                                                                                                                                            >  3:69591/1‑181 (MQ=255)
                                                                                                 gCTATTAGAAATTCCTATGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAACGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTCATGACTGGGGCGTTTGTTTTTGGGCGAAAGAAAAGATCCGTAATGCCTGATGTGCTATGTTTATCGGGCCAACGGTAGAATTTGAATCTATTGAATTTa                                                                                                                                                                              >  4:203677/1‑274 (MQ=255)
                                                                                                    aTTAGAAATTCCTATGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTggg                                                                                                                                                                                                                                                                        <  1:69591/181‑1 (MQ=255)
                                                                                                         aaaTTCCCTTGCAAAAAATGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAAAAAAAAAGAGTAAAGTTAGCCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGa                                                                                                                                                                     <  4:107842/275‑1 (MQ=255)
                                                                                                                    caacaaCTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATc                                                                                                                                                          <  2:283215/275‑1 (MQ=255)
                                                                                                                         aattaaaaaaaaATACAAAAAGTGATTTCCGAACTGAACAAAAAAGAGGAAAGTTGGTCGCGTAGGGTACAGAGGTAATATGTTCTAGCTTTCAGACCGTTTACTTCACGTAATCTGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGca                                                                                                                                                     <  4:236230/271‑1 (MQ=255)
                                                                                                                          cTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCcagtcagt                                                                                                                                                                                                                                                                                  <  5:643421/149‑1 (MQ=255)
                                                                                                                          cTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCcagtcagt                                                                                                                                                                                                                                                                                  >  8:835223/1‑149 (MQ=255)
                                                                                                                            gAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTaa                                                                                                                                                                                                                                                                                <  7:727896/149‑1 (MQ=255)
                                                                                                                            gAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTaa                                                                                                                                                                                                                                                                                >  7:605322/1‑149 (MQ=255)
                                                                                                                            gAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTaa                                                                                                                                                                                                                                                                                <  8:11550/149‑1 (MQ=255)
                                                                                                                                   aaTTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTgg                                                                                                                                                                                                                                                                         <  5:251234/149‑1 (MQ=255)
                                                                                                                                   aaTTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTgg                                                                                                                                                                                                                                                                         <  5:575635/149‑1 (MQ=255)
                                                                                                                                    aTTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTggg                                                                                                                                                                                                                                                                        <  7:387400/149‑1 (MQ=255)
                                                                                                                                    aTTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTggg                                                                                                                                                                                                                                                                        <  6:1321535/149‑1 (MQ=255)
                                                                                                                                         aaaaaGTGCATTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGTCGCCCCAGTCAGTAATGACTGGGGCGtc                                                                                                                                                                                                                                                                   >  8:1256095/1‑148 (MQ=255)
                                                                                                                                           aaaGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGtttt                                                                                                                                                                                                                                                                 >  7:1354295/1‑149 (MQ=255)
                                                                                                                                            aaGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGttttt                                                                                                                                                                                                                                                                <  6:360812/149‑1 (MQ=255)
                                                                                                                                                gCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTAtt                                                                                                                                                                                                                                                            <  6:701602/149‑1 (MQ=255)
                                                                                                                                                  tttCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTgg                                                                                                                                                                                                                                                          <  7:865116/149‑1 (MQ=255)
                                                                                                                                                  tttCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTgg                                                                                                                                                                                                                                                          <  8:570815/149‑1 (MQ=255)
                                                                                                                                                        aactgaacAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCATTACTGACTGGGGc                                                                                                                                                                                                                                                                      <  6:1429745/131‑1 (MQ=255)
                                                                                                                                                         actgaacAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCATTACTGACTGGGGc                                                                                                                                                                                                                                                                      <  8:139196/130‑1 (MQ=255)
                                                                                                                                                           tgaacAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCgaaagaa                                                                                                                                                                                                                                                 >  5:1347379/1‑149 (MQ=255)
                                                                                                                                                             aacAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGaaagaaaa                                                                                                                                                                                                                                               <  6:423959/149‑1 (MQ=255)
                                                                                                                                                             aacAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCa                                                                                                                                                                                                               >  2:350124/1‑181 (MQ=255)
                                                                                                                                                             aacAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCgg                                                                                                                                                       >  1:33300/1‑237 (MQ=255)
                                                                                                                                                             aacAAAAAAGAGTAAAGTGAGTCGCGTAGGGTAGAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGaaagaaaa                                                                                                                                                                                                                                               <  7:666025/149‑1 (MQ=255)
                                                                                                                                                              acAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCgaaa                                                                                                                                                                                                                                                    <  6:1027359/143‑1 (MQ=255)
                                                                                                                                                                aaaaaaGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGAt                                                                                                                                                                                                                                            <  7:1406837/149‑1 (MQ=255)
                                                                                                                                                                 aaaaaGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGcc                                                                                                                                                                                                           <  4:350124/181‑1 (MQ=255)
                                                                                                                                                                         tAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGc                                                                                                                                                                                                                                   >  5:164254/1‑149 (MQ=255)
                                                                                                                                                                         tAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGc                                                                                                                                                                                                                                   >  8:1182595/1‑149 (MQ=255)
                                                                                                                                                                                aGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCa                                                                                              >  2:326559/1‑275 (MQ=255)
                                                                                                                                                                                     cgTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTAt                                                                                                                                                                                                                       <  5:1113387/149‑1 (MQ=255)
                                                                                                                                                                                         gggTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGtttttt                                                                                                                                                                                                                                                               >  7:1103865/1‑105 (MQ=255)
                                                                                                                                                                                         gggTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGtttttt                                                                                                                                                                                                                                                               <  8:1103865/105‑1 (MQ=255)
                                                                                                                                                                                            tACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATc                                                                                                                                                                                                                >  5:679933/1‑149 (MQ=255)
                                                                                                                                                                                             aCAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCa                                                                                                                                                                                                               >  8:1419930/1‑149 (MQ=255)
                                                                                                                                                                                              cAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAg                                                                                                                                                                                                              <  8:1064068/149‑1 (MQ=255)
                                                                                                                                                                                              cAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAg                                                                                                                                                                                                              <  8:1426798/149‑1 (MQ=255)
                                                                                                                                                                                               agagGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGtttt                                                                                                                                                                                                                                                                 <  7:467094/97‑1 (MQ=255)
                                                                                                                                                                                               agagGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGtttt                                                                                                                                                                                                                                                                 >  8:467094/1‑97 (MQ=255)
                                                                                                                                                                                               agagGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAgg                                                                                                                                                                                                             >  6:1212444/1‑149 (MQ=255)
                                                                                                                                                                                                  ggTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTAt                                                                            <  4:326559/275‑1 (MQ=255)
                                                                                                                                                                                                   gTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCaa                                                                                                                                                                                                         >  8:1129458/1‑149 (MQ=255)
                                                                                                                                                                                                     aaGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACg                                                                                                                                                                                                       >  5:736577/1‑149 (MQ=255)
                                                                                                                                                                                                           ttCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGaa                                                                                                                                                                                                 <  6:1144444/149‑1 (MQ=255)
                                                                                                                                                                                                           ttCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGaa                                                                                                                                                                                                 >  8:362082/1‑149 (MQ=255)
                                                                                                                                                                                                           ttCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGaa                                                                                                                                                                                                 >  8:597781/1‑149 (MQ=255)
                                                                                                                                                                                                             cTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAAtt                                                                                                                                                                                               <  7:597781/149‑1 (MQ=255)
                                                                                                                                                                                                             cTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCCCTAAACCAGGTATCCAGTTTCTGc                                                                 >  4:41641/1‑275 (MQ=255)
                                                                                                                                                                                                              tATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATg                                                                                                                                                                                                                             <  8:883070/118‑1 (MQ=255)
                                                                                                                                                                                                              tATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATg                                                                                                                                                                                                                             >  7:883070/1‑118 (MQ=255)
                                                                                                                                                                                                              tATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTg                                                                                                                                                                                              >  8:1452556/1‑149 (MQ=255)
                                                                                                                                                                                                                  tttCAGACCTTTTACTTCACGTAATCGGTTTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAAt                                                                                                                                                                                          >  7:1046828/1‑149 (MQ=255)
                                                                                                                                                                                                                  tttCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTa                                                                                                                                                                                                    <  6:991860/139‑1 (MQ=255)
                                                                                                                                                                                                                  tttCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAAt                                                                                                                                                                                          >  7:821734/1‑149 (MQ=255)
                                                                                                                                                                                                                  tttCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAAt                                                                                                                                                                                          >  6:746705/1‑149 (MQ=255)
                                                                                                                                                                                                                   ttCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCg                                                                                                                                                        >  2:382171/1‑182 (MQ=255)
                                                                                                                                                                                                                     cAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCAtct                                                                                           >  4:53431/1‑241 (MQ=255)
                                                                                                                                                                                                                       gACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGcac                                                                                                                                                    <  4:382170/182‑1 (MQ=255)
                                                                                                                                                                                                                        aCCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTg                                                                                                                                                                                    >  5:383652/1‑149 (MQ=255)
                                                                                                                                                                                                                           ttttACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAAt                                                                                                                                                                                 >  5:281998/1‑149 (MQ=255)
                                                                                                                                                                                                                                 ttCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACgg                                                                                                                                                                           >  7:572763/1‑149 (MQ=255)
                                                                                                                                                                                                                                 ttCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACgg                                                                                                                                                                           >  8:302278/1‑149 (MQ=255)
                                                                                                                                                                                                                                 ttCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACgg                                                                                                                                                                           <  5:1385947/149‑1 (MQ=255)
                                                                                                                                                                                                                                 ttCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACgg                                                                                                                                                                           <  6:281998/149‑1 (MQ=255)
                                                                                                                                                                                                                                     cGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCg                                                                                                                                                                       >  5:1000543/1‑149 (MQ=255)
                                                                                                                                                                                                                                          tCGGATTTGGCTGAATATTTTAGCCTCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTAtt                                                                                                                                                                                     >  5:79545/1‑130 (MQ=255)
                                                                                                                                                                                                                                          tCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGTCGTTTTTTATTGGGCGAAAGAAAAGATCCGTa                                                                                                                                                                                                                                       >  6:103160/1‑80 (MQ=255)
                                                                                                                                                                                                                                          tCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTa                                                                                                                                                                                                                                       <  5:103160/80‑1 (MQ=255)
                                                                                                                                                                                                                                          tCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTAtt                                                                                                                                                                                     <  6:79545/130‑1 (MQ=255)
                                                                                                                                                                                                                                          tCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATAc                                                                                                                                                                  >  8:699714/1‑149 (MQ=255)
                                                                                                                                                                                                                                          tCGGATTTGGCTGAATATTTTAGCCGCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACg                                                                                                                                                                 >  7:331295/1‑149 (MQ=255)
                                                                                                                                                                                                                                                   gCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATcc                                                                                                                                                         <  6:1131382/149‑1 (MQ=255)
                                                                                                                                                                                                                                                    cTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCg                                                                                                                                                        >  6:94427/1‑149 (MQ=255)
                                                                                                                                                                                                                                                    cTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCg                                                                                                                                                        >  7:1184792/1‑149 (MQ=255)
                                                                                                                                                                                                                                                     tGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCgg                                                                                                                                                       >  8:160978/1‑149 (MQ=255)
                                                                                                                                                                                                                                                      gAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGc                                                                                                                                                      >  7:853652/1‑149 (MQ=255)
                                                                                                                                                                                                                                                          atctttgccgccCCCGTCAGTTATGACTGGGGTGTTTTTTATTGGGCAAAAGAAAAGGTCCGGTATGCTTGAGGCGCGAGGTTTATCAGGCCAACGGTCGAATTGTAATCAATGGAATTTACCGGGCGGGTAAGCCACAGCCGGCACAAGCATTAATGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGTCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCGGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCa                    <  4:336926/269‑1 (MQ=255)
                                                                                                                                                                                                                                                          atTTTAGTGGCCCCAGTTAGGAAGGACTGGGGCTGTTTTTATTGGGCGAAAGAAAAGTTCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTc                     <  2:33944/274‑1 (MQ=255)
                                                                                                                                                                                                                                                          atTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACaa                                                                                                                                                  <  7:699714/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                gccgccCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCg                                                                                                                                                        >  5:130099/1‑137 (MQ=255)
                                                                                                                                                                                                                                                                gccgccCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCg                                                                                                                                                        <  6:130099/137‑1 (MQ=255)
                                                                                                                                                                                                                                                                 ccgccCCCGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGTAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTaa                                                                                                                                           <  7:1256095/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                 ccgccCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTaa                                                                                                                                           <  6:250954/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                  cgccCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAg                                                                                                                                          <  8:962551/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                    ccccAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGc                                                                                                                                        >  6:176210/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                    ccccAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGc                                                                                                                                        >  5:255514/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                     cccAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCa                                                                                                                                       >  8:1170120/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                     cccAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCa                                                                                                                                       >  5:5120/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                     cccAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCaa                                                                                                                                                                                                         >  3:53525/1‑83 (MQ=255)
                                                                                                                                                                                                                                                                      ccAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCaa                                                                                                                                      >  6:315636/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                      ccAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCaa                                                                                                                                      >  7:1162220/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                         gtcagtAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGt                                                                                                                                                                                                     <  1:53525/83‑1 (MQ=255)
                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                 
GTGGGGAGTGGCGAAAACTGGGATCAGGCGCACTAAGATTCGCCTGAACATGCCTTTTTTCGTAAGTAAGCAACATAAGCTGTCACGTTTTGTGATGGCTATTAGAAATTCCTATGCAACAACTGAAAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTA  >  NC_000913/4049717‑4050259

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