Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,133,454 T→C 7.0% intergenic (‑57/+36) wcaA ← / ← wzc putative colanic acid biosynthesis glycosyl transferase/protein‑tyrosine kinase Wzc

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,133,4540TC 7.0% 299.0 / 9.3 100intergenic (‑57/+36)wcaA/wzcputative colanic acid biosynthesis glycosyl transferase/protein‑tyrosine kinase Wzc
Reads supporting (aligned to +/- strand):  ref base T (52/41);  new base C (2/5);  total (54/46)
Fisher's exact test for biased strand distribution p-value = 2.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

CTGGCGCAAAACCGATTTTATTGCCCGAATCGCCAGTTGTTGGCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGC  >  NC_000913/2133308‑2133601
                                                                                                                                                  |                                                                                                                                                   
cTGGCGCAAAACCGATTTTATTGCCCGAATCGCCAGTTGTTGGCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATc                                                                                                                                                  >  1:8255/1‑150 (MQ=255)
  ggCGCAAAACCGATTTTATTGCCCGAATCGCCAGTTGTTGGCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCg                                                                                                                                                <  2:1091432/150‑1 (MQ=255)
    cgcAAAACCGATTTTATTGCCCGAATCGCCAGTTGTTGGCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                              <  1:941267/150‑1 (MQ=255)
     gcAAAACCGATTTTATTGCCCGAATCGCCAGTTGTTGGCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGGCGGATAAGGCGTTCAGGCCGCATCCggga                                                                                                                                             <  2:703352/150‑3 (MQ=255)
     gcAAAACCGATTTTATTGCCCGAATCGCCAGTTGTTGGCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                             >  1:229011/1‑150 (MQ=255)
     gcAAAACCGATTTTATTGCCCGAATCGCCAGTTGTTGGCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                             <  2:65287/150‑1 (MQ=255)
     gcAAAACCGATTTTATTGCCCGAATCGCCAGTTGTTGGCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATAAGGCGTTCACGCCGCATCCGGca                                                                                                                                             <  1:96601/150‑2 (MQ=255)
      cAAAACCGATTTTATTGCCCGAATCGCCAGTTGTTGGCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATAAGGCGTTCACGCCGCATCCGGcaa                                                                                                                                            <  1:1172620/150‑3 (MQ=255)
      cAAAACCGATTTTATTGCCCGAATCGCCAGTTGTTGGCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATAAGGCGTTCACGCCGCATCCGGcaa                                                                                                                                            <  1:154807/150‑3 (MQ=255)
       aaaaCCGATTTTATTGCCCGAATCGCCAGTTGTTGGCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTa                                                                                                                                           <  1:1132543/150‑1 (MQ=255)
         aaCCGATTTTATTGCCCGAATCGCCAGTTGTTGGCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTa                                                                                                                                                    <  1:411088/139‑1 (MQ=255)
         aaCCGATTTTATTGCCCGAATCGCCAGTTGTTGGCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTa                                                                                                                                                    >  2:737340/1‑139 (MQ=255)
         aaCCGATTTTATTGCCCGAATCGCCAGTTGTTGGCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTa                                                                                                                                                    >  2:411088/1‑139 (MQ=255)
         aaCCGATTTTATTGCCCGAATCGCCAGTTGTTGGCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTa                                                                                                                                                    <  1:737340/139‑1 (MQ=255)
             gATTTTATTGCCCGAATCGCCAGTTGTTGGCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGc                                                                                                                                     <  2:1050909/150‑1 (MQ=255)
               ttttATTGCCCGAATCGCCAGTTGTTGGCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCgg                                                                                                                                   <  1:1528987/150‑1 (MQ=255)
                 ttATTGCCCGAATCGCCAGTTGTTGGCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCATACGGGGCGGTg                                                                                                                                 >  1:1077895/1‑150 (MQ=255)
                  tATTGCCCGAATCGCCAGTTGTTGGCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGc                                                                                                                                <  2:727211/150‑1 (MQ=255)
                      gCCCGAATCGCCAGTTGTTGGCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAAt                                                                                                                            >  2:1087276/1‑150 (MQ=255)
                       cccGAATCGCCAGTTGTTGGCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATg                                                                                                                           >  1:504315/1‑150 (MQ=255)
                        ccGAATCGCCAGTTGTTGGCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGc                                                                                                                          <  1:1073256/150‑1 (MQ=255)
                         cGAATCGCCAGTTGTTGGCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCa                                                                                                                         >  1:89988/1‑150 (MQ=255)
                          gAATCGCCAGTTGTTGGCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAg                                                                                                                        <  1:817165/150‑1 (MQ=255)
                          gAATCGCCAGTTGTTGGCGGTTCCAGGTCGGCATATAGATCGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAg                                                                                                                        <  1:1149420/150‑1 (MQ=255)
                           aaTCGCCAGTTGTTGGCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGGTGTTCTCCCGGCTTATGAGCGGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAgg                                                                                                                       <  1:440533/150‑1 (MQ=255)
                              cGCCAGTTGTTGGCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGcgt                                                                                                                    >  2:321844/1‑150 (MQ=255)
                              cGCCAGTTGTTGGCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGcgt                                                                                                                    >  2:553339/1‑150 (MQ=255)
                               gCCAGTTGTTGGCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGcgtc                                                                                                                   >  1:777718/1‑150 (MQ=255)
                               gCCAGTTGTTGGCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATAAGGCGTTCACGCCGCATCCGGcaag                                                                                                                                           >  1:778304/1‑123 (MQ=38)
                               gCCAGTTGTTGGCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATAAGGCGTTCACGCCGCATCCGGcaag                                                                                                                                           <  2:778304/126‑4 (MQ=38)
                                ccAGTTGTTGGCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACCGGGCGGTGCGAATGCAGGcgtcg                                                                                                                  >  2:1518268/1‑150 (MQ=255)
                                 cAGTTGTTGGCGGTTCCAGGTCGGCGTATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGcgtcgt                                                                                                                 >  1:539663/1‑150 (MQ=255)
                                 cAGTTGTTGGCGGTTCCAGGTCGGCGTATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGcgtcgt                                                                                                                 >  1:534414/1‑150 (MQ=255)
                                       ttgGCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTc                                                                                                           <  2:997395/150‑1 (MQ=255)
                                        tgGCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGTGGTGCGAATGCAGGCGTCGTTATTTCg                                                                                                          <  1:1338152/150‑1 (MQ=255)
                                        tgGCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTTGTTATTTCg                                                                                                          <  2:1340652/150‑1 (MQ=255)
                                        tgGCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTTGTTATTTCg                                                                                                          <  2:1335274/150‑1 (MQ=255)
                                        tgGCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCg                                                                                                          <  2:912785/150‑1 (MQ=255)
                                             gTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATcc                                                                                                     >  1:338716/1‑150 (MQ=255)
                                                 cAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACt                                                                                                 >  2:29015/1‑150 (MQ=255)
                                                   ggTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTta                                                                                               <  1:510193/150‑1 (MQ=255)
                                                      cGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTatat                                                                                            >  2:387920/1‑150 (MQ=255)
                                                        gCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTAtattc                                                                                          >  1:845257/1‑150 (MQ=255)
                                                         cATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTAtattcg                                                                                         >  2:172711/1‑150 (MQ=255)
                                                         cATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTAtattcg                                                                                         >  2:308454/1‑150 (MQ=255)
                                                             taGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTAtattcgtatt                                                                                     >  1:1378181/1‑150 (MQ=255)
                                                               gATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTAtattcgtattcg                                                                                   >  2:822751/1‑150 (MQ=255)
                                                                   gagaTCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGGTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGGCGTTATTTCGCATCCGACTTATATTCGTATTCGTAAt                                                                               <  1:238698/150‑1 (MQ=255)
                                                                     gaTCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAg                                                                             >  2:329668/1‑150 (MQ=255)
                                                                         aGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCAt                                                                         <  2:89988/150‑1 (MQ=255)
                                                                            ggATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAAt                                                                      >  2:425157/1‑150 (MQ=255)
                                                                            ggATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAAt                                                                      <  1:1217696/150‑1 (MQ=255)
                                                                                 gTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGa                                                                 >  1:209733/1‑150 (MQ=255)
                                                                                    tttCATAGGTCGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATAc                                                              >  1:730474/1‑150 (MQ=255)
                                                                                     ttCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATAcg                                                             <  2:1378181/150‑1 (MQ=255)
                                                                                       cATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGCATTATATTCGTATTAGTAATAGCCATAATCATGATAcgcg                                                           >  2:679957/1‑150 (MQ=255)
                                                                                       cATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATAcgcg                                                           <  1:1221521/150‑1 (MQ=255)
                                                                                         tAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCt                                                         <  2:1507567/150‑1 (MQ=255)
                                                                                         tAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCt                                                         >  2:1483338/1‑150 (MQ=255)
                                                                                           ggTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGACTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCTCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTgg                                                       >  2:1292060/1‑150 (MQ=255)
                                                                                              tgttCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGgcgc                                                    >  2:276329/1‑150 (MQ=255)
                                                                                                ttCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCgcgg                                                  >  2:1074655/1‑150 (MQ=255)
                                                                                                 tctcCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCGGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCgcggc                                                 >  2:1336969/1‑150 (MQ=255)
                                                                                                  ctcCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCgcggcg                                                >  2:1117713/1‑150 (MQ=255)
                                                                                                  ctcCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCgcggcg                                                >  2:471321/1‑150 (MQ=255)
                                                                                                    cccGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGa                                              >  1:1001968/1‑150 (MQ=255)
                                                                                                    cccGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGa                                              >  1:883130/1‑150 (MQ=255)
                                                                                                         cTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATg                                         >  1:278146/1‑150 (MQ=255)
                                                                                                          ttATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATgg                                        <  1:469254/150‑1 (MQ=255)
                                                                                                          ttATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATgg                                        <  1:387920/150‑1 (MQ=255)
                                                                                                          ttATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATgg                                        <  1:534503/150‑1 (MQ=255)
                                                                                                          ttATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATgg                                        <  2:1042676/150‑1 (MQ=255)
                                                                                                            aTGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAg                                      <  2:229011/150‑1 (MQ=255)
                                                                                                             tGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGt                                     >  2:1090641/1‑150 (MQ=255)
                                                                                                              gAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGtt                                    >  2:276613/1‑150 (MQ=255)
                                                                                                                  tGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGa                                <  2:1001968/150‑1 (MQ=255)
                                                                                                                   ggTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGaa                               >  2:1006447/1‑150 (MQ=255)
                                                                                                                       tgcttgcCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCAc                           <  1:1090641/150‑1 (MQ=255)
                                                                                                                        gcttgcCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCAcc                          <  1:1087276/150‑1 (MQ=255)
                                                                                                                        gcttgcCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCAcc                          >  1:482403/1‑150 (MQ=255)
                                                                                                                        gcttgcCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCAcc                          >  2:578968/1‑150 (MQ=255)
                                                                                                                        gcttgcCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTGATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGCAGCAcc                          >  1:118046/1‑150 (MQ=255)
                                                                                                                         cttgcCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCAccc                         <  2:777718/150‑1 (MQ=255)
                                                                                                                           tgcCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGaa                                             <  1:1080067/128‑1 (MQ=255)
                                                                                                                           tgcCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGaa                                             >  2:1080067/1‑128 (MQ=255)
                                                                                                                            gcCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCtt                      >  2:415979/1‑150 (MQ=255)
                                                                                                                            gcCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCtt                      >  2:410798/1‑150 (MQ=255)
                                                                                                                            gcCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCtt                      >  2:917278/1‑150 (MQ=255)
                                                                                                                              cGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTAGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTc                    <  2:845257/150‑1 (MQ=255)
                                                                                                                               ggataaggcgTTCACGCCGCATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCa                   >  1:176214/7‑150 (MQ=255)
                                                                                                                                 aTGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCAcc                 <  1:1117713/150‑1 (MQ=255)
                                                                                                                                  tGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCg                >  1:1424902/1‑150 (MQ=255)
                                                                                                                                       cgTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATa           >  2:471970/1‑150 (MQ=255)
                                                                                                                                            aCGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGtt      >  1:1313040/1‑150 (MQ=255)
                                                                                                                                              gCCTTATCCGTCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGtttt    <  1:275720/150‑1 (MQ=255)
                                                                                                                                               ccTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGcgcg                                                   >  2:624316/1‑102 (MQ=255)
                                                                                                                                               ccTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGcgcg                                                   <  1:624316/102‑1 (MQ=255)
                                                                                                                                               ccTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGCAATACCGTTTTg   >  2:1050541/1‑150 (MQ=255)
                                                                                                                                                cTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGc  <  1:276329/150‑1 (MQ=255)
                                                                                                                                                cTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGc  <  1:615440/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CTGGCGCAAAACCGATTTTATTGCCCGAATCGCCAGTTGTTGGCGGTTCCAGGTCGGCATATAGATTGAGATCAGCGGATTGTTTTTCATAGGTTGTTCTCCCGGCTTATGAGCTGGTTTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGCGGTGCGAATGCAGGCGTCGTTATTTCGCATCCGACTTATATTCGTATTCGTAATAGCCATAATCCTGATACGCGCTGGCGCGGCGGAAGATGGAGTTCAGAATCACCCCTTTCACCGGAATACCGTTTTGC  >  NC_000913/2133308‑2133601

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: