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% 327.5 / 13.5 93intergenic (‑29/+26)wcaJ/cpsGcolanic biosynthesis UDP‑glucose lipid carrier transferase/phosphomannomutase
Reads supporting (aligned to +/- strand):  ref base G (45/42);  new base . (3/3);  total (48/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGCCA  >  NC_000913/2121436‑2121731
                                                                                                                                                   |                                                                                                                                                     
cTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCg                                                                                                                                                       >  2:122756/1‑148 (MQ=255)
cTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCg                                                                                                                                                       >  2:860676/1‑148 (MQ=255)
 tGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCGTTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCa                                                                                                                                                     >  3:548586/1‑149 (MQ=255)
 tGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCa                                                                                                                                                     >  4:441141/1‑149 (MQ=255)
  gACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAg                                                                                                                                                    <  3:422844/149‑1 (MQ=255)
     ttCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTg‑g                                                                                                                                                >  4:313426/1‑149 (MQ=255)
      tCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTg‑g                                                                                                                                                <  2:529655/148‑1 (MQ=255)
      tCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTg‑g                                                                                                                                                >  2:114446/1‑148 (MQ=255)
           aaCCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAg                                                                                                                                          >  3:475027/1‑149 (MQ=255)
            aCCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAgg                                                                                                                                         <  3:701513/149‑1 (MQ=255)
            aCCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAgg                                                                                                                                         >  3:310739/1‑149 (MQ=255)
             ccAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCC‑CACTGCGGGTAAggg                                                                                                                                        <  3:759390/149‑1 (MQ=255)
              cAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCC‑CACTGCGGGTAAggg                                                                                                                                        <  1:252239/148‑1 (MQ=255)
                 ccATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTg‑g                                                                                                                                                <  4:135441/137‑1 (MQ=255)
                 ccATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTg‑g                                                                                                                                                >  3:135441/1‑137 (MQ=255)
                 ccATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAg                                                                                                                                     <  1:778649/148‑1 (MQ=255)
                 ccATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGa                                                                                                                                    >  4:652347/1‑149 (MQ=255)
                  cATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGAt                                                                                                                                   >  3:566029/1‑149 (MQ=255)
                   aTAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATc                                                                                                                                  <  3:326597/149‑1 (MQ=255)
                      gTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCTCTAGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCgt                                                                                                                               <  4:784704/149‑2 (MQ=255)
                       tCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGAc                                                                                                                              >  4:852811/1‑149 (MQ=255)
                          ggCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACAtt                                                                                                                           <  4:132191/149‑1 (MQ=255)
                            cAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACGTTa                                                                                                                          <  2:222970/148‑1 (MQ=255)
                              aaCATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACt                                                                                                                        <  1:679623/148‑1 (MQ=255)
                              aaCATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACt                                                                                                                        <  2:859818/148‑1 (MQ=255)
                                 atgatgGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGtt                                                                                                                    <  3:740308/149‑1 (MQ=255)
                                 atgatgGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGt                                                                                                                     <  2:250143/148‑1 (MQ=255)
                                   gatgGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTc                                                                                                                   >  2:446175/1‑148 (MQ=255)
                                       gTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCaa                                                                                                              <  4:144013/149‑1 (MQ=255)
                                             tCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTTTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGGAGATCCGACATTACTCGTTCATCAACGTcag                                                                                                        >  4:535206/1‑149 (MQ=255)
                                             tCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCC‑TACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTcag                                                                                                        >  2:110015/1‑148 (MQ=255)
                                                gAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAg                                                                                                                                          >  1:263211/1‑112 (MQ=255)
                                                gAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAg                                                                                                                                          <  2:263211/112‑1 (MQ=255)
                                                gAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGaa                                                                                                                                      <  3:468194/116‑1 (MQ=255)
                                                gAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGaa                                                                                                                                      >  4:468194/1‑116 (MQ=255)
                                                gAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATAACTCGTTCAGCAACGTcagcag                                                                                                     >  3:778284/1‑149 (MQ=255)
                                                 aaaaGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCaga                                                                                                    <  4:564930/149‑1 (MQ=255)
                                                     gCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAgttc                                                                                                <  3:471496/149‑1 (MQ=255)
                                                                  gagaTTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTc                                                                                   <  3:9669/149‑1 (MQ=255)
                                                                      ttAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATc                                                                               <  3:62153/149‑1 (MQ=255)
                                                                       tAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCa                                                                              <  4:566029/149‑1 (MQ=255)
                                                                             tGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGc                                                                         <  2:457945/148‑1 (MQ=255)
                                                                                aTTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACAt                                                                     >  3:734103/1‑149 (MQ=255)
                                                                                aTTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACAt                                                                     >  3:386352/1‑149 (MQ=255)
                                                                                 ttGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAgttcgcgttcgc                                                                                        <  3:888030/129‑1 (MQ=255)
                                                                                 ttGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAgttcgcgttcgc                                                                                        >  4:888030/1‑129 (MQ=255)
                                                                                 ttGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACAt                                                                     >  2:259132/1‑148 (MQ=255)
                                                                                 ttGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATc                                                                    >  3:653989/1‑149 (MQ=255)
                                                                                     ttttCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCTATCAGCGGCACATCACcg                                                                >  3:341429/1‑149 (MQ=255)
                                                                                          gctcgctcgcGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACcgcgcg                                                            <  1:446175/148‑1 (MQ=255)
                                                                                          gctcgctcgcGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGa                                                           >  4:606709/1‑149 (MQ=255)
                                                                                           ctcgctcgcGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGAt                                                          >  3:568346/1‑149 (MQ=255)
                                                                                              gctcgcGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTc                                                        <  2:270262/148‑1 (MQ=255)
                                                                                                 cgcgCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCcaca                                                    <  3:852811/149‑1 (MQ=255)
                                                                                                    gcTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTCCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTc                                                 <  3:487070/149‑1 (MQ=255)
                                                                                                     cTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCC‑TACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCa                                                <  1:110015/148‑1 (MQ=255)
                                                                                                      ttttttAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAg                                               <  3:313426/149‑1 (MQ=255)
                                                                                                       tttttAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAg                                               >  1:790059/1‑148 (MQ=255)
                                                                                                          ttAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGgcgc                                           >  3:758271/1‑149 (MQ=255)
                                                                                                             gATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACAt                                                                     >  1:550715/1‑120 (MQ=255)
                                                                                                             gATTTGTCATCGTTGTTCCTGTGATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACAt                                                                     <  2:550715/120‑1 (MQ=255)
                                                                                                                ttGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGcaccac                                      >  2:428737/1‑148 (MQ=255)
                                                                                                                     aTCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGtt                                 <  1:860676/148‑1 (MQ=255)
                                                                                                                     aTCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTc                                >  3:500671/1‑149 (MQ=255)
                                                                                                                     aTCGTTGTTCCTGTTATTAGCCCCTTACCC‑CACTGCGGGTAAgg                                                                                                                                         >  3:334976/1‑44 (MQ=17)
                                                                                                                     aTCGTTGTTCCTGTTATTAGCCCCTTACCC‑CACTGCGGGTAAgg                                                                                                                                         <  4:334976/44‑1 (MQ=17)
                                                                                                                      tCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTc                                >  2:404295/1‑148 (MQ=255)
                                                                                                                      tCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCg                               >  3:95973/1‑149 (MQ=255)
                                                                                                                       cGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGGACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCgg                              <  4:707987/149‑1 (MQ=255)
                                                                                                                       cGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCgg                              <  4:386352/149‑1 (MQ=255)
                                                                                                                       cGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCg                               >  2:575324/1‑148 (MQ=255)
                                                                                                                           gttCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTAt                           <  2:790059/148‑1 (MQ=255)
                                                                                                                            ttCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTg                         >  3:5321/1‑149 (MQ=255)
                                                                                                                             tCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTg                         >  1:435047/1‑148 (MQ=255)
                                                                                                                                tGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAgg                     <  3:606709/149‑1 (MQ=255)
                                                                                                                                 gTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCa                                                >  1:745639/1‑121 (MQ=255)
                                                                                                                                 gTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCa                                                <  2:745639/121‑1 (MQ=255)
                                                                                                                                  ttattaGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTgg                        <  1:675067/144‑1 (MQ=255)
                                                                                                                                  ttattaGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTgg                        >  2:675067/1‑144 (MQ=255)
                                                                                                                                   tattaGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTg                   >  2:761184/1‑148 (MQ=255)
                                                                                                                                         tcccTTACCC‑CACTGCGGGTAAGGGGAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAgg             >  2:144718/2‑148 (MQ=255)
                                                                                                                                         tcccTTACCC‑CACTGCGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGt            >  3:603080/2‑149 (MQ=255)
                                                                                                                                         ccccTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAgg             >  2:819731/1‑148 (MQ=255)
                                                                                                                                         ccccTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAg              <  3:864866/147‑1 (MQ=255)
                                                                                                                                         ccccTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAg              >  4:864866/1‑147 (MQ=255)
                                                                                                                                         ccccTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGt            >  3:549206/1‑149 (MQ=255)
                                                                                                                                         ccccTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGt            >  4:722715/1‑149 (MQ=255)
                                                                                                                                           ccTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCGGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTcc                                                       <  4:126478/104‑1 (MQ=255)
                                                                                                                                           ccTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCGGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTcc                                                       >  3:126478/1‑104 (MQ=255)
                                                                                                                                           ccTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTa          >  3:379995/1‑149 (MQ=255)
                                                                                                                                            cTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTa          >  1:1918/1‑148 (MQ=255)
                                                                                                                                             ttACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTaaa        <  4:475027/149‑1 (MQ=255)
                                                                                                                                               aCCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAgc      <  4:653989/149‑1 (MQ=255)
                                                                                                                                               aCCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGAGTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAgc      <  4:548586/149‑1 (MQ=255)
                                                                                                                                                 ccGCAGTG‑GGGAAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACGTTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAgcgc    >  4:81992/1‑149 (MQ=255)
                                                                                                                                                   gCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGCCa  <  4:310739/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTG‑GGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGCCA  >  NC_000913/2121436‑2121731

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: ‑