Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,133,601 C→T 5.5% E684E (GAG→GAA)  wzc ← colanic acid production tyrosine‑protein kinase; autokinase; Ugd phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,133,6010CT 5.5% 234.6 / 2.2 73E684E (GAG→GAA) wzccolanic acid production tyrosine‑protein kinase; autokinase; Ugd phosphorylase
Reads supporting (aligned to +/- strand):  ref base C (32/37);  new base T (2/2);  total (34/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.30e-01

ATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCACTTCTTTCAATGTGTTGACCGCATAACGCGCCACCATTAACGTGGTTCCGACATGACGACCAACAATTGCCGCATCGGTCACTGCCAGAATCGGCGGCGTATCAATCAACACCAGGTCATAGT  >  NC_000913/2133455‑2133749
                                                                                                                                                  |                                                                                                                                                    
aTCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGCt                                                                                                                                                     <  4:210403/148‑1 (MQ=255)
 tCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGc                                                                                                                                                      >  4:339796/1‑146 (MQ=255)
 tCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGc                                                                                                                                                      <  3:339796/146‑1 (MQ=255)
   cGGCCTACGGGGCGGTGCGAATGCAGGCATCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGCTCaaa                                                                                                                                                 <  2:316426/149‑1 (MQ=255)
                  tgagAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGt                                                                                                                                  <  8:418285/146‑1 (MQ=255)
                    cGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGttt                                                                                                                                >  7:295458/1‑149 (MQ=255)
                    cGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGtt                                                                                                                                 <  3:7393/148‑1 (MQ=255)
                     gAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTc                                                                                                                               <  8:116699/149‑1 (MQ=255)
                          cAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCa                                                                                                                             >  4:251478/1‑146 (MQ=255)
                          cAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCa                                                                                                                             <  3:251478/146‑1 (MQ=255)
                           aGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCActtc                                                                                                                         >  6:129821/1‑149 (MQ=255)
                           aGGCGGAGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCAGCGGAATACGGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCActtc                                                                                                                         <  6:90151/149‑1 (MQ=255)
                            ggCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCActtct                                                                                                                        <  5:129821/149‑1 (MQ=255)
                            ggCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCActtct                                                                                                                        <  6:205227/149‑1 (MQ=255)
                            ggCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGACATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCActtct                                                                                                                        >  5:205227/1‑149 (MQ=255)
                             gCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCActtctt                                                                                                                       >  2:103260/1‑149 (MQ=255)
                             gCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCActtc                                                                                                                         <  4:20141/147‑1 (MQ=255)
                              cgtcgtTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCTCTTCttt                                                                                                                      <  2:166972/149‑1 (MQ=255)
                               gtcgtTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCACTTCttt                                                                                                                      >  5:35696/1‑148 (MQ=255)
                                tcgtTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCACTTCTTTc                                                                                                                     >  2:62061/1‑148 (MQ=255)
                                tcgtTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCACTTCTTTCa                                                                                                                    <  1:62061/149‑1 (MQ=255)
                                  gtTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCACTTCTTTCaa                                                                                                                   <  6:35696/148‑1 (MQ=255)
                                  gtTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCACTTCTTTCAAt                                                                                                                  >  7:395714/1‑149 (MQ=255)
                                  gtTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCACTTCTTTCAAt                                                                                                                  <  8:295458/149‑1 (MQ=255)
                                         cGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCACTTCTTTCAATGTGTTGa                                                                                                           >  3:130177/1‑149 (MQ=255)
                                            aTCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTc                                                                                                                               >  7:413745/1‑126 (MQ=255)
                                            aTCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTc                                                                                                                               <  8:413745/126‑1 (MQ=255)
                                                     atatTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCACTTCTTTCAATGTGTTGACCGCATAAcgc                                                                                                <  1:282670/148‑1 (MQ=255)
                                                      tattcgtattcgTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCACTTCTTTCAATGTGTTGACCGCATAACgcgc                                                                                              >  1:410759/1‑149 (MQ=255)
                                                      tattcgtattcgTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCACTTCTTTCAATGTGTTGACCGCATAACgcgc                                                                                              >  5:61058/1‑149 (MQ=255)
                                                      tattcgtattcgTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCACTTCTTTCAATGTGTTGACCGCATAACgcgc                                                                                              <  3:275545/149‑1 (MQ=255)
                                                       attcgtattcgtAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCACTTCTTTCAATGTGTTGACCGCATAACgcgc                                                                                              >  2:236610/1‑148 (MQ=255)
                                                        ttcgtattcgtaATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCACTTCTTTCAATGTGTTGACCGCATAACGCGcca                                                                                            >  1:101506/1‑149 (MQ=255)
                                                         tcgtattcgtaATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCACTTCTTTCAATGTGTTGACCGCATAACGCGccac                                                                                           >  1:46626/1‑149 (MQ=255)
                                                             attcgtaATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCActtct                                                                                                                        <  4:172956/116‑1 (MQ=255)
                                                             attcgtaATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCActtct                                                                                                                        >  3:172956/1‑116 (MQ=255)
                                                                 gtaATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCACTTCTTTCAATGTGTTGACCGCATAACGCGCCACCATTAAc                                                                                    <  4:130177/148‑1 (MQ=255)
                                                                   aaTAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGTAATAACGTTTTGCTCAAAGCGGCTCAGACTGCTTTCCACTTCTTTCAATGTGTTGACCGCATAACGCCCCACCATTAACGTg                                                                                 >  6:29754/1‑149 (MQ=255)
                                                                          aTAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCACTTCTTTCAATGTGTTGACCGCATAACGCGCCACCATTAACGTGGTTCCGa                                                                          <  8:389950/149‑1 (MQ=255)
                                                                          aTAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCACTTCTTTCAATGTGTTGACCGCATAACGCGCCACCATTAACGTGGTTCCGa                                                                          <  7:155926/149‑1 (MQ=255)
                                                                          aTAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCACTTCTTTCAATGTGTTGACCGCATAACGCGCCACCATTAACGTGGTTCCGa                                                                          >  7:306104/1‑149 (MQ=255)
                                                                           tAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCACTTCTTTCAATGTGTTGACCGCATAACGCGCCACCATTAACGTGGTTCCGAc                                                                         <  1:293271/149‑1 (MQ=255)
                                                                                cTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCACTACTTTCAATGTGTTGCCCGCATAACGCGCCACCATTAACGTCGTTCCGACAAgac                                                                    >  8:248618/1‑149 (MQ=255)
                                                                                 tGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCACTTCTTTCAATGTGTTGACCGCATAACGCGCCACCATTAACGTGGTTCCGACATgacg                                                                   >  4:280836/1‑149 (MQ=255)
                                                                                  gATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCACTTCTTTCAATGTGTTGACCGCATAACGCGCCCCCATTAACGTGGTTCCGACATgacga                                                                  >  6:153936/1‑149 (MQ=255)
                                                                                   aTACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCACTTCTTTCAATGTGTTGACCGCATAACGCGCCACCATTAACGTGGTTCCGACATgacgac                                                                 <  2:46626/149‑1 (MQ=255)
                                                                                      cgcgCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCACTTCTTTCAATGTGTTGACCGCATAACGCGCCACCATTAACGTGGTTCCGACATGACGACcaa                                                              <  5:215341/149‑1 (MQ=255)
                                                                                      cgcgCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCACTTCTTTCAATGTGTTGACCGCATAACGCGCCACCATTAACGTGGTTCCGACATGACGACcaa                                                              <  6:61058/149‑1 (MQ=255)
                                                                                         gcTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCACTTCTTTCAATGTGTTGACCGCATAACGCGCCACCATTAACGTGGTTCCGACATGACGACcaacaa                                                           <  1:236610/149‑1 (MQ=255)
                                                                                             gcgcGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCACTTCTTTCAATGTGTTGACCGCATAACGCGCCACCATTAACGTGGTTCCGACATGACGACCAACAAtt                                                         <  3:9976/147‑1 (MQ=255)
                                                                                             gcgcGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCACTTCTTTCAATGTGTTGACCGCATAACGCGCCACCATTAACGTGGTTCCGACATGACGACCAACAATTg                                                        <  8:306104/148‑1 (MQ=255)
                                                                                              cgcgGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCACTTCTTTCAATGTGTTGACCGCATAACGCGCCACCATTAACGTGGTTCCGACATGACGACCAACAATTGcc                                                      <  3:280836/149‑1 (MQ=255)
                                                                                                 ggcggAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCACTTCTTTCAATGTGTTGACCGCATAACGCGCCACCATTAACGTTGTTCCGACATGACGACCAACAATTGCCGCa                                                   <  2:346720/149‑1 (MQ=255)
                                                                                                       aGATGGAGTTCAGCATCACCCCTTACACCGGAATACCGATTTGCTCAAAGCGGCTCAGACTCGTTTCAACTTCTTTCAATGTGTTGACCGCATAACGCGACACCATTAACGTGGTTCCGACATGACGACCAACAATTGCAGCAGCGGt                                              >  5:252765/1‑148 (MQ=255)
                                                                                                        gATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCACTTCTTTCAATGTGTTGACCGCATAACGCGCCACCATTAACGTGGTTCCGACATGACGACCAACAATTGCCGCATCGGTCa                                            <  2:410759/149‑1 (MQ=255)
                                                                                                         aTGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGTTCAAAGCGGCTCAGACTGGTTTCCACTTCTTTCAATGTGTTGAcc                                                                                                         >  3:390708/1‑87 (MQ=255)
                                                                                                         aTGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGTTCAAAGCGGCTCAGACTGGTTTCCACTTCTTTCAATGTGTTGAcc                                                                                                         >  3:146582/1‑87 (MQ=255)
                                                                                                         aTGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGTTCAAAGCGGCTCAGACTGGTTTCCACTTCTTTCAATGTGTTGAcc                                                                                                         <  4:146582/87‑1 (MQ=255)
                                                                                                         aTGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGTTCAAAGCGGCTCAGACTGGTTTCCACTTCTTTCAATGTGTTGAcc                                                                                                         <  4:390708/87‑1 (MQ=255)
                                                                                                             aGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCACTTCTTTCAATGTGTTGCCCGCATAACGCGCCACCATTAACGTGGTTCCGACATGACGACCAACAATTGCCGCATCGGTCACTGc                                        >  2:35958/1‑148 (MQ=255)
                                                                                                              gTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCACTTCTTTCACTGTGTTGACCGCATAACGCGCCACCATTAACGTGGTTCCGACATGACGACCAACAATTGCCGCATCGGTCACTGCCa                                      <  3:235147/149‑1 (MQ=255)
                                                                                                              gTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCACTTCTTTCAATGTGTTGACCGCATAACGCGCCACCATTAACGTGGTTCCGACATGACGACCAACAATTGCCGCATCGGTCACTGCCa                                      <  7:63794/149‑1 (MQ=255)
                                                                                                               ttCAGAATCACCCCTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCACTTCTTTCAATGTGTTGACCGCATAACGCGCCACCATTAACGTGGTTCCGACATGACGACCAACAATTGCCGCATCGGTCACTGc                                        >  3:70743/1‑146 (MQ=255)
                                                                                                                 cAGAATCACCCCTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCACTTCTTTCAATGTGTTGACCGCATAACGCGCCACCATTAACGTGGTTCCGACATGACGACCAACAATTGCCGCATCGGTCACTGCCAGaa                                   >  8:221258/1‑149 (MQ=255)
                                                                                                                   gAATCACCCCTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCACTTCTTTCAATGTGTTGACCGCATAACGCGCCACCATTAACGTGGTTCCGACATGACGACCAACAATTGCCGCATCGGTCACTGCCAGAATc                                 >  8:401580/1‑149 (MQ=255)
                                                                                                                      tCACCCCTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCACTTCTTTCAATGTGTTGACCGCATAACGCGCCACCATTAACGTGGTTCCGACATGACGACCAACAATTGCCGCATCGGTCACTGCCAGAATcggc                              <  7:401580/149‑1 (MQ=255)
                                                                                                                       cACCCCTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCACTTCTTTCAATGTGTTGACCGCATAACGCGCCACCATTAACGTGGTTCCGACATGACGACCAACAATTGCCGCATCGGTCACTGCCAGAATcggcg                             >  8:234655/1‑149 (MQ=255)
                                                                                                                        aCCCCTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCACTTCTTTCAATGTGTTGACCGCATAACGCGCCACCATTAACGTGGTTCCGACATGACGACCAACAATTGCCGCATCGGTCACTGCCAGAATcggcgg                            <  7:221258/149‑1 (MQ=255)
                                                                                                                          cccTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCACTTCTTTCAATGTGTTGACCGCATAACGCGCCACCATTAACGTGGTTCCGACATGACGACCAACAATTGCCGCATCGGTCACTGCCAGAATCGgcggcg                          >  1:89475/1‑149 (MQ=255)
                                                                                                                            cTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCACTTCTTTCAATGTGTTGACCGCATAACGCGCCACCATTAACGTGGTTCCGACATGACGACCAACAATTGCCGCATCGGTCACTGCCAGAATCGGCGGCGTa                        >  4:367713/1‑149 (MQ=255)
                                                                                                                             tttCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCACTTCTTTCAATGTGTTGACCGCATAACGCGCCACCATTAACGTGGTTCCGACATGACGACCAACAATTGCCGCATCGGTCACTGCCAGAATCGGCGGCGTat                       >  3:79901/1‑149 (MQ=255)
                                                                                                                                 aCCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCACTTCTTTCAATGTGTTGACCGCATAACGCGCCACCATTAACGTGGTTCCGACATGACGACCAACAATTGCCGCATCGGTCACTGCCAGAATCGGCGGCGTatcaat                   <  7:248618/149‑1 (MQ=255)
                                                                                                                                        tACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCACTTCTTTCAATGTGTTGACCGCATAACGCGCCACCATTAACGTGGTTCCGACATGACGACCAACAATTGCCGCATCGGTCACTGCCAGAATCGGCGGCGTATCAATCAacac             >  1:169330/1‑148 (MQ=255)
                                                                                                                                         aCCGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCACTTCTTTCAATGTGTTGACCGCATAACGCGCCACCATTAACGTGGTTCCGACATGACGACCAACAATTGCCGCATCGGTCACTGCCAGAATCGGCGGCGTATCAATCAacac             >  1:86949/1‑147 (MQ=255)
                                                                                                                                                 gCTCAAAGCGGCTCAGACTGGTTTCCACTTCTTTCAATGTGTTGACCGCATAACGCGCCACCATTAACGTGGTTCCGACATGACGACCAACAATTGCCGCATCGGTCACTGCCAGAATCGGCGGCGTATCAATCAACACCAGGTCATa    <  4:70743/148‑1 (MQ=255)
                                                                                                                                                  cTCAAAGCGGCTCAGACTGGTTTCCACTTCTTTCAATGTGTTGACCGCATAACGCGCCACCATTAACGTGGTTCCGACATGACGACCAACAATTGCCGCATCGGTCACTGCCAGAATCGGCGGCGTATCAATCAACACCAGGTCATAGt  >  1:165447/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
ATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGCTCAAAGCGGCTCAGACTGGTTTCCACTTCTTTCAATGTGTTGACCGCATAACGCGCCACCATTAACGTGGTTCCGACATGACGACCAACAATTGCCGCATCGGTCACTGCCAGAATCGGCGGCGTATCAATCAACACCAGGTCATAGT  >  NC_000913/2133455‑2133749

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