Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,121,586 G→C 5.6% 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 5.6% 247.6 / 5.6 71intergenic (‑32/+23)wcaJ/cpsGcolanic biosynthesis UDP‑glucose lipid carrier transferase/phosphomannomutase
Reads supporting (aligned to +/- strand):  ref base G (37/30);  new base C (2/2);  total (39/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-01

GACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAA‑GGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGCC  >  NC_000913/2121438‑2121730
                                                                                                                                                    |                                                                                                                                                  
gACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAg                                                                                                                                                    >  3:51608/1‑149 (MQ=255)
   ttCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTg‑g                                                                                                                                                >  4:226646/1‑149 (MQ=255)
   ttCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTg‑g                                                                                                                                                <  3:397701/149‑1 (MQ=255)
   ttCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTg‑g                                                                                                                                                >  3:206652/1‑149 (MQ=255)
   ttCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTg‑g                                                                                                                                                >  2:249148/1‑149 (MQ=255)
        aaaCCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTCCACCATAGAGATGAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTaa                                                                                                                                           >  4:455799/1‑149 (MQ=255)
         aaCCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAA‑g                                                                                                                                         <  3:609489/149‑1 (MQ=255)
          aCCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCC‑CACTGCGGGTAA‑gg                                                                                                                                        <  4:259674/149‑1 (MQ=255)
          aCCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCC‑CACTGCGGGTAA‑gg                                                                                                                                        <  2:609940/149‑1 (MQ=255)
           ccAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAA‑ggg                                                                                                                                       >  4:178461/1‑149 (MQ=255)
            cAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAA‑GGGa                                                                                                                                      >  1:491865/1‑149 (MQ=255)
             aGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAA‑GGGaa                                                                                                                                     <  4:6179/149‑1 (MQ=255)
                cATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAA‑GGGAAGAt                                                                                                                                  >  4:140955/1‑149 (MQ=255)
                   aGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAA‑GGGAAGATCCg                                                                                                                               >  1:119690/1‑149 (MQ=255)
                    gTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAA‑GGGAAGATCCGa                                                                                                                              <  1:425118/149‑1 (MQ=255)
                     tCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAA‑GGGAAGATCCGAc                                                                                                                             >  2:478064/1‑149 (MQ=255)
                       cGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAA‑GGGAAGATCCGACAt                                                                                                                           <  3:528407/149‑1 (MQ=255)
                         gCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAA‑GGGAAGATCCGACATTa                                                                                                                         <  3:281961/149‑1 (MQ=255)
                            aaCATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAA‑GGGAAGATCCGACATTACTc                                                                                                                      <  2:119690/149‑1 (MQ=255)
                                tgatgGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAA‑GGGAAGATCCGACATTACTCGTTc                                                                                                                  >  3:581886/1‑149 (MQ=255)
                                      tGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAA‑GGGAAGATCCGACATTACTCGTTCAGCAAc                                                                                                            <  3:215519/149‑1 (MQ=255)
                                      tGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAA‑GGGAAGATCCGACATTACTCGTTCAGCAAc                                                                                                            <  1:191884/149‑1 (MQ=255)
                                        atatCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAA‑GGGAAGATCCGACATTACTCGTTCAGCAACGt                                                                                                          >  1:433812/1‑149 (MQ=255)
                                           tCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAA‑GGGAAGATCCGACATTACTCGTTCAGCAACGTcag                                                                                                       >  1:464963/1‑149 (MQ=255)
                                              gAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAA‑GGGAAGATCCGACATTACTCGTTCAGCAACGTcagcag                                                                                                    >  4:355953/1‑149 (MQ=255)
                                               aaaaGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAA‑GGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCaga                                                                                                   >  2:33622/1‑149 (MQ=255)
                                                     gTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAA‑GGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAgttcgc                                                                                             >  1:275436/1‑149 (MQ=255)
                                                        gCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAA‑GGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAgttcgcgtt                                                                                          >  4:556927/1‑149 (MQ=255)
                                                         cACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAA‑GGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAgttcgcgttc                                                                                         <  3:178461/149‑1 (MQ=255)
                                                                gagaTTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAA‑GGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTc                                                                                  >  1:148991/1‑149 (MQ=255)
                                                                gagaTTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAA‑GGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTc                                                                                  <  2:168714/149‑1 (MQ=255)
                                                                      aaCGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAA‑GGGAAGATCCGACATTACTCGTTCAGCAACGt                                                                                                          <  2:486752/119‑1 (MQ=255)
                                                                      aaCGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAA‑GGGAAGATCCGACATTACTCGTTCAGCAACGt                                                                                                          >  1:486752/1‑119 (MQ=255)
                                                                       acgctacATTGGTTTTCGCTCGCTCGCCCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAA‑GGGAAGATCCGACATTACTCGTTCAGCGACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGc                                                                           <  4:536574/143‑1 (MQ=255)
                                                                          aTGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAA‑GGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGc                                                                        <  3:355953/149‑1 (MQ=255)
                                                                              aTTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAA‑GGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACAt                                                                    >  3:104442/1‑149 (MQ=255)
                                                                                  gTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAA‑GGGAAGATCCGACATTACTCGTTCAGCAAc                                                                                                            >  1:119383/1‑105 (MQ=255)
                                                                                  gTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAA‑GGGAAGATCCGACATTACTCGTTCAGCAAc                                                                                                            <  2:119383/105‑1 (MQ=255)
                                                                                  gTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAA‑GGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCAcc                                                                <  3:226646/149‑1 (MQ=255)
                                                                                   ttttCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAA‑GGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACcg                                                               >  2:121181/1‑149 (MQ=255)
                                                                                        gctcgctcgcGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAA‑GGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGa                                                          <  1:33622/149‑1 (MQ=255)
                                                                                          tcgctcgcGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAA‑GGTAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGAtt                                                        <  4:206652/149‑1 (MQ=255)
                                                                                           cgctcgcGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAA‑GGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTc                                                       <  3:140955/149‑1 (MQ=255)
                                                                                              tcgcGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAA‑GGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCcac                                                    >  3:257871/1‑149 (MQ=255)
                                                                                               cgcgCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAA‑GGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCcaca                                                   <  3:196505/149‑1 (MQ=255)
                                                                                                gcgcTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAA‑GGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACAt                                                  >  3:242523/1‑149 (MQ=255)
                                                                                                 cgcTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAA‑GGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACAtt                                                 >  3:559217/1‑149 (MQ=255)
                                                                                                        ttAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAA‑GGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGgcgc                                          <  2:483724/149‑1 (MQ=255)
                                                                                                           gATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAA‑GGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCTGGCGcacc                                       >  2:508281/1‑149 (MQ=255)
                                                                                                           gATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAA‑GGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGATCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGcacc                                       >  1:268676/1‑149 (MQ=255)
                                                                                                                     cGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAA‑GGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCgg                             <  2:275436/149‑1 (MQ=255)
                                                                                                                     cGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAA‑GGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCgg                             <  1:478064/149‑1 (MQ=255)
                                                                                                                        tgttCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAA‑GGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGATCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTAt                          <  2:268676/149‑1 (MQ=255)
                                                                                                                          ttCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAA‑GGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTg                        >  3:205436/1‑149 (MQ=255)
                                                                                                                            ccTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAA‑GGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGa                      <  1:121181/149‑1 (MQ=255)
                                                                                                                             cTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAA‑GGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAg                     >  2:86425/1‑149 (MQ=255)
                                                                                                                             cTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAA‑GGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAg                     <  2:464963/149‑1 (MQ=255)
                                                                                                                             cTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAA‑GGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAg                     >  3:361427/1‑149 (MQ=255)
                                                                                                                              tGTTATTAGCCCCTTACCCGCAGTG‑GGGTAA‑GGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAgg                    >  4:328636/1‑149 (MQ=255)
                                                                                                                                  gatcttcccTTACCC‑CACTGCGGGTAA‑GGGGAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTgcg                >  4:274023/7‑149 (MQ=255)
                                                                                                                                     cttcccTTACCC‑CACTGCGGGTAAGGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGTGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCa              >  1:504612/4‑149 (MQ=255)
                                                                                                                                      gCCCCTTACCCGCAGTG‑GGGTAA‑GGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAgg            <  3:455799/149‑1 (MQ=255)
                                                                                                                                      gCCCCTTACCCGCAGTG‑GGGTAA‑GGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAgg            <  2:189416/149‑1 (MQ=255)
                                                                                                                                        cccTTACCCGCAGTG‑GGGTAA‑GGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGtt                                <  2:344795/127‑1 (MQ=255)
                                                                                                                                        cccTTACCCGCAGTG‑GGGTAA‑GGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGtt                                >  1:344795/1‑127 (MQ=255)
                                                                                                                                        cccTTACCCGCAGTG‑GGGTAA‑GGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGtt          >  2:233239/1‑149 (MQ=255)
                                                                                                                                         ccTTACCCGCAGTG‑GGGTAC‑GGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTa         >  3:513897/1‑149 (MQ=255)
                                                                                                                                         ccTTACCCGCAGTG‑GGGTAA‑GGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTa         >  1:155271/1‑149 (MQ=255)
                                                                                                                                           ttACCCGCAGTG‑GGGTAA‑GGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCTCATCACCGCGCGAGTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTaaa       >  2:513289/1‑149 (MQ=255)
                                                                                                                                               ccGCAGTG‑GGGTAA‑GGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAgcgc   >  1:297963/1‑149 (MQ=255)
                                                                                                                                                cGCAGTG‑GGGTAA‑GGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGcc  <  2:491865/149‑1 (MQ=255)
                                                                                                                                                cGCAGTG‑GGGTAA‑GGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGcc  >  2:82190/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAA‑GGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGCC  >  NC_000913/2121438‑2121730

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