Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,121,583 Δ1 bp 6.4% intergenic (‑29/+26) 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,5830G. 6.4% 375.6 / 18.2 109intergenic (‑29/+26)wcaJ/cpsGcolanic biosynthesis UDP‑glucose lipid carrier transferase/phosphomannomutase
Reads supporting (aligned to +/- strand):  ref base G (51/51);  new base . (5/2);  total (56/53)
Fisher's exact test for biased strand distribution p-value = 4.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 109 total aligned reads displayed.

TGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGC  >  NC_000913/2121437‑2121729
                                                                                                                                                  |                                                                                                                                                   
tGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGc                                                                                                                                                    >  3:698004/1‑148 (MQ=255)
   cTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGt                                                                                                                                                 >  3:270172/1‑148 (MQ=255)
    ttCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTg                                                                                                                                                <  4:432206/148‑1 (MQ=255)
     tCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTg‑g                                                                                                                                              <  3:825704/148‑1 (MQ=255)
      cGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTg‑gg                                                                                                                                             >  4:7120/1‑148 (MQ=255)
       gCAAACCAACCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGt                                                                                                                                           <  1:9409/149‑1 (MQ=255)
         aaaCCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTa                                                                                                                                          <  4:123368/148‑1 (MQ=255)
          aaCCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTaa                                                                                                                                         >  4:937264/1‑148 (MQ=255)
          aaCCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAg                                                                                                                                        <  1:253806/149‑1 (MQ=255)
             cAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAggg                                                                                                                                      <  3:937264/148‑1 (MQ=255)
             cAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAggg                                                                                                                                      <  3:222222/148‑1 (MQ=255)
             cAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGa                                                                                                                                     >  1:74824/1‑149 (MQ=255)
              aGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGa                                                                                                                                     >  4:194394/1‑148 (MQ=255)
               gCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGaa                                                                                                                                    >  3:768270/1‑148 (MQ=255)
               gCCATACTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAg                                                                                                                                   >  2:525622/1‑149 (MQ=255)
                ccATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAg                                                                                                                                   <  4:270172/148‑1 (MQ=255)
                 cATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGAt                                                                                                                                 <  2:149259/149‑1 (MQ=255)
                 cATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGAt                                                                                                                                 >  1:633445/1‑149 (MQ=255)
                  aTAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATc                                                                                                                                <  1:577594/149‑1 (MQ=255)
                     gTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCg                                                                                                                              <  3:60524/148‑1 (MQ=255)
                      tCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGa                                                                                                                             >  3:1041515/1‑148 (MQ=255)
                          gCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCg                                                                                                                              <  4:447388/143‑1 (MQ=255)
                           cAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTAc                                                                                                                       <  2:619228/149‑1 (MQ=255)
                             aaCATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCg                                                                                                                              >  3:577742/1‑140 (MQ=255)
                             aaCATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACt                                                                                                                      <  3:530210/148‑1 (MQ=255)
                             aaCATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTc                                                                                                                     >  1:905043/1‑149 (MQ=255)
                                 tgatgGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGtt                                                                                                                  <  4:768270/148‑1 (MQ=255)
                                  gatgGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCa                                                                                                                <  1:525622/149‑1 (MQ=255)
                                  gatgGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCa                                                                                                                <  2:748690/149‑1 (MQ=255)
                                     ggTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCa                                                                                                             >  1:1007964/1‑149 (MQ=255)
                                       tGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGATATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAAc                                                                                                           <  2:894368/149‑1 (MQ=255)
                                           atCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTca                                                                                                       >  2:365855/1‑149 (MQ=255)
                                            tCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGTCATTACTCGTTCAGCCACGTcag                                                                                                      <  2:858472/149‑1 (MQ=255)
                                            tCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTcag                                                                                                      <  2:905043/149‑1 (MQ=255)
                                               gAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGa                                                                                                                             <  3:965617/123‑1 (MQ=255)
                                               gAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGa                                                                                                                             >  4:965617/1‑123 (MQ=255)
                                               gAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCC‑CACTGCGGGTAAggg                                                                                                                                      >  4:144297/1‑114 (MQ=255)
                                               gAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCC‑CACTGCGGGTAAggg                                                                                                                                      <  3:144297/114‑1 (MQ=255)
                                                aaaaGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCATACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCaga                                                                                                  <  2:791377/149‑1 (MQ=255)
                                                 aaaGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCagag                                                                                                 >  2:596970/1‑149 (MQ=255)
                                                   aGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAgt                                                                                                >  3:732171/1‑148 (MQ=255)
                                                        tGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAgttcgcgt                                                                                          <  2:962636/149‑1 (MQ=255)
                                                        tGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAgttcgcg                                                                                           <  4:14521/148‑1 (MQ=255)
                                                            ccATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAgttcgcgttcgc                                                                                      >  1:522618/1‑149 (MQ=255)
                                                             cATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAgttcgcgttcgc                                                                                      >  3:1004120/1‑148 (MQ=255)
                                                               tAGAGATTAACGCTGCATTGGCTTTCGCTCGCTCGCGCTTTTTTAGATTTATCATCGTTGTTCGTGTTATGAGCCCCTTACCCGCAGTG‑GGGTAAGCGAAGATCCGACATCACTCGTTCAGCAACGTCAGCAGATTACGCGTTCGCGCt                                                                                   >  2:81113/1‑149 (MQ=255)
                                                               tAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCTTCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCgcgc                                                                                    >  3:249609/1‑148 (MQ=255)
                                                                 gagaTTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTc                                                                                 <  2:18771/149‑1 (MQ=255)
                                                                   gaTTAACGGTGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTcc                                                                                <  3:656468/148‑1 (MQ=255)
                                                                    aTTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCa                                                                               <  4:1004120/148‑1 (MQ=255)
                                                                      tAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATc                                                                             <  4:391825/148‑1 (MQ=255)
                                                                            tGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGc                                                                       >  4:865298/1‑148 (MQ=255)
                                                                               aTTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAgttcgcgttc                                                                                        <  1:894747/128‑1 (MQ=255)
                                                                               aTTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAgttcgcgttc                                                                                        >  2:894747/1‑128 (MQ=255)
                                                                                ttGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACAt                                                                   >  4:880673/1‑148 (MQ=255)
                                                                                 tGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATc                                                                  >  3:36330/1‑148 (MQ=255)
                                                                                   gTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCAcc                                                               >  1:1037360/1‑149 (MQ=255)
                                                                                         gctcgctcgcGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACcgcgcg                                                          >  4:916209/1‑148 (MQ=255)
                                                                                         gctcgctcgcGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGa                                                         <  2:672166/149‑1 (MQ=255)
                                                                                             gctcgcGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTc                                                      <  4:261899/148‑1 (MQ=255)
                                                                                             gctcgcGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTc                                                      <  4:758454/148‑1 (MQ=255)
                                                                                               tcgcGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAgt                                                                                                >  2:949353/1‑104 (MQ=255)
                                                                                               tcgcGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAgt                                                                                                <  1:949353/104‑1 (MQ=255)
                                                                                                cgcgCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCC‑CACTGCGGGTAAgg                                                                                                                                       >  3:443769/1‑64 (MQ=25)
                                                                                                cgcgCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCC‑CACTGCGGGTAAgg                                                                                                                                       <  4:443769/64‑1 (MQ=25)
                                                                                                  cgcTTTTTAAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACAtt                                                <  2:927049/149‑1 (MQ=255)
                                                                                                   gcTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACAtt                                                <  4:36330/148‑1 (MQ=255)
                                                                                                     ttttttAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAg                                             >  2:410387/1‑149 (MQ=255)
                                                                                                             aTTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGcacc                                      >  4:391811/1‑148 (MQ=255)
                                                                                                              tttGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGcacca                                     >  3:751316/1‑148 (MQ=255)
                                                                                                              tttGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGcacca                                     >  3:735991/1‑148 (MQ=255)
                                                                                                                   cATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCATAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCATGCGCACCACCGGtt                               <  1:807407/149‑1 (MQ=255)
                                                                                                                    aTCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTc                              >  2:938405/1‑149 (MQ=255)
                                                                                                                    aTCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTc                              <  2:522618/149‑1 (MQ=255)
                                                                                                                       gttgttCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGt                           >  2:587997/1‑149 (MQ=255)
                                                                                                                         tgttCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTa                          >  4:599703/1‑148 (MQ=255)
                                                                                                                           ttCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTg                       <  2:404048/149‑1 (MQ=255)
                                                                                                                            tCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTgg                      <  1:155531/149‑1 (MQ=255)
                                                                                                                            tCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTgg                      <  1:496304/149‑1 (MQ=255)
                                                                                                                              cTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCTCCGGTTCGGTATTGGa                     <  4:403715/148‑1 (MQ=255)
                                                                                                                               tGTTATTGGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCagag                                                                                                 >  2:574575/1‑71 (MQ=255)
                                                                                                                               tGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAgg                   >  1:1029287/1‑149 (MQ=255)
                                                                                                                               tGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAgg                   >  1:49880/1‑149 (MQ=255)
                                                                                                                                gTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGt                  <  1:48019/149‑1 (MQ=255)
                                                                                                                                 ttattaGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTg                 >  2:320093/1‑149 (MQ=255)
                                                                                                                                  tattaGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTg                 >  4:301455/1‑148 (MQ=255)
                                                                                                                                       ttcccTTACCC‑CACTGCGGGTAAGGGGAGGTCCGACTTTACTCGTTCAGCAACGTCAGCTGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGGGCGATTCCAGATTCAGGCGCACCAACGGTTCGGTATTGGAGGTGCGCAgg           >  1:360635/3‑149 (MQ=255)
                                                                                                                                        tcccTTACCC‑CACTGCGGGTAAGGGGAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAgg           >  4:934455/2‑148 (MQ=255)
                                                                                                                                        tcccTTACCC‑CACTGCGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAgg           >  4:11759/2‑148 (MQ=255)
                                                                                                                                            ttACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTaa       <  3:891198/148‑1 (MQ=255)
                                                                                                                                             tACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCg                                                                         <  3:961787/81‑1 (MQ=255)
                                                                                                                                             tACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCg                                                                         >  4:961787/1‑81 (MQ=255)
                                                                                                                                             tACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGTGCAGGTTAAAg     <  1:560876/149‑1 (MQ=255)
                                                                                                                                             tACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTaaa      >  3:1030803/1‑148 (MQ=255)
                                                                                                                                              aCCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGGGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGt                  >  1:681429/1‑135 (MQ=255)
                                                                                                                                              aCCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGGGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGt                  <  2:681429/135‑1 (MQ=255)
                                                                                                                                              aCCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAgc    <  1:938405/149‑1 (MQ=255)
                                                                                                                                                 cGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAgcgc  >  4:1050244/1‑148 (MQ=255)
                                                                                                                                                 cGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAgcgc  >  4:545010/1‑148 (MQ=255)
                                                                                                                                                 cGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAgcgc  <  3:542024/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGC  >  NC_000913/2121437‑2121729

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