Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,124,534 T→C 5.2% I405V (ATT→GTT)  wcaI ← putative colanic biosynthesis glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,124,5340TC 5.2% 234.0 / 2.6 76I405V (ATT→GTT) wcaIputative colanic biosynthesis glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base T (31/41);  new base C (2/2);  total (33/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.38e-01

GGTGGTTTGCAGCATGGTGAGATCGCCTTTCAGGCATAAAAACTGCTTGGGATAAAGTACGCGGGAAAGCGGCCATAAGCGGCTACCGGAGCCACCTGCCATCACAACTGGATAGAGTTTCGACTGCGCCATAATTATCCCCGAATATCATTTATAAATTGACGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCGCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCG  >  NC_000913/2124389‑2124680
                                                                                                                                                 |                                                                                                                                                  
ggtggtTTGCAGCATGGTGAGATCGCCTTTCATGCATAAAAACTGCTTGGGATAAAGTACGCGGGAAAGCGGCCATAAGCGGCTACCGGAGCCACCTGCCATCTCAACTGGATAGAGTTGCGACTGCGCCATAATTATCCCCGAATATCa                                                                                                                                                <  2:168815/150‑1 (MQ=255)
  gggtTTGCAGCATGGTGAGATCGCCTTTCAGGCATAAAAACTGCTTGGGATAAAGTACGCGGGAAAGCGGCCATAAGCGGCTACCGGAGCCACCTGCCATCACAACTGGATAGAGTTTCGACTGCGCCATAATTATCCCCGAATATCAtt                                                                                                                                              <  2:1621051/149‑1 (MQ=255)
   ggtTTGCAGCATGGTGAGATCGCCTTTCAGGCATAAAAACTGCTTGGGATAAAGTACGCGGGAAAGCGGCCATAAGCGGCTACCGGAGCCACCTGCCATCACAACTGGATAGAGTTTCGACTGCGCCATAATTATCCCCGAATATCAttt                                                                                                                                             <  1:446836/150‑1 (MQ=255)
    gtTTGCAGCATGGTGAGATCGCCTTTCAGGCATAAAAACTGCTTGGGATAAAGTACGCGGGAAAGCGGCCATAAGCGGCTACCGGAGCCACCTGCCATCACAACTGGATAGAGTTTCGACTGCGCATAATTATCCCCGAATATCATTtat                                                                                                                                           <  2:656778/150‑1 (MQ=255)
     tttGCAGCATGGTGAGATCGCCTTTCAGGCATAAAAACTGCTTGGGATAAAGTACGCGGGAAAGCGGCCATAAGCGGCTACCGGAGCCACCTGCCATCACAACTGGATAGAGTTTCGACTGCGCCATAATTATCCCCGAATATCATTtat                                                                                                                                           <  1:71498/150‑1 (MQ=255)
     tttGCAGCATGGTGAGATCGCCTTTCAGGCATAAAAACTGCTTGGGATAAAGTACGCGGGAAAGCGGCCATAAGCGGCTACCGGAGCCACCTGCCATCACAACTGGATAGAGTTTCGACTGCGCCATAATTATCCCCGAATATCATTtat                                                                                                                                           >  1:1135119/1‑150 (MQ=255)
     tttGCAGCATGGTGAGATCGCCTTTCAGGCATAAAAACTGCTTGGGATAAAGTACGCGGGAAAGCGGCCATAAGCGGCTACCGGAGCCACCTGCCATCACAACTGGATAGAGTTTCGACTGCGCCATAATTATCCCCGAATATCATTtat                                                                                                                                           <  1:1077190/150‑1 (MQ=255)
        gcagcaTGGTGAGATCGCCTTTCAGGCATAAAAACTGCTTGGGATAAAGTACGCGGGAAAGCGGTCATAAGCGGCTACCGGAGCCACCTGCCATCACAACTGGATAGAGTTTCGACTGCGCCATAATTATCCCCGAATATCATTTATaaa                                                                                                                                        >  1:1256629/1‑150 (MQ=255)
         cagcaTGGTGAGATCGCCTTTCAGGCATAAAAACTGCTTGGGATAAAGTACGCGGGAAAGCGGCCATAAGCGGCTACCGGAGCCACCTGCCATCACAACTGGATAGAGTTTCGACTGCGCCATAATTATCCCCGAATATCATTTATAAAt                                                                                                                                       >  1:674692/1‑150 (MQ=255)
                gTGAGATCGCCTTTCAGGCATAAAAACTGCTTGGGATAAAGTACGCGGGAAAGCGGCCATAAGCGGCTACCGGAGCCACCTGCCATCACAACTGGATAGAGTTTCGACTGCGCCATAATTATCCCCGAATATCATTTATAAATTGACGTa                                                                                                                                >  2:1454132/1‑150 (MQ=255)
                   agaTCGCCTTTCAGGCATAAAAACTGCTTGGGATAAAGTACGCGGGAAAGCGGTCATAAGCGGCTACCGGAGCCACCTGCCATCACAACTGGATAGAGTTTCGACTGCGCCATAATTATCCCCGAATATCATTTATAAATTGACGTAaca                                                                                                                             <  1:1867991/150‑1 (MQ=255)
                   agaTCGCCTTTCAGGCATAAAAACTGCTTGGGATAAAGTACGCGGGAAAGCGGCCATAAGCGGCTACCGGAGCCACCTGCCATCACAACTGGATAGAGTTTCGACTGCGCCATAATTATCCCCGAATATCATTTATAAATTGACGTAcca                                                                                                                             >  1:1305383/1‑150 (MQ=255)
                     aTCGCCTTTCAGGCATAAAAACTGCTTGGGATAAAGTACGCGGGAAAGCGGCCATAAGCGGCTACCGGAGCCACCTGCCATCACAACTGGATAGAGTTTCGACTGCGCCATAATTATCCCCGAATATCATTTATAAATTGACGTAACACg                                                                                                                           <  1:1398843/150‑1 (MQ=255)
                           tttCAGGCATAAAAACTGCTTGGGATAAAGTACGCGGGAAAGCGGCCATAAGCGGCTACCGGAGCCACCTGCCATCACAACTGGATAGAGTTTCGACTGCGCCATAATTATCCCCGAATATCATTTATAAATTGACGTAACACGTTctct                                                                                                                     <  2:1967652/150‑1 (MQ=255)
                             tCAGGCATAAAAACTGCTTGGGATAAAGTACGCGGGAAAGCGGCCATAAGCGGCTACCGGAGCCACCTGCCATCACAACTGGATAGAGTTTCGACTGCGCCATAATTATCCCCGAATATCATTTATAAATTGACGTAACACGTTCTCttt                                                                                                                   >  2:680444/1‑150 (MQ=255)
                              cAGGCATAAAAACTGCTTGGGATAAAGTACGCGGGAAAGCGGCCATAAGCGGCTACCGGAGCCACCTGCCATCACAACTGGATAGAGTTTCGACTGCGCCATAATTATCCCCGAATATCATTTATAAATTGACGTAACACGTTCTCTTTa                                                                                                                  >  2:263112/1‑150 (MQ=255)
                                ggCATAAAAACTGCTTGGGATAAAGTACGCGGGAAAGCGGCCATAAGCGGCTACCGGAGCCACCTGCCATCACAACTGGATAGAGTTTCGACTGCGCCATAATTATCCCCGAATATCATTTATAAATTGACGTAACACGTTCTCTTTATc                                                                                                                >  2:736702/1‑150 (MQ=255)
                                  cATAAAAACTGCTTGGGATAAAGTACGCGGGAAAGCGGCCATAAGCGGCTACCGGAGCCACCTGCCATCACAACTGGATAGAGTTTCGACTGCGCCATAATTATCCCCGAATATCATTTATAAATTGACGTAACACGTTCTCTTTATCGa                                                                                                              >  2:1126929/1‑150 (MQ=255)
                                    tAAAAACTGCTTGGGATAAAGTACGCGGGAAAGCGGCCATAAGCGGCTACCGGAGCCACCTGCCATCACAACTGGATAGAGTTTCGACTGCGCCATAATTATCCCCGAATATCATTTATAAATTGACGTAACACGTTCTCTTTATCGAgc                                                                                                            >  2:279940/1‑150 (MQ=255)
                                         aCTGCTTGGGATAAAGTACGCGGGAAAGCGGCCATAAGCGGCTACCGGAGCCACCTGCCATCACAACTGGATAGAGTTTCGACTGCGCCATAATTATCCCCGAATATCATTTATAAATTGACGTAACACGTTCTCTTTATCGAgcgtgcg                                                                                                       <  2:444164/150‑1 (MQ=255)
                                            tcTTGTGATAAAGTACGCGGGAAAGCGGCCATATTCGGCTACAGGAGCCACCTGCCATCACAACTGGATAGAGTTTCGACTGCTCCATAATTATCCCCGAATATCATTTATAAATTGACGTAACACGTTCTCTTTATCGAGCGTGCGTTc                                                                                                    <  2:669281/149‑1 (MQ=255)
                                              ttGGGATAAAGTACGCGGGAAAGCGGCCATAAGCGGCTACCGGAGCCACCTGCCATCACAACTGGATAGAGTTTCGACTGCGCCTTAATTATCCCCGAATATCATTTATAAATTGACGTAACACGTTCTCTTTATCGAGCGTGCGTTCgg                                                                                                  <  2:1965591/150‑1 (MQ=255)
                                              ttGGGATAAAGTACGCGGGAAAGCGGCCATAAGCGGCTACCGGAGCCACCTGCCATCACAACTGGATAGAGTTTCGACTGCGCCATAATTATCCCCGAATATCATTTATAAATTGACGTAACACGTTCTCTTTATCGAGCGTGCGTTCgg                                                                                                  <  2:1458768/150‑1 (MQ=255)
                                               tGGGATAAAGTACGCGGGAAAGCGGCCATAAGCGGCTACCGGAGCCACCTGCCATCACAACTGGATAGAGTTTCGACTGCGCCATAATTATCCCCGAATATCATTTATAAATTGACGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGc                                                                                                 <  2:947999/150‑1 (MQ=255)
                                                gggATAAAGTACGCGGGAAAGCGGCCATAAGCGGCTACCGGAGCCACCTGCCATCACAACTGGATAGAGTTTCGACTGCGCCATAATTATCCCCGAATATCATTTATAAATTGACGTAAGACGTTCTCTTTATCGAGCGTGCGTTCGGCa                                                                                                <  1:1219914/150‑1 (MQ=255)
                                                  gATAAAGTACGCGGGAAAGCGGCCATAAGCGGCTACCGGAGCCACCTGCCATCACAACTGGATAGAGTTTCGACTGCGCCATAATTATCCCCGAATATCATTTATAAATTGACGTAACACGTTCTCTTTATCGAgcgtg                                                                                                         >  2:1891814/1‑139 (MQ=255)
                                                  gATAAAGTACGCGGGAAAGCGGCCATAAGCGGCTACCGGAGCCACCTGCCATCACAACTGGATAGAGTTTCGACTGCGCCATAATTATCCCCGAATATCATTTATAAATTGACGTAACACGTTCTCTTTATCGAgcgtg                                                                                                         <  1:1891814/139‑1 (MQ=255)
                                                  gATAAAGTACGCGGGAAAGCGGCCATAAGCGGCTACCGGAGCCACCTGCCATCACAACTGGATAGAGTTTCGACTGCGCCATAATTATCCCCGAATATCATTTATAAATTGACGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCata                                                                                              <  2:953793/150‑1 (MQ=255)
                                                   aTAAAGTACGCGGGAAAGCGGCCATAAGCGGCTACCGGAGCCACCTGCCATCACAACTGGATAGAGTTTCGACTGCGCCATAATTATCCCCGAATATCATTTATAAATTGACGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCatat                                                                                             <  1:1964783/150‑1 (MQ=255)
                                                    tAAAGTACGCGGGAAAGCGGCCATAAGCGGCTACCGGAGCCACCTGCCATCACAACTGGATAGAGTTTCGACTGCGCCATAATTATCCCCGAATATCATTTATAAATTGACGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATAtt                                                                                            <  1:1383356/150‑1 (MQ=255)
                                                                gAAAGCGGCCATAAGCGGCTACCGGAGCCACCTGCCATCACAACTGGATAGAGTTTCGACTGCGCCATAATTATCCCCGAATATCATTTATAAATTGACGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCg                                                                                <  1:1126929/150‑1 (MQ=255)
                                                                   aGCGGCCATAAGCGGCTACCGGAGCCACCTGCCATCACAACTGTATAGAGTTTCGACTGCGCCATAATTATCCCCGAATATCATTTATAAATTGACGTAACACGTTCTCTTTATCGATCGTGCGTTCGGCATATTCACGTTCCACCgtgt                                                                             <  2:1063511/150‑1 (MQ=255)
                                                                    gCGGCCATAAGCGGCTACCGGAGCCACCTGCCATCACAACTGGATAGAGTTTCGACTGCGCCATAATTATCCCCGAATATCATTTATAAATTGACGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCgtgtt                                                                            >  2:515589/1‑150 (MQ=255)
                                                                     cGGCCATAAGCGGCTACCGGAGCCACCTGCCATCACAACTGGATAGAGTTTCGACTGCGCCATAATTATCCCCGAATATCATTTATAAATTGACGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCgtgttg                                                                           <  2:1511069/150‑1 (MQ=255)
                                                                      ggCCATAAGCGGCTACCGGAGCCACCTGCCATCACAACTGGATAGAGTTTCGACTGCGCCATAATTATCCCCGAATATCATTTATAAATTGACGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCgtgttgt                                                                          <  2:330729/150‑1 (MQ=255)
                                                                           tAAGCGGCTACCGGAGCCACCTGCCATCACAACTGGATAGAGTTTCGACTGCGCCATAATTATCCCCGAATATCATTTATAAATTGACGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTg                                                                     >  1:327798/1‑150 (MQ=255)
                                                                            aaGCGGCTACCGGAGCCACCTGCCATCACAACTGGATAGAGTTTCGACTGCGCCATAATTATCCCCGAATATCATTTATAAATTGACGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTgg                                                                    <  2:490321/150‑1 (MQ=255)
                                                                             aGCGGCTACCGGAGCCACCTGCCATCACAACTGGATAGAGTTTCGACTGCGCCATAATTATCCCCGAATATCATTTATAAATTGACGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTggg                                                                   <  1:736702/150‑1 (MQ=255)
                                                                             aGCGGCTACCGGAGCCACCTGCCATCACAACTGGATAGAGTTTCGACTGCGCCATAATTATCCCCGAATATCATTTATAAATTGACGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTggg                                                                   >  1:1946825/1‑150 (MQ=255)
                                                                              gCGGCTACCGGAGCCACCTGCCATCACAACTGGATAGAGTTTCGACTGCGCCATAATTATCCCCGAATATCATTTATAAATTGACGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGgc                                                                  <  1:1024308/150‑1 (MQ=255)
                                                                              gCGGCTACCGGAGCCACCTGCCATCACAACTGGATAGAGAATCGACTGCGCCATAATTATCCCCGAATATCATTTATAAATTGACGTAACACGTTCTCTTTATCGAGCGTGCGATCGGCATAAACACGTGCCACCGTGTTGTGTTTGGgc                                                                  >  1:649722/1‑150 (MQ=255)
                                                                                   tACCGGAGCCACCTGCCATCACAACTGGATAGAGTTTCGACTGCGCCATAATTATCCCCGAATATCATTTATAAATTGACGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGcagcag                                                             >  1:1737951/1‑150 (MQ=255)
                                                                                    aCCGGAGCCACCTGCCATCACAACTGGATAGAGTTTCGACTGCGCCATAATTATCCCCGAATATCATTTATAAATTGACGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAgg                                                            >  1:858302/1‑150 (MQ=255)
                                                                                          gCCACCTGCCATCACAACTGGATAGAGTTTCGACTGCGCCATAATTATCCCCGAATATCATTTATAAATTGACGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCGcc                                                      >  2:63591/1‑150 (MQ=255)
                                                                                            cACCTGCCATCACAACTGGATAGAGTTTCGACTGCGCCATAATTATCCCCGAATATCATTTATAAATTGACGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCgtgttgt                                                                          >  2:1255916/1‑128 (MQ=255)
                                                                                            cACCTGCCATCACAACTGGATAGAGTTTCGACTGCGCCATAATTATCCCCGAATATCATTTATAAATTGACGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCgtgttgt                                                                          <  1:1255916/128‑1 (MQ=255)
                                                                                            cACCTGCCATCACAACTGGATAGAGTTTCGACTGCGCCATAATTATCCCCGAACATCATTTATAAATTGACGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCGCCTg                                                    >  2:419105/1‑150 (MQ=255)
                                                                                            cACCTGCCATCACAACTGGATAGAGTTTCGACTGCGCCATAATTATCCCCGAACATCATTTATAAATTGACGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCGCCTg                                                    >  2:411824/1‑150 (MQ=255)
                                                                                             aCCTGCCATCACAACTGGATAGAGTTTCGACTGCGCCATAATTATCCCCGAATATCATTTATAAATTGACGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCGCCTGa                                                   >  1:84126/1‑150 (MQ=255)
                                                                                               cTGCCATCACAACTGGATAGAGTTTCGACTGCGCCATAATTATCCCCGAATATCATTTATAAATTGACGTAACACGTTCTCTTTATCGTGCGTGCGTTCGGCATATTCACGTGCCAGCGTGTTGTGTTTGGGCAGCAGGAGCGCCTGACg                                                 <  1:384409/150‑1 (MQ=255)
                                                                                                tGCAATCACAACTGGATAGAGTTTCGACTGCGCCATAATTATCCCCGAATATCATTTATAAATTGACGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCGCCTGACgg                                                <  1:508954/150‑1 (MQ=255)
                                                                                                 gCCATCACAACTGGATAGAGTTTCGACTGCGCCATAATTATCCCCGAATATCATTTATAAATTGACGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCGCCTGACGGa                                               >  2:527522/1‑150 (MQ=255)
                                                                                                      cacaACTGGATAGAGTTTCGACTGCGCCATAATTATCCCCGAATATCATTTATAAATTGACGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCGCCTGACGGATcccc                                          >  2:1482779/1‑150 (MQ=255)
                                                                                                            tGGATAGAGTTTCGACTGCGCCATAATTATCCCCGAACATCATTTATAAATTGACGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCGCCTGACGGATCCCCGccacc                                    <  1:411824/150‑1 (MQ=255)
                                                                                                            tGGATAGAGTTTCGACTGCGCCATAATTATCCCCGAACATCATTTATAAATTGACGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCGCCTGACGGATCCCCGccacc                                    <  1:419105/150‑1 (MQ=255)
                                                                                                               aTAGAGTTTCGACTGCGCCATAATTATCCCCGAATATCATTTATAAATTGACGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCGCCTGACGGATCCCCGCCACCAgc                                 <  1:515589/150‑1 (MQ=255)
                                                                                                                 agagTTTCGACTGCGCCATAATTATCCCCGAATATCATTTATAAATTGACGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCGCCTGACGGATCCCCGCCACCAgcgc                               >  1:832216/1‑150 (MQ=255)
                                                                                                                   agTTTCGACTGCGCCATAATTATCCCCGAATATCATTTATAAATTGACGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCGCCTGACGGATCCCCGCCACCAGCGCCt                             <  1:279940/150‑1 (MQ=255)
                                                                                                                    gTTTCGACTGCGCCATAATTATCCCCGAATATCATTTATAAATTGACGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCGCCTGACGGATCCCCGCCACCAGCGCCTc                            <  1:680444/150‑1 (MQ=255)
                                                                                                                      ttCGACTGCGCCATAATTATCCCCGAATATCATTTATAAATTGACGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCGCCTGACGGATCCCCGCCACCAGCGCCTCGa                          >  1:694925/1‑150 (MQ=255)
                                                                                                                         gACTGCGCCATAATTATCCCCGAATATCATTTATAAATTGACGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCGCCTGACGGATCCCCGCCACCAGCGCCTCGACCg                       <  1:1702289/150‑1 (MQ=255)
                                                                                                                         gACTGCGCCATAATTATCCCCGAATATCATTTATAAATTGACGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCGCCTGACGGATCCCCGCCACCAGCGCCTCGACCg                       >  2:1056528/1‑150 (MQ=255)
                                                                                                                             gcgcCATAATTATCCCCGAATATCATTTATAAATTGACGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCACCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTc                   >  1:1480589/1‑150 (MQ=255)
                                                                                                                              cgcCATAATTATCCCCGAATATCATTTATAAATTGACGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCGCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTcc                  <  2:1135119/150‑1 (MQ=255)
                                                                                                                               gcCATAATTATCCCCGAATATCATTTATAAATTGACGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCGCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCg                 >  2:610885/1‑150 (MQ=255)
                                                                                                                                ccATAATTATCCCCGAATATCATTTATAAATTGACGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCGCCTGGCGGATCCCCGCCACCAGCGCCTCGACCGATTCCgg                >  2:26088/1‑150 (MQ=255)
                                                                                                                                ccATAATTATCCCCGAATATCATTTATAAATTGACGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCGCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCgg                <  2:283654/150‑1 (MQ=255)
                                                                                                                                 cATAATTATCCCCGAATATCATTTATAAATTGACGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCGCCTGACGGATCCCCGCCACCAGCGCCTCg                           <  2:1237272/138‑1 (MQ=255)
                                                                                                                                 cATAATTATCCCCGAATATCATTTATAAATTGACGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCGCCTGACGGATCCCCGCCACCAGCGCCTCg                           >  1:1237272/1‑138 (MQ=255)
                                                                                                                                      ttATCCCCGAATATCATTTATAAATTGACGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCGCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAAc          <  1:527522/150‑1 (MQ=255)
                                                                                                                                        aTCCCCGAATATCATTTATAAATTGACGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCGCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGc        <  2:674692/150‑1 (MQ=255)
                                                                                                                                          ccccGAATATCATTTATAAATTGACGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCGCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCaa      >  2:1718031/1‑150 (MQ=255)
                                                                                                                                           cccGAATATCATTTATAAATTGACGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGTAGCAGGAGCGCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCaaa     >  2:993502/1‑150 (MQ=255)
                                                                                                                                            ccGAATATCATTTATAAATTGACGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCGCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAAc    <  2:1946825/150‑1 (MQ=255)
                                                                                                                                             cGAATATCATTTATAAATTGACGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCGCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAAcc   <  2:257299/150‑1 (MQ=255)
                                                                                                                                              gAATATCATTTATAAATTGACGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCGCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCg  <  1:263112/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GGTGGTTTGCAGCATGGTGAGATCGCCTTTCAGGCATAAAAACTGCTTGGGATAAAGTACGCGGGAAAGCGGCCATAAGCGGCTACCGGAGCCACCTGCCATCACAACTGGATAGAGTTTCGACTGCGCCATAATTATCCCCGAATATCATTTATAAATTGACGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCGCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCG  >  NC_000913/2124389‑2124680

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