Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,121,586 G→C 8.3% intergenic (‑32/+23) wcaJ ← / ← cpsG colanic biosynthesis UDP‑glucose lipid carrier transferase/phosphomannomutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,121,5860GC 8.3% 364.3 / 19.3 109intergenic (‑32/+23)wcaJ/cpsGcolanic biosynthesis UDP‑glucose lipid carrier transferase/phosphomannomutase
Reads supporting (aligned to +/- strand):  ref base G (50/50);  new base C (2/7);  total (52/57)
Fisher's exact test for biased strand distribution p-value = 1.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.88e-01

Only 100 of 109 total aligned reads displayed.

TTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGCC  >  NC_000913/2121441‑2121730
                                                                                                                                                 |                                                                                                                                                 
ttCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTg‑g                                                                                                                                               <  4:52325/149‑1 (MQ=255)
gttgCAACCCAGCCATAGACAGTCAAACATGATGGTGCTATCTGGAAAGCGTTGCACCATAGAGATTAACGATGCATAGGTTGTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTg‑g                                                                                                                                               >  2:613964/4‑149 (MQ=255)
   gCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCTTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGt                                                                                                                                            <  4:384863/149‑1 (MQ=255)
    cAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTa                                                                                                                                           >  1:238329/1‑149 (MQ=255)
      aaCCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAg                                                                                                                                         <  2:389899/149‑1 (MQ=255)
      aaCCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAg                                                                                                                                         <  1:654568/149‑1 (MQ=255)
       aCCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAgg                                                                                                                                        <  2:319931/149‑1 (MQ=255)
       aCCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAgg                                                                                                                                        <  2:482666/149‑1 (MQ=255)
        ccAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCC‑CACTGCGGGTAAggg                                                                                                                                       <  3:734603/149‑1 (MQ=255)
        ccAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCC‑CACTGCGGGTAAggg                                                                                                                                       <  2:650969/149‑1 (MQ=255)
        ccAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCC‑CACTGCGGGTAAggg                                                                                                                                       <  3:200770/149‑1 (MQ=255)
        ccAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCC‑CACTGCGGGTAAggg                                                                                                                                       <  3:658505/149‑1 (MQ=255)
        ccAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCC‑CACTGCGGGTAAggg                                                                                                                                       <  4:387924/149‑1 (MQ=255)
         cAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCC‑CACTGCGGGTAAgggg                                                                                                                                      <  2:829598/149‑2 (MQ=255)
          aGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGaa                                                                                                                                     <  3:490874/149‑1 (MQ=255)
             cATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGAt                                                                                                                                  >  1:88362/1‑149 (MQ=255)
             cATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGAt                                                                                                                                  >  2:113413/1‑149 (MQ=255)
                 gTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCTTTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGCAAGGGAAGATCCGa                                                                                                                              <  4:346978/149‑1 (MQ=255)
                      gtaaaCATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGATTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTa                                                                                                                         <  4:757563/147‑1 (MQ=255)
                      gCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTg‑gg                                                                                                                                              >  4:228854/1‑128 (MQ=255)
                      gCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTg‑gg                                                                                                                                              <  3:228854/128‑1 (MQ=255)
                       cAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTAc                                                                                                                        <  1:277332/149‑1 (MQ=255)
                         aaCATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTc                                                                                                                      >  1:180411/1‑149 (MQ=255)
                             tgatgGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTc                                                                                                                  >  3:349015/1‑149 (MQ=255)
                              gatgGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCa                                                                                                                 <  4:155849/149‑1 (MQ=255)
                               atgGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTATGGGAAGATCCGACATTACTCGTTCAg                                                                                                                <  4:654048/149‑1 (MQ=255)
                               atgGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCC‑CACTGCGGGTAAgg                                                                                                                                        >  1:129713/1‑125 (MQ=255)
                               atgGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCC‑CACTGCGGGTAAgg                                                                                                                                        <  2:129713/125‑1 (MQ=255)
                                  gTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCaa                                                                                                             >  3:220843/1‑149 (MQ=255)
                                   tGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAAc                                                                                                            <  1:769216/149‑1 (MQ=255)
                                   tGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAAc                                                                                                            <  4:349015/149‑1 (MQ=255)
                                   tGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAAc                                                                                                            <  2:835281/149‑1 (MQ=255)
                                        tCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTcag                                                                                                       >  4:632994/1‑149 (MQ=255)
                                         cTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGt                                                                                                          >  3:337761/1‑145 (MQ=255)
                                         cTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGt                                                                                                          <  4:337761/145‑1 (MQ=255)
                                          tGAAAAGCGTTGCACCATAGATATTAACGATGCATTGGTTTTCGCTCGCTCACGCTTTTTTACATTTGTCATCGTTGTTCCTGTTATTAGCTCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACAGTACTCGTTCAGCAACGTcagca                                                                                                     <  3:569036/149‑1 (MQ=255)
                                          tGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAAc                                                                                                            <  4:495303/142‑1 (MQ=255)
                                          tGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAAc                                                                                                            >  3:495303/1‑142 (MQ=255)
                                          tGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTcagca                                                                                                     >  3:359352/1‑149 (MQ=255)
                                          tGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTcagca                                                                                                     >  4:125639/1‑149 (MQ=255)
                                            aaaaGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCaga                                                                                                   >  3:332713/1‑149 (MQ=255)
                                               aGCGTTGCACCATAGAGATTCACGATGCATTGGTTTTCGCTCGCGGGCGCTATGTTAGATTTGTCATTGTTGTTCCTGTTATTAGCCCCTTACCCGCCGCG‑GGGTCCGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCagagct                                                                                                >  2:6038/1‑147 (MQ=255)
                                                  gTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAgttcgc                                                                                             <  4:6000/149‑1 (MQ=255)
                                                  gTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAgttcgc                                                                                             >  1:873094/1‑149 (MQ=255)
                                                          aTAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACTCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCgcgc                                                                                     >  1:86270/1‑149 (MQ=255)
                                                          aTAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCgcgc                                                                                     >  4:298221/1‑149 (MQ=255)
                                                          aTAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCgcgc                                                                                     <  2:254395/149‑1 (MQ=255)
                                                            agagATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCtt                                                                                   <  4:149371/149‑1 (MQ=255)
                                                             gagaTTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTc                                                                                  <  2:88362/149‑1 (MQ=255)
                                                              agaTTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTcc                                                                                 <  3:822057/149‑1 (MQ=255)
                                                                aTTAACGATGCATTGGTTCTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCAt                                                                               >  3:118284/1‑149 (MQ=255)
                                                                  tAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCa                                                                             <  2:120937/149‑1 (MQ=255)
                                                                    aCGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGc                                                                           <  2:733729/149‑1 (MQ=255)
                                                                      gATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCgg                                                                         >  4:835664/1‑149 (MQ=255)
                                                                       aTGCATTGGTTTTCGCTCCGTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATGACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGc                                                                        <  2:436779/149‑1 (MQ=255)
                                                                         gCATTGGTTTTCGCTCGCTAGCGCTTTTTTAGATTTGTCATCGTTGTTGCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACAGTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGcac                                                                      <  4:769011/149‑1 (MQ=255)
                                                                            ttGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAgg                                                                                                                                        >  3:848298/1‑80 (MQ=255)
                                                                            ttGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAgg                                                                                                                                        <  4:848298/80‑1 (MQ=255)
                                                                                       tcgctcgcGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGAtt                                                        <  1:113413/149‑1 (MQ=255)
                                                                                         gctcgcGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTcc                                                      <  3:423159/149‑1 (MQ=255)
                                                                                         gctcgcGCTTGTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTcc                                                      >  2:677477/1‑149 (MQ=255)
                                                                                           tcgcGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCa                                                                  >  4:262845/1‑135 (MQ=255)
                                                                                           tcgcGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCa                                                                  <  3:262845/135‑1 (MQ=255)
                                                                                             gcgcTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTCCTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACAt                                                  >  3:532842/1‑149 (MQ=255)
                                                                                               gcTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTc                                                <  3:579150/149‑1 (MQ=255)
                                                                                                cTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCa                                               <  4:744062/149‑1 (MQ=255)
                                                                                                 ttttttAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAg                                              >  1:140649/1‑149 (MQ=255)
                                                                                                  tttttAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAgg                                             <  4:754874/149‑1 (MQ=255)
                                                                                                       aGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGcac                                        >  4:758918/1‑149 (MQ=255)
                                                                                                              tCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGt                                 >  3:469122/1‑149 (MQ=255)
                                                                                                              tCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGt                                 >  3:352840/1‑149 (MQ=255)
                                                                                                                 tCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCg                              >  1:326609/1‑149 (MQ=255)
                                                                                                                    ttgttCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGACGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTa                           <  2:140649/149‑1 (MQ=255)
                                                                                                                      gttCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTAtt                         >  3:29856/1‑149 (MQ=255)
                                                                                                                          cTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAg                     >  3:69631/1‑149 (MQ=255)
                                                                                                                          cTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAg                     <  3:141475/149‑1 (MQ=255)
                                                                                                                          cTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAg                     >  3:186994/1‑149 (MQ=255)
                                                                                                                           tGTTATTAGCCTCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCATCAACGTcag                                                                                                       >  4:606655/1‑66 (MQ=255)
                                                                                                                           tGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTATTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAgg                    >  4:647138/1‑149 (MQ=255)
                                                                                                                           tGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTcag                                                                                                       <  3:606655/66‑1 (MQ=255)
                                                                                                                           tGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAgg                    >  1:450828/1‑149 (MQ=255)
                                                                                                                           tGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAgg                    >  3:99565/1‑149 (MQ=255)
                                                                                                                            gTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGt                   <  1:231742/149‑1 (MQ=255)
                                                                                                                             ttattaGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTg                  >  1:280646/1‑149 (MQ=255)
                                                                                                                              tattaGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTgc                 <  1:6996/149‑1 (MQ=255)
                                                                                                                              tattaGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGAgc                 <  4:469122/149‑1 (MQ=255)
                                                                                                                                ttaGCCCCTTTCCCGCAGTG‑GGGTAAGGGAAGATCCGCCATGACTCAGTCACCAACGTCAGCAGAGTTCTCGGTCGCGCTGCCACCAGGGGCACATCACCCCGCGATTCTACGTTCGGCAGCACCACCGGTTCGGTATTGGAGGTgcgc               <  4:270063/149‑1 (MQ=255)
                                                                                                                                ttaGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTgcgc               >  4:871576/1‑149 (MQ=255)
                                                                                                                                  aGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAg             <  2:326609/149‑1 (MQ=255)
                                                                                                                                   ttcccTTACCC‑CACTGCGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAgg            >  4:153656/3‑149 (MQ=255)
                                                                                                                                    ccccTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGt           >  4:498551/1‑149 (MQ=255)
                                                                                                                                     cccTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTa                           >  4:173843/1‑132 (MQ=255)
                                                                                                                                     cccTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTa                           <  3:173843/132‑1 (MQ=255)
                                                                                                                                     cccTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGtt          >  1:40472/1‑149 (MQ=255)
                                                                                                                                      ccTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTa         >  3:330943/1‑149 (MQ=255)
                                                                                                                                      ccTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTa         >  2:171755/1‑149 (MQ=255)
                                                                                                                                       cTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTaa        >  1:434001/1‑149 (MQ=255)
                                                                                                                                          aCCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTACAgc     <  4:29856/149‑1 (MQ=255)
                                                                                                                                            ccGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAgcgc   >  2:128213/1‑149 (MQ=255)
                                                                                                                                             cGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGcc  <  3:417271/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
TTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGCC  >  NC_000913/2121441‑2121730

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