Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,134,397 A→G V419A (GTG→GCG)  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,134,3970AG98.5% 484.3 / ‑4.1 134V419A (GTG→GCG) wzccolanic acid production tyrosine‑protein kinase; autokinase; Ugd phosphorylase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (56/76);  total (57/77)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.73e-01

Only 100 of 136 total aligned reads displayed.

GGCATAGACGCTAATACCGTGTTCTTCCAGTACTTGCGGGCTTTCAATGCCGCGATTAAACAACGAGCGCAGCAGCACACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCACGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGTTGCATATAGACCTGCTGACCAGACTCGACATCGCGGGTCAGACGGACAATCTC  >  NC_000913/2134249‑2134542
                                                                                                                                                    |                                                                                                                                                 
ggCATAGACGCTAATACCGTGTTCTTCCAGTACTTGCGGGCTTTCAATGCCGCGATTAAACAACGAGCGCAGCAGCACACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGcgcg                                                                                                                                                 <  1:395248/151‑1 (MQ=255)
 gCATAGACGCTAATACCGTGTTCTTCCAGTACTTGCGGGCTTTCAATGCCGCGATTAAACAACGAGCGCAGCAGCACACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGcgcga                                                                                                                                                <  1:25605/151‑2 (MQ=255)
    tAGACGCTAATACCGTGTTCTTCCAGTACTTGCGGGCTTTCAATGCCGCGATTAAACAACGAGCGCAGCAGCACACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGcc                                                                                                                                               >  2:854006/1‑149 (MQ=255)
      gACGCTAATACCGTGTTCTTCCAGTACTTGCGGGCTTTCAATGCCGCGATTAAACAACGAGCGCAGCAGCACACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGc                                                                                                                                           <  1:149994/151‑1 (MQ=255)
         gCTAATACCGTGTTCTTCCAGTACTTGCGGGCTTTCAATGCCGCGATTAAACAACGAGCGCAGCAGCACACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCt                                                                                                                                          <  1:6946/149‑1 (MQ=255)
            aaTACCGTGTTCTTCCAGTACTTGCGGGCTTTCAATGCCGCGATTAAACAACGAGCGCAGCAGCACACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGCGGATAAGCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAg                                                                                                                                       <  4:434830/149‑1 (MQ=255)
            aaTACCGTGTTCTTCCAGTACTTGCGGGCTTTCAATGCCGCGATTAAACAACGAGCGCAGCAGCACACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGTGCGCCAGGCTGAGt                                                                                                                                      <  2:1141392/150‑1 (MQ=255)
            aaTACCGTGTTCTTCCAGTACTTGCGGGCTTTCAATGCCGCGATTAAACAACGAGCGCAGCAGCACACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAg                                                                                                                                       >  4:520832/1‑149 (MQ=255)
            aaTACCGTGTTCTTCCAGTACTTGCGGGCTTTCAATGCCGCGATTAAACAACGAGCGCAGCAGCACACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTg                                                                                                                                     >  2:811799/1‑151 (MQ=255)
             aTACCGTGTTCTTCCAGTACTTGCGGGCTTTCAATGCCGCGATTAAACAACGAGCGCAGCAGCACACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGa                                                                                                                                    <  3:450334/151‑1 (MQ=255)
             aTACCGTGTTCTTCCAGTACTTGCGGGCTTTCAATGCCGCGATTAAACAACGAGCGCAGCAGCACACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAACCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGt                                                                                                                                      <  2:649099/149‑1 (MQ=255)
               aCCGTGTTCTTCCAGTACTTGCGGGCTTTCAATGCCGCGATTAAACAACGAGCGCAGCAGCACACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGAtt                                                                                                                                  <  3:307430/151‑1 (MQ=255)
                ccGTGTTCTTCCAGTACTTGCGGGCTTTCAATGCCGCGATTAAACAACGAGCGCAGCAGCACACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGAtt                                                                                                                                  <  2:144251/150‑1 (MQ=255)
                 cGTGTTCTTCCAGTACTTGCGGGCTTTCAATGCCGCGATTAAACAACGAGCGCAGCAGCACACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGc                                                                                                                                >  2:495170/1‑151 (MQ=255)
                 cGTGTTCTTCCAGTACTTGCGGGCTTTCAATGCCGCGATTAAACAACGAGCGCAGCAGCACACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGc                                                                                                                                <  2:977324/151‑1 (MQ=255)
                      tcttcCAGTACTTGCGGGCTTTCAATGCCGCGATTAAACAACGAGCGCAGCAGCACACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGt                                                                                                                           >  2:1401124/1‑151 (MQ=255)
                      tcttcCAGAACTTGCGGGCTTTCAATGCCGCGATTAAACAACGAGCGCAGCAGCCCACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGCTTCAGCGCGCCAGGCTGAGTGATTGCCGGGt                                                                                                                           >  1:889587/1‑151 (MQ=255)
                          ccAGTACTTGCGGGCTTTCAATGCCGCGATTAAACAACGAGCGCAGCAGCACACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTcaa                                                                                                                        <  2:1577256/150‑1 (MQ=255)
                           cAGTACTTGCGGGCTTTCAATGCCGCGATTAAACAACGAGCGCAGCAGCACACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTcaaca                                                                                                                      <  1:1735734/151‑1 (MQ=255)
                           cAGTACTTGCGGGCTTTCAATGCCGCGATTAAACAACGAGCGCAGCAGCACACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTcaaca                                                                                                                      <  2:1183565/151‑1 (MQ=255)
                               aCTTGCGGGCTTCCAATGCCGCGATTAAACAACGAGCGCAGCAGCACACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATgc                                                                                                                  >  2:1665924/1‑151 (MQ=255)
                                  tGCGGGCTTTCAATGCCGCGATTAAACAACGAGCGCAGCAGCACACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATgcgc                                                                                                                >  1:764850/1‑150 (MQ=255)
                                    cGGGCTTTCAATGCCGCGATTAAACAACGAGCGCAGCAGCACACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGcaca                                                                                                             <  2:1086147/151‑1 (MQ=255)
                                     gggCTTTCAATGCCGCGATTAAACAACGAGCGCAGCAGCACACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACAt                                                                                                            <  2:616216/151‑1 (MQ=255)
                                      ggCTTTCAATGCCGCGATTAAACAACGAGCGCAGCAGCACACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATc                                                                                                           <  1:53978/151‑1 (MQ=255)
                                         tttCAATGCCGCGATTAAACAACGAGCGCAGCAGCACACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGcc                                                                                                        >  1:595518/1‑151 (MQ=255)
                                          ttCAATGCCGCGATTAAACAACGAGCGCAGCAGCACACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGc                                                                                                         <  1:1499547/149‑1 (MQ=255)
                                            cAATGCCGCGATTAAACAACGAGCGCAGCAGCACACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGAc                                                                                                     <  1:89305/151‑1 (MQ=255)
                                             aaTGCCGCGATTAAACAACGAGCGCAGCAGCACACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACg                                                                                                    <  1:1147544/151‑1 (MQ=255)
                                               tGCCGCGATTAAACAACGCGCGCAGCAGCACACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACgg                                                                                                   >  1:973382/1‑150 (MQ=255)
                                                gCCGCGATTAAACAACGAGCGCAGCAGCACACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTg                                                                                                 <  2:726551/151‑1 (MQ=255)
                                                 ccGCGATTAAACAACGAGCGCAGCAGCACACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGc                                                                                                >  4:465533/1‑151 (MQ=255)
                                                  cgcgATTAAACAACGAGCGCAGCAGCACACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCt                                                                                               <  4:387671/151‑1 (MQ=255)
                                                  cgcgATTAAACAACGAGCGCAGCAGCACACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCt                                                                                               <  2:1268400/151‑1 (MQ=255)
                                                   gcgATTAAACAACGAGCGCAGCAGCACACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTg                                                                                              <  4:388642/151‑1 (MQ=255)
                                                    cgATTAAACAACGAGCGCAGCAGCACACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTgg                                                                                             >  3:562967/1‑151 (MQ=255)
                                                    cgATTAAACAACGAGCGCAGCAGCACACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTgg                                                                                             >  2:988184/1‑151 (MQ=255)
                                                     gATTAAACAACGAGCGCAGCAGCACACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGc                                                                                            <  3:204024/151‑1 (MQ=255)
                                                       ttAAACAACGAGCGCAGCAGCACACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCt                                                                                          >  3:299434/1‑151 (MQ=255)
                                                        tAAACAACGAGCGCAGCAGCACACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCt                                                                                          <  2:393737/150‑1 (MQ=255)
                                                        tAAACAACGAGCGCAGCAGCACACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTc                                                                                         >  1:215306/1‑151 (MQ=255)
                                                        tAAACAACGAGCGCAGCAGCACACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTc                                                                                         >  1:60528/1‑151 (MQ=255)
                                                          aacaacGAGCGCAGCAGCACACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCgg                                                                                       >  1:1129780/1‑151 (MQ=255)
                                                             aacGAGCGCAGCAGCACACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTg                                                                                     >  2:796679/1‑150 (MQ=255)
                                                               cGAGCGCAGCAGCACACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGAtt                                                                                  <  1:681896/151‑1 (MQ=255)
                                                               cGAGCGCAGCAGCACACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGAtt                                                                                  >  1:401155/1‑151 (MQ=255)
                                                                gAGCGCAGCAGCACACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGAttt                                                                                 >  2:672990/1‑151 (MQ=255)
                                                                gAGCGCAGCAGCACACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGAttt                                                                                 <  2:581057/151‑1 (MQ=255)
                                                                   cgcAGCAGCACACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTc                                                                               >  4:326213/1‑150 (MQ=255)
                                                                     cagcagcaCACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGc                                                                            <  1:479164/151‑1 (MQ=255)
                                                                      agcagcaCACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCt                                                                           <  4:394124/151‑1 (MQ=255)
                                                                       gcagcaCACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTc                                                                          >  2:772861/1‑151 (MQ=255)
                                                                         agcaCACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCt                                                                        <  1:504837/151‑1 (MQ=255)
                                                                             cacCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGctc                                                                    <  2:1285422/151‑1 (MQ=255)
                                                                               cccACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCtt                                                                  <  1:1409043/151‑1 (MQ=255)
                                                                                  aCGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTAt                                                               >  2:314248/1‑151 (MQ=255)
                                                                                   cGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTAtt                                                              >  3:295744/1‑151 (MQ=255)
                                                                                         aaaGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTcagcag                                                        <  1:208898/151‑1 (MQ=255)
                                                                                         aaaGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCACGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTcagc                                                          >  2:1509997/1‑149 (MQ=255)
                                                                                          aaGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCTGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGt                                                       <  2:1193522/151‑1 (MQ=255)
                                                                                           aGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCCCATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGtt                                                      <  3:119484/151‑1 (MQ=255)
                                                                                            gCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGTTCTTTATTCAGCAGTTg                                                     <  1:598481/151‑1 (MQ=255)
                                                                                            gCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGTTg                                                     <  1:1304651/151‑1 (MQ=255)
                                                                                            gCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGTTg                                                     <  2:451301/151‑1 (MQ=255)
                                                                                            gCATCAGACCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGTTg                                                     >  2:208734/1‑151 (MQ=255)
                                                                                             cATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGTTGc                                                    >  2:851798/1‑151 (MQ=255)
                                                                                                 aGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGTTGCatat                                                >  2:1715972/1‑151 (MQ=255)
                                                                                                   gCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGTTGCATATAg                                              <  1:1685480/151‑1 (MQ=255)
                                                                                                    ccAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGTTGCATATAGa                                             <  1:1234221/151‑1 (MQ=255)
                                                                                                     cAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGTTGCATATAGAc                                            >  2:335828/1‑151 (MQ=255)
                                                                                                      aaGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGTTGCATATAGAcc                                           >  2:1193996/1‑151 (MQ=255)
                                                                                                       aGGATAATCGCTCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCACGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGTTGCATATAGa                                             <  2:1355167/148‑1 (MQ=255)
                                                                                                          aTAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGTTGCATATAGACctgct                                       <  1:1248962/151‑1 (MQ=255)
                                                                                                          aTAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGTTGCATATAGACctgct                                       <  1:546294/151‑1 (MQ=255)
                                                                                                               cgcgCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGTTGCATATAGACCTGCTGAcc                                   <  4:493072/150‑1 (MQ=255)
                                                                                                               cgcgCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGTTGCATATAGACCTGCTGAcc                                   <  1:1353128/150‑1 (MQ=255)
                                                                                                                 cgcCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGTTGCATATAGACCTGCTGAcc                                   >  1:823701/1‑148 (MQ=255)
                                                                                                                   ccGAGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGTTGCATATAGACCTGCTGACCAGACt                              >  4:336316/1‑151 (MQ=255)
                                                                                                                      aGGATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGTTGCATATAGACCTGCTGACCAGACTCg                            >  3:341765/1‑150 (MQ=255)
                                                                                                                       ggATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGTTGCATATAGACCTGCTGACCAGACTCGAc                          <  2:373302/151‑1 (MQ=255)
                                                                                                                        gATAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGTTGCATATAGACCTGCTGACCAGACTCGACa                         <  3:23113/151‑1 (MQ=255)
                                                                                                                         aTAATCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGTTGCATATAGACCTGCTGACCAGACTCGACAt                        >  1:1424397/1‑151 (MQ=255)
                                                                                                                           aaTCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGTTGCATATAGACCTGCTGACCAGACTCGACATcg                      <  3:515848/151‑1 (MQ=255)
                                                                                                                           aaTCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGTTGCATATAGACCTGCTGACCAGACTCGACATcg                      >  2:38688/1‑151 (MQ=255)
                                                                                                                            aTCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGTTGCATATAGACCTGCTGACCAGACTCGACATcgc                     <  4:534295/151‑1 (MQ=255)
                                                                                                                            aTCAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGTTGCATATAGACCTGCTGACCAGACTCGACATcgc                     >  2:1219545/1‑151 (MQ=255)
                                                                                                                              cAGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGTTGCATATAGACCTGCTGACCAGACTCGACATCGCgg                   <  1:495092/151‑1 (MQ=255)
                                                                                                                               aGCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGTTGCATATAGACCTGCTGACCAGACTCGACATCGCggg                  <  3:207549/151‑1 (MQ=255)
                                                                                                                                gCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGTTGCATATAGACCTGCTGACCAGACTCGACATCGCGGGt                 >  4:120537/1‑151 (MQ=255)
                                                                                                                                gCCCTTTCTTCGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGTTGCATATAGACCTGCTGACCAGACTCGACATCGCGGGt                 <  2:139468/151‑1 (MQ=255)
                                                                                                                                     ttcttcGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGTTGCATATAGACCTGCTGACCAGACTCGACATCGCGGGTCAg              <  2:100846/149‑1 (MQ=255)
                                                                                                                                      tcttcGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGTTGCATATAGACCTGCTGACCAGACTCGACATCGCGGGTCAGAc            <  2:168576/150‑1 (MQ=255)
                                                                                                                                        ttcGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGTTGCATATAGACCTGCTGACCAGACTCGACATCGCGGGTCAGACgg          >  3:449776/1‑150 (MQ=255)
                                                                                                                                        ttcGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGTTGCATATAGACCTGCTGACCAGACTCGACATCGCGGGTCAGACgg          >  1:398998/1‑150 (MQ=255)
                                                                                                                                          cGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGTTGCATATAGACCTGCTGACCAGACTCGACATCGCGGGTCAGACGGAc        <  2:1264121/150‑1 (MQ=255)
                                                                                                                                          cGGTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGTTGCATATAGACCTGCTGACCAGACTCGACATCGCGGGTCAGACGGACa       <  1:159998/151‑1 (MQ=255)
                                                                                                                                           ggTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGTTGCATATAGACCTGCTGACCAGACTCGACATCGCGGGTCAGACGGACa       >  2:1368386/1‑150 (MQ=255)
                                                                                                                                           ggTTTCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGTTGCATATAGACCTGCTGACCAGACTCGACATCGCGGGTCAGACGGACa       >  2:1360075/1‑150 (MQ=255)
                                                                                                                                             tttCAGCGCGCCAGGATGAGTAATTACCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCATCGGAGATTTTCAGCTCCTGGTCTTTATTCAGCAGTTGCATAAAGACCTGCTGACCAGACTCGACATCGCGGGTCAGACGGACAAtc    >  1:1161633/1‑151 (MQ=255)
                                                                                                                                               tCAGCGCGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGTTGCATATAGACCTGCTGACCAGACTCGACATCGCGGGTCAGACGGACAAtctc  >  2:550115/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GGCATAGACGCTAATACCGTGTTCTTCCAGTACTTGCGGGCTTTCAATGCCGCGATTAAACAACGAGCGCAGCAGCACACCCACGATAGAAAGCATCAGGCCAAGGATAATCGCGCCGAGGATAATCAGCCCTTTCTTCGGTTTCAGCACGCCAGGCTGAGTGATTGCCGGGTCAACAATGCGCACATCGCCGACGGTGCTGGCCTCGGTGATTTTCAGCTCCTGCTCTTTATTCAGCAGTTGCATATAGACCTGCTGACCAGACTCGACATCGCGGGTCAGACGGACAATCTC  >  NC_000913/2134249‑2134542

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