Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,121,586 G→C 6.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 6.3% 366.8 / 8.5 109intergenic (‑32/+23)wcaJ/cpsGcolanic biosynthesis UDP‑glucose lipid carrier transferase/phosphomannomutase
Reads supporting (aligned to +/- strand):  ref base G (51/51);  new base C (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.

CTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGCCA  >  NC_000913/2121440‑2121731
                                                                                                                                                  |                                                                                                                                                  
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:222222/148‑1 (MQ=255)
          cAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAggg                                                                                                                                        <  3:937264/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                                                                                                                                   >  1:633445/1‑149 (MQ=255)
              cATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGAt                                                                                                                                   <  2:149259/149‑1 (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‑GGGTAAGGGAAGATCCg                                                                                                                                <  4:577742/140‑1 (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)
                                                                                        tcgctcgcGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGAt                                                          >  4:197856/1‑148 (MQ=255)
                                                                                          gctcgcGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTc                                                        <  4:758454/148‑1 (MQ=255)
                                                                                            tcgcGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAgt                                                                                                  <  1:949353/104‑1 (MQ=255)
                                                                                            tcgcGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAgt                                                                                                  >  2:949353/1‑104 (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)
                                                                                                                        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:49880/1‑149 (MQ=255)
                                                                                                                            tGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAgg                     >  1:1029287/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)
                                                                                                                                   aGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCa               <  3:865298/148‑1 (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)
                                                                                                                                          tACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCg                                                                           >  4:961787/1‑81 (MQ=255)
                                                                                                                                          tACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCg                                                                           <  3:961787/81‑1 (MQ=255)
                                                                                                                                          tACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGTGCAGGTTAAAg       <  1:560876/149‑1 (MQ=255)
                                                                                                                                          tACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTaaa        >  3:1030803/1‑148 (MQ=255)
                                                                                                                                           aCCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGGGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGt                    <  2:681429/135‑1 (MQ=255)
                                                                                                                                           aCCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGGGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGt                    >  1:681429/1‑135 (MQ=255)
                                                                                                                                           aCCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAgc      <  1:938405/149‑1 (MQ=255)
                                                                                                                                              cGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAgcgc    >  4:545010/1‑148 (MQ=255)
                                                                                                                                              cGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAgcgc    >  4:1050244/1‑148 (MQ=255)
                                                                                                                                              cGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAgcgc    <  3:542024/148‑1 (MQ=255)
                                                                                                                                                cAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGCCa  >  4:686130/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGCCA  >  NC_000913/2121440‑2121731

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